• <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)的部分可觀時(shí)間Petri網(wǎng)故障診斷

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

    劉久富, 張治國(guó), 鄭 銳, 劉海陽(yáng), 楊 忠, 王志勝

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

    0 引 言

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

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

    1 部分可觀時(shí)間PN的故障檢測(cè)問題

    1.1 部分可觀時(shí)間PN

    定義1PN定義為一個(gè)四元數(shù)組N=(P,T,Pre,Post),P是m個(gè)庫(kù)所的集合,T是n個(gè)變遷的集合,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描述前一個(gè)變遷對(duì)后一個(gè)變遷的影響,代表系統(tǒng)的實(shí)時(shí)狀態(tài)。

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

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

    以τ0=0對(duì)應(yīng)標(biāo)識(shí)(M0,τ0)為初始點(diǎn)的TPNNd記為TPN系統(tǒng),簡(jiǎn)記為〈Nd,M0〉。一個(gè)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∪{ε}為每個(gè)變遷分配一個(gè)標(biāo)簽L或空符號(hào)ε。不可觀標(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 故障檢測(cè)問題

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

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

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

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

    (1)

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

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

    (2)

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

    (3)

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

    (4)

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

    2.1 構(gòu)建MSCG圖

    在SCG方法[11]的基礎(chǔ)上添加時(shí)間變量和時(shí)間約束,構(gòu)建基于標(biāo)簽TPN的MSCG。MSCG繼承于SCG,區(qū)別主要是MSCG引入了時(shí)間變量和時(shí)間約束的信息,可以估計(jì)狀態(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é)點(diǎn)。MSCT算法具體如下:

    算法1構(gòu)建MSCG

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

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

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

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

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

    for任意ti∈A(Mk)do

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

    for?tr∈A(Mq)do

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

    Pre(:,tr)then

    else

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

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

    else

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

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

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

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

    2.1.2 MSCG節(jié)點(diǎn)分析

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

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

    證畢

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

    π=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是可觀測(cè)變遷tiq∈log(δo)最后觸發(fā)時(shí)間。給定路徑π不滿足約束條件,則該路徑不可能觸發(fā)。

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

    (7)

    (8)

    (9)

    證畢

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

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

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

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

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

    for所有π∈Ⅱ路徑do

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

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

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

    for所有π∈Ⅱτdo

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

    forπ∈Ⅱdo

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

    for所有π∈Ⅱτdo

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

    步驟5故障診斷。

    fori=1,…,rdo

    elseif?π∈Ⅱ,π′∈Ⅱ

    elseif?π∈Ⅱ,?

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

    2.3 復(fù)雜度分析

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

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

    3 實(shí)例分析與驗(yàn)證

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

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

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

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

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

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

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

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

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

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

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

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

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

    可觀序列δo/ms截止時(shí)間τ/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可知,第一次試驗(yàn),可觀測(cè)變遷時(shí)間序列為(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計(jì)算所有滿足時(shí)間約束的TTS集合σ共4條,其中正常變遷σ為4條,故障變遷σ為0條,系統(tǒng)故障發(fā)生的置信度為0,同步距離sd(t1,t2)為2表明系統(tǒng)目前運(yùn)行正常,工作在理想狀態(tài)。可觀測(cè)變遷(t8,24.9)觸發(fā)后,由算法2計(jì)算所有滿足時(shí)間約束的TTS集合σ共6條,其中正常變遷σ為4條,故障變遷σ為2條,系統(tǒng)故障發(fā)生的置信度為0.33,同步距離sd(t1,t2)為2表明氧化劑加注閥目前可能發(fā)生故障。

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

    3.3 算法比較

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

    表4 診斷性能比較

    4 結(jié) 論

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

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

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

    未來(lái)研究工作中,將對(duì)算法的可診斷性開展研究,構(gòu)建MSCG的復(fù)雜性依賴于系統(tǒng)結(jié)構(gòu)和變遷數(shù)目,當(dāng)遇到高復(fù)雜性系統(tǒng)時(shí),所提方法還需進(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é)報(bào), 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é)報(bào), 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ā)動(dòng)機(jī)氣路部件健康診斷方法[J].航空學(xué)報(bào),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] 陳金豹, 翟國(guó)富, 王淑娟, 等.航天電子設(shè)備多余物檢測(cè)信號(hào)特征的影響因素分析[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無(wú)人機(jī)系統(tǒng)
    ZC系列無(wú)人機(jī)遙感系統(tǒng)
    故障一點(diǎn)通
    連通與提升系統(tǒng)的最后一塊拼圖 Audiolab 傲立 M-DAC mini
    奔馳R320車ABS、ESP故障燈異常點(diǎn)亮
    因果圖定性分析法及其在故障診斷中的應(yīng)用
    故障一點(diǎn)通
    江淮車故障3例
    基于LCD和排列熵的滾動(dòng)軸承故障診斷
    欧美色欧美亚洲另类二区| 99视频精品全部免费 在线| 色av中文字幕| 亚洲欧美精品综合久久99| 国产真实乱freesex| 亚洲欧美清纯卡通| 女人十人毛片免费观看3o分钟| 欧美性猛交╳xxx乱大交人| 国产精品一及| 久久亚洲精品不卡| 精品久久久久久久久亚洲 | 高清毛片免费观看视频网站| 日本a在线网址| 狂野欧美激情性xxxx在线观看| 嫩草影院精品99| 亚洲国产欧美人成| 欧美一区二区国产精品久久精品| 88av欧美| 亚州av有码| 亚洲一级一片aⅴ在线观看| 久久婷婷人人爽人人干人人爱| avwww免费| 男女做爰动态图高潮gif福利片| 色在线成人网| 亚洲成人中文字幕在线播放| 人人妻人人澡欧美一区二区| 嫩草影视91久久| 大型黄色视频在线免费观看| 日韩人妻高清精品专区| 国产乱人视频| 久久久午夜欧美精品| 亚洲真实伦在线观看| 日韩欧美精品v在线| 91av网一区二区| 亚洲国产精品sss在线观看| 日本色播在线视频| 婷婷精品国产亚洲av在线| av天堂在线播放| 国产三级在线视频| 99国产极品粉嫩在线观看| 午夜免费成人在线视频| 日韩人妻高清精品专区| 大型黄色视频在线免费观看| 欧美日本亚洲视频在线播放| 日本色播在线视频| 国产一区二区在线观看日韩| 亚洲美女黄片视频| 亚洲真实伦在线观看| 人人妻人人澡欧美一区二区| 午夜福利18| 婷婷精品国产亚洲av在线| 日本色播在线视频| 国产蜜桃级精品一区二区三区| 香蕉av资源在线| 99热只有精品国产| 亚洲av免费在线观看| 色哟哟·www| 国产精品女同一区二区软件 | 午夜爱爱视频在线播放| 成人午夜高清在线视频| 极品教师在线免费播放| 国产黄色小视频在线观看| 校园人妻丝袜中文字幕| 欧美xxxx黑人xx丫x性爽| 精品欧美国产一区二区三| 国产亚洲av嫩草精品影院| 亚洲av免费高清在线观看| 国产真实乱freesex| 干丝袜人妻中文字幕| 成人二区视频| 日本三级黄在线观看| 国产色婷婷99| 国产伦一二天堂av在线观看| 国语自产精品视频在线第100页| 又紧又爽又黄一区二区| 欧美三级亚洲精品| 国内少妇人妻偷人精品xxx网站| 2021天堂中文幕一二区在线观| 五月伊人婷婷丁香| 日韩欧美一区二区三区在线观看| 亚洲精品色激情综合| 久久国产乱子免费精品| 在线天堂最新版资源| 天堂动漫精品| 99久久中文字幕三级久久日本| 亚洲熟妇中文字幕五十中出| 哪里可以看免费的av片| 床上黄色一级片| а√天堂www在线а√下载| 久久久久九九精品影院| 国产精品美女特级片免费视频播放器| 国产真实乱freesex| 99热这里只有是精品在线观看| 欧美另类亚洲清纯唯美| 久久久久精品国产欧美久久久| 窝窝影院91人妻| 久久精品国产亚洲av涩爱 | 亚洲精华国产精华精| 精品久久久久久久末码| 免费观看的影片在线观看| av女优亚洲男人天堂| 久久草成人影院| 高清在线国产一区| 免费一级毛片在线播放高清视频| 两性午夜刺激爽爽歪歪视频在线观看| 国产老妇女一区| 国产精品久久久久久亚洲av鲁大| 国产主播在线观看一区二区| 在线免费观看的www视频| 国产一区二区三区在线臀色熟女| 狠狠狠狠99中文字幕| 禁无遮挡网站| 别揉我奶头~嗯~啊~动态视频| 成年女人毛片免费观看观看9| 国产 一区 欧美 日韩| netflix在线观看网站| 在线观看66精品国产| 亚洲 国产 在线| 欧美日韩精品成人综合77777| АⅤ资源中文在线天堂| 69人妻影院| 午夜福利视频1000在线观看| 精品人妻熟女av久视频| 成人亚洲精品av一区二区| 老女人水多毛片| 中文资源天堂在线| 波野结衣二区三区在线| ponron亚洲| 国内揄拍国产精品人妻在线| АⅤ资源中文在线天堂| 校园人妻丝袜中文字幕| 男人和女人高潮做爰伦理| aaaaa片日本免费| 精品久久久久久久久久久久久| 日韩在线高清观看一区二区三区 | 久久精品夜夜夜夜夜久久蜜豆| 国产成年人精品一区二区| 欧美性感艳星| 国产亚洲av嫩草精品影院| 极品教师在线视频| 久久人妻av系列| 久久久久免费精品人妻一区二区| 又紧又爽又黄一区二区| 神马国产精品三级电影在线观看| 国产精品综合久久久久久久免费| 久久精品国产亚洲av香蕉五月| 精品久久久久久久久av| а√天堂www在线а√下载| 国内精品宾馆在线| 成人鲁丝片一二三区免费| 亚洲av成人精品一区久久| 国产免费男女视频| 麻豆久久精品国产亚洲av| 成人国产麻豆网| 欧美性猛交黑人性爽| 99riav亚洲国产免费| 国产精华一区二区三区| 老司机午夜福利在线观看视频| 国产综合懂色| netflix在线观看网站| 亚洲成人免费电影在线观看| 人人妻,人人澡人人爽秒播| 亚洲av.av天堂| 一个人免费在线观看电影| 久久久色成人| 免费黄网站久久成人精品| 国产在线男女| 成人特级黄色片久久久久久久| 国产女主播在线喷水免费视频网站 | xxxwww97欧美| 亚洲人成伊人成综合网2020| 日韩 亚洲 欧美在线| 成年女人毛片免费观看观看9| 欧美最黄视频在线播放免费| 欧美精品啪啪一区二区三区| 亚洲一区二区三区色噜噜| 男人舔女人下体高潮全视频| 乱人视频在线观看| 午夜老司机福利剧场| 国产视频内射| 精品一区二区三区av网在线观看| 亚洲,欧美,日韩| 国产亚洲精品久久久久久毛片| 久久热精品热| 我的女老师完整版在线观看| 日日干狠狠操夜夜爽| 美女 人体艺术 gogo| 国内精品久久久久精免费| 亚洲最大成人中文| 最好的美女福利视频网| 日韩欧美三级三区| 欧美黑人巨大hd| xxxwww97欧美| 看片在线看免费视频| 久久香蕉精品热| 无遮挡黄片免费观看| 亚洲国产欧美人成| 桃红色精品国产亚洲av| 在线免费观看不下载黄p国产 | 亚洲精品成人久久久久久| 在线免费观看的www视频| 亚洲精品影视一区二区三区av| 成年免费大片在线观看| 免费一级毛片在线播放高清视频| 在线a可以看的网站| 日本五十路高清| 久久精品影院6| 成人无遮挡网站| 国产精品久久久久久久久免| 成人欧美大片| 成人特级黄色片久久久久久久| 男女做爰动态图高潮gif福利片| 日本a在线网址| aaaaa片日本免费| 赤兔流量卡办理| 欧美日韩黄片免| 免费电影在线观看免费观看| 搡老妇女老女人老熟妇| 99精品久久久久人妻精品| 国产精品美女特级片免费视频播放器| 一个人观看的视频www高清免费观看| 日韩中文字幕欧美一区二区| 黄色视频,在线免费观看| 久久久久久久亚洲中文字幕| 免费大片18禁| 亚洲av.av天堂| 国产蜜桃级精品一区二区三区| 亚洲内射少妇av| 国产精华一区二区三区| 99热这里只有精品一区| 看免费成人av毛片| 搡老岳熟女国产| 嫩草影院精品99| 可以在线观看的亚洲视频| 观看免费一级毛片| 国产精品国产高清国产av| 亚洲精品在线观看二区| 88av欧美| 亚洲欧美清纯卡通| 精品99又大又爽又粗少妇毛片 | 有码 亚洲区| 真人一进一出gif抽搐免费| 一卡2卡三卡四卡精品乱码亚洲| 午夜免费男女啪啪视频观看 | 国产爱豆传媒在线观看| 美女高潮喷水抽搐中文字幕| 99国产精品一区二区蜜桃av| 97超级碰碰碰精品色视频在线观看| 亚洲在线观看片| 亚洲av二区三区四区| 婷婷色综合大香蕉| 在线观看美女被高潮喷水网站| 热99re8久久精品国产| 啦啦啦观看免费观看视频高清| 久久精品91蜜桃| 我的女老师完整版在线观看| 嫩草影院入口| 亚洲美女搞黄在线观看 | 国产精品99久久久久久久久| 人人妻,人人澡人人爽秒播| 精品久久国产蜜桃| 国产亚洲av嫩草精品影院| 真实男女啪啪啪动态图| 久久久久国产精品人妻aⅴ院| 国产午夜精品久久久久久一区二区三区 | 女的被弄到高潮叫床怎么办 | 制服丝袜大香蕉在线| 欧美3d第一页| 搡老熟女国产l中国老女人| 99国产极品粉嫩在线观看| 久久久久久大精品| 最近在线观看免费完整版| 亚洲欧美日韩高清在线视频| ponron亚洲| 综合色av麻豆| 午夜激情欧美在线| 免费观看精品视频网站| 全区人妻精品视频| 国产精品国产三级国产av玫瑰| 国模一区二区三区四区视频| 亚洲人成网站在线播| 黄色日韩在线| 久久99热这里只有精品18| 亚洲七黄色美女视频| 99久国产av精品| 久久久久免费精品人妻一区二区| 2021天堂中文幕一二区在线观| 乱人视频在线观看| 久久久国产成人精品二区| 在线观看舔阴道视频| 神马国产精品三级电影在线观看| 我要看日韩黄色一级片| 国产美女午夜福利| 成人欧美大片| 97碰自拍视频| or卡值多少钱| 天堂动漫精品| 国产精华一区二区三区| 中文字幕人妻熟人妻熟丝袜美| 国语自产精品视频在线第100页| 国产免费av片在线观看野外av| 中国美女看黄片| 又爽又黄a免费视频| 夜夜看夜夜爽夜夜摸| 美女高潮喷水抽搐中文字幕| 国产淫片久久久久久久久| 国产在线精品亚洲第一网站| 亚洲中文字幕日韩| 国产精品人妻久久久影院| 久久天躁狠狠躁夜夜2o2o| 蜜桃久久精品国产亚洲av| 国产伦一二天堂av在线观看| 国产一区二区三区av在线 | 久久精品国产亚洲av涩爱 | 欧美日韩亚洲国产一区二区在线观看| 国产精品一区二区三区四区免费观看 | 国产探花在线观看一区二区| 国产熟女欧美一区二区| 久久久久久久久久黄片| 1024手机看黄色片| 国产精品一区二区免费欧美| 22中文网久久字幕| or卡值多少钱| 久久精品国产99精品国产亚洲性色| 亚洲最大成人中文| 男女下面进入的视频免费午夜| 亚洲精品久久国产高清桃花| 国产成人a区在线观看| 久久久久久九九精品二区国产| 国产精品精品国产色婷婷| 中文字幕免费在线视频6| 99视频精品全部免费 在线| 九九热线精品视视频播放| 亚洲图色成人| 成人av一区二区三区在线看| 日本爱情动作片www.在线观看 | 欧美日韩亚洲国产一区二区在线观看| 精品一区二区三区视频在线观看免费| 久久久久久久久大av| 最近在线观看免费完整版| 日韩一区二区视频免费看| av中文乱码字幕在线| 少妇裸体淫交视频免费看高清| 蜜桃亚洲精品一区二区三区| 亚洲av二区三区四区| 国产精品一区二区三区四区久久| 美女 人体艺术 gogo| 两性午夜刺激爽爽歪歪视频在线观看| 日本 av在线| 高清毛片免费观看视频网站| 看片在线看免费视频| 少妇人妻一区二区三区视频| 亚洲av第一区精品v没综合| 欧美国产日韩亚洲一区| 永久网站在线| 俄罗斯特黄特色一大片| 欧美高清性xxxxhd video| 网址你懂的国产日韩在线| 波多野结衣高清作品| 亚洲成人久久爱视频| 亚洲,欧美,日韩| 色综合亚洲欧美另类图片| 久9热在线精品视频| 99热精品在线国产| 日韩人妻高清精品专区| 嫩草影院新地址| 级片在线观看| 一区二区三区高清视频在线| 综合色av麻豆| 婷婷丁香在线五月| 久久香蕉精品热| 九色成人免费人妻av| 久久久国产成人精品二区| 免费av不卡在线播放| 国国产精品蜜臀av免费| 很黄的视频免费| h日本视频在线播放| 久久人人爽人人爽人人片va| 国产视频内射| 九九在线视频观看精品| 舔av片在线| 色尼玛亚洲综合影院| 亚洲av中文av极速乱 | 中文字幕av在线有码专区| 精品久久久久久久久亚洲 | 国产精品1区2区在线观看.| 亚洲综合色惰| 别揉我奶头~嗯~啊~动态视频| 中文资源天堂在线| 国产精品av视频在线免费观看| 亚洲欧美日韩高清专用| 搡老熟女国产l中国老女人| 他把我摸到了高潮在线观看| 桃红色精品国产亚洲av| 人人妻人人看人人澡| 久久久成人免费电影| 精品一区二区三区视频在线| 成人精品一区二区免费| 国产精品av视频在线免费观看| 熟女人妻精品中文字幕| 久久久久国产精品人妻aⅴ院| 精品久久国产蜜桃| 亚洲欧美精品综合久久99| 国产精品一区二区性色av| 搡老妇女老女人老熟妇| 老师上课跳d突然被开到最大视频| 欧美日本亚洲视频在线播放| 男插女下体视频免费在线播放| 搡老熟女国产l中国老女人| 亚洲综合色惰| 午夜激情欧美在线| 免费av不卡在线播放| 久久人人爽人人爽人人片va| 亚洲精品456在线播放app | 免费人成视频x8x8入口观看| av专区在线播放| 深夜a级毛片| 在现免费观看毛片| 欧美绝顶高潮抽搐喷水| 日韩欧美三级三区| 在线免费观看不下载黄p国产 | 中国美女看黄片| 欧美高清成人免费视频www| 国产伦精品一区二区三区视频9| 可以在线观看的亚洲视频| 国产免费男女视频| 欧美精品国产亚洲| 老司机福利观看| 12—13女人毛片做爰片一| 丝袜美腿在线中文| 窝窝影院91人妻| 日本色播在线视频| 久久人人爽人人爽人人片va| 很黄的视频免费| 日本-黄色视频高清免费观看| 亚洲中文日韩欧美视频| 久久久久国产精品人妻aⅴ院| 长腿黑丝高跟| 校园春色视频在线观看| 亚洲一区二区三区色噜噜| 精品免费久久久久久久清纯| 黄色欧美视频在线观看| 亚洲精品成人久久久久久| 亚洲人与动物交配视频| 丰满乱子伦码专区| 午夜福利成人在线免费观看| 三级男女做爰猛烈吃奶摸视频| 丰满乱子伦码专区| 中文字幕人妻熟人妻熟丝袜美| 国产一区二区三区av在线 | 中文亚洲av片在线观看爽| 日本a在线网址| 国产精品电影一区二区三区| 精品久久久久久久久亚洲 | 国产高清视频在线播放一区| 夜夜爽天天搞| 欧美日韩国产亚洲二区| 国产av不卡久久| 国产精品永久免费网站| 自拍偷自拍亚洲精品老妇| 美女cb高潮喷水在线观看| h日本视频在线播放| 听说在线观看完整版免费高清| av黄色大香蕉| 18禁黄网站禁片免费观看直播| 亚洲人与动物交配视频| 99视频精品全部免费 在线| 最近中文字幕高清免费大全6 | 日韩欧美 国产精品| 69av精品久久久久久| 亚洲午夜理论影院| 精品午夜福利在线看| 成人亚洲精品av一区二区| 最好的美女福利视频网| 国产欧美日韩一区二区精品| 免费av不卡在线播放| 日本免费a在线| 精品久久久噜噜| 亚洲av.av天堂| 特大巨黑吊av在线直播| 欧美一区二区精品小视频在线| 18禁裸乳无遮挡免费网站照片| 又黄又爽又免费观看的视频| 黄色视频,在线免费观看| 欧美一级a爱片免费观看看| 深夜a级毛片| 亚洲人与动物交配视频| 亚洲第一区二区三区不卡| 亚洲av一区综合| a级一级毛片免费在线观看| 久久精品国产清高在天天线| 啦啦啦啦在线视频资源| 国产爱豆传媒在线观看| 成年人黄色毛片网站| 免费av不卡在线播放| 在线观看美女被高潮喷水网站| 亚洲三级黄色毛片| 在线观看一区二区三区| 国产成人av教育| 99久久久亚洲精品蜜臀av| 午夜精品在线福利| 免费av不卡在线播放| 岛国在线免费视频观看| 亚洲,欧美,日韩| 成人特级av手机在线观看| 亚洲在线自拍视频| 黄色女人牲交| 成熟少妇高潮喷水视频| 亚洲精品日韩av片在线观看| 99热这里只有是精品在线观看| 老熟妇仑乱视频hdxx| 午夜福利视频1000在线观看| 伦理电影大哥的女人| 精品免费久久久久久久清纯| 久久久久免费精品人妻一区二区| 女生性感内裤真人,穿戴方法视频| 国产 一区 欧美 日韩| 成人国产麻豆网| 99热这里只有是精品在线观看| 老熟妇仑乱视频hdxx| 在线免费十八禁| 免费在线观看影片大全网站| 久久久久九九精品影院| h日本视频在线播放| 精品免费久久久久久久清纯| 日韩国内少妇激情av| 99在线人妻在线中文字幕| 99热这里只有精品一区| 国产午夜精品久久久久久一区二区三区 | 最近在线观看免费完整版| 亚洲成人免费电影在线观看| 一个人观看的视频www高清免费观看| 久久人妻av系列| 国产男靠女视频免费网站| 精品人妻熟女av久视频| 看十八女毛片水多多多| 国产一区二区在线av高清观看| 尤物成人国产欧美一区二区三区| 国产精品av视频在线免费观看| 欧美另类亚洲清纯唯美| 少妇被粗大猛烈的视频| 成人亚洲精品av一区二区| 3wmmmm亚洲av在线观看| 国产人妻一区二区三区在| 身体一侧抽搐| 桃红色精品国产亚洲av| 欧美最新免费一区二区三区| 99国产精品一区二区蜜桃av| 日韩,欧美,国产一区二区三区 | 亚洲人成网站在线播| av中文乱码字幕在线| 非洲黑人性xxxx精品又粗又长| 亚洲专区中文字幕在线| 国产成人av教育| 此物有八面人人有两片| 成人av在线播放网站| 一区二区三区高清视频在线| 简卡轻食公司| 身体一侧抽搐| 乱系列少妇在线播放| 夜夜看夜夜爽夜夜摸| 免费观看的影片在线观看| 国产午夜福利久久久久久| 免费观看的影片在线观看| av.在线天堂| 美女cb高潮喷水在线观看| 亚洲中文字幕一区二区三区有码在线看| 俺也久久电影网| 又紧又爽又黄一区二区| 俺也久久电影网| 全区人妻精品视频| 男人和女人高潮做爰伦理| 91久久精品电影网| 久99久视频精品免费| 国产高潮美女av| 一级黄片播放器| 午夜福利18| 免费看av在线观看网站| 国产一区二区三区av在线 | 久久99热这里只有精品18| 级片在线观看| 亚洲精品色激情综合| 联通29元200g的流量卡| 九九热线精品视视频播放| 国产亚洲精品av在线| 最近视频中文字幕2019在线8| 久久久久久久久中文| 我要看日韩黄色一级片| 国产午夜精品论理片| 国产黄色小视频在线观看| 国产精品一区二区三区四区免费观看 | 成人国产一区最新在线观看| 简卡轻食公司| 欧美又色又爽又黄视频| 久久久国产成人精品二区| avwww免费| 99久久成人亚洲精品观看| 久久久久久伊人网av| 十八禁国产超污无遮挡网站| 深夜精品福利| 国产探花在线观看一区二区| 成人特级av手机在线观看| 尾随美女入室| 亚洲av熟女| 麻豆成人av在线观看| 真人一进一出gif抽搐免费| 男女啪啪激烈高潮av片| 99热这里只有是精品在线观看| 国产午夜精品久久久久久一区二区三区 | 日本免费一区二区三区高清不卡| 国产成人av教育| 亚洲乱码一区二区免费版| 极品教师在线视频| 九九久久精品国产亚洲av麻豆| 99热这里只有是精品50|