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

    一種補(bǔ)償t-OPT噪聲的冗余捷聯(lián)慣組故障檢測方法*

    2023-07-05 13:56:26王永剛毛博年
    航天控制 2023年3期
    關(guān)鍵詞:慣組奇偶捷聯(lián)

    王永剛,毛博年,高 東

    1. 中國科學(xué)院國家空間科學(xué)中心復(fù)雜航天系統(tǒng)綜合電子與信息技術(shù)重點(diǎn)實(shí)驗(yàn)室,北京 100190 2. 中國科學(xué)院大學(xué)計(jì)算機(jī)科學(xué)與技術(shù)學(xué)院,北京 100049

    0 引言

    在對(duì)系統(tǒng)的可靠性和安全性有極高要求的航空航天領(lǐng)域,如果捷聯(lián)慣組出現(xiàn)故障,可能會(huì)直接導(dǎo)致任務(wù)的失敗,所以捷聯(lián)慣組的可靠性對(duì)于提高載體的生存能力具有重要意義。與提高慣性傳感器的元器件質(zhì)量和制造工藝等方式相比,冗余設(shè)計(jì)通過增加器件的數(shù)量大幅提升了系統(tǒng)的可靠性[1],同時(shí)為故障檢測方法提供了硬件基礎(chǔ)。

    目前對(duì)于冗余捷聯(lián)慣組的故障檢測,常采用基于等價(jià)空間原理的故障檢測方法[2],等價(jià)空間法對(duì)冗余測量信息解耦,得到與系統(tǒng)狀態(tài)無關(guān)的奇偶?xì)埐?通過檢測殘差量的變化判斷系統(tǒng)故障信息,可以得到很好的檢測效果。許多文獻(xiàn)對(duì)等價(jià)空間法進(jìn)行了相關(guān)的研究和擴(kuò)展,文獻(xiàn)[3-5]解決了等價(jià)空間法對(duì)于四軸冗余慣性傳感器無法實(shí)現(xiàn)故障隔離的問題;文獻(xiàn)[6-7]提出了檢測冗余捷聯(lián)慣組多故障的方法;文獻(xiàn)[8-9]中解耦矩陣的構(gòu)造方式可以檢測出系統(tǒng)冗余捷聯(lián)慣組特定軸的故障。

    根據(jù)具體實(shí)現(xiàn)方式的不同,可將等價(jià)空間法分為廣義似然比法[10]、最優(yōu)奇偶向量法(Optimal Parity Vector Technology,OPT)[11-12]和奇異值分解法[13],其中,OPT方法提高了系統(tǒng)解耦的魯棒性,但在構(gòu)造故障檢測函數(shù)時(shí),未考慮奇偶?xì)埐罱y(tǒng)計(jì)特性的未知性,降低了故障檢測方法的適用性。

    冗余捷聯(lián)慣組載體在機(jī)動(dòng)時(shí),慣性傳感器的刻度系數(shù)誤差、安裝誤差等誤差項(xiàng)會(huì)使奇偶?xì)埐钇?導(dǎo)致故障檢測函數(shù)在系統(tǒng)未發(fā)生故障時(shí)就超過故障閾值,因此需要提前采取措施將測量誤差標(biāo)定[14]或者使用卡爾曼濾波[15-16]等方法將誤差項(xiàng)補(bǔ)償?shù)?。文獻(xiàn)[17]針對(duì)濾波方法會(huì)抑制奇偶?xì)埐钪泄收享?xiàng)的問題,提出使用交互式多模型方法補(bǔ)償奇偶?xì)埐?但這些補(bǔ)償方法均沒有補(bǔ)償隨機(jī)噪聲,導(dǎo)致檢測方法的虛警率較高,而且無法有效檢測低故障幅值的常值漂移故障。

    本文依據(jù)基于等價(jià)空間原理的OPT方法,提出了一種補(bǔ)償故障檢測函數(shù)噪聲的t檢驗(yàn)最優(yōu)奇偶向量法(t-OPT),該方法的特點(diǎn)在于:使用奇偶?xì)埐顦颖镜慕y(tǒng)計(jì)量構(gòu)造了新的故障檢測函數(shù),使用卡爾曼濾波算法補(bǔ)償了故障檢測函數(shù)中的隨機(jī)噪聲,降低了隨機(jī)噪聲對(duì)故障判決的干擾,提高了故障檢測方法在冗余捷聯(lián)慣組故障檢測中的適用性和魯棒性,最后通過仿真實(shí)驗(yàn)證實(shí)了該方法的有效性。

    1 等價(jià)空間原理

    等價(jià)空間法的基本原理是根據(jù)冗余捷聯(lián)慣組的測量方程構(gòu)造解耦矩陣,利用解耦矩陣將系統(tǒng)狀態(tài)與系統(tǒng)故障分離,然后根據(jù)解耦后的等價(jià)方程進(jìn)行系統(tǒng)的故障檢測。

    1.1 冗余捷聯(lián)慣組測量模型

    捷聯(lián)慣組的冗余方式包括系統(tǒng)級(jí)冗余(如3套捷聯(lián)慣組同軸共基座安裝)和器件級(jí)冗余(如5個(gè)陀螺沿圓錐面均勻分布安裝,6個(gè)陀螺沿正十二面體的6個(gè)法向安裝),以三套捷聯(lián)慣組同軸共基座安裝方式為例,冗余捷聯(lián)慣組的傳感器安裝矩陣H為[9]:

    (1)

    假設(shè)冗余捷聯(lián)慣組有n(n>3)個(gè)陀螺或者加速度計(jì),傳感器測量誤差只包含隨機(jī)噪聲,可得冗余捷聯(lián)慣組的測量方程為[4]:

    Z=HX+bf+ε

    (2)

    式中:Z∈Rn為n個(gè)慣性傳感器的輸出值,n為慣性傳感器的個(gè)數(shù),H∈Rn×3為冗余捷聯(lián)慣組的傳感器安裝矩陣,X∈R3為三維狀態(tài)信息(加速度或者角速度信息),bf∈Rn為故障向量,其對(duì)應(yīng)于有故障慣性傳感器的元不為0,其他元均為0,ε∈Rn是均值為0,方差為σ2In的隨機(jī)噪聲。

    1.2 系統(tǒng)解耦

    定義矩陣V∈R(n-3)×n為式(2)所示測量方程的解耦矩陣,解耦矩陣的約束條件為:

    (3)

    式中:H∈Rn×3為冗余捷聯(lián)慣組的傳感器安裝矩陣,n為慣性傳感器的個(gè)數(shù)。

    對(duì)于解耦矩陣V的構(gòu)造,文獻(xiàn)[9]提出以傳感器安裝矩陣H的正交投影陣W=I-H(HTH)-1HT為基礎(chǔ),選取正交投影陣行向量組中的一個(gè)極大無關(guān)組,并進(jìn)行施密特正交化,從而得到冗余測量方程的解耦矩陣V。

    利用解耦矩陣V,根據(jù)式(2)可得到如下所示的等價(jià)方程:

    P=VZ=Vbf+Vε

    (4)

    式中:P∈Rn-3稱為等價(jià)向量。

    由式(4)可見,等價(jià)向量僅與隨機(jī)噪聲和可能出現(xiàn)的故障相關(guān),如果慣性傳感器發(fā)生故障,等價(jià)向量的均值會(huì)發(fā)生變化,由等價(jià)向量在無故障和有故障情況下的不一致性,為故障檢測提供了參考。但僅根據(jù)等價(jià)向量無法實(shí)現(xiàn)對(duì)故障器件的定位,需要采取措施對(duì)等價(jià)空間法進(jìn)行具體的實(shí)現(xiàn),所以本文在OPT方法的基礎(chǔ)上提出了t-OPT方法。

    2 t-OPT方法

    等價(jià)空間法中的解耦矩陣V由冗余捷聯(lián)慣組的傳感器安裝矩陣H確定,為了降低傳感器安裝矩陣編排方式對(duì)故障檢測性能的影響,OPT方法構(gòu)造了對(duì)特定慣性傳感器敏感的最優(yōu)奇偶向量(又稱最優(yōu)解耦向量),提高了系統(tǒng)解耦的魯棒性[11]。但OPT方法在構(gòu)造故障檢測函數(shù)時(shí)引入了未知參數(shù),為提高OPT方法在冗余捷聯(lián)慣組故障檢測中的適用性,本文提出了t-OPT方法,該方法在OPT方法的基礎(chǔ)上構(gòu)造了新的故障檢測函數(shù)。

    2.1 OPT方法

    OPT方法針對(duì)冗余系統(tǒng)中的每個(gè)慣性傳感器,定義了與各個(gè)慣性傳感器對(duì)應(yīng)的最優(yōu)解耦向量vi,考慮解耦約束條件viTH=0,最優(yōu)解耦向量可以表示為式(3)中解耦矩陣V各行的線性組合,然后通過最大化性能指標(biāo)函數(shù)Si,使得最優(yōu)解耦向量vi對(duì)第i個(gè)慣性傳感器敏感,而對(duì)其他慣性傳感器和隨機(jī)噪聲不敏感,最優(yōu)解耦向量的約束條件可總結(jié)為[12]:

    (5)

    式中:vi∈Rn為第i個(gè)慣性傳感器的最優(yōu)解耦向量,n為慣性傳感器的個(gè)數(shù),H∈Rn×3為冗余捷聯(lián)慣組的傳感器安裝矩陣,V∈R(n-3)×n為解耦矩陣,ci為解耦矩陣V線性組合的坐標(biāo)向量,ei為n階單位矩陣In的第i個(gè)列向量,viTei和viTej分別表示第i個(gè)慣性傳感器的最優(yōu)解耦向量對(duì)第i個(gè)和第j個(gè)慣性傳感器故障的敏感度,‖viT‖表示對(duì)隨機(jī)噪聲的敏感度。

    根據(jù)式(5)所示的約束條件可求得各個(gè)慣性傳感器的最優(yōu)解耦向量:

    vi=VT(V(2In-eieiT)VT)-1Vei

    (6)

    利用最優(yōu)解耦向量vi,根據(jù)式(2)可得到各個(gè)慣性傳感器的奇偶?xì)埐?

    ri=viTZ=viTbf+viTε

    (7)

    式中:ri為第i個(gè)慣性傳感器的奇偶?xì)埐?bf∈Rn為故障向量,ε∈Rn是均值為0,方差為σ2In的隨機(jī)噪聲。

    OPT方法為構(gòu)造故障檢測函數(shù),以正態(tài)檢驗(yàn)法為基礎(chǔ),假設(shè)慣性傳感器無故障,由式(7)可得奇偶?xì)埐顁i~N(0,σ2‖viT‖2),對(duì)奇偶?xì)埐钸M(jìn)行標(biāo)準(zhǔn)化處理,得到各個(gè)慣性傳感器的故障檢測函數(shù):

    (8)

    對(duì)給定的顯著性水平α,由標(biāo)準(zhǔn)正態(tài)分布的分位點(diǎn)得到OPT方法的故障閾值:

    Td=z1-α

    (9)

    根據(jù)式(8)~(9)和正態(tài)檢驗(yàn)法,OPT方法的故障判決策略為:比較n個(gè)慣性傳感器的故障檢測函數(shù),如果Dm為Di中的最大值,當(dāng)Dm>Td時(shí),則判定第m個(gè)慣性傳感器發(fā)生故障,否則認(rèn)為所有慣性傳感器均無故障。

    2.2 OPT方法的改進(jìn)

    由式(8)可見,OPT方法的故障檢測函數(shù)中含有參數(shù)σ,在實(shí)際情況中,隨機(jī)噪聲的標(biāo)準(zhǔn)差σ是未知的,所以O(shè)PT方法不能依據(jù)正態(tài)檢驗(yàn)法構(gòu)造故障檢測函數(shù)。為了提高OPT方法的適用性,本文提出了t-OPT方法,該方法依據(jù)t檢驗(yàn)使用奇偶?xì)埐顦颖镜慕y(tǒng)計(jì)量構(gòu)造了新的故障檢測函數(shù),同時(shí)為提高樣本的代表性,增加了奇偶?xì)埐畹臉颖救萘俊?/p>

    根據(jù)式(7)中奇偶?xì)埐钤跓o故障和有故障情況下統(tǒng)計(jì)特性的差異,針對(duì)第i(i=1,2,…,n)個(gè)慣性傳感器,提出如下無故障假設(shè)H0與有故障假設(shè)H1:

    H0:ri~N(0,σ2‖viT‖2)

    (10)

    H1:ri~N(viTbf,σ2‖viT‖2)

    (11)

    (12)

    (13)

    由式(12)和(13)可見,樣本均值和樣本方差的概率密度函數(shù)中均含有未知參數(shù)σ,為消除未知參數(shù),根據(jù)t分布的定義[18],使用樣本均值和樣本方差構(gòu)造如下統(tǒng)計(jì)量:

    (14)

    化簡后得到服從t分布的檢驗(yàn)統(tǒng)計(jì)量:

    (15)

    根據(jù)式(15),構(gòu)造各個(gè)慣性傳感器的故障檢測函數(shù):

    (16)

    對(duì)給定的顯著性水平α,根據(jù)t分布的分位點(diǎn)得到t-OPT方法的故障閾值:

    Td*=t1-α(l-1)

    (17)

    式中:l為樣本個(gè)數(shù)。

    由式(16)可見,故障檢測函數(shù)中的所有參數(shù)均可由奇偶?xì)埐顦颖镜慕y(tǒng)計(jì)量求得,所以t-OPT方法相較于OPT方法具有較好的適用性。

    3 故障檢測函數(shù)噪聲的補(bǔ)償

    由式(7)可見,各個(gè)慣性傳感器的奇偶?xì)埐顑H與傳感器的故障和噪聲相關(guān),通過降低噪聲的方差,可以使檢測方法具有較低的虛警率和較好的小故障檢測能力。所以為進(jìn)一步提高t-OPT方法的適用性和魯棒性,本文采用卡爾曼濾波算法補(bǔ)償了t-OPT方法中故障檢測函數(shù)的隨機(jī)噪聲。

    以冗余捷聯(lián)慣組中的一個(gè)慣性傳感器為例,由式(16)得到該慣性傳感器的故障檢測函數(shù)g:

    (18)

    將故障檢測函數(shù)建模為隨機(jī)常數(shù)f作為離散狀態(tài)方程,將式(18)所示的方程作為離散測量方程,并且定義測量方程中的噪聲為ε′,得到慣性傳感器故障檢測函數(shù)計(jì)算模型:

    fk=fk-1+βk-1

    (19)

    (20)

    (21)

    Jk/k-1=Jk-1+Qk-1

    (22)

    Kk=Jk/k-1/(Jk/k-1+Rk)

    (23)

    (24)

    Jk=(1-Kk)Jk/k-1

    (25)

    因此,補(bǔ)償故障檢測函數(shù)的噪聲后,t-OPT方法的故障檢測流程可總結(jié)為:

    1)根據(jù)式(3)計(jì)算冗余捷聯(lián)慣組測量方程的解耦矩陣V;

    2)根據(jù)式(6)計(jì)算每個(gè)慣性傳感器的最優(yōu)解耦向量vi;

    3)根據(jù)式(7),由最優(yōu)解耦向量vi得到各個(gè)慣性傳感器的奇偶?xì)埐?

    5)對(duì)于給定的顯著性水平α,根據(jù)式(17)得到故障閾值Td*;

    4 仿真校驗(yàn)

    為了驗(yàn)證方法的有效性,建立仿真平臺(tái),設(shè)定故障檢測仿真條件如下:

    1)三套捷聯(lián)慣組為同軸共基座冗余配置方式,陀螺序號(hào)與式(1)中傳感器安裝矩陣的行號(hào)一致;

    2)冗余捷聯(lián)慣組中陀螺的噪聲標(biāo)準(zhǔn)差σ=0.1(°)/h;

    3)設(shè)定顯著性水平α=0.95,選取樣本個(gè)數(shù)l=16,則OPT方法的故障閾值為Td=z1-0.95=1.645,t-OPT方法的故障閾值為Td*=t1-0.95(16-1)=1.7531;

    4)在慣性傳感器的刻度系數(shù)誤差、安裝誤差等誤差項(xiàng)已被標(biāo)定或補(bǔ)償?shù)臈l件下,設(shè)定載體運(yùn)動(dòng)場景1為載體勻速直線運(yùn)動(dòng),運(yùn)動(dòng)場景2為載體機(jī)動(dòng)運(yùn)行(三軸角運(yùn)動(dòng)設(shè)為X:50sin(8πt)(°)/s,Y:50cos(πt)(°)/s,Z:50t(°)/s),分別在兩種運(yùn)動(dòng)場景下選取冗余捷聯(lián)慣組載體運(yùn)動(dòng)時(shí)長為10 s的陀螺輸出數(shù)據(jù)作為方法性能驗(yàn)證數(shù)據(jù),采樣頻率為100 Hz;

    5)故障建模:常值漂移和線性漂移是慣性傳感器常見且不易察覺的故障,在仿真中,將常值漂移建模為幅值可調(diào)的階躍函數(shù),將線性漂移建模為斜率可調(diào)的斜坡函數(shù);

    6)虛警率:冗余系統(tǒng)無故障時(shí),所有陀螺故障檢測函數(shù)中的最大值超過故障閾值的概率;漏檢率:冗余系統(tǒng)有故障時(shí),故障陀螺的故障檢測函數(shù)不是最大值或者是最大值但低于故障閾值的概率。

    4.1 常值漂移故障檢測仿真

    當(dāng)載體分別處于運(yùn)動(dòng)場景1和2時(shí),以設(shè)定的仿真條件為基礎(chǔ),在第4秒對(duì)冗余捷聯(lián)慣組的陀螺1注入故障幅值為5σ的常值漂移故障,使用補(bǔ)償噪聲后的t-OPT方法對(duì)冗余捷聯(lián)慣組進(jìn)行故障檢測,得到兩種運(yùn)動(dòng)場景下所有陀螺的常值漂移故障檢測曲線如圖1~2所示。

    圖1 運(yùn)動(dòng)場景1下所有陀螺的常值漂移故障檢測曲線

    圖2 運(yùn)動(dòng)場景2下所有陀螺的常值漂移故障檢測曲線

    為比較不同故障檢測方法的故障檢測能力,當(dāng)載體處于運(yùn)動(dòng)場景1時(shí),在第4秒對(duì)冗余捷聯(lián)慣組的陀螺1注入故障幅值為5σ的常值漂移故障,分別使用OPT方法、t-OPT方法和補(bǔ)償噪聲后的t-OPT方法對(duì)陀螺1進(jìn)行故障檢測,得到陀螺1的常值漂移故障檢測曲線如圖3所示。

    圖3 運(yùn)動(dòng)場景1下陀螺1的常值漂移故障檢測曲線

    為分析對(duì)比故障檢測方法的小故障檢測能力,當(dāng)載體處于運(yùn)動(dòng)場景1時(shí),在第4秒對(duì)陀螺1注入不同故障幅值的常值漂移故障,進(jìn)行1000次隨機(jī)實(shí)驗(yàn)并記錄仿真結(jié)果的虛警率和漏檢率,實(shí)驗(yàn)結(jié)果如表1和2所示。

    表1 運(yùn)動(dòng)場景1下常值漂移故障檢測虛警率

    表2 運(yùn)動(dòng)場景1下常值漂移故障檢測漏檢率

    分析常值漂移故障檢測仿真結(jié)果:

    1)如圖1~2所示,注入故障后,所有陀螺的故障檢測函數(shù)中,陀螺1的故障檢測函數(shù)值最大且超出故障閾值,與t-OPT方法的故障判決策略一致,而且故障檢測結(jié)果與載體運(yùn)動(dòng)場景無關(guān)。

    2)如圖3所示,注入故障后,3種方法的故障檢測函數(shù)在第4秒后均明顯變大且超出故障閾值。但OPT方法使用了仿真條件中的噪聲標(biāo)準(zhǔn)差σ,在實(shí)際情況中噪聲標(biāo)準(zhǔn)差是未知的。而t-OPT方法使用奇偶?xì)埐顦颖镜慕y(tǒng)計(jì)量就可以判斷出故障陀螺,提高了方法的適用性。

    3)如表1所示,OPT方法和t-OPT方法的故障檢測虛警率均較高,導(dǎo)致陀螺在無故障時(shí)經(jīng)常被判定為出現(xiàn)故障。使用卡爾曼濾波算法補(bǔ)償t-OPT方法中故障檢測函數(shù)的噪聲后,故障檢測虛警率降至0.55%,進(jìn)一步提高了方法的適用性。

    4)如表2所示,當(dāng)故障幅值低于7σ時(shí),OPT方法的故障檢測漏檢率較高,而t-OPT方法的故障檢測漏檢率較低,補(bǔ)償噪聲后,t-OPT方法的漏檢率得到了進(jìn)一步的降低。

    5)結(jié)合表1和2,當(dāng)陀螺故障幅值較低時(shí),通過補(bǔ)償t-OPT方法中故障檢測函數(shù)的噪聲,可以明顯降低故障檢測的虛警率和漏檢率,使檢測方法可以有效檢測到低故障幅值的常值漂移故障,說明補(bǔ)償噪聲可以提高t-OPT方法的魯棒性。

    4.2 線性漂移故障檢測仿真

    當(dāng)載體處于運(yùn)動(dòng)場景1時(shí),在第4秒對(duì)陀螺1注入斜率為0.001的線性漂移故障,分別使用OPT方法、t-OPT方法和補(bǔ)償噪聲后的t-OPT方法對(duì)冗余捷聯(lián)慣組進(jìn)行故障檢測,得到陀螺1的線性漂移故障檢測曲線如圖4所示。

    圖4 運(yùn)動(dòng)場景1下陀螺1的線性漂移故障檢測曲線

    如圖4的線性漂移故障檢測曲線所示, OPT方法的故障檢測時(shí)延約為4 s,而t-OPT方法的故障檢測時(shí)延僅為0.2 s左右,但虛警現(xiàn)象仍較嚴(yán)重,通過補(bǔ)償t-OPT方法中故障檢測函數(shù)的噪聲,有效降低了故障檢測的虛警率。

    由于線性漂移故障的漸變特性,陀螺出現(xiàn)故障初期,故障幅值極低,所以線性漂移的故障檢測有一定的延時(shí)性。由表1~2可知,補(bǔ)償噪聲后的t-OPT方法可以有效檢測到低故障幅值的常值漂移故障,所以在線性漂移故障檢測中表現(xiàn)為故障檢測時(shí)延較低。

    5 結(jié)論

    對(duì)于冗余捷聯(lián)慣組的故障檢測,在OPT方法的基礎(chǔ)上,提出了一種補(bǔ)償故障檢測函數(shù)噪聲的t-OPT方法。該方法使用奇偶?xì)埐顦颖镜慕y(tǒng)計(jì)量構(gòu)造了新的故障檢測函數(shù),補(bǔ)償了t-OPT方法中故障檢測函數(shù)的隨機(jī)噪聲,提高了故障檢測方法的適用性和魯棒性。仿真結(jié)果表明,該方法相較于OPT方法具有較低的虛警率,可以準(zhǔn)確檢測到低故障幅值的常值漂移,有效降低了線性漂移故障的檢測時(shí)延。

    猜你喜歡
    慣組奇偶捷聯(lián)
    三招求解“奇偶項(xiàng)交織”遞推數(shù)列問題
    談?wù)勂媾己瘮?shù)的應(yīng)用
    n分奇偶時(shí),如何求數(shù)列的通項(xiàng)
    活用奇偶函數(shù)的性質(zhì)妙解題
    基于LabVIEW的光纖慣組自動(dòng)測試軟件的設(shè)計(jì)與實(shí)現(xiàn)
    彈道導(dǎo)彈的捷聯(lián)慣性/天文組合導(dǎo)航方法
    捷聯(lián)慣性/天文/雷達(dá)高度表組合導(dǎo)航
    激光捷聯(lián)慣組減振系統(tǒng)設(shè)計(jì)計(jì)算及試驗(yàn)研究
    半捷聯(lián)雷達(dá)導(dǎo)引頭視線角速度提取
    帶斜置元件的光纖陀螺捷聯(lián)慣組標(biāo)定方法
    亚洲无线观看免费| 在线观看免费视频日本深夜| 久久久久久久亚洲中文字幕| av女优亚洲男人天堂| 美女xxoo啪啪120秒动态图| 成人三级黄色视频| 成人特级黄色片久久久久久久| 国产精品一区二区免费欧美| 91麻豆av在线| 99riav亚洲国产免费| 天堂√8在线中文| 深夜a级毛片| 国产爱豆传媒在线观看| 九九在线视频观看精品| 少妇的逼水好多| 亚洲天堂国产精品一区在线| 免费看av在线观看网站| 成人性生交大片免费视频hd| 91在线精品国自产拍蜜月| 日韩欧美一区二区三区在线观看| 久久久午夜欧美精品| 亚洲人成网站在线播放欧美日韩| 搞女人的毛片| 亚洲美女搞黄在线观看 | 亚洲成人中文字幕在线播放| 国产精品久久久久久久电影| 欧美绝顶高潮抽搐喷水| 老女人水多毛片| 久久草成人影院| 桃红色精品国产亚洲av| 有码 亚洲区| 亚洲美女视频黄频| 嫩草影视91久久| 人妻制服诱惑在线中文字幕| 亚洲在线观看片| 91麻豆av在线| 男人和女人高潮做爰伦理| 欧美在线一区亚洲| 琪琪午夜伦伦电影理论片6080| 日韩欧美一区二区三区在线观看| 我的老师免费观看完整版| 国产伦在线观看视频一区| 亚洲av成人精品一区久久| 成人av在线播放网站| 国产精品99久久久久久久久| 国产高潮美女av| avwww免费| 午夜福利18| 欧美又色又爽又黄视频| 国产高潮美女av| h日本视频在线播放| 看片在线看免费视频| 日韩欧美国产在线观看| 老司机午夜福利在线观看视频| 午夜老司机福利剧场| 亚洲七黄色美女视频| 亚洲内射少妇av| 婷婷六月久久综合丁香| 成人国产一区最新在线观看| 在线观看66精品国产| 欧美日韩国产亚洲二区| 国产探花在线观看一区二区| 免费黄网站久久成人精品| 久久精品国产亚洲网站| 亚洲av免费高清在线观看| 午夜福利在线观看吧| 一夜夜www| 国产国拍精品亚洲av在线观看| 国产高清激情床上av| 欧美日本亚洲视频在线播放| 国产一级毛片七仙女欲春2| 女人被狂操c到高潮| 乱系列少妇在线播放| 性色avwww在线观看| 搡老妇女老女人老熟妇| 狠狠狠狠99中文字幕| 国产欧美日韩精品亚洲av| 午夜a级毛片| 午夜影院日韩av| 久久热精品热| 久久香蕉精品热| 日韩精品中文字幕看吧| 午夜免费男女啪啪视频观看 | 成年免费大片在线观看| 亚洲av免费高清在线观看| 亚洲欧美日韩卡通动漫| 午夜亚洲福利在线播放| 欧美xxxx黑人xx丫x性爽| 听说在线观看完整版免费高清| 深夜精品福利| 国产精品久久久久久精品电影| 亚洲国产精品sss在线观看| 国产精品综合久久久久久久免费| 免费av不卡在线播放| 老司机深夜福利视频在线观看| 免费在线观看日本一区| 国产又黄又爽又无遮挡在线| 亚洲美女搞黄在线观看 | 欧美色视频一区免费| 久久久久国内视频| 亚洲aⅴ乱码一区二区在线播放| 夜夜爽天天搞| 亚洲黑人精品在线| 少妇猛男粗大的猛烈进出视频 | 亚洲四区av| 亚洲最大成人手机在线| av女优亚洲男人天堂| 99精品久久久久人妻精品| 国产欧美日韩一区二区精品| 久久久色成人| 国产美女午夜福利| 在线看三级毛片| 日本黄色片子视频| 精品无人区乱码1区二区| 国产69精品久久久久777片| 免费看av在线观看网站| 91麻豆av在线| 亚洲精品乱码久久久v下载方式| 在线播放国产精品三级| 国产 一区 欧美 日韩| 免费观看精品视频网站| av专区在线播放| 波多野结衣巨乳人妻| 国产亚洲精品综合一区在线观看| 亚洲,欧美,日韩| 色综合色国产| 我的老师免费观看完整版| 老熟妇仑乱视频hdxx| 人妻夜夜爽99麻豆av| 一区二区三区高清视频在线| 老熟妇仑乱视频hdxx| 99久久久亚洲精品蜜臀av| 国产乱人视频| 1000部很黄的大片| 午夜福利高清视频| 久久久久久九九精品二区国产| 在线观看av片永久免费下载| 亚洲精品亚洲一区二区| 日韩精品青青久久久久久| 久久精品91蜜桃| 身体一侧抽搐| 欧美一区二区精品小视频在线| 国产爱豆传媒在线观看| 国产乱人视频| 日日啪夜夜撸| 免费看a级黄色片| 亚洲欧美日韩无卡精品| av国产免费在线观看| 18禁黄网站禁片免费观看直播| 久久欧美精品欧美久久欧美| 不卡一级毛片| 国产精品综合久久久久久久免费| 一进一出抽搐gif免费好疼| 悠悠久久av| 欧美+亚洲+日韩+国产| 欧美日韩黄片免| 波野结衣二区三区在线| 日本成人三级电影网站| 亚洲成人免费电影在线观看| 成年版毛片免费区| h日本视频在线播放| 国产精品国产高清国产av| 白带黄色成豆腐渣| 婷婷亚洲欧美| 亚洲经典国产精华液单| 午夜福利成人在线免费观看| 高清毛片免费观看视频网站| 中亚洲国语对白在线视频| 国产精品久久视频播放| 欧美成人a在线观看| 一区二区三区免费毛片| 欧美xxxx黑人xx丫x性爽| 久久天躁狠狠躁夜夜2o2o| 色精品久久人妻99蜜桃| 亚洲五月天丁香| 美女cb高潮喷水在线观看| 欧美xxxx黑人xx丫x性爽| 久久久久久久精品吃奶| 男女视频在线观看网站免费| 精品久久久久久久末码| 天天躁日日操中文字幕| 亚洲精华国产精华精| 成年人黄色毛片网站| 又爽又黄a免费视频| 欧美日韩国产亚洲二区| 国产精品一区二区三区四区免费观看 | 热99在线观看视频| 亚洲av免费在线观看| 成年女人永久免费观看视频| 超碰av人人做人人爽久久| 天天躁日日操中文字幕| 99精品久久久久人妻精品| 亚洲一区高清亚洲精品| 亚洲美女搞黄在线观看 | 91精品国产九色| 国产综合懂色| 中文资源天堂在线| av天堂中文字幕网| 国产高清不卡午夜福利| 中文字幕av成人在线电影| 午夜影院日韩av| www.色视频.com| 久久久久久久午夜电影| 岛国在线免费视频观看| 国产亚洲精品久久久com| 日韩 亚洲 欧美在线| 国产精品久久电影中文字幕| 国内精品久久久久精免费| 欧美日本视频| 波野结衣二区三区在线| 精品人妻偷拍中文字幕| 热99在线观看视频| 久久久久久久午夜电影| 禁无遮挡网站| 美女高潮的动态| 午夜久久久久精精品| 在线国产一区二区在线| 欧美激情久久久久久爽电影| 国产一区二区三区视频了| 精品久久久久久,| 久久婷婷人人爽人人干人人爱| 国产成人福利小说| 两个人视频免费观看高清| 1024手机看黄色片| 一级av片app| 精品一区二区三区av网在线观看| 国产伦一二天堂av在线观看| 成年女人看的毛片在线观看| 精品久久久久久,| 99热只有精品国产| 日韩欧美免费精品| 男人和女人高潮做爰伦理| 老司机福利观看| 亚洲三级黄色毛片| 少妇熟女aⅴ在线视频| 熟女电影av网| 2021天堂中文幕一二区在线观| 丰满的人妻完整版| 深爱激情五月婷婷| 少妇猛男粗大的猛烈进出视频 | 久久久色成人| 麻豆av噜噜一区二区三区| 亚洲第一电影网av| 国产精品人妻久久久久久| 国产精品98久久久久久宅男小说| 久久欧美精品欧美久久欧美| 国产不卡一卡二| 亚洲熟妇熟女久久| 欧美在线一区亚洲| 日本一二三区视频观看| 男人狂女人下面高潮的视频| 精品久久久久久久久久免费视频| 成人永久免费在线观看视频| 婷婷精品国产亚洲av| 欧美+亚洲+日韩+国产| 精品一区二区免费观看| 亚洲狠狠婷婷综合久久图片| 国模一区二区三区四区视频| 成人性生交大片免费视频hd| 直男gayav资源| 日韩在线高清观看一区二区三区 | 大型黄色视频在线免费观看| 国产伦精品一区二区三区视频9| 免费高清视频大片| 日本a在线网址| 真人一进一出gif抽搐免费| 欧美丝袜亚洲另类 | 久久亚洲真实| 欧美成人a在线观看| 超碰av人人做人人爽久久| 亚洲国产精品sss在线观看| av天堂中文字幕网| 亚洲av电影不卡..在线观看| 欧美性感艳星| 露出奶头的视频| 两个人视频免费观看高清| 乱系列少妇在线播放| 高清在线国产一区| 国产激情偷乱视频一区二区| 日本成人三级电影网站| 无遮挡黄片免费观看| 成年女人永久免费观看视频| 99riav亚洲国产免费| 老司机福利观看| 亚洲国产日韩欧美精品在线观看| 免费观看在线日韩| 亚洲精品粉嫩美女一区| www.www免费av| 99热这里只有是精品在线观看| 国产精品一区二区免费欧美| 又黄又爽又免费观看的视频| 人人妻人人澡欧美一区二区| 国产精品,欧美在线| 亚洲av第一区精品v没综合| 成人性生交大片免费视频hd| 中文在线观看免费www的网站| 精品福利观看| 欧美又色又爽又黄视频| 哪里可以看免费的av片| 国产精华一区二区三区| 国产精品久久电影中文字幕| 久久久久免费精品人妻一区二区| 亚洲电影在线观看av| 国产精品av视频在线免费观看| 黄色视频,在线免费观看| 日本欧美国产在线视频| 欧美另类亚洲清纯唯美| 久久久久国产精品人妻aⅴ院| 欧美又色又爽又黄视频| 国产一区二区在线av高清观看| 久久午夜福利片| 极品教师在线视频| 韩国av一区二区三区四区| 国产高潮美女av| 男女做爰动态图高潮gif福利片| 亚洲av日韩精品久久久久久密| 嫩草影院入口| 91精品国产九色| 欧美日本视频| 日韩人妻高清精品专区| 美女高潮喷水抽搐中文字幕| 亚洲黑人精品在线| 性欧美人与动物交配| 国产伦精品一区二区三区四那| 成人特级黄色片久久久久久久| 日韩中字成人| 国产精品女同一区二区软件 | av天堂中文字幕网| 亚洲成人中文字幕在线播放| 中文字幕熟女人妻在线| 午夜精品久久久久久毛片777| 十八禁国产超污无遮挡网站| 国产黄片美女视频| 色在线成人网| 亚洲狠狠婷婷综合久久图片| 可以在线观看的亚洲视频| 国产成人a区在线观看| 一区二区三区激情视频| 波多野结衣高清作品| 午夜免费激情av| 日本在线视频免费播放| 网址你懂的国产日韩在线| 午夜a级毛片| 91久久精品国产一区二区成人| 国产v大片淫在线免费观看| 内射极品少妇av片p| av在线观看视频网站免费| 最近视频中文字幕2019在线8| 可以在线观看的亚洲视频| 99精品在免费线老司机午夜| 成人av一区二区三区在线看| 日韩欧美一区二区三区在线观看| 毛片一级片免费看久久久久 | 日韩欧美在线乱码| 久久精品人妻少妇| 久久久久久久午夜电影| 很黄的视频免费| 亚洲无线观看免费| 国产av在哪里看| 一级毛片久久久久久久久女| 久久精品国产亚洲网站| 亚洲精华国产精华液的使用体验 | 91久久精品国产一区二区成人| 精品日产1卡2卡| av天堂中文字幕网| 国产91精品成人一区二区三区| 欧美黑人巨大hd| 欧美国产日韩亚洲一区| 国语自产精品视频在线第100页| 在线免费观看的www视频| 精品欧美国产一区二区三| 人妻丰满熟妇av一区二区三区| 一区二区三区高清视频在线| 亚洲精品国产成人久久av| 在线看三级毛片| 亚洲欧美日韩东京热| 搡女人真爽免费视频火全软件 | 精品久久久久久久久亚洲 | 成人永久免费在线观看视频| 一级a爱片免费观看的视频| 欧美成人性av电影在线观看| 亚洲美女搞黄在线观看 | 国产毛片a区久久久久| 国产精品爽爽va在线观看网站| 可以在线观看的亚洲视频| 搡老熟女国产l中国老女人| 3wmmmm亚洲av在线观看| 日韩欧美精品免费久久| 国产探花在线观看一区二区| 午夜福利在线在线| 搡老妇女老女人老熟妇| 国产精品永久免费网站| 日本黄色视频三级网站网址| 日韩欧美三级三区| 国产精品一区www在线观看 | 免费高清视频大片| 美女 人体艺术 gogo| 日韩国内少妇激情av| 久久久成人免费电影| 国产精品嫩草影院av在线观看 | 欧美一区二区精品小视频在线| 我要看日韩黄色一级片| 亚洲av电影不卡..在线观看| 国产熟女欧美一区二区| 国产一区二区在线av高清观看| 久久久久性生活片| 精品不卡国产一区二区三区| АⅤ资源中文在线天堂| 蜜桃亚洲精品一区二区三区| 老熟妇仑乱视频hdxx| 成人午夜高清在线视频| 亚洲国产精品久久男人天堂| 国产熟女欧美一区二区| 国产精品日韩av在线免费观看| 亚洲精品456在线播放app | 一夜夜www| 校园春色视频在线观看| 日本黄大片高清| 国产精品自产拍在线观看55亚洲| 国产精品久久久久久久久免| 日本成人三级电影网站| 国产一区二区激情短视频| 成人二区视频| 精品不卡国产一区二区三区| 床上黄色一级片| 干丝袜人妻中文字幕| 国产高清不卡午夜福利| 少妇的逼好多水| 久久草成人影院| 国产精品久久久久久久电影| 综合色av麻豆| 亚洲国产精品久久男人天堂| 成人av一区二区三区在线看| а√天堂www在线а√下载| 十八禁国产超污无遮挡网站| 日日啪夜夜撸| 99九九线精品视频在线观看视频| 99久久成人亚洲精品观看| 国产精品综合久久久久久久免费| 制服丝袜大香蕉在线| 亚洲人成网站高清观看| 亚洲乱码一区二区免费版| 性插视频无遮挡在线免费观看| 99热精品在线国产| 婷婷色综合大香蕉| 久久九九热精品免费| 国产探花极品一区二区| 亚洲18禁久久av| 校园春色视频在线观看| 日本-黄色视频高清免费观看| 国产精品日韩av在线免费观看| 啦啦啦啦在线视频资源| 欧美日韩综合久久久久久 | 嫁个100分男人电影在线观看| 亚洲中文日韩欧美视频| 搡女人真爽免费视频火全软件 | 国产蜜桃级精品一区二区三区| 久久精品夜夜夜夜夜久久蜜豆| 69av精品久久久久久| 免费观看精品视频网站| 亚洲乱码一区二区免费版| 97热精品久久久久久| 成人精品一区二区免费| 午夜福利欧美成人| 99在线视频只有这里精品首页| 国产探花极品一区二区| 亚洲自偷自拍三级| 蜜桃亚洲精品一区二区三区| 国产极品精品免费视频能看的| 亚洲av一区综合| 色尼玛亚洲综合影院| 亚洲专区中文字幕在线| 91久久精品国产一区二区成人| 日日摸夜夜添夜夜添小说| 国产精品av视频在线免费观看| 国产中年淑女户外野战色| 又紧又爽又黄一区二区| 免费观看人在逋| 美女被艹到高潮喷水动态| 最新中文字幕久久久久| x7x7x7水蜜桃| 九九在线视频观看精品| 日韩强制内射视频| 国产aⅴ精品一区二区三区波| 亚洲自拍偷在线| 麻豆成人午夜福利视频| 亚洲av五月六月丁香网| 亚洲狠狠婷婷综合久久图片| 直男gayav资源| 亚洲精品一区av在线观看| 禁无遮挡网站| 成人国产综合亚洲| 久久久久久久午夜电影| 国产综合懂色| 国产免费一级a男人的天堂| 国国产精品蜜臀av免费| 99久久成人亚洲精品观看| 日韩中文字幕欧美一区二区| 精品久久国产蜜桃| 尤物成人国产欧美一区二区三区| 久久人人精品亚洲av| 久9热在线精品视频| 变态另类成人亚洲欧美熟女| 制服丝袜大香蕉在线| 真人做人爱边吃奶动态| 97超视频在线观看视频| 欧美高清性xxxxhd video| 国产三级在线视频| 欧美xxxx性猛交bbbb| 日韩中文字幕欧美一区二区| 亚洲国产精品成人综合色| 国产成人一区二区在线| 91久久精品国产一区二区三区| 在线播放国产精品三级| 91av网一区二区| 在线看三级毛片| 一区二区三区高清视频在线| 国产成人影院久久av| 99久久九九国产精品国产免费| 日韩强制内射视频| 欧美最黄视频在线播放免费| 美女黄网站色视频| 国产成人一区二区在线| 琪琪午夜伦伦电影理论片6080| 免费在线观看日本一区| 久久婷婷人人爽人人干人人爱| 久久精品国产99精品国产亚洲性色| 免费在线观看影片大全网站| 亚洲精品影视一区二区三区av| 亚洲国产日韩欧美精品在线观看| 长腿黑丝高跟| 乱码一卡2卡4卡精品| 久久精品久久久久久噜噜老黄 | 97超级碰碰碰精品色视频在线观看| 99视频精品全部免费 在线| 深爱激情五月婷婷| 久久久精品大字幕| 在线观看一区二区三区| 99热6这里只有精品| 女同久久另类99精品国产91| 俄罗斯特黄特色一大片| 成人美女网站在线观看视频| 在线观看一区二区三区| 99热6这里只有精品| 人妻夜夜爽99麻豆av| 免费人成视频x8x8入口观看| 看黄色毛片网站| 国产成年人精品一区二区| 欧美另类亚洲清纯唯美| 国内少妇人妻偷人精品xxx网站| 国内精品宾馆在线| 亚洲人成网站在线播放欧美日韩| 久久精品综合一区二区三区| 又爽又黄a免费视频| 国产男靠女视频免费网站| 久久久国产成人免费| 麻豆成人av在线观看| or卡值多少钱| 欧美潮喷喷水| 国产不卡一卡二| 色视频www国产| 国产爱豆传媒在线观看| 18禁裸乳无遮挡免费网站照片| 国产成人一区二区在线| 国产精品三级大全| 亚洲天堂国产精品一区在线| 五月伊人婷婷丁香| x7x7x7水蜜桃| 色哟哟哟哟哟哟| 亚洲国产精品成人综合色| 中文资源天堂在线| 少妇的逼好多水| 日本在线视频免费播放| 欧美日本亚洲视频在线播放| 日韩欧美精品免费久久| 99精品在免费线老司机午夜| 自拍偷自拍亚洲精品老妇| 免费大片18禁| 老司机午夜福利在线观看视频| 91在线观看av| 在线观看av片永久免费下载| 免费大片18禁| 男女那种视频在线观看| 久久99热这里只有精品18| 亚洲男人的天堂狠狠| 久久99热这里只有精品18| 午夜久久久久精精品| 亚州av有码| 人妻丰满熟妇av一区二区三区| 哪里可以看免费的av片| 中文字幕高清在线视频| 亚洲最大成人av| 免费看光身美女| 男女视频在线观看网站免费| 小说图片视频综合网站| 国产亚洲精品久久久com| 免费在线观看影片大全网站| 国产美女午夜福利| 亚洲精品一区av在线观看| 欧美中文日本在线观看视频| 日韩中字成人| 亚洲美女搞黄在线观看 | 国产免费av片在线观看野外av| 久久久久精品国产欧美久久久| 亚洲av免费在线观看| 啦啦啦韩国在线观看视频| 国产精品爽爽va在线观看网站| 国产午夜精品论理片| 欧美成人性av电影在线观看| 网址你懂的国产日韩在线| 亚洲欧美清纯卡通| 日韩欧美 国产精品| 日韩 亚洲 欧美在线| 欧美又色又爽又黄视频| 日本免费a在线|