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

    航天器質(zhì)量特性在線辨識方法及地面驗證

    2022-11-11 06:07:08張桀睿馬廣程夏紅偉
    中國慣性技術(shù)學(xué)報 2022年4期
    關(guān)鍵詞:轉(zhuǎn)動慣量航天器特性

    徐 旭,張桀睿,王 冠,馬廣程,夏紅偉

    (哈爾濱工業(yè)大學(xué)空間控制與慣性技術(shù)研究中心,哈爾濱150001)

    隨著航天科技的發(fā)展,航天器結(jié)構(gòu)變得愈加復(fù)雜,新型航天器上往往搭載多種活動載荷,例如機械臂、繩網(wǎng)捕獲機構(gòu)以及視覺敏感裝置等。在活動載荷工作時,航天器質(zhì)心以及轉(zhuǎn)動慣量等質(zhì)量特性參數(shù)會產(chǎn)生一定的改變,給航天器的姿態(tài)控制帶來一定影響。為保證航天器姿態(tài)控制的精度、快速性和穩(wěn)定性,實時獲得航天器的質(zhì)量特性參數(shù)是十分必要的。

    在工程實際中,航天器的質(zhì)量特性往往僅能夠通過理論分析或者間接測量獲得,例如通過質(zhì)量特性測量裝置測量、質(zhì)量分解計算以及有限元分析等方法。但這些方法測量精度不高且對參數(shù)時變系統(tǒng)無能為力,這就導(dǎo)致了在航天器執(zhí)行空間任務(wù)時,實時變化的質(zhì)量特性參數(shù)會對航天器產(chǎn)生未知擾動,嚴重影響航天器的性能[1-3]。

    對于質(zhì)量特性變化產(chǎn)生的擾動,可以使用魯棒性強的控制律來抑制,如滑模控制、模糊控制以及自適應(yīng)控制。但此類控制律均是被動抑制方法,而質(zhì)量特性辨識方法則是從根源上補償了該類擾動,利用慣組等測量元件和飛輪、噴氣發(fā)動機等執(zhí)行元件,在線動態(tài)辨識出航天器質(zhì)量特性參數(shù)及其改變情況,對控制器參數(shù)設(shè)計提供依據(jù)[4-6]。

    Bergman、Williams以及Wilson等人相繼提出了高斯二階濾波辨識方法、基于指數(shù)加權(quán)遞歸最小二乘法等算法[7],Wilson提出了基于指數(shù)加權(quán)遞歸最小二乘辨識衛(wèi)星的質(zhì)量特性參數(shù)的方法[8]?;谝陨涎芯浚谧钚《朔ǖ谋孀R手段在各類工程應(yīng)用中愈加廣泛[9-11]。此外,王書廷、侯振東等人提出了采用推力器組合為激勵的衛(wèi)星質(zhì)量特性辨識方法,該方法將質(zhì)心和慣量矩陣進行解耦辨識,減小了耦合帶來的辨識誤差[12];何驍?shù)柔槍教炱鞔蠼嵌葯C動問題,基于卡爾曼濾波的撓性附件振動模態(tài)狀態(tài)估計,提出了一種帶撓性航天器轉(zhuǎn)動慣量的辨識方法[13];倪智宇等提出了改進的基于遞歸預(yù)測器的子空間識別(Recursive Predictor-based Subspace Identification,RPBSID)方法并應(yīng)用于衛(wèi)星時變模態(tài)參數(shù)辨識,減小了辨識過程中的數(shù)據(jù)量,顯著節(jié)約了在軌辨識過程的計算時間[14]。但以上方法均為離線或半在線辨識,無法滿足實時、在線辨識任務(wù)的需求,在質(zhì)量特性動態(tài)改變上缺少或具有較差的跟蹤性。除此之外,以上方法大都是基于理論和數(shù)字仿真層面,缺少全物理試驗環(huán)節(jié)和對實際應(yīng)用的指導(dǎo)性。

    由于航天器工作環(huán)境的特殊性,空間任務(wù)需要在地面上進行仿真試驗驗證。目前,國內(nèi)外許多科研單位已對航天器地面仿真技術(shù)進行了相應(yīng)的研究,且大多數(shù)基于數(shù)字仿真,而為了充分進行空間任務(wù)的地面驗證,有必要開展地面半物理仿真以及全物理仿真試驗[15]。常用的方法有懸掛法,落體法,液浮法和氣浮法,其中氣浮法的應(yīng)用最為廣泛[16,17]。我國對基于氣浮仿真的航天器全物理仿真系統(tǒng)研究起步較晚,北京控制工程研究所開啟了國內(nèi)三軸氣浮臺技術(shù)應(yīng)用的先河。隨后哈爾濱工業(yè)大學(xué)和上海衛(wèi)星工程研究所等單位開展了一系列有關(guān)三軸氣浮臺的研究[18]。

    在這種背景下,針對上述航天器質(zhì)量特性的可辨識性及辨識問題,本文分析了影響辨識模型可辨識性的因素,根據(jù)此類辨識模型并基于漸消記憶遞推最小二乘辨識理論,提出了一種改進遞推最小二乘法的高效率在線辨識方法,與傳統(tǒng)方法相比,該方法可解決航天器質(zhì)量特性在軌實時、在線辨識的任務(wù)需求,為高精度姿態(tài)控制提供充分的先驗信息,解決了以往傳統(tǒng)離線或半在線辨識方法實時跟蹤性能差,實用性不高的缺點,且具有收斂速度快,跟蹤效率高以及辨識誤差小的優(yōu)點;同時為了驗證該算法在實際工程任務(wù)中的可行性,本文提出了一種基于三軸氣浮臺的航天器地面全物理仿真系統(tǒng),通過數(shù)字仿真和全物理仿真試驗的方式驗證了上述算法的可行性,解決了質(zhì)量特性辨識方法常常停留在理論分析及數(shù)學(xué)仿真層面上的問題,彌補了質(zhì)量特性辨識在全物理仿真試驗上的空缺。

    本文首先針對航天器質(zhì)量特性辨識問題給出了一類航天器質(zhì)量特性辨識模型以及其在線辨識算法;然后,針對該類模型分析了其可辨識性的影響因素,提出了一種基于三軸氣浮臺的地面全物理仿真系統(tǒng);最后,基于該全物理仿真系統(tǒng),通過數(shù)字仿真與全物理試驗兩種方法分析了在線辨識算法的收斂速度與辨識誤差,驗證了該算法的正確性及可行性。

    1 數(shù)學(xué)模型

    1.1 坐標系定義

    在研究航天器軌道動力學(xué)模型中,需對不同的坐標系進行定義,以便于航天器姿態(tài)描述。本文中所需要的坐標系定義如下:

    (1)地心慣性坐標系

    地心慣性坐標系Ox y z原點為地球質(zhì)心,x軸與地球赤道平面交于春分點,z軸垂直于赤道平面指向地球自轉(zhuǎn)角速度方向。

    (2)軌道坐標系

    軌道坐標系Ox0y0z0原點為航天器質(zhì)心,z0軸從原點指向地心,x0軸在軌道面內(nèi)與z0軸垂直,并指向航天器速度方向。y0軸在軌道平面負法線方向,形成右手正交坐標系。

    (3)本體坐標系

    本體坐標系Oxbybzb原點為航天器質(zhì)心。初始狀態(tài)下,航天器的本體坐標系與軌道坐標系重合,當(dāng)航天器發(fā)生姿態(tài)運動時,本體坐標系開始旋轉(zhuǎn),不再與軌道坐標系重合。

    1.2 動力學(xué)分析

    僅有飛輪作為航天器的執(zhí)行機構(gòu)時,航天器的動力學(xué)方程為:

    其中,J為航天器的轉(zhuǎn)動慣量矩陣;Jf為飛輪的轉(zhuǎn)動慣量;ω為航天器本體轉(zhuǎn)動在慣性坐標系下投影的角速度;ωf為飛輪角速度。

    整理式(1)可得動力學(xué)差分方程如下:

    其中

    其中ωx、ωy和ωz分別為航天器各個軸的角速度;Tx、Ty和Tz分別為航天器各個軸的輸入力矩,該力矩由飛輪提供;dt為系統(tǒng)的采樣時間;ax、ay和az分別為航天器三個軸上飛輪的加速度。

    1.3 一種改進遞推最小二乘法質(zhì)量特性在線辨識算法

    上一小節(jié)給出了航天器的動力學(xué)模型,本小節(jié)根據(jù)該模型給出一種改進遞推最小二乘法的質(zhì)量特性在線辨識算法。

    定義:

    對于實對稱陣P,定義:

    根據(jù)系統(tǒng)辨識理論,該航天器質(zhì)量特性辨識模型的最優(yōu)解為:

    其中

    由于基本最小二乘法不適用于具有參數(shù)時變特性的動態(tài)系統(tǒng),為了在線辨識出航天器質(zhì)量特性參數(shù)及其實時變化情況,采用改進遞推最小二乘法具有更好的效果。

    依據(jù)式(2)和式(3),可建立以ε1、ε2和ε3為待辨識參數(shù)的遞推最小二乘公式:

    其中,E為單位矩陣;ξ∈(0,1)為收斂因子,用于調(diào)節(jié)辨識收斂速度。需指出,當(dāng)選取較大的收斂因子時,迭代運算前期辨識結(jié)果振蕩幅度越小,迭代過程辨識結(jié)果收斂速度越慢;當(dāng)選取較小的收斂因子時,迭代過程辨識結(jié)果收斂速度較快,但前期振蕩明顯,很容易導(dǎo)致辨識結(jié)果發(fā)散。因此,在該因子的選取上,需兼顧考慮上述因素,以獲得迭代辨識過程的最佳收斂效果。λ∈(0,1)為遺忘因子,通過對舊數(shù)據(jù)按指數(shù)加權(quán),一定程度上衰減舊數(shù)據(jù)的作用,進而達到辨識數(shù)據(jù)更新的效果。通常該因子越小,舊數(shù)據(jù)的衰減效果越明顯。

    分別對x、y和z軸進行遞推運算,可得:

    該算法流程如表1所示。

    表1 改進遞推最小二乘法的在線辨識算法Tab.1 Online identification method based on the improved asymptotic least squares

    2 可辨識性分析

    在工程應(yīng)用中,必須要考慮參數(shù)的可辨識性,其代表著模型中的每一個參數(shù)是否能夠通過輸入輸出的觀測數(shù)據(jù)來確定。如果上述說法對于某一模型成立,則代表著模型是參數(shù)可辨識的。

    模型的參數(shù)可辨識性是模型本身的固有屬性,依賴于模型的觀測數(shù)據(jù),不依賴于模型參數(shù)的辨識算法。一種方法是分析辨識模型來得到模型的參數(shù)可辨識性。該方法適用于理論證明,但在工程應(yīng)用中缺少應(yīng)用性。而工程實際中常用的方法是設(shè)計一些激勵信號來保證系統(tǒng)是可辨識的。

    2.1 模型的參數(shù)可辨識性分析

    由式(10)知,系數(shù)矩陣A由系統(tǒng)的姿態(tài)數(shù)據(jù)構(gòu)成,故只需姿態(tài)數(shù)據(jù)保證T0>A A成立,即可滿足參數(shù)可辨識條件2。

    同時滿足條件1和條件2,該模型參數(shù)可辨識。

    2.2 激勵信號分析

    系數(shù)矩陣A由航天器輸出的姿態(tài)數(shù)據(jù)組成,所以在實際應(yīng)用中,分析系數(shù)矩陣A確定系統(tǒng)參數(shù)可辨識性的方法比較被動,我們難以主動地保證系統(tǒng)是可辨識的。為了解決上述問題,我們要主動選擇并設(shè)計一些激勵信號,將模型的參數(shù)可辨識性從被動證明轉(zhuǎn)為主動滿足。

    從系統(tǒng)辨識理論來看,可辨識性與模型的輸入密切相關(guān),其依賴于模型的輸入輸出觀測數(shù)據(jù)。為了保證模型是可辨識的,激勵信號要在整個觀測周期上持續(xù)激勵出系統(tǒng)的所有模態(tài)。

    當(dāng)系統(tǒng)參數(shù)無法通過觀測數(shù)據(jù)辨識出來時,一個原因可能是輸入輸出的觀測數(shù)據(jù)不合格。換言之很有可能是激勵信號無法有效激勵出系統(tǒng)的各個模態(tài),不能滿足參數(shù)辨識的要求,導(dǎo)致系統(tǒng)參數(shù)不能正確地辨識出來。

    引理2[19]:滿足需求的激勵信號能夠激勵出系統(tǒng)的全部模態(tài),使得系統(tǒng)的輸出數(shù)據(jù)包含系統(tǒng)的全部信息,此種信號就被定義為是持續(xù)激勵信號。

    其數(shù)學(xué)表述為:

    對于激勵信號r(t)∈R,定義一個n維激勵信號列向量:

    下式中I為單位陣,若存在整數(shù)N≥n和常數(shù)α> 0使得以下公式成立:

    其中

    t≥N+n-2

    則稱激勵信號r(t)是n階持續(xù)激勵信號。該條件被稱為是持續(xù)激勵條件。

    其中正常數(shù)α在不同的激勵信號背景下是有不同取值的。若存在某一待辨識系統(tǒng)有n個參數(shù),則要求其激勵信號至少要滿足n階條件。除此之外,對于確定性系統(tǒng)或模型,其激勵信號不需要對每一個持續(xù)激勵條件中的t都成立,只需存在一個t≥N+n-2使得該條件成立即可。

    3 地面全物理仿真系統(tǒng)設(shè)計方案

    為了提高算法的部署效率,本文研究了基于三軸氣浮臺的航天器地面全物理仿真系統(tǒng)的設(shè)計與實現(xiàn),并將本文所提算法在該系統(tǒng)中進行試驗、分析與驗證。

    3.1 系統(tǒng)組成

    基于三軸氣浮臺的航天器地面全物理仿真系統(tǒng)分為臺上綜合試驗系統(tǒng)和臺下監(jiān)控系統(tǒng)。

    臺上試驗系統(tǒng)由三軸氣浮臺及儀表平臺、臺上管理控制單元、電源模塊、調(diào)平衡模塊、數(shù)據(jù)傳輸功能模塊、轉(zhuǎn)動慣量測量模塊等組成,如圖2-3所示。儀表平臺主要為臺上試驗系統(tǒng)各組成部分提供硬件安裝接口;臺上管理控制單元主要完成以下功能:

    圖2 三軸氣浮臺實物圖Fig.2 Physical drawing of the three-axis air bearing testbed

    (1)對臺上的調(diào)平衡模塊的控制;

    (2)對姿態(tài)測量信息的采集和管理;

    (3)對臺上的電源模塊進行綜合管理和控制;

    (4)對臺上仿真試驗參試設(shè)備的控制和管理;

    (5)對臺上工作狀態(tài)的記錄系統(tǒng)控制和管理。

    電源模塊用來給臺上試驗系統(tǒng)所有電氣設(shè)備供電;調(diào)平衡模塊用于三軸氣浮臺于試驗前的質(zhì)心調(diào)節(jié)工作;數(shù)據(jù)傳輸功能模塊用于收集臺上所有實驗數(shù)據(jù)以及各設(shè)備的工作狀態(tài)并發(fā)送給臺下監(jiān)控系統(tǒng);轉(zhuǎn)動慣量測量模塊負責(zé)辨識分析三軸氣浮臺的轉(zhuǎn)動慣量參數(shù)。

    臺下監(jiān)控系統(tǒng)由數(shù)據(jù)傳輸功能模塊、安全輔助功能模塊以及臺下監(jiān)控軟件組成。數(shù)據(jù)傳輸功能模塊負責(zé)接收由臺上數(shù)據(jù)傳輸功能模塊發(fā)送來的實驗數(shù)據(jù);安全輔助功能模塊包括傘狀支撐和千斤頂,分別用于三軸氣浮臺運動限位以及支撐;臺下監(jiān)控軟件用于顯示、分析、存儲和回放該實驗數(shù)據(jù),并向臺上綜合試驗系統(tǒng)下達控制指令。

    圖3 地面全物理仿真系統(tǒng)組成示意圖Fig.3 Ground full physicssimulation system compositi on diagram

    3.2 工作原理

    三軸氣浮臺通過供氣系統(tǒng)向氣浮球和球軸承中提供壓縮空氣形成氣膜來使其儀表平臺浮起,來模擬衛(wèi)星在軌時微干擾力微干擾力矩的動力學(xué)環(huán)境。

    通過在三軸氣浮臺上安裝星上載荷、星載計算機、慣性測量單元以及角動量交換裝置等設(shè)備來搭建三軸氣浮臺的測控系統(tǒng)。再配置好其微干擾力矩的工作環(huán)境,包括降低環(huán)境風(fēng)干擾,調(diào)節(jié)三軸氣浮臺質(zhì)心至其旋轉(zhuǎn)中心位置,減小重力偏心矩的影響,來使其承載的試驗系統(tǒng)滿足試驗要求,以完成衛(wèi)星在軌動力學(xué)的各項地面試驗工作。

    4 數(shù)字仿真與地面全物理試驗分析

    第一節(jié)給出了航天器在軌動力學(xué)模型以及質(zhì)量特性在線辨識算法,第二節(jié)給出了模型的參數(shù)可辨識性,本節(jié)通過數(shù)學(xué)仿真以及地面試驗驗證該辨識算法的正確性與可行性。

    4.1 數(shù)字仿真分析

    設(shè)該仿真模型中,航天器的轉(zhuǎn)動慣量矩陣J(單位為kg·m2)為:

    設(shè)該模型中三個軸上的飛輪輸出幅值為±100mNm的力矩方波信號,信號周期分別為14s,16 s和18s,采樣時間為0.02s;初始化辨識信息θ(0)=0,P(0)=αI,α=106~1010;設(shè)置遺忘因子λ=0.99,收斂因子ξ=0.15。得到的辨識結(jié)果如圖4-6所示。

    圖4 Jx辨識仿真結(jié)果對比圖Fig.4 Jxidentification simulation result comparisonchart

    圖5 Jy辨識仿真結(jié)果對比圖Fig.5 Jyidentification simulat ion result comparisonchart

    圖6 Jz 辨識仿真結(jié)果對比圖Fig.6 Jzidentification simulation result comparisonchart

    改進遞推最小二乘法和傳統(tǒng)方法仿真結(jié)果精度一致,分別為Jx= 386.4 kg·m2,Jy=401.7kg·m2以及Jz= 385.6kg·m2,誤差率分別為1.2%,0.7%和1.6%。如圖7所示,以每5s數(shù)據(jù)變化率小于0.5%作為收斂指標,改進遞推最小二乘法較傳統(tǒng)方法具有更快的收斂速度,在實際應(yīng)用中更具優(yōu)勢。

    圖7 辨識仿真結(jié)果對比圖Fig.7 Identification simulation resultscomparison chart

    4.2 地面全物理仿真試驗分析

    下面通過地面全物理仿真試驗來驗證本文算法的正確性與可行性。在試驗過程中,三個軸上的飛輪力矩為該系統(tǒng)的力矩輸入,通過該力矩信息以及慣性部件采集到的姿態(tài)數(shù)據(jù)即可辨識出轉(zhuǎn)動慣量矩陣。用于該項地面試驗的三軸氣浮臺的真實轉(zhuǎn)動慣量矩陣(單位為kg·m2)為:

    令三軸氣浮臺x,y和z軸上飛輪輸出與仿真分析中相同的力矩方波信號,即信號幅值為±100mNm,信號周期分別為14s,16s和18s。采樣時間為0.02s,辨識得到三軸氣浮臺各個軸轉(zhuǎn)動慣量如圖8-10所示。

    圖8 Jx辨識結(jié)果圖Fig.8 Jxidentification result c hart

    圖9 Jy辨識結(jié)果圖Fig.9 Jyidentification result chart

    辨識得到的三軸氣浮臺各個軸的轉(zhuǎn)動慣量分別為Jx= 402.4 kg·m2,Jy=408.1kg·m2以及Jz= 406.7kg·m2。誤差率分別為2.8%,2.4%以及3.8%。

    試驗時還分析了數(shù)據(jù)采樣時間對辨識精度的影響。分別選取數(shù)據(jù)采樣時間為0.01s,0.02s,0.04s,0.05s和0.1s,得到了不同收斂速度下的辨識結(jié)果如圖11-13所示。

    圖1 坐標系定義圖Fig.1 Coordinate system definition diagram

    圖10 Jz辨識結(jié)果圖Fig.10 Jzidentificati on result chart

    圖11 Jx在不同數(shù)據(jù)采樣時間下的辨識結(jié)果Fig.11 Identification resultsof Jxat different sampling times

    圖12 Jy在不同數(shù)據(jù)采樣時間下的辨識結(jié)果Fig.12 Identification resultsof Jyat different sam pling times

    圖13 Jz在不同數(shù)據(jù)采樣時間下的辨識結(jié)果Fig.13 Identification results of Jzat different sampling times

    地面全物理仿真試驗結(jié)果討論:

    (1)由于地面試驗中慣性測量單元采集到的姿態(tài)信息會不可避免地存在環(huán)境干擾及噪聲,所以在不同的數(shù)據(jù)采樣時間下,會得到不同的轉(zhuǎn)動慣量辨識收斂速度和結(jié)果。

    (2)在該次地面試驗的平臺及環(huán)境條件下,采樣時間為0.02 s下的辨識收斂速度最快,辨識誤差最小,為2.8%。除采樣時間0.01 s以外,采樣時間越長,辨識收斂速度越慢,辨識誤差越大。

    (3)采樣時間為0.01 s時放大了測量噪聲對辨識的影響,導(dǎo)致其辨識誤差及收斂速度不如采樣時間為0.02 s的。

    (4)五種數(shù)據(jù)采樣時間下的辨識誤差率均在4.5%以內(nèi),滿足工程應(yīng)用需求。

    通過本次仿真分析以及地面試驗,充分驗證了本文給出的轉(zhuǎn)動慣量在線辨識算法的可行性與有效性。

    5 結(jié)論

    本文針對航天器質(zhì)量特性在線辨識問題,提出一種高效率在線辨識方法和地面全物理仿真試驗系統(tǒng)。通過建立適用于航天器轉(zhuǎn)動慣量辨識的數(shù)學(xué)模型,分析了其可辨識性和數(shù)據(jù)采樣時間的影響機理,基于漸消記憶遞推最小二乘辨識理論,提出了一種改進遞推最小二乘法的在線辨識方法。為了驗證在線辨識方法的正確性與可行性,本文設(shè)計了一種基于三軸氣浮臺的地面全物理仿真試驗系統(tǒng)。數(shù)學(xué)分析表明,改進遞推最小二乘法的在線辨識方法收斂速度優(yōu)于傳統(tǒng)辨識方法;全物理試驗表明,該算法的辨識精度優(yōu)于96.2%,驗證了本文提出的航天器質(zhì)量特性在線辨識方法的可行性與高效性,突出了該方法具有實用性強、精度高、收斂速度快的優(yōu)點,對航天器質(zhì)量特性在線辨識的算法設(shè)計和地面驗證具有一定的參考應(yīng)用價值。

    猜你喜歡
    轉(zhuǎn)動慣量航天器特性
    2022 年第二季度航天器發(fā)射統(tǒng)計
    國際太空(2022年7期)2022-08-16 09:52:50
    谷稗的生物學(xué)特性和栽培技術(shù)
    差值法巧求剛體轉(zhuǎn)動慣量
    三線擺測剛體轉(zhuǎn)動慣量誤差分析及改進
    色彩特性
    流行色(2020年9期)2020-07-16 08:08:54
    2019 年第二季度航天器發(fā)射統(tǒng)計
    國際太空(2019年9期)2019-10-23 01:55:34
    2018 年第三季度航天器發(fā)射統(tǒng)計
    國際太空(2018年12期)2019-01-28 12:53:20
    進一步凸顯定制安裝特性的優(yōu)勢 Integra DRX-5.2
    2018年第二季度航天器發(fā)射統(tǒng)計
    國際太空(2018年9期)2018-10-18 08:51:32
    Quick Charge 4:什么是新的?
    CHIP新電腦(2017年6期)2017-06-19 09:41:44
    一进一出好大好爽视频| 亚洲最大成人手机在线| 香蕉丝袜av| 18禁美女被吸乳视频| www.色视频.com| 日本黄大片高清| 在线观看舔阴道视频| 免费看光身美女| 熟妇人妻久久中文字幕3abv| 观看免费一级毛片| 日韩欧美在线二视频| 99久国产av精品| 色尼玛亚洲综合影院| 真实男女啪啪啪动态图| 日韩欧美在线乱码| 国产精品一区二区三区四区免费观看 | av欧美777| 午夜a级毛片| 午夜a级毛片| 亚洲熟妇中文字幕五十中出| 久久久精品欧美日韩精品| 欧美性猛交╳xxx乱大交人| 最近最新中文字幕大全免费视频| 成人午夜高清在线视频| 久久九九热精品免费| 中国美女看黄片| 99在线人妻在线中文字幕| 极品教师在线免费播放| 麻豆国产av国片精品| 1024手机看黄色片| 色视频www国产| 精品久久久久久成人av| 精品久久久久久久久久久久久| 成人一区二区视频在线观看| 国产精品亚洲一级av第二区| 99久久精品热视频| 久久久久九九精品影院| 丰满人妻熟妇乱又伦精品不卡| 天天一区二区日本电影三级| 精品久久久久久,| 舔av片在线| 久久久国产成人免费| 精品人妻一区二区三区麻豆 | 岛国视频午夜一区免费看| 色精品久久人妻99蜜桃| 久久久久精品国产欧美久久久| 一个人看视频在线观看www免费 | 中文字幕人妻熟人妻熟丝袜美 | 亚洲国产色片| 亚洲在线自拍视频| 三级国产精品欧美在线观看| 一区福利在线观看| 美女cb高潮喷水在线观看| 亚洲精品乱码久久久v下载方式 | 国产麻豆成人av免费视频| 日日摸夜夜添夜夜添小说| av在线蜜桃| 免费高清视频大片| 亚洲人与动物交配视频| 国产午夜精品论理片| 久久久久亚洲av毛片大全| 久久久久久久午夜电影| 熟女人妻精品中文字幕| 国产欧美日韩精品一区二区| 国产精品久久久久久亚洲av鲁大| 国产精品精品国产色婷婷| 成人av在线播放网站| 丁香六月欧美| 91字幕亚洲| a级毛片a级免费在线| 国产麻豆成人av免费视频| 看黄色毛片网站| 伊人久久大香线蕉亚洲五| 亚洲熟妇熟女久久| 久99久视频精品免费| 亚洲av成人精品一区久久| 日本黄大片高清| 变态另类丝袜制服| 国产中年淑女户外野战色| 中文字幕高清在线视频| 美女大奶头视频| 69av精品久久久久久| 老司机福利观看| 成人国产综合亚洲| 亚洲精品色激情综合| 日本一二三区视频观看| 免费av毛片视频| 亚洲av熟女| 欧洲精品卡2卡3卡4卡5卡区| 小说图片视频综合网站| 久久香蕉国产精品| av片东京热男人的天堂| 一个人看的www免费观看视频| 男女下面进入的视频免费午夜| 亚洲aⅴ乱码一区二区在线播放| 亚洲男人的天堂狠狠| 在线观看66精品国产| 国产伦精品一区二区三区视频9 | 国产乱人伦免费视频| 亚洲激情在线av| 久久精品91无色码中文字幕| 九色国产91popny在线| 国产一级毛片七仙女欲春2| 国产高清视频在线观看网站| bbb黄色大片| 国产熟女xx| 女人高潮潮喷娇喘18禁视频| 成人午夜高清在线视频| 国内精品久久久久精免费| 欧美日韩黄片免| 亚洲精品美女久久久久99蜜臀| 嫁个100分男人电影在线观看| 日韩大尺度精品在线看网址| 欧美一区二区精品小视频在线| 久久久国产精品麻豆| 99热6这里只有精品| 窝窝影院91人妻| 女人高潮潮喷娇喘18禁视频| 亚洲国产精品成人综合色| 成人亚洲精品av一区二区| 亚洲色图av天堂| 一区福利在线观看| 成人性生交大片免费视频hd| 国产亚洲精品综合一区在线观看| 久久精品国产自在天天线| 免费搜索国产男女视频| 国产精品电影一区二区三区| 国产男靠女视频免费网站| 亚洲七黄色美女视频| 精品日产1卡2卡| 精品福利观看| 国产欧美日韩一区二区三| 男人和女人高潮做爰伦理| 怎么达到女性高潮| 婷婷精品国产亚洲av在线| 国内毛片毛片毛片毛片毛片| 在线看三级毛片| 日韩国内少妇激情av| 99久久精品热视频| 国产精品美女特级片免费视频播放器| 丁香六月欧美| 婷婷六月久久综合丁香| 男插女下体视频免费在线播放| 丁香六月欧美| 国产精品av视频在线免费观看| 99国产精品一区二区蜜桃av| 黑人欧美特级aaaaaa片| 久久久久久九九精品二区国产| 国产精品日韩av在线免费观看| 国产视频内射| 女同久久另类99精品国产91| www.999成人在线观看| 亚洲成av人片免费观看| 真人做人爱边吃奶动态| 91av网一区二区| 午夜福利在线观看吧| 啦啦啦观看免费观看视频高清| 亚洲av成人不卡在线观看播放网| 亚洲欧美一区二区三区黑人| 性欧美人与动物交配| 香蕉久久夜色| 一二三四社区在线视频社区8| 国产亚洲欧美98| 久久久久国产精品人妻aⅴ院| 久久久久久久久大av| 熟女人妻精品中文字幕| 小说图片视频综合网站| 免费大片18禁| 国产成人av激情在线播放| av中文乱码字幕在线| 欧美激情在线99| 日本一本二区三区精品| 好男人电影高清在线观看| 69人妻影院| 最新美女视频免费是黄的| 日韩免费av在线播放| 免费大片18禁| 国产色婷婷99| 中文字幕人成人乱码亚洲影| 亚洲人成电影免费在线| 国语自产精品视频在线第100页| 久久久成人免费电影| 国产欧美日韩一区二区精品| 国产精品国产高清国产av| www.999成人在线观看| 最新在线观看一区二区三区| 波野结衣二区三区在线 | 国产 一区 欧美 日韩| 欧美一区二区精品小视频在线| 一夜夜www| 国产aⅴ精品一区二区三区波| 一进一出抽搐gif免费好疼| 91在线精品国自产拍蜜月 | 综合色av麻豆| 亚洲在线自拍视频| 日本 av在线| 在线观看日韩欧美| 日韩国内少妇激情av| 美女高潮喷水抽搐中文字幕| 欧美中文日本在线观看视频| 国产一区二区三区视频了| 免费观看人在逋| 国产 一区 欧美 日韩| 99国产极品粉嫩在线观看| 国产精品久久久久久人妻精品电影| 中文字幕人妻熟人妻熟丝袜美 | 国内久久婷婷六月综合欲色啪| 成人av一区二区三区在线看| 亚洲国产精品成人综合色| 啦啦啦韩国在线观看视频| 色在线成人网| 国产成人av教育| 一个人看视频在线观看www免费 | 婷婷亚洲欧美| 国内毛片毛片毛片毛片毛片| 看免费av毛片| 日韩欧美精品v在线| 性色avwww在线观看| 人人妻人人看人人澡| 中文字幕人成人乱码亚洲影| 欧美色视频一区免费| 成人欧美大片| 国产乱人伦免费视频| 国产精品美女特级片免费视频播放器| 日日摸夜夜添夜夜添小说| 亚洲欧美日韩无卡精品| 噜噜噜噜噜久久久久久91| 久久久国产成人精品二区| 99久久九九国产精品国产免费| 久久精品国产亚洲av香蕉五月| 一本久久中文字幕| 亚洲中文日韩欧美视频| 亚洲精品日韩av片在线观看 | 国产久久久一区二区三区| 久久久久久久久中文| 中文字幕av成人在线电影| 美女免费视频网站| 免费看十八禁软件| 日本在线视频免费播放| 成年女人看的毛片在线观看| 成年女人永久免费观看视频| aaaaa片日本免费| netflix在线观看网站| 每晚都被弄得嗷嗷叫到高潮| 国产成人福利小说| 成人鲁丝片一二三区免费| 在线播放无遮挡| 久久欧美精品欧美久久欧美| 国产精品精品国产色婷婷| 老司机午夜福利在线观看视频| 国产精品乱码一区二三区的特点| 3wmmmm亚洲av在线观看| 欧美xxxx黑人xx丫x性爽| 国产在线精品亚洲第一网站| 不卡一级毛片| 午夜福利欧美成人| 欧美又色又爽又黄视频| 国产亚洲欧美在线一区二区| e午夜精品久久久久久久| 亚洲成av人片在线播放无| avwww免费| 青草久久国产| 免费一级毛片在线播放高清视频| 嫩草影视91久久| 国产精品女同一区二区软件 | 日本一二三区视频观看| 波多野结衣巨乳人妻| 少妇人妻一区二区三区视频| 午夜免费激情av| 亚洲无线观看免费| 国产精品av视频在线免费观看| 别揉我奶头~嗯~啊~动态视频| 天堂av国产一区二区熟女人妻| 人人妻人人澡欧美一区二区| 欧美激情在线99| 久久人妻av系列| 一进一出抽搐gif免费好疼| 极品教师在线免费播放| 亚洲男人的天堂狠狠| 国产午夜福利久久久久久| 色综合站精品国产| 免费看日本二区| 黄色成人免费大全| 老汉色av国产亚洲站长工具| 在线看三级毛片| 成年版毛片免费区| 一二三四社区在线视频社区8| x7x7x7水蜜桃| 亚洲色图av天堂| 长腿黑丝高跟| 夜夜夜夜夜久久久久| 欧美日韩精品网址| 日韩欧美三级三区| 很黄的视频免费| 欧美激情久久久久久爽电影| 18禁国产床啪视频网站| 51国产日韩欧美| 亚洲精品一卡2卡三卡4卡5卡| 亚洲内射少妇av| 欧美黄色淫秽网站| 日韩精品中文字幕看吧| 午夜免费观看网址| av片东京热男人的天堂| www日本在线高清视频| av欧美777| 免费看十八禁软件| 99精品在免费线老司机午夜| 嫩草影视91久久| 久久6这里有精品| 久久人妻av系列| 国产精品影院久久| 最近最新免费中文字幕在线| 99热精品在线国产| 少妇人妻一区二区三区视频| www.www免费av| 美女被艹到高潮喷水动态| 日韩成人在线观看一区二区三区| 精品一区二区三区人妻视频| 亚洲国产欧洲综合997久久,| 国产精品久久电影中文字幕| 国产主播在线观看一区二区| 日本在线视频免费播放| 国产精品久久久人人做人人爽| 老司机午夜福利在线观看视频| 最好的美女福利视频网| 欧美3d第一页| 午夜影院日韩av| 在线免费观看不下载黄p国产 | 在线十欧美十亚洲十日本专区| 久久久久久国产a免费观看| 男女那种视频在线观看| 黄片小视频在线播放| 亚洲国产精品久久男人天堂| 国产蜜桃级精品一区二区三区| 午夜福利成人在线免费观看| 91av网一区二区| 首页视频小说图片口味搜索| 男人舔奶头视频| 国产美女午夜福利| 老汉色∧v一级毛片| www.www免费av| 在线观看66精品国产| 夜夜爽天天搞| 国产精品亚洲美女久久久| 国产精品精品国产色婷婷| 18+在线观看网站| 九九热线精品视视频播放| 国产伦精品一区二区三区四那| 亚洲18禁久久av| 欧洲精品卡2卡3卡4卡5卡区| 免费av毛片视频| 久久久久精品国产欧美久久久| 黄色女人牲交| 欧洲精品卡2卡3卡4卡5卡区| 久久久久久久久大av| 国产高清videossex| 久久99热这里只有精品18| svipshipincom国产片| 搡女人真爽免费视频火全软件 | 精品人妻1区二区| 别揉我奶头~嗯~啊~动态视频| 黄色日韩在线| 丰满乱子伦码专区| 午夜福利在线观看吧| 黄片大片在线免费观看| 成人欧美大片| 亚洲精品成人久久久久久| 国产精华一区二区三区| 国产精品亚洲美女久久久| xxx96com| 亚洲美女黄片视频| 在线观看免费午夜福利视频| av福利片在线观看| 欧美性感艳星| 最近最新中文字幕大全电影3| 免费看a级黄色片| 黄片大片在线免费观看| 天天一区二区日本电影三级| 此物有八面人人有两片| 久久久久久大精品| 校园春色视频在线观看| 国内精品美女久久久久久| 国产黄色小视频在线观看| 99精品在免费线老司机午夜| 男女做爰动态图高潮gif福利片| 国产精品嫩草影院av在线观看 | 日韩欧美一区二区三区在线观看| aaaaa片日本免费| 人妻夜夜爽99麻豆av| 性色avwww在线观看| 久久精品国产亚洲av涩爱 | 丁香欧美五月| 成人av在线播放网站| 亚洲专区国产一区二区| 最近视频中文字幕2019在线8| www.熟女人妻精品国产| 国产精品影院久久| 真人一进一出gif抽搐免费| 免费无遮挡裸体视频| 国产精品日韩av在线免费观看| 国产aⅴ精品一区二区三区波| 日韩欧美精品v在线| 最近在线观看免费完整版| 两个人的视频大全免费| 免费在线观看成人毛片| 久久久久久九九精品二区国产| 美女cb高潮喷水在线观看| 婷婷精品国产亚洲av| 宅男免费午夜| 色噜噜av男人的天堂激情| 一区二区三区免费毛片| 亚洲午夜理论影院| 少妇丰满av| 成年免费大片在线观看| 特级一级黄色大片| 亚洲国产欧美人成| 最新美女视频免费是黄的| 亚洲精品成人久久久久久| 在线免费观看不下载黄p国产 | 亚洲人成电影免费在线| 久久国产精品影院| 亚洲国产精品成人综合色| 国产精品99久久久久久久久| 日韩欧美一区二区三区在线观看| 亚洲精品国产精品久久久不卡| 欧美另类亚洲清纯唯美| 91字幕亚洲| 免费搜索国产男女视频| 99久久99久久久精品蜜桃| 老汉色av国产亚洲站长工具| 蜜桃久久精品国产亚洲av| 久久久久久九九精品二区国产| 俺也久久电影网| 精品福利观看| 搡老妇女老女人老熟妇| 色综合婷婷激情| 长腿黑丝高跟| xxx96com| 99久久精品热视频| 91麻豆av在线| 久久久久亚洲av毛片大全| 精品一区二区三区视频在线 | 丰满乱子伦码专区| 夜夜躁狠狠躁天天躁| 色噜噜av男人的天堂激情| 免费看日本二区| 国产野战对白在线观看| 成人性生交大片免费视频hd| 色哟哟哟哟哟哟| 国产亚洲av嫩草精品影院| 波野结衣二区三区在线 | 久久久久精品国产欧美久久久| 国产不卡一卡二| 狂野欧美激情性xxxx| 午夜影院日韩av| 性欧美人与动物交配| 五月伊人婷婷丁香| 老司机午夜福利在线观看视频| eeuss影院久久| 欧美另类亚洲清纯唯美| 国产真人三级小视频在线观看| 亚洲人成伊人成综合网2020| 久久国产乱子伦精品免费另类| 法律面前人人平等表现在哪些方面| 日日摸夜夜添夜夜添小说| 最近视频中文字幕2019在线8| 少妇熟女aⅴ在线视频| 男女视频在线观看网站免费| 亚洲国产高清在线一区二区三| 操出白浆在线播放| 国产真人三级小视频在线观看| 国产又黄又爽又无遮挡在线| 亚洲内射少妇av| 欧美乱色亚洲激情| 最后的刺客免费高清国语| 久久久久久人人人人人| 亚洲欧美精品综合久久99| 国产精品久久久人人做人人爽| 日日夜夜操网爽| 一个人免费在线观看电影| 欧美日韩精品网址| 一级作爱视频免费观看| av国产免费在线观看| 国产色爽女视频免费观看| 久久久久久大精品| 亚洲色图av天堂| 国产欧美日韩一区二区精品| 在线视频色国产色| xxx96com| 国内毛片毛片毛片毛片毛片| 国产老妇女一区| 国产成人av教育| 一a级毛片在线观看| 欧美国产日韩亚洲一区| 久久久久久久午夜电影| 免费无遮挡裸体视频| 最近在线观看免费完整版| 精品久久久久久,| 国产综合懂色| av在线天堂中文字幕| 免费无遮挡裸体视频| 国产淫片久久久久久久久 | 国产精品,欧美在线| 美女被艹到高潮喷水动态| 51午夜福利影视在线观看| 国产老妇女一区| 美女 人体艺术 gogo| av国产免费在线观看| 天堂网av新在线| 极品教师在线免费播放| 在线观看免费视频日本深夜| 久久久久性生活片| 婷婷六月久久综合丁香| 午夜免费男女啪啪视频观看 | 欧美+亚洲+日韩+国产| 亚洲国产中文字幕在线视频| 床上黄色一级片| 久久久久国产精品人妻aⅴ院| tocl精华| 亚洲av熟女| 不卡一级毛片| 色综合欧美亚洲国产小说| 他把我摸到了高潮在线观看| 亚洲成av人片在线播放无| 欧美成人一区二区免费高清观看| 大型黄色视频在线免费观看| 亚洲国产精品久久男人天堂| 亚洲无线在线观看| 亚洲va日本ⅴa欧美va伊人久久| 国产精华一区二区三区| 在线视频色国产色| 国产色婷婷99| 欧美zozozo另类| 老熟妇仑乱视频hdxx| 熟女少妇亚洲综合色aaa.| 麻豆国产97在线/欧美| 亚洲av五月六月丁香网| 51国产日韩欧美| 国产免费av片在线观看野外av| a级一级毛片免费在线观看| 男女视频在线观看网站免费| 此物有八面人人有两片| 国产成人欧美在线观看| 中文字幕人成人乱码亚洲影| 亚洲avbb在线观看| 免费在线观看日本一区| 丰满的人妻完整版| 午夜久久久久精精品| 最新中文字幕久久久久| 国产精品久久电影中文字幕| 五月伊人婷婷丁香| 一卡2卡三卡四卡精品乱码亚洲| 在线视频色国产色| 特大巨黑吊av在线直播| 一区福利在线观看| 亚洲avbb在线观看| 一进一出好大好爽视频| 国产成人a区在线观看| 国产三级黄色录像| 在线观看舔阴道视频| 两个人视频免费观看高清| 免费av不卡在线播放| 日本一二三区视频观看| 12—13女人毛片做爰片一| 国内精品久久久久久久电影| 国产成人系列免费观看| 亚洲国产中文字幕在线视频| 免费电影在线观看免费观看| 亚洲av成人精品一区久久| 成年女人永久免费观看视频| 国产熟女xx| 欧美在线一区亚洲| 午夜两性在线视频| 国产真实乱freesex| 午夜免费成人在线视频| 成年女人毛片免费观看观看9| 国产精品久久久人人做人人爽| 黑人欧美特级aaaaaa片| 亚洲国产欧美人成| 欧美+日韩+精品| 国产精品久久久久久人妻精品电影| 国产午夜福利久久久久久| 色综合欧美亚洲国产小说| 久久久久久国产a免费观看| av天堂中文字幕网| 97超级碰碰碰精品色视频在线观看| 国产成年人精品一区二区| 婷婷亚洲欧美| 在线看三级毛片| www.999成人在线观看| 国语自产精品视频在线第100页| 国产日本99.免费观看| 国产高清视频在线观看网站| 老司机在亚洲福利影院| 桃色一区二区三区在线观看| 天堂动漫精品| 黄色成人免费大全| 又爽又黄无遮挡网站| 国产亚洲欧美在线一区二区| 国产亚洲精品久久久久久毛片| 欧美一级毛片孕妇| 亚洲国产欧美人成| 亚洲成人免费电影在线观看| 国产v大片淫在线免费观看| 国产精品免费一区二区三区在线| 麻豆一二三区av精品| 亚洲aⅴ乱码一区二区在线播放| 日韩精品中文字幕看吧| 国产成人系列免费观看| 日本黄色片子视频| 久久这里只有精品中国| 精品久久久久久久久久免费视频| 亚洲av日韩精品久久久久久密| 欧美绝顶高潮抽搐喷水| 日本 欧美在线|