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

    雙組元推進(jìn)系統(tǒng)的部分可觀時間Petri網(wǎng)故障診斷

    2018-06-07 08:05:49劉久富張治國劉海陽王志勝
    關(guān)鍵詞:故障診斷故障系統(tǒng)

    劉久富, 張治國, 鄭 銳, 劉海陽, 楊 忠, 王志勝

    (1. 南京航空航天大學(xué)自動化學(xué)院, 江蘇 南京 210016;2. 東南大學(xué)電子信息工程學(xué)院, 江蘇 南京 210096)

    0 引 言

    航天器故障診斷具有故障樣本獲取困難、診斷經(jīng)驗不足、知識獲取存在瓶頸等問題,對實時性、準(zhǔn)確性和完整性提出了更高的要求[1-2]。為解決知識獲取的瓶頸問題,滿足在軌診斷的實時性要求,國外診斷系統(tǒng)一般采用基于定性模型的診斷方法,根據(jù)系統(tǒng)組件之間的連接關(guān)系(或參數(shù)依賴關(guān)系)建立故障診斷模型,根據(jù)系統(tǒng)模型和系統(tǒng)輸入,通過定性推理技術(shù),推導(dǎo)出系統(tǒng)正常情況下的預(yù)期行為并與實際輸出相比較,從而判斷系統(tǒng)是否存在故障。這種方法的優(yōu)點是實現(xiàn)簡單,實時性較好,能識別某些未預(yù)知故障等,缺點是分辨率低,對潛在的早期故障診斷能力不足[3]。為克服這種方法的不足,文獻(xiàn)[4]提出在故障傳播有向圖中加入故障傳播概率和時間,來辨識候選故障集合。文獻(xiàn)[5-6]提出基于符號有向圖(signed directed graph, SDG)模型和模糊推理相結(jié)合的故障診斷模型,通過部件故障診斷概率和故障傳播權(quán)重,對候選故障源進(jìn)行故障可能性排序,提高了診斷的分辨率。文獻(xiàn)[7]提出利用貝葉斯網(wǎng)絡(luò)推理算法和區(qū)間灰數(shù)運算規(guī)則,給出系統(tǒng)處于不同故障狀態(tài)的概率和部件重要度等可靠性特征信息,為不確定條件下的多態(tài)系統(tǒng)可靠性分析提供了完整的解決方案。

    基于定性模型的故障診斷方法主要是對可觀節(jié)點的故障診斷,本文提出一種基于部分可觀時間Petri網(wǎng)(Petri nets,PN)的故障診斷方法[8-15],利用可觀測節(jié)點的時序信息對不可觀節(jié)點進(jìn)行故障診斷。時間間隔比固定延時具有更強(qiáng)的表達(dá)能力,間隔既能反映事件發(fā)生的不確定性,又能刻畫事件發(fā)生的時限性。時間PN可以描述事件間異步、并發(fā)等邏輯關(guān)系,還可以反映事件時間層次的信息,是基本PN理論的功能擴(kuò)充[16-27]。本文在部分可觀時間PN和狀態(tài)類圖(state class graph,SCG)[28]的基礎(chǔ)上結(jié)合變遷屬性為節(jié)點分配標(biāo)簽,將其分為可觀和不可觀變遷,構(gòu)建包含時間信息的修正SCG(modified SCG,MSCG),把復(fù)雜的并行事件轉(zhuǎn)變成尋找滿足時序約束的路徑集合,能很好地解決系統(tǒng)的故障診斷問題。

    1 部分可觀時間PN的故障檢測問題

    1.1 部分可觀時間PN

    定義1PN定義為一個四元數(shù)組N=(P,T,Pre,Post),P是m個庫所的集合,T是n個變遷的集合,Pre:P×T→N和Post:P×T→N分別是輸入弧和輸出弧的關(guān)聯(lián)函數(shù)。定義矩陣C=Post-Pre為N的關(guān)聯(lián)矩陣,關(guān)聯(lián)矩陣C描述前一個變遷對后一個變遷的影響,代表系統(tǒng)的實時狀態(tài)。

    庫所P的標(biāo)識M:P→N為每個庫所分配一個非負(fù)整數(shù),每個庫所具有的資源個數(shù)記為M(p)。當(dāng)且僅當(dāng)M≥Pre(·,t) 變遷t邏輯使能,變遷t觸發(fā)后標(biāo)識M′=M+C(·,t)。標(biāo)識M代表系統(tǒng)資源的實時分布情況,是關(guān)聯(lián)矩陣C的一個重要子集。M處觸發(fā)使能的變遷集合為:A(M)={t∈T|M≥Pre(·,t)}。

    定義2時間PN(time PN, TPN)定義為Nd=(N,Q),其中,N=(P,T,Pre,Post)是PN的四元組,Q:T→Q×(Q∪{∞})是變遷初始時間約束相關(guān)聯(lián)的閉區(qū)間集合,變遷ti∈T對應(yīng)的Q由兩個自然數(shù)組成Q(ti)=(li,ui),li≥0,ui≥li,li≠∞當(dāng)且僅當(dāng)在時間區(qū)間[li,ui]內(nèi),變遷ti才能觸發(fā)。

    以τ0=0對應(yīng)標(biāo)識(M0,τ0)為初始點的TPNNd記為TPN系統(tǒng),簡記為〈Nd,M0〉。一個TPN系統(tǒng)中所有滿足M0[σ〉M的TTS集合σ記為Rt(Nd,M0)。TPN系統(tǒng)〈Nd,M0〉為有界系統(tǒng),?M∈Rt(Nd,M0), ?k∈N滿足M(p)≤k。

    給定TPN系統(tǒng)〈Nd,M0〉,M0[σ〉Mk且t1,t2∈σ滿足t1,t2∈A(Mk), 定義t1,t2的同步距離為sd(t1,t2),同步距離是并發(fā)變遷的最大觸發(fā)間隔。

    定義3標(biāo)簽PN(label PN, LPN)定義為〈N,M0,L〉,標(biāo)簽函數(shù)L:T→L∪{ε}為每個變遷分配一個標(biāo)簽L或空符號ε。不可觀標(biāo)簽和靜止標(biāo)簽集合為Tu,Tu={t∈T|L(t)=ε},可觀標(biāo)簽集合為To,To={t∈T|L(t)=L}。

    擴(kuò)展標(biāo)簽函數(shù)L:T*→L*滿足如下關(guān)系:

    tj∈To,則L(tj)=γ,γ∈L。

    tj∈Tu,則L(tj)=ε。

    σ∈T*∧tj∈T,則L(σtj)=L(σ)L(tj)。

    λ為空序列,則L(λ)=ε。

    1.2 故障檢測問題

    假設(shè)1任意變遷ti∈T在時間約束區(qū)間Q(ti)=(li,ui)(l,u∈Q)內(nèi)只發(fā)生一次且兩次觸發(fā)時間間隔不為零。

    假設(shè)2同一標(biāo)簽γ∈L可以關(guān)聯(lián)多個變遷,若 L(t1)=L(t2)=γ,變遷t1,t2相互獨立。

    L(σ)=δo,tl(σ)=τh,τh≤τ,

    t∈Tu:M[t〉,rt(M0,σ)≤τ-τh}

    (1)

    式中,rt(M0,σ)表示以M0為初始標(biāo)識經(jīng)過變遷時間序列σ后對應(yīng)變遷t觸發(fā)后剩余時間;τh是可觀測變遷最后觸發(fā)時間。式(1)保證在時刻τ-τh,標(biāo)識M可以使能變遷,但必須為可觀測變遷且變遷觸發(fā)時刻τeh滿足τeh≤τ。

    C(δo,τ)={M∈Nm|M0[σ〉M,σ∈∑(δo,τ)

    (2)

    S(δo,τ)={Mk,Θk|Mk∈C(δo,τ)

    (3)

    觀測序列中故障集出現(xiàn)的可信度稱為置信度,用Fb表示。給定故障變遷tf對應(yīng)的時間序列集合為

    (4)

    2 基于標(biāo)簽TPN的故障診斷方法

    2.1 構(gòu)建MSCG圖

    在SCG方法[11]的基礎(chǔ)上添加時間變量和時間約束,構(gòu)建基于標(biāo)簽TPN的MSCG。MSCG繼承于SCG,區(qū)別主要是MSCG引入了時間變量和時間約束的信息,可以估計狀態(tài)集合S(δo,τ)的信息。

    2.1.1 構(gòu)建MSCG算法

    構(gòu)建MSCG圖首先構(gòu)建修正狀態(tài)類樹(modified state class tree,MSCT),然后復(fù)制MSCT中標(biāo)記為B(C)=Z的節(jié)點。MSCT算法具體如下:

    算法1構(gòu)建MSCG

    輸入標(biāo)簽時間Petri網(wǎng)系統(tǒng)。

    輸出對應(yīng)的修正狀態(tài)類圖MSCG。

    步驟2構(gòu)建新節(jié)點Cq的過程。

    while存在標(biāo)記B(C)=E的節(jié)點do

    選擇一新節(jié)點Ck標(biāo)記B(C)=E

    for任意ti∈A(Mk)do

    ti在標(biāo)識Mk觸發(fā),Mq=Mk+C(·,ti)

    for?tr∈A(Mq)do

    iftr∈A(Mk)且Mk-Pre(:,ti)≥

    Pre(:,tr)then

    else

    if存和Cq相同節(jié)點then

    Cq節(jié)點標(biāo)記B(C)=Z

    else

    Cq節(jié)點標(biāo)記B(C)=E

    步驟3在MSCT中復(fù)制標(biāo)記為B(C)=Z的節(jié)點,構(gòu)建MSCG圖。

    算法1中步驟1初始化節(jié)點C0及其對應(yīng)變遷ti的初始時間約束Θ0,該時間約束Θ0僅是靜態(tài)時間約束,不同變遷序列中需根據(jù)步驟2修正時間約束。使用標(biāo)記B:C→{E,Z}避免重復(fù)構(gòu)建相同MSCT節(jié)點,在標(biāo)記B(C)=E節(jié)點構(gòu)建新的MSCT節(jié)點。

    步驟3將MSCT中標(biāo)記為B(C)=Z節(jié)點復(fù)制構(gòu)建對應(yīng)的MSCG圖。

    2.1.2 MSCG節(jié)點分析

    MSCG是有向圖,每個節(jié)點包含標(biāo)識M∈Rt(Nd,M0)和標(biāo)識M使能變遷的時間約束不等式Θk={li≤θk≤ui}:tk∈A(Mk)。節(jié)點之間邊界記為(t,γ,Δ∈[l*,u*]),變遷t∈T是源節(jié)點到目標(biāo)節(jié)點觸發(fā)的變遷,γ=L(t)是變遷t對應(yīng)的標(biāo)簽屬性,Δ∈[l*,u*]是源節(jié)點到目標(biāo)節(jié)點的時間約束。兩個節(jié)點有相同的標(biāo)識M和相同的時間約束Θ,則認(rèn)為其是等價的。

    定理1給定TPN系統(tǒng)〈Nd,M0〉,根據(jù)算法1構(gòu)建的MSCG圖節(jié)點是有限的。

    證畢

    2.2 基于TPN的狀態(tài)估計

    π=Cq0(Mq0,Θq0),ti1,L(ti1)

    Cq1(Mq1,Θq1),ti2,L(ti2)

    ?

    Cqk-1(Mqk-1,Θqk-1),tik,L(tik-1)

    (5)

    式中,out(Cq)為狀態(tài)集合Cq的下邊界;τh是可觀測變遷tiq∈log(δo)最后觸發(fā)時間。給定路徑π不滿足約束條件,則該路徑不可能觸發(fā)。

    狀態(tài)集合(Mk,Θk)是最后節(jié)點Cqk對應(yīng)狀態(tài),其中

    (7)

    (8)

    (9)

    證畢

    基于MSCG的狀態(tài)估計與診斷算法具體如下:

    算法2基于MSCG的狀態(tài)估計與診斷算法

    輸入TPN系統(tǒng)及對應(yīng)MSCG圖,TLS集合δo=(γi1,τ1)…(γih,τh),時間τ。

    步驟1初始化。δo=ε,S(δo,0)=φ。

    步驟2遍歷τ=0,δo=ε對應(yīng)的路徑π。

    for所有π∈Ⅱ路徑do

    S(ε,0)=S(ε,0)∪{(M,Θ)},(M,Θ)是路徑π最后節(jié)點。

    步驟3遍歷τ>0,δo=ε對應(yīng)的路徑π。

    if給定S(ε,τ),τ>0then

    for所有π∈Ⅱτdo

    S(ε,τ)=S(ε,τ)∪{(M,Θ′)},Θ′是根據(jù)式(6)、式(8)和式(9)得出的最后節(jié)點時間約束。

    forπ∈Ⅱdo

    S(δo,τ)=S(δo,τ)∪{(M,Θ)}

    for所有π∈Ⅱτdo

    S(δo,τ)=S(δo,τ)∪{(M,Θ′)}

    步驟5故障診斷。

    fori=1,…,rdo

    elseif?π∈Ⅱ,π′∈Ⅱ

    elseif?π∈Ⅱ,?

    給定路徑π,定義obs(π)為其對應(yīng)的可觀序列δo,式(5)中,obs(π)=L(ti1,ti2,…,tik)。

    2.3 復(fù)雜度分析

    算法1構(gòu)建MSCG圖,步驟2根據(jù)變遷ti∈A(Mq)在標(biāo)識Mk使能情況,修正其對應(yīng)時間約束,建立新節(jié)點Cq,最壞情況下變遷ti每觸發(fā)一次,所有變遷都使能,對應(yīng)計算復(fù)雜度為O(|T|2)。由假設(shè)1可知,給定截止時間τ,滿足δ0的路徑個數(shù)有限,定義lmax為MSCG中滿足δ0的最長路徑,所有滿足δ0的σ數(shù)最多為|T|lmax,即步驟2執(zhí)行次數(shù)最多為|T|lmax,所以算法1的計算復(fù)雜度為O(|T|lmax)。

    算法2中步驟2~步驟4判斷所有滿足obs(π)=ε路徑及obs(π)=δ0的路徑是否滿足式(6)~式(9),修正最后節(jié)點時間約束并將路徑放入相應(yīng)集合,對應(yīng)計算復(fù)雜度為O(|T|lmax)。步驟2~步驟4每次循環(huán)需根據(jù)式(6)判斷路徑是否滿足時間約束,其中式(6.1)需計算次數(shù)為1,式(6.2)建立不等式個數(shù)為Cqk使能變遷數(shù),所需計算次數(shù)最大為|T|。為保證變遷時間間隔Δ(l)滿足修正時間約束Θ′和δ0,式(6.3)和式(6.4)所需計算次數(shù)最多為3lmax|T|,式(6.5)建立不等式個數(shù)與σ中包含可觀變遷數(shù)量有關(guān),最壞情況整個路徑變遷都是可觀,所需計算次數(shù)為lmax,因此式(6)的計算復(fù)雜度為O(lmax|T|lmax),綜上,算法2的計算復(fù)雜度為O(lmax|T|lmax)。

    3 實例分析與驗證

    推進(jìn)系統(tǒng)的作用是實現(xiàn)軌道控制、姿態(tài)控制等。由于工作要求和零部件設(shè)計的結(jié)構(gòu)差異,推進(jìn)系統(tǒng)工作時各零部件存在時序偏差,若零部件出現(xiàn)故障便會導(dǎo)致其工作觸發(fā)時間超過設(shè)計要求。以某型號雙組元推進(jìn)系統(tǒng)為例驗證本文提出故障診斷算法[29-30]。

    3.1 推進(jìn)系統(tǒng)模型構(gòu)建

    圖1以推進(jìn)系統(tǒng)工作過程中關(guān)鍵節(jié)點為庫所,關(guān)鍵動作為變遷建立PN模型, 模擬推進(jìn)系統(tǒng)工作過程,各庫所、變遷含義如表1和表2所示。

    圖1 推進(jìn)系統(tǒng)PN模型 Fig1 PN model of propulsion system

    庫所庫所含義庫所庫所含義P0系統(tǒng)準(zhǔn)備就緒P9加排閥工作使能P11支路氦氣接通P10加排閥準(zhǔn)備就緒P22支路氦氣接通P11氧化劑管路測溫使能P31支路氣閥穩(wěn)定工作P12燃料管路測溫使能P42支路氣閥穩(wěn)定工作P13氧化劑管路溫度測試正常P5氧化劑儲箱壓力達(dá)標(biāo)P14燃料管路溫度測試正常P6燃料儲箱壓力達(dá)標(biāo)P15氧化劑閥門工作P7氧化劑閥門穩(wěn)定工作P16燃料準(zhǔn)備加注P8燃料劑閥門穩(wěn)定工作P17吹除使能

    表2 圖1中各變遷物理含義、可觀測性及設(shè)計動作時間

    雙組元推進(jìn)系統(tǒng)工作過程為:液路管路預(yù)冷完成后地面控制站發(fā)出動作命令,氣路電爆閥起爆,減壓閥和單向閥受氣體壓力影響開啟,儲箱開始增壓,管路測溫完成后液路閥門相繼打開,最后推進(jìn)器閥門打開,發(fā)動機(jī)工作。

    3.2 推進(jìn)系統(tǒng)故障診斷

    3.2.1 建立推進(jìn)系統(tǒng)MSCG

    構(gòu)建MSCG如圖2所示,每個節(jié)點由狀態(tài)集合S(δo,τ)組成,包括標(biāo)識M及對應(yīng)時間約束Θ,所有可能路徑為{t0t1t2t3t4t10t11t8,t0t1t2t3t4t11t10t8,t0t1t2t3t10t4t11t8,t0t1t2t4t3t10t11t8,t0t1t2t4t3t11t10t8,t0t1t2t4t11t3t10t8}等。圖2是部分序列觸發(fā)對應(yīng)的MSCG,初始節(jié)點C0對應(yīng)標(biāo)識M0和時間約束Θ0,可觀測變遷t0觸發(fā)使能;變遷t0在C0觸發(fā)得到節(jié)點C1,對應(yīng)標(biāo)識M1[011000000010000000],可觀測變遷t1,t2觸發(fā)使能Q(t1)=[2.5,6.7],Q(t2)=[4.0,7.9];變遷t1在C1觸發(fā)得到節(jié)點C2,對應(yīng)標(biāo)識M2[001100000010000000],為保證t2觸發(fā)必須滿足Δ1≤min{u1,u2},修正t2的時間約束為Q(t2)=[4.0-Δ1,7.9-Δ1]使能變遷為t2,t3;變遷t2在節(jié)點C2觸發(fā)得到節(jié)點C3,對應(yīng)標(biāo)識M3[000110000010000000],為保證可觀測變遷t3觸發(fā),必須滿足Δ2≤min{u2,u3},修正t3的時間約束為Q(t3)=[2.2-Δ2,6.5-Δ2]使能變遷為t3,t4,同理構(gòu)建節(jié)點C4~C20。

    3.2.2 基于MSCG推進(jìn)系統(tǒng)故障診斷結(jié)果

    雙組元推進(jìn)系統(tǒng)中加注閥故障會導(dǎo)致推進(jìn)系統(tǒng)故障,影響系統(tǒng)正常工作。本文PN模型中,無故障情況下氧化劑加注閥(變遷t5)會在氧化劑管路測溫完成后3.1~5.2 ms內(nèi)打開,若閥門故障或管路泄露便會導(dǎo)致閥門打開時間推遲到5.2~12.3 ms內(nèi)打開。根據(jù)可觀測變遷觸發(fā)時間和約束式(6)計算每個σ觸發(fā)可能性,利用算法2判斷系統(tǒng)的故障狀況。表3是根據(jù)實驗仿真數(shù)據(jù)及本文提出的算法進(jìn)行故障診斷得出的結(jié)果,其中,τ是采集有效數(shù)據(jù)的時間上限;N表示系統(tǒng)無故障;U表示系統(tǒng)可能故障;F表示系統(tǒng)故障。

    圖2 推進(jìn)系統(tǒng)MSCG Fig2 MSCG of the propulsion system

    可觀序列δo/ms截止時間τ/ms故障變遷Tif變遷log(σ)同步距離sd(t1,t2)置信度Fb診斷結(jié)果Γ(δo,τ,Tif)(t0,0)(t3,7.6)15.5t5{t0t1t2t3,t0t2t1t3}10N(t0,0)(t3,7.6)(t10,8.6)16.0t5{t0t1t2t3t10,t0t2t1t3t10}20N(t0,0)(t3,7.6)(t10,8.6)(t4,11.7)18.3t5{t0t1t2t3t10t4,t0t2t1t3t10t4}20N(t0,0)(t3,7.6)(t10,8.6)(t4,11.7)(t11,13.1)44.3t5{t0t1t2t3t10t4t11,t0t2t1t3t10t4t11,t0t1t2t3t10t4t5t11,t0t2t1t3t10t4t5t11}20N(t0,0)(t3,7.6)(t10,8.6)(t4,11.7)(t11,13.1)(t8,24.9)64.3t5{t0t1t2t3t10t4t11t5t7t6t8,t0t2t1t3t10t4t11t5t7t6t8,t0t1t2t3t10t4t5t11t7t6t8,t0t2t1t3t10t4t5t11t7t6t8,t0t1t2t3t10t4t11tft7t6t8,t0t2t1t3t10t4t11tft7t6t8}20.33U(t0,0)(t3,7.3)15.5t5{t0t1t2t3,t0t2t1t3}10N(t0,0)(t3,7.3)(t4,8.3)15.5t5{t0t1t2t3t4,t0t2t1t3t4}10N(t0,0)(t3,7.3)(t4,8.3)(t10,8.5)18.3t5{t0t1t2t3t4t10,t0t2t1t3t4t10}10N(t0,0)(t3,7.3)(t4,8.3)(t10,8.5)(t11,9.8)44.3t5{t0t1t2t3t4t10t11,t0t2t1t3t4t10t11}10N(t0,0)(t3,7.3)(t4,8.3)(t10,8.5)(t11,9.8)(t8,31.2)64.3t5{t0t1t2t3t4t10t11tft7t6t8,t0t2t1t3t4t10t11tft7t6t8}11F

    由表3可知,第一次試驗,可觀測變遷時間序列為(t0,0)(t3,7.6)(t10,8.6)(t4,11.7)(t11,13.1)(t8,24.9),其中,(t0,0)(t3,7.6)(t10,8.6)(t4,11.7)(t11,13.1)觸發(fā)后,由算法2計算所有滿足時間約束的TTS集合σ共4條,其中正常變遷σ為4條,故障變遷σ為0條,系統(tǒng)故障發(fā)生的置信度為0,同步距離sd(t1,t2)為2表明系統(tǒng)目前運行正常,工作在理想狀態(tài)??捎^測變遷(t8,24.9)觸發(fā)后,由算法2計算所有滿足時間約束的TTS集合σ共6條,其中正常變遷σ為4條,故障變遷σ為2條,系統(tǒng)故障發(fā)生的置信度為0.33,同步距離sd(t1,t2)為2表明氧化劑加注閥目前可能發(fā)生故障。

    第二次試驗,可觀測時間序列分別為(t0,0)(t3,7.3)(t4,8.3)(t10,8.5)(t11,9.8)(t8,31.2)}其中,可觀測變遷(t0,0)(t3,7.3)(t4,8.3)(t10,8.5)(t11,9.8)觸發(fā)后,由算法2計算所有滿足時間約束的TTS集合σ共2條,其中,正常變遷σ為2條,故障變遷σ為0條,系統(tǒng)故障發(fā)生的置信度為0,同步距離sd(t1,t2)為1,系統(tǒng)運行正常,工作在理想狀態(tài)??捎^測變遷(t8,31.2)觸發(fā)后,由算法2計算所有滿足時間約束的TTS集合σ共2條,其中,正常變遷σ為0條,故障變遷σ為2條,系統(tǒng)故障發(fā)生的置信度為1,同步距離sd(t1,t2)為1,表明氧化劑加注閥目前發(fā)生故障。

    3.3 算法比較

    表4是文獻(xiàn)[31]與本文算法性能比較。采用文獻(xiàn)[31]故障診斷算法至少需要8處傳感器,以保證任意兩個不可觀變遷不連接同時滿足算法可診斷性要求。為確定故障診斷候選集合并對系統(tǒng)進(jìn)行故障診斷,文獻(xiàn)[31]計算復(fù)雜度為3.93×108,其中,確定故障候選集合并對故障可能性排序的計算復(fù)雜度為2.61×107,根據(jù)候選集合對系統(tǒng)進(jìn)行故障診斷的計算復(fù)雜度為3.67×108。本文故障診斷算法在構(gòu)建MSCG過程中,不斷修正時間約束,使得系統(tǒng)對可觀變遷配置要求降低,只需6處傳感器便可進(jìn)行故障診斷。修正時間約束的應(yīng)用,使算法2的計算復(fù)雜度降低,構(gòu)建MSCT的計算復(fù)雜度為1.68×107,在MSCG中尋找滿足時間約束的路徑σ和判斷系統(tǒng)是否發(fā)生故障對應(yīng)計算復(fù)雜度為1.34×108,因此使用本文方法診斷系統(tǒng)的計算復(fù)雜度為1.51×108。修正時間約束的應(yīng)用,在實現(xiàn)故障診斷的同時,又減少了故障診斷的計算復(fù)雜度。

    表4 診斷性能比較

    4 結(jié) 論

    (1) 針對PN模型中不可觀節(jié)點故障診斷問題,本文基于SCG故障診斷方法和變遷時間約束信息,并根據(jù)變遷觸發(fā)關(guān)系修正時間約束,提出一種構(gòu)建MSCG的故障診斷方法,根據(jù)可觀序列集合S(δo,τ)和算法2在MSCG中尋找所有滿足時間約束的序列σ并判斷是否包含故障變遷,確定系統(tǒng)是否發(fā)生故障。

    (2) 建立雙組元推進(jìn)系統(tǒng)的部分可觀TPN故障診斷模型,并根據(jù)組件之間觸發(fā)關(guān)系修正相應(yīng)變遷時間約束,構(gòu)建對應(yīng)的MSCG,通過本文提出的算法診斷系統(tǒng)的故障狀態(tài)。仿真結(jié)果證明本文提出的故障診斷算法能夠滿足實際應(yīng)用要求。

    (3) 對于有嚴(yán)格時序要求且無法及時有效獲取內(nèi)部信息的系統(tǒng),本文提出的故障診斷算法能充分利用系統(tǒng)的時序信息,限制滿足可觀序列的無效路徑數(shù),減少故障診斷計算量,提高故障診斷效率。

    未來研究工作中,將對算法的可診斷性開展研究,構(gòu)建MSCG的復(fù)雜性依賴于系統(tǒng)結(jié)構(gòu)和變遷數(shù)目,當(dāng)遇到高復(fù)雜性系統(tǒng)時,所提方法還需進(jìn)一步改進(jìn)。

    參考文獻(xiàn):

    [1] YANG H, JIANG B, COCQUEMPOT V, et al. Spacecraft formation stabilization and fault tolerance: a state-varying switched system approach[J]. Systems & Control Letters, 2013, 62(9): 715-722.

    [2] ZHANG X, TANG L, DECASTRO J. Robust fault diagnosis of aircraft engines: a nonlinear adaptive estimation-based approach[J]. IEEE Trans.on Control Systems Technology,2013,21(3):861-868.

    [3] 姜連祥, 李華旺, 楊根慶, 等. 航天器自主故障診斷技術(shù)研究進(jìn)展[J].宇航學(xué)報, 2009, 30(4): 28-34.

    JIANG L X, LI H W, YANG G Q. A survey of spacecraft autonomous fault diagnosis research[J]. Journal of Astronautics, 2009, 30(4): 28-34.

    [4] SEMMEL G S, DAVIS S R, LEUCHT K W, et al. Space shuttle ground processing with monitoring agents[J]. IEEE Intelligent Systems, 2006, 21(1): 68-73.

    [5] NACI Z, GREG R L. Stability of gas pressure regulators[J]. Applied Mathematical Modeling, 2008, 32(1): 61-82.

    [6] 宋其江,徐敏強(qiáng),王日新.基于分層有向圖的航天器故障診斷[J].航空學(xué)報, 2009, 30(6): 1058-1062.

    SONG Q J, XU M Q, WANG R X. Spacecraft fault diagnosis based on hierarchical digraphs[J]. Acta Aeronautica et Astronautica Sinica, 2009, 30(6): 1058-1062.

    [7] HWANG I, KIM S, KIM Y, et al. A survey of fault detection, isolation, and reconfiguration methods[J]. IEEE Trans.on Control Systems Technology, 2010, 18(3): 636-653.

    [8] CABASINO M P, GIUA A, SEATZU C. Diagnosability of discrete event systems using labeled Petri nets[J]. IEEE Trans.on Automation Science and Engineering, 2014, 11(1): 144-153.

    [9] DECLERCK P, BONHOMMME P. State estimation of timed labeled Petri nets with unobservable transitions[J]. IEEE Trans.on Automation Science and Engineering, 2014, 11(1): 103-110.

    [10] BASILE F, CORDONE R, PIRODDI L. A branch and bound approach for the design of decentralized supervisors in Petri net models[J]. Automatica, 2015, 52(C):322-333.

    [11] SHEN Q, QIU J, LIU G, et al. Intermittent fault’s parameter framework and stochastic petri net based formalization model[J]. Eksploatacja I Niezawodnosc-Maintenance and Reliability, 2016, 18(2): 210-217.

    [12] WANG L, CHEN Q, GAO Z, et al. Knowledge representation and general Petri net models for power grid fault diagnosis[J]. IET Generation, Transmission & Distribution, 2015, 9(9): 866-873.

    [13] CABRAL F G, MOREIRA M V, DIENE O, et al. A Petri net diagnoser for discrete event systems modeled by finite state automata[J]. IEEE Trans.on Automatic Control, 2015, 60(1): 59-71.

    [14] BONHOMME P. Marking estimation of P-time Petri nets with unobservable transitions[J]. IEEE Trans.on Systems, Man, and Cybernetics: Systems, 2015, 45(3): 508-518.

    [15] SIMONA B, JAVIER C, JOSE M. Timing-failure risk assessment of UML design using time Petri net bound techniques[J]. IEEE Trans.on Industrial Informatics, 2011, 7(1): 90-104.

    [16] GIUA A, SEATZU C. Observability of place/transition nets[J]. IEEE Trans.on Automatic Control, 2002, 47(9):1424-1437.

    [17] CABASINO M P, GIUA A, SEATZU C. Fault detection for discrete event system using Petri nets with unobserved transitions[J]. Automatica, 2010, 46(9):1531-1539.

    [18] CORDONE R, PIRODDI L. Parsimonious monitor control of Petri net models of flexible manufacturing systems[J]. IEEE Trans.on System, Man, and Cybernetics: Systems, 2013,43(1): 215-221.

    [19] CHEN Y F, LI Z W, Barkaoui K, et al. New Petri nets structure and its application to optimal supervisory control: interval inhibitor arcs[J]. IEEE Trans.on Systems, Man, and Cybernetics: Systems, 2014, 44(10), 1384-1400.

    [20] CABASINO M P, GIUA A, SEATZU C, et al. Fault diagnosis of an ABS system using Petri nets[C]∥Proc.of the IEEE International Conference on Automation Science and Engineering, 2011: 24-27.

    [21] LEFEBVRE D. On-line fault diagnosis with partially observed Petri nets[J]. IEEE Trans.on Automatic Control, 2014, 59(7): 1919-1924.

    [22] BASILE F, CHIACCHIO P, TOMMASI D G. On diagnosability of Petri nets via integer linear programming[J]. Automatica, 2012, 48(9):2047-2058.

    [23] DARBY R. The dynamic response of pressure relief valves in vapor or gas service[J]. Journal of Loss Prevention in the Process Industries, 2013, 26(6): 1262-1268.

    [24] LI L, HADJICOSTIS C N. Minimun initial marking estimation in labeled Petri nets[J]. IEEE Trans.on Automatic Control, 2013, 58(1):198-203.

    [25] 王亮,呂衛(wèi)民,滕克難,等.基于Petri網(wǎng)的復(fù)雜設(shè)備健康狀態(tài)退化分析[J].系統(tǒng)工程與電子技術(shù),2014,36(10):1973-1981.

    WANG L, Lü W M, TENG K N, et al. Health degradation analysis of complex equipment based on Petri nets[J]. Systems Engineering and Electronics, 2014, 36(10): 1973-1981.

    [26] KO S. Performance comparison of covariance-assignment state estimators with intermittent observations[J]. International Journal of Control Automation and Systems, 2015, 13(6): 1391-1401.

    [27] CABASINO M P, LAFORTUNE S, SEATZU C. Optimal sensor selection for ensuring diagnosability in labeled petri nets[J]. Automatica, 2013, 49(8): 2373-2383.

    [28] LIME D, ROUX O H. Model checking of time Petri nets using the state class timed automaton[J]. Discrete Event Dynamic Systems: Theory and Application, 2006, 16(2): 179-205.

    [29] 魯峰,黃金泉,呂怡秋,等.基于非線性自適應(yīng)濾波的發(fā)動機(jī)氣路部件健康診斷方法[J].航空學(xué)報,2013,34(11):2529-2538.

    LU F, HUANG J Q, Lü Y Q, et al. Aircraft engine gas-path components health diagnosis based on nonlinear adaptive filters[J]. Acta Aeronautica et Astronautica Sinica, 2013, 34(11): 2529-2538.

    [30] 陳金豹, 翟國富, 王淑娟, 等.航天電子設(shè)備多余物檢測信號特征的影響因素分析[J]. 系統(tǒng)工程與電子技術(shù), 2013, 35(4): 889-894.

    CHEN J B, ZHAI G F, WANG S J, et al. Factors affecting characteristics of acoustic signals in particle impact noise detection for aerospace devices[J]. Systems Engineering and Electronics, 2013, 35(4): 889-894.

    [31] WANG X, CRISTIAN M, MANUEL S. Diagnosis of time Petri nets using fault diagnosis graph[J]. IEEE Trans.on Automatic Control, 2015, 60(9): 2321-2335.

    猜你喜歡
    故障診斷故障系統(tǒng)
    Smartflower POP 一體式光伏系統(tǒng)
    WJ-700無人機(jī)系統(tǒng)
    ZC系列無人機(jī)遙感系統(tǒng)
    北京測繪(2020年12期)2020-12-29 01:33:58
    故障一點通
    連通與提升系統(tǒng)的最后一塊拼圖 Audiolab 傲立 M-DAC mini
    奔馳R320車ABS、ESP故障燈異常點亮
    因果圖定性分析法及其在故障診斷中的應(yīng)用
    故障一點通
    江淮車故障3例
    基于LCD和排列熵的滾動軸承故障診斷
    欧美xxxx性猛交bbbb| 亚洲国产日韩欧美精品在线观看| 国产成人91sexporn| 精品不卡国产一区二区三区| 免费av观看视频| 亚洲精品第二区| 免费大片18禁| 男女边吃奶边做爰视频| 高清日韩中文字幕在线| h日本视频在线播放| videossex国产| 亚洲欧美成人精品一区二区| 久久人人爽人人爽人人片va| 中文字幕av在线有码专区| 亚洲精品色激情综合| 嘟嘟电影网在线观看| 欧美bdsm另类| 亚洲欧美成人精品一区二区| 国产久久久一区二区三区| 久久精品久久久久久久性| 婷婷色av中文字幕| 久久久久久久久久久丰满| 国产在线男女| 日韩欧美国产在线观看| 2021少妇久久久久久久久久久| 美女国产视频在线观看| 日本熟妇午夜| 久久久久九九精品影院| 搡老乐熟女国产| 亚洲熟女精品中文字幕| 97人妻精品一区二区三区麻豆| 午夜视频国产福利| 欧美丝袜亚洲另类| 天天躁日日操中文字幕| 嫩草影院精品99| 久久久亚洲精品成人影院| 波多野结衣巨乳人妻| 国内精品宾馆在线| 免费电影在线观看免费观看| 午夜精品一区二区三区免费看| 亚洲av.av天堂| videos熟女内射| 色综合亚洲欧美另类图片| 美女被艹到高潮喷水动态| 简卡轻食公司| 日韩伦理黄色片| 中文字幕人妻熟人妻熟丝袜美| 亚洲高清免费不卡视频| 亚洲av成人精品一二三区| 99热这里只有是精品50| 国产黄色小视频在线观看| 亚洲熟女精品中文字幕| 成人性生交大片免费视频hd| 国产高清国产精品国产三级 | 少妇熟女aⅴ在线视频| 亚洲欧美清纯卡通| 日韩欧美精品v在线| 亚洲精品中文字幕在线视频 | 99久国产av精品国产电影| 伦理电影大哥的女人| 亚洲精品456在线播放app| 日本一二三区视频观看| 99久久中文字幕三级久久日本| 不卡视频在线观看欧美| 久久精品国产亚洲网站| 精华霜和精华液先用哪个| 国产毛片a区久久久久| 免费观看的影片在线观看| 日日撸夜夜添| 欧美日韩一区二区视频在线观看视频在线 | 免费观看a级毛片全部| 国产极品天堂在线| 精品酒店卫生间| 看免费成人av毛片| 99久久九九国产精品国产免费| 白带黄色成豆腐渣| 免费av不卡在线播放| av福利片在线观看| 看十八女毛片水多多多| 国产白丝娇喘喷水9色精品| 成人综合一区亚洲| 在线免费观看的www视频| 网址你懂的国产日韩在线| 久久精品熟女亚洲av麻豆精品 | 国精品久久久久久国模美| 一区二区三区高清视频在线| 午夜精品一区二区三区免费看| 777米奇影视久久| 三级经典国产精品| 日韩中字成人| 亚洲精品乱码久久久久久按摩| 国产白丝娇喘喷水9色精品| 丝瓜视频免费看黄片| 久久人人爽人人片av| 国产美女午夜福利| 能在线免费观看的黄片| 69av精品久久久久久| 久久久a久久爽久久v久久| 国产黄a三级三级三级人| 久久久久久伊人网av| 蜜臀久久99精品久久宅男| 欧美区成人在线视频| 美女脱内裤让男人舔精品视频| 欧美变态另类bdsm刘玥| 一级爰片在线观看| 青春草视频在线免费观看| 人人妻人人看人人澡| 国产亚洲av片在线观看秒播厂 | 最近的中文字幕免费完整| 久久久久九九精品影院| 亚洲av免费高清在线观看| 亚洲va在线va天堂va国产| 国产在视频线精品| 综合色丁香网| 免费看av在线观看网站| 久久久精品94久久精品| 久久久国产一区二区| 亚洲国产欧美在线一区| 亚洲精品乱久久久久久| 国产精品国产三级国产av玫瑰| 色综合亚洲欧美另类图片| 久久久久精品久久久久真实原创| 肉色欧美久久久久久久蜜桃 | av黄色大香蕉| 久久精品国产亚洲av天美| 欧美成人精品欧美一级黄| 天堂网av新在线| 久久精品夜色国产| 国产69精品久久久久777片| 天堂俺去俺来也www色官网 | 久久精品国产亚洲av天美| 丝袜喷水一区| 成人av在线播放网站| 久久久久精品性色| 国产中年淑女户外野战色| av女优亚洲男人天堂| 亚洲高清免费不卡视频| 久久99热这里只频精品6学生| 狂野欧美激情性xxxx在线观看| 国产欧美另类精品又又久久亚洲欧美| 国产精品av视频在线免费观看| 又大又黄又爽视频免费| 精品久久国产蜜桃| 国产成人91sexporn| 综合色av麻豆| 噜噜噜噜噜久久久久久91| 女人十人毛片免费观看3o分钟| 一级毛片我不卡| 超碰97精品在线观看| 狠狠精品人妻久久久久久综合| 欧美成人精品欧美一级黄| 久久久久久久亚洲中文字幕| 久久精品熟女亚洲av麻豆精品 | 国产免费一级a男人的天堂| 午夜福利网站1000一区二区三区| 亚洲激情五月婷婷啪啪| 国产三级在线视频| 亚洲美女视频黄频| 亚洲欧洲日产国产| 日日摸夜夜添夜夜爱| 午夜爱爱视频在线播放| 色视频www国产| 高清毛片免费看| 日韩欧美一区视频在线观看 | videos熟女内射| 一本久久精品| 国语对白做爰xxxⅹ性视频网站| 欧美成人一区二区免费高清观看| 五月天丁香电影| 国产精品久久久久久精品电影小说 | 嫩草影院入口| 嫩草影院入口| 99热6这里只有精品| 日本黄色片子视频| 国产v大片淫在线免费观看| 成年av动漫网址| 色尼玛亚洲综合影院| 成人午夜精彩视频在线观看| 少妇的逼水好多| 日韩国内少妇激情av| 国产精品.久久久| .国产精品久久| 91午夜精品亚洲一区二区三区| 亚洲精品,欧美精品| 2022亚洲国产成人精品| 只有这里有精品99| 哪个播放器可以免费观看大片| 热99在线观看视频| 久久精品国产自在天天线| 婷婷色综合大香蕉| 欧美精品国产亚洲| 国精品久久久久久国模美| 亚洲国产精品成人综合色| 久久久久九九精品影院| 国产淫片久久久久久久久| 国产成人aa在线观看| a级毛色黄片| 日韩视频在线欧美| 赤兔流量卡办理| 99视频精品全部免费 在线| av国产免费在线观看| 欧美日韩在线观看h| av又黄又爽大尺度在线免费看| 99九九线精品视频在线观看视频| 美女内射精品一级片tv| 啦啦啦啦在线视频资源| 国产午夜精品一二区理论片| 久久久精品免费免费高清| 久久精品国产亚洲av涩爱| 精品一区二区三区人妻视频| eeuss影院久久| 搡女人真爽免费视频火全软件| 2018国产大陆天天弄谢| 久久97久久精品| 亚洲四区av| 麻豆乱淫一区二区| 建设人人有责人人尽责人人享有的 | 日本午夜av视频| 日本午夜av视频| 精品国产三级普通话版| 美女xxoo啪啪120秒动态图| 午夜老司机福利剧场| 国产精品久久久久久精品电影| freevideosex欧美| 免费播放大片免费观看视频在线观看| 日韩av免费高清视频| 日韩欧美国产在线观看| 成人亚洲精品av一区二区| 高清在线视频一区二区三区| 免费观看精品视频网站| 18+在线观看网站| 免费大片黄手机在线观看| 男人舔奶头视频| 亚洲av.av天堂| 亚洲欧美日韩卡通动漫| 一级黄片播放器| 国产 一区 欧美 日韩| 亚洲国产精品国产精品| 久久久久久久久久久丰满| 精品人妻一区二区三区麻豆| 秋霞伦理黄片| 精品人妻偷拍中文字幕| 国产亚洲91精品色在线| 美女黄网站色视频| av免费在线看不卡| 青青草视频在线视频观看| 亚洲aⅴ乱码一区二区在线播放| 亚洲综合精品二区| 在线 av 中文字幕| 精品亚洲乱码少妇综合久久| 亚洲,欧美,日韩| 全区人妻精品视频| 国产一区二区在线观看日韩| 又粗又硬又长又爽又黄的视频| 嫩草影院入口| 亚洲国产精品国产精品| 免费无遮挡裸体视频| 亚洲精品乱码久久久久久按摩| 最新中文字幕久久久久| 床上黄色一级片| 99久久九九国产精品国产免费| 国产伦理片在线播放av一区| 99视频精品全部免费 在线| 亚洲精品,欧美精品| 亚洲精品色激情综合| 少妇猛男粗大的猛烈进出视频 | 五月伊人婷婷丁香| 国产欧美另类精品又又久久亚洲欧美| 我要看日韩黄色一级片| 麻豆成人午夜福利视频| 水蜜桃什么品种好| 亚洲久久久久久中文字幕| a级一级毛片免费在线观看| 大陆偷拍与自拍| 亚洲精品成人久久久久久| 一级片'在线观看视频| 国产成人a∨麻豆精品| 免费黄网站久久成人精品| 午夜免费观看性视频| 亚洲在久久综合| 亚洲精品亚洲一区二区| 成人特级av手机在线观看| 午夜免费激情av| 亚洲国产精品国产精品| 看十八女毛片水多多多| 日日啪夜夜撸| 国产亚洲午夜精品一区二区久久 | 听说在线观看完整版免费高清| 好男人在线观看高清免费视频| 久久久久久久久久久丰满| 欧美不卡视频在线免费观看| 国产不卡一卡二| 国产综合精华液| 最近手机中文字幕大全| 国产麻豆成人av免费视频| 欧美一级a爱片免费观看看| 久久久久久久久久人人人人人人| 日产精品乱码卡一卡2卡三| 男女边吃奶边做爰视频| 乱系列少妇在线播放| 欧美高清性xxxxhd video| 亚洲av不卡在线观看| 久久精品国产鲁丝片午夜精品| 老师上课跳d突然被开到最大视频| 亚洲天堂国产精品一区在线| 中文字幕久久专区| 亚洲av二区三区四区| 菩萨蛮人人尽说江南好唐韦庄| 国产精品爽爽va在线观看网站| 久99久视频精品免费| 久久久a久久爽久久v久久| 欧美不卡视频在线免费观看| 亚洲精品亚洲一区二区| 国产av不卡久久| 国产精品嫩草影院av在线观看| 国产在线一区二区三区精| 国产精品不卡视频一区二区| 久久热精品热| 九色成人免费人妻av| 国产精品熟女久久久久浪| 99久久精品国产国产毛片| 九九爱精品视频在线观看| 神马国产精品三级电影在线观看| 日韩电影二区| freevideosex欧美| 久久97久久精品| 如何舔出高潮| 国产老妇女一区| 日日啪夜夜爽| 99久久九九国产精品国产免费| 久久精品人妻少妇| 亚洲美女视频黄频| 亚洲aⅴ乱码一区二区在线播放| 午夜激情福利司机影院| 国产免费又黄又爽又色| 少妇高潮的动态图| 亚洲欧美日韩无卡精品| 乱人视频在线观看| 欧美3d第一页| 青春草视频在线免费观看| 99视频精品全部免费 在线| 国产 亚洲一区二区三区 | 亚洲欧洲日产国产| 秋霞在线观看毛片| av黄色大香蕉| 久久精品国产亚洲av涩爱| 婷婷色麻豆天堂久久| 深爱激情五月婷婷| 成人二区视频| 在线天堂最新版资源| 床上黄色一级片| 成人性生交大片免费视频hd| 可以在线观看毛片的网站| 日韩欧美精品v在线| 亚洲精品影视一区二区三区av| 九九爱精品视频在线观看| 看十八女毛片水多多多| 精品久久久久久久久亚洲| 嫩草影院新地址| 久久久久久久国产电影| 亚洲av二区三区四区| 2022亚洲国产成人精品| 狂野欧美激情性xxxx在线观看| 99热6这里只有精品| 中文资源天堂在线| 久久国产乱子免费精品| 嫩草影院精品99| 建设人人有责人人尽责人人享有的 | 久久97久久精品| 女的被弄到高潮叫床怎么办| 97超碰精品成人国产| 汤姆久久久久久久影院中文字幕 | 亚洲欧美清纯卡通| 在线观看人妻少妇| 国产免费福利视频在线观看| 亚洲av免费在线观看| 乱系列少妇在线播放| 人妻制服诱惑在线中文字幕| 成人国产麻豆网| 午夜福利视频1000在线观看| 午夜精品一区二区三区免费看| ponron亚洲| 中国国产av一级| 成人亚洲精品av一区二区| 精品国产一区二区三区久久久樱花 | av国产免费在线观看| 久久99热这里只频精品6学生| 国产成人a区在线观看| 成人毛片a级毛片在线播放| 精华霜和精华液先用哪个| 欧美成人一区二区免费高清观看| 久久精品国产鲁丝片午夜精品| 欧美激情久久久久久爽电影| 人人妻人人澡欧美一区二区| 听说在线观看完整版免费高清| 久久精品国产亚洲网站| 国产成人一区二区在线| 男人舔女人下体高潮全视频| 国模一区二区三区四区视频| 精品人妻一区二区三区麻豆| 精品熟女少妇av免费看| 国产精品三级大全| 国产精品一二三区在线看| 日韩国内少妇激情av| 国产精品av视频在线免费观看| 一本一本综合久久| 国产成人freesex在线| 草草在线视频免费看| 亚洲精品国产av成人精品| 一级爰片在线观看| 欧美日韩一区二区视频在线观看视频在线 | 嫩草影院新地址| 成人二区视频| 国产男女超爽视频在线观看| 亚洲精品乱久久久久久| 国产精品嫩草影院av在线观看| 小蜜桃在线观看免费完整版高清| 91aial.com中文字幕在线观看| 久久6这里有精品| 观看免费一级毛片| 午夜福利网站1000一区二区三区| av一本久久久久| 欧美成人精品欧美一级黄| 欧美最新免费一区二区三区| 五月玫瑰六月丁香| 亚洲欧美一区二区三区国产| 99re6热这里在线精品视频| 26uuu在线亚洲综合色| 美女内射精品一级片tv| 成人一区二区视频在线观看| 国产午夜精品久久久久久一区二区三区| 欧美成人a在线观看| 欧美变态另类bdsm刘玥| 亚洲怡红院男人天堂| 非洲黑人性xxxx精品又粗又长| 黄色配什么色好看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 深夜a级毛片| 日本免费a在线| 黄色欧美视频在线观看| 亚洲精品视频女| 久久午夜福利片| 国产精品不卡视频一区二区| 亚洲欧美成人综合另类久久久| 精品欧美国产一区二区三| 精品少妇黑人巨大在线播放| 日日撸夜夜添| 国产精品伦人一区二区| 久久久久久久久久久免费av| 亚洲在线观看片| 久久99精品国语久久久| 性色avwww在线观看| 国产午夜精品论理片| 日日撸夜夜添| 三级国产精品片| 全区人妻精品视频| 91在线精品国自产拍蜜月| 天堂网av新在线| 欧美日韩在线观看h| 欧美zozozo另类| 国产免费又黄又爽又色| 精品人妻视频免费看| 精品久久久久久久末码| 国产又色又爽无遮挡免| 亚洲国产av新网站| 亚洲精品第二区| 久久鲁丝午夜福利片| 国产亚洲91精品色在线| 久久精品久久久久久久性| 美女脱内裤让男人舔精品视频| 亚洲精品国产成人久久av| 久久精品国产鲁丝片午夜精品| 22中文网久久字幕| 一区二区三区免费毛片| 亚洲国产日韩欧美精品在线观看| 99久久中文字幕三级久久日本| 一级av片app| 成年女人看的毛片在线观看| 日日干狠狠操夜夜爽| 国产色婷婷99| 精华霜和精华液先用哪个| 成年人午夜在线观看视频 | 精品少妇黑人巨大在线播放| 亚洲自偷自拍三级| 精品99又大又爽又粗少妇毛片| 亚洲精品视频女| 春色校园在线视频观看| 亚洲性久久影院| 日韩av免费高清视频| 亚洲美女搞黄在线观看| 国产av在哪里看| 亚洲国产精品成人久久小说| 伊人久久精品亚洲午夜| 久久久久久国产a免费观看| 高清在线视频一区二区三区| 免费高清在线观看视频在线观看| 2021少妇久久久久久久久久久| 欧美97在线视频| 国产不卡一卡二| av女优亚洲男人天堂| 老司机影院毛片| 国产一区二区亚洲精品在线观看| 精品久久久久久久久久久久久| 免费观看a级毛片全部| 毛片女人毛片| 欧美激情久久久久久爽电影| 国产精品三级大全| 国产大屁股一区二区在线视频| 亚洲成人精品中文字幕电影| 好男人在线观看高清免费视频| 岛国毛片在线播放| 精品一区在线观看国产| 99久久人妻综合| 日日啪夜夜撸| 久久久久网色| 色综合色国产| 亚洲精品久久午夜乱码| 国产激情偷乱视频一区二区| 色吧在线观看| 久久久久久久大尺度免费视频| 男人舔奶头视频| 国产美女午夜福利| 亚洲av日韩在线播放| 卡戴珊不雅视频在线播放| 国产成人午夜福利电影在线观看| 一级毛片 在线播放| 亚洲av不卡在线观看| 哪个播放器可以免费观看大片| 亚洲精品日本国产第一区| 少妇高潮的动态图| 欧美性猛交╳xxx乱大交人| 午夜精品一区二区三区免费看| 亚洲av成人精品一二三区| 日韩一区二区三区影片| 人妻夜夜爽99麻豆av| 国产免费视频播放在线视频 | 日本色播在线视频| 啦啦啦啦在线视频资源| 国产精品无大码| 国产精品国产三级国产专区5o| 精品一区二区三区人妻视频| 一级毛片我不卡| 精品酒店卫生间| 国产综合精华液| 精品久久久久久久人妻蜜臀av| 国产黄a三级三级三级人| 国产激情偷乱视频一区二区| 亚洲国产欧美在线一区| 蜜桃久久精品国产亚洲av| 街头女战士在线观看网站| 国产人妻一区二区三区在| 国产有黄有色有爽视频| 嫩草影院入口| 国产毛片a区久久久久| av在线天堂中文字幕| 日韩欧美一区视频在线观看 | 在线观看一区二区三区| 免费看av在线观看网站| 国产精品国产三级专区第一集| 成年人午夜在线观看视频 | 97人妻精品一区二区三区麻豆| 久久精品夜色国产| 国产永久视频网站| 在线免费十八禁| 日本黄大片高清| 亚洲av免费高清在线观看| 免费观看在线日韩| 午夜老司机福利剧场| 九草在线视频观看| 国产黄片视频在线免费观看| 欧美人与善性xxx| 偷拍熟女少妇极品色| 久久亚洲国产成人精品v| 色综合站精品国产| 亚洲av成人精品一区久久| 在线观看一区二区三区| a级毛色黄片| ponron亚洲| 亚洲国产精品成人久久小说| 三级国产精品欧美在线观看| 国产亚洲精品久久久com| 你懂的网址亚洲精品在线观看| 性插视频无遮挡在线免费观看| 亚洲在久久综合| 久久99热这里只频精品6学生| 中国美白少妇内射xxxbb| 久久久精品欧美日韩精品| 精品少妇黑人巨大在线播放| 国产免费视频播放在线视频 | 欧美另类一区| 中文资源天堂在线| 天美传媒精品一区二区| 少妇熟女aⅴ在线视频| 淫秽高清视频在线观看| 日韩在线高清观看一区二区三区| 久久精品综合一区二区三区| 欧美潮喷喷水| 黑人高潮一二区| 青春草视频在线免费观看| 午夜激情久久久久久久| 中文字幕久久专区| 麻豆乱淫一区二区| 国产爱豆传媒在线观看| 国产黄色视频一区二区在线观看| 亚洲国产高清在线一区二区三| 性色avwww在线观看| 老司机影院毛片| 国产综合精华液| 成人午夜精彩视频在线观看| 99九九线精品视频在线观看视频| 亚洲性久久影院| 国精品久久久久久国模美| 国产男女超爽视频在线观看| 日本一二三区视频观看| 永久网站在线|