• <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
    看十八女毛片水多多多| 精品亚洲乱码少妇综合久久| 午夜福利高清视频| 亚洲在久久综合| 国产乱人偷精品视频| 亚洲人与动物交配视频| 欧美zozozo另类| 一区二区av电影网| 人妻一区二区av| 国产日韩欧美在线精品| 成人二区视频| 亚洲av不卡在线观看| 狂野欧美激情性xxxx在线观看| 极品教师在线视频| 在线观看免费高清a一片| 久久影院123| 精品久久久久久久人妻蜜臀av| 成人鲁丝片一二三区免费| 简卡轻食公司| 久久久久网色| 91久久精品国产一区二区三区| 26uuu在线亚洲综合色| 亚洲成人av在线免费| 1000部很黄的大片| 91久久精品电影网| 久久精品人妻少妇| 久久99精品国语久久久| 少妇人妻一区二区三区视频| 亚洲成人精品中文字幕电影| 搡女人真爽免费视频火全软件| 黄色怎么调成土黄色| 亚洲精品成人av观看孕妇| 精品视频人人做人人爽| 日韩av在线免费看完整版不卡| 国产国拍精品亚洲av在线观看| 日本一二三区视频观看| av专区在线播放| 亚洲av日韩在线播放| 我的女老师完整版在线观看| 中文乱码字字幕精品一区二区三区| 老女人水多毛片| 国产老妇女一区| 久久久久久久大尺度免费视频| 欧美日韩一区二区视频在线观看视频在线 | 久久久久久久久大av| 欧美xxxx黑人xx丫x性爽| 人妻夜夜爽99麻豆av| 免费高清在线观看视频在线观看| 亚洲av电影在线观看一区二区三区 | 日本爱情动作片www.在线观看| 91久久精品国产一区二区三区| 国产片特级美女逼逼视频| 国产精品麻豆人妻色哟哟久久| 色网站视频免费| 美女国产视频在线观看| 日韩亚洲欧美综合| 男人添女人高潮全过程视频| 国产女主播在线喷水免费视频网站| 中国三级夫妇交换| 久久国内精品自在自线图片| 欧美潮喷喷水| 自拍偷自拍亚洲精品老妇| av免费观看日本| 色视频在线一区二区三区| 成人鲁丝片一二三区免费| 在线精品无人区一区二区三 | 在线观看国产h片| 少妇人妻精品综合一区二区| 在线精品无人区一区二区三 | 免费观看无遮挡的男女| 精品人妻一区二区三区麻豆| 美女脱内裤让男人舔精品视频| 哪个播放器可以免费观看大片| 午夜亚洲福利在线播放| 丝瓜视频免费看黄片| 亚洲成色77777| 亚洲国产av新网站| 亚洲国产精品成人久久小说| 欧美老熟妇乱子伦牲交| 男女下面进入的视频免费午夜| 午夜免费鲁丝| 伊人久久国产一区二区| 边亲边吃奶的免费视频| 久久久久久国产a免费观看| 最近的中文字幕免费完整| 欧美变态另类bdsm刘玥| 美女主播在线视频| 最近手机中文字幕大全| 久久久欧美国产精品| 久久精品人妻少妇| 一二三四中文在线观看免费高清| 内地一区二区视频在线| 日韩在线高清观看一区二区三区| 国产精品99久久99久久久不卡 | 丝袜脚勾引网站| 久久97久久精品| 日本猛色少妇xxxxx猛交久久| 99久久精品热视频| 少妇丰满av| 精品午夜福利在线看| 激情 狠狠 欧美| 亚洲精品视频女| 免费av观看视频| 久久这里有精品视频免费| 国产伦精品一区二区三区四那| 亚洲av二区三区四区| 2021天堂中文幕一二区在线观| 精品久久久久久久久亚洲| 九色成人免费人妻av| 日韩av免费高清视频| 七月丁香在线播放| 国产黄a三级三级三级人| 一级二级三级毛片免费看| 国产亚洲午夜精品一区二区久久 | 亚洲性久久影院| 久久国产乱子免费精品| 婷婷色综合www| 国产免费又黄又爽又色| 亚洲成人精品中文字幕电影| 插阴视频在线观看视频| 在线观看一区二区三区激情| 国产毛片a区久久久久| 国产免费福利视频在线观看| 日韩av在线免费看完整版不卡| 婷婷色麻豆天堂久久| 欧美bdsm另类| 亚洲怡红院男人天堂| 中文欧美无线码| 久久精品久久久久久噜噜老黄| 黑人高潮一二区| 美女xxoo啪啪120秒动态图| 国产免费一区二区三区四区乱码| 国产 一区 欧美 日韩| 国产午夜福利久久久久久| 精品午夜福利在线看| 欧美日韩国产mv在线观看视频 | 国产精品国产三级国产av玫瑰| 亚洲怡红院男人天堂| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 91久久精品电影网| 日日啪夜夜撸| 国产免费一区二区三区四区乱码| 亚洲精品乱码久久久v下载方式| 亚洲丝袜综合中文字幕| 国产色婷婷99| 亚洲人成网站在线观看播放| 国产午夜精品一二区理论片| 国产欧美亚洲国产| 日韩av在线免费看完整版不卡| 一区二区三区精品91| www.色视频.com| 亚洲av二区三区四区| 国产一区有黄有色的免费视频| 在线精品无人区一区二区三 | 亚洲av中文字字幕乱码综合| 免费看光身美女| 国产极品天堂在线| 在线观看三级黄色| 一区二区三区乱码不卡18| 久久精品久久久久久久性| 在线观看一区二区三区激情| av.在线天堂| 99久久中文字幕三级久久日本| 成人国产av品久久久| 免费播放大片免费观看视频在线观看| 国产av国产精品国产| 亚洲国产成人一精品久久久| 女人久久www免费人成看片| 久久国产乱子免费精品| 亚洲精品成人久久久久久| 少妇猛男粗大的猛烈进出视频 | 国产精品久久久久久精品电影| 日韩欧美精品免费久久| 国产高清三级在线| 精品国产露脸久久av麻豆| 高清日韩中文字幕在线| 日韩免费高清中文字幕av| 欧美一级a爱片免费观看看| 亚洲av国产av综合av卡| 欧美性猛交╳xxx乱大交人| 中文字幕免费在线视频6| 老司机影院毛片| 国产亚洲91精品色在线| 性色av一级| 欧美激情久久久久久爽电影| 国产一区二区在线观看日韩| 日本av手机在线免费观看| 蜜桃久久精品国产亚洲av| 美女脱内裤让男人舔精品视频| 菩萨蛮人人尽说江南好唐韦庄| 色哟哟·www| 亚洲国产色片| 涩涩av久久男人的天堂| 欧美精品人与动牲交sv欧美| 夫妻午夜视频| 久久韩国三级中文字幕| 亚洲久久久久久中文字幕| 免费看光身美女| 亚洲精品日韩av片在线观看| 99九九线精品视频在线观看视频| 在线免费观看不下载黄p国产| 免费人成在线观看视频色| 欧美精品一区二区大全| 少妇人妻 视频| 丝袜美腿在线中文| 国产乱人偷精品视频| 波野结衣二区三区在线| 亚洲国产av新网站| 王馨瑶露胸无遮挡在线观看| 亚洲精品456在线播放app| 高清视频免费观看一区二区| 亚洲精品视频女| 男女啪啪激烈高潮av片| eeuss影院久久| 色哟哟·www| av线在线观看网站| 熟女av电影| 亚洲精品aⅴ在线观看| 2018国产大陆天天弄谢| 一级毛片我不卡| 国产一级毛片在线| 免费少妇av软件| 亚洲精品自拍成人| 少妇裸体淫交视频免费看高清| 久久精品国产亚洲av天美| 国产熟女欧美一区二区| 18禁裸乳无遮挡动漫免费视频 | 日韩电影二区| 国产免费一级a男人的天堂| 久久精品国产自在天天线| 禁无遮挡网站| 免费人成在线观看视频色| 国产精品成人在线| 观看美女的网站| 99热这里只有精品一区| 国产高清不卡午夜福利| 黑人高潮一二区| 国产亚洲精品久久久com| 国产一区有黄有色的免费视频| 边亲边吃奶的免费视频| 欧美另类一区| 国产v大片淫在线免费观看| 欧美xxⅹ黑人| 亚洲av电影在线观看一区二区三区 | 欧美xxxx性猛交bbbb| 精品久久久久久久人妻蜜臀av| 久久99蜜桃精品久久| 精品少妇久久久久久888优播| 在线观看一区二区三区激情| 22中文网久久字幕| 男女边吃奶边做爰视频| 免费看av在线观看网站| 五月开心婷婷网| 亚洲自拍偷在线| 男男h啪啪无遮挡| 国产久久久一区二区三区| 国产精品国产av在线观看| 久久99精品国语久久久| 亚洲精品亚洲一区二区| 一区二区av电影网| 毛片女人毛片| kizo精华| 成人国产麻豆网| 久久久色成人| 欧美人与善性xxx| 日韩制服骚丝袜av| 日韩欧美精品免费久久| 一级爰片在线观看| 久久久国产一区二区| 春色校园在线视频观看| 国产91av在线免费观看| av在线蜜桃| 久久久精品94久久精品| 一级av片app| 老司机影院毛片| 少妇的逼好多水| 麻豆乱淫一区二区| 免费大片黄手机在线观看| 99久久精品热视频| 成人无遮挡网站| videos熟女内射| 边亲边吃奶的免费视频| 久久影院123| 亚洲人成网站在线播| 欧美xxⅹ黑人| a级毛片免费高清观看在线播放| 久久久久网色| av在线亚洲专区| 亚洲av福利一区| 又爽又黄无遮挡网站| 国产成人精品久久久久久| 97在线视频观看| 晚上一个人看的免费电影| 日韩在线高清观看一区二区三区| 欧美bdsm另类| 禁无遮挡网站| 免费av观看视频| 日韩中字成人| 国产精品嫩草影院av在线观看| 91精品一卡2卡3卡4卡| 日本色播在线视频| 美女内射精品一级片tv| 亚洲不卡免费看| 亚洲av二区三区四区| 中国三级夫妇交换| av女优亚洲男人天堂| 亚洲精华国产精华液的使用体验| 国内揄拍国产精品人妻在线| 下体分泌物呈黄色| 亚洲性久久影院| 日日啪夜夜撸| 春色校园在线视频观看| 在线播放无遮挡| 少妇丰满av| 日韩欧美精品免费久久| 免费av观看视频| 在线观看一区二区三区| 亚洲最大成人手机在线| 国产精品人妻久久久久久| 日韩欧美精品v在线| 亚洲av一区综合| 午夜精品一区二区三区免费看| 亚洲天堂av无毛| 午夜福利视频1000在线观看| 久久亚洲国产成人精品v| 免费av毛片视频| 午夜激情福利司机影院| 免费大片黄手机在线观看| 高清日韩中文字幕在线| 免费观看的影片在线观看| 亚洲天堂国产精品一区在线| av天堂中文字幕网| 中文精品一卡2卡3卡4更新| 性色av一级| 免费大片黄手机在线观看| 少妇人妻久久综合中文| 亚洲人与动物交配视频| 国产精品.久久久| 少妇的逼好多水| 欧美3d第一页| 永久免费av网站大全| 午夜激情久久久久久久| 国产精品蜜桃在线观看| 精华霜和精华液先用哪个| a级毛色黄片| 久久久午夜欧美精品| 成年人午夜在线观看视频| 天堂中文最新版在线下载 | 黄色怎么调成土黄色| 亚洲精华国产精华液的使用体验| 国产精品一及| 亚洲欧美中文字幕日韩二区| 亚洲精品aⅴ在线观看| 久久久久网色| 日韩av不卡免费在线播放| 国产精品一区www在线观看| 国产免费一级a男人的天堂| 久久久色成人| 国产永久视频网站| 少妇 在线观看| 男女边摸边吃奶| 午夜爱爱视频在线播放| 日韩av在线免费看完整版不卡| 国产美女午夜福利| 欧美亚洲 丝袜 人妻 在线| 狂野欧美白嫩少妇大欣赏| 国产成人精品久久久久久| 国内精品宾馆在线| 国内少妇人妻偷人精品xxx网站| 老司机影院成人| 制服丝袜香蕉在线| 丝袜脚勾引网站| 91久久精品电影网| 麻豆乱淫一区二区| 免费观看a级毛片全部| 美女视频免费永久观看网站| 成人国产麻豆网| 亚洲不卡免费看| 久久99热这里只有精品18| 高清毛片免费看| 菩萨蛮人人尽说江南好唐韦庄| 亚洲精品乱码久久久v下载方式| 搡女人真爽免费视频火全软件| 日本三级黄在线观看| 少妇人妻久久综合中文| 成人美女网站在线观看视频| 国产黄a三级三级三级人| 亚洲丝袜综合中文字幕| 舔av片在线| 国产成人一区二区在线| 嫩草影院新地址| 免费大片黄手机在线观看| 午夜福利视频精品| 在线a可以看的网站| 成人美女网站在线观看视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | av播播在线观看一区| 女人被狂操c到高潮| 亚洲天堂av无毛| 久久久久久久国产电影| 我的女老师完整版在线观看| 日韩视频在线欧美| 狂野欧美激情性xxxx在线观看| 在线播放无遮挡| 在线观看美女被高潮喷水网站| 日韩国内少妇激情av| 国产免费视频播放在线视频| 在线天堂最新版资源| 最近手机中文字幕大全| 欧美精品人与动牲交sv欧美| 国产精品一二三区在线看| 国产一区二区三区av在线| 国产精品久久久久久精品古装| 少妇人妻一区二区三区视频| 日本三级黄在线观看| 一级毛片 在线播放| 国产成人a区在线观看| 各种免费的搞黄视频| 亚洲精品乱码久久久v下载方式| 日产精品乱码卡一卡2卡三| 水蜜桃什么品种好| 精品一区二区三区视频在线| 男人爽女人下面视频在线观看| 美女脱内裤让男人舔精品视频| 噜噜噜噜噜久久久久久91| 国产 一区 欧美 日韩| 日韩精品有码人妻一区| 性色av一级| 亚洲精品久久久久久婷婷小说| 国产精品久久久久久精品电影小说 | 高清视频免费观看一区二区| 久久久久久国产a免费观看| 高清在线视频一区二区三区| 婷婷色综合大香蕉| 在线免费观看不下载黄p国产| 精品熟女少妇av免费看| 啦啦啦中文免费视频观看日本| av在线亚洲专区| 亚洲aⅴ乱码一区二区在线播放| 亚洲成人av在线免费| 国产精品嫩草影院av在线观看| 国内精品美女久久久久久| 一级毛片我不卡| 国产大屁股一区二区在线视频| 日韩视频在线欧美| 国产午夜精品久久久久久一区二区三区| 一本久久精品| 99久久九九国产精品国产免费| 色视频在线一区二区三区| 99久久精品一区二区三区| 各种免费的搞黄视频| 久久精品国产亚洲av天美| av免费观看日本| 国产69精品久久久久777片| 欧美性猛交╳xxx乱大交人| 九九在线视频观看精品| 成人特级av手机在线观看| 久久久久久久久久久丰满| 亚洲欧美日韩无卡精品| 中文乱码字字幕精品一区二区三区| 人人妻人人澡人人爽人人夜夜| 97在线视频观看| 中文在线观看免费www的网站| 欧美日韩在线观看h| 成年人午夜在线观看视频| 熟女电影av网| 一区二区三区免费毛片| 麻豆久久精品国产亚洲av| 国产乱来视频区| 日韩精品有码人妻一区| 一级爰片在线观看| 深爱激情五月婷婷| 国产精品一区二区在线观看99| 2021少妇久久久久久久久久久| 久久热精品热| 国产欧美日韩一区二区三区在线 | www.色视频.com| 天堂俺去俺来也www色官网| 国产老妇女一区| 我要看日韩黄色一级片| 乱码一卡2卡4卡精品| 亚洲美女视频黄频| 又爽又黄无遮挡网站| 身体一侧抽搐| av在线天堂中文字幕| 美女视频免费永久观看网站| 国内精品美女久久久久久| 免费观看性生交大片5| 少妇的逼水好多| 爱豆传媒免费全集在线观看| 九九爱精品视频在线观看| 免费观看a级毛片全部| a级毛色黄片| 精品国产三级普通话版| 亚洲三级黄色毛片| 校园人妻丝袜中文字幕| 精品熟女少妇av免费看| 啦啦啦啦在线视频资源| 18禁裸乳无遮挡动漫免费视频 | 深爱激情五月婷婷| 美女脱内裤让男人舔精品视频| 51国产日韩欧美| 欧美日韩精品成人综合77777| 日本一本二区三区精品| 在线观看美女被高潮喷水网站| 人体艺术视频欧美日本| 亚洲精品第二区| eeuss影院久久| 夜夜爽夜夜爽视频| 日本av手机在线免费观看| 肉色欧美久久久久久久蜜桃 | 深夜a级毛片| 美女视频免费永久观看网站| 少妇丰满av| 日韩电影二区| 日本av手机在线免费观看| 亚洲电影在线观看av| 午夜爱爱视频在线播放| 麻豆久久精品国产亚洲av| 噜噜噜噜噜久久久久久91| av在线老鸭窝| 草草在线视频免费看| 性色avwww在线观看| 80岁老熟妇乱子伦牲交| av又黄又爽大尺度在线免费看| 亚洲久久久久久中文字幕| 天天躁夜夜躁狠狠久久av| 一个人观看的视频www高清免费观看| 免费av观看视频| 有码 亚洲区| 丝袜喷水一区| 亚洲成人精品中文字幕电影| 能在线免费看毛片的网站| 91狼人影院| 听说在线观看完整版免费高清| 久久这里有精品视频免费| 国产一区亚洲一区在线观看| videos熟女内射| 免费av毛片视频| 好男人视频免费观看在线| 日韩国内少妇激情av| 欧美激情在线99| 精品酒店卫生间| 女人十人毛片免费观看3o分钟| 国产精品不卡视频一区二区| 午夜激情久久久久久久| 国产一区有黄有色的免费视频| 七月丁香在线播放| 久热久热在线精品观看| 国产黄频视频在线观看| 国产亚洲一区二区精品| 日韩欧美精品v在线| 日本黄色片子视频| 国产精品国产三级专区第一集| 啦啦啦啦在线视频资源| 久久99蜜桃精品久久| 美女视频免费永久观看网站| 少妇丰满av| 99久久精品一区二区三区| 国产精品久久久久久av不卡| 亚洲av国产av综合av卡| 在线观看人妻少妇| 色视频在线一区二区三区| 国产精品一区二区性色av| 国产欧美日韩一区二区三区在线 | 免费人成在线观看视频色| 色播亚洲综合网| 久久久久久九九精品二区国产| 亚洲不卡免费看| 久久ye,这里只有精品| 一级毛片黄色毛片免费观看视频| 成人鲁丝片一二三区免费| 国产男人的电影天堂91| 亚洲国产日韩一区二区| 国产成人aa在线观看| 免费av毛片视频| 人人妻人人看人人澡| 精品少妇久久久久久888优播| 精品酒店卫生间| 欧美xxxx黑人xx丫x性爽| 女人久久www免费人成看片| 69av精品久久久久久| 国产精品99久久久久久久久| 国产亚洲午夜精品一区二区久久 | av国产免费在线观看| 中文字幕亚洲精品专区| av国产免费在线观看| 午夜精品国产一区二区电影 | 69av精品久久久久久| h日本视频在线播放| 秋霞在线观看毛片| 各种免费的搞黄视频| 黄色配什么色好看| 欧美日韩一区二区视频在线观看视频在线 | 一级爰片在线观看| av在线播放精品| 久久久久久九九精品二区国产| 尤物成人国产欧美一区二区三区| 青青草视频在线视频观看| 边亲边吃奶的免费视频| 日日啪夜夜爽| 深爱激情五月婷婷| 欧美精品一区二区大全| 欧美成人午夜免费资源| 在线观看国产h片| 国产免费一级a男人的天堂| 99久久人妻综合| 青青草视频在线视频观看| 亚洲av欧美aⅴ国产| videos熟女内射| 一区二区三区四区激情视频| 国产永久视频网站| 丝袜喷水一区| 少妇猛男粗大的猛烈进出视频 | 97在线人人人人妻|