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

    基于北斗/iNEMO慣性模塊的機(jī)器人運(yùn)動姿態(tài)解算方法

    2017-03-21 12:53:20王廣龍薄洪波
    計算機(jī)測量與控制 2017年2期
    關(guān)鍵詞:磁強(qiáng)計對準(zhǔn)陀螺儀

    陳 凱,王廣龍,王 成,薄洪波

    (1.66008部隊,天津 300250;2.軍械工程學(xué)院 納米技術(shù)與微系統(tǒng)實(shí)驗室,石家莊 050003)

    基于北斗/iNEMO慣性模塊的機(jī)器人運(yùn)動姿態(tài)解算方法

    陳 凱1,王廣龍2,王 成2,薄洪波1

    (1.66008部隊,天津 300250;2.軍械工程學(xué)院 納米技術(shù)與微系統(tǒng)實(shí)驗室,石家莊 050003)

    基于北斗/iNEMO慣性模塊研究了一種組合式機(jī)器人運(yùn)動姿態(tài)測量系統(tǒng)和解算方法;設(shè)計了以ARM為核心的嵌入式組合姿態(tài)解算平臺,通過四元數(shù)法和卡爾曼濾波技術(shù),對iNEMO慣性組件測量數(shù)據(jù)進(jìn)行融合,進(jìn)而在高動態(tài)環(huán)境中實(shí)時解算出機(jī)器人的運(yùn)動姿態(tài);同時,以北斗接收機(jī)輸出的1PPS上升沿脈沖作為數(shù)據(jù)融合的同步信號,并利用其輸出的導(dǎo)航信息輔助iNEMO慣性模塊實(shí)現(xiàn)精對準(zhǔn);測試結(jié)果驗證了設(shè)計方案的可行性和正確性,為機(jī)器人姿態(tài)解算提供了一種高精度、高穩(wěn)定性、小體積、低功耗的解決方案。

    姿態(tài)解算;機(jī)器人;慣性組件;北斗系統(tǒng);初始對準(zhǔn)

    0 引言

    移動機(jī)器人是一個集環(huán)境感知、動態(tài)決策與規(guī)劃、行為控制與執(zhí)行等多功能于一體的綜合系統(tǒng),是目前科學(xué)技術(shù)發(fā)展最活躍的領(lǐng)域之一[1]。它可以協(xié)助或取代人類的工作,已經(jīng)廣泛地應(yīng)用于地質(zhì)勘查、軍事偵察、災(zāi)難救援等領(lǐng)域[2]。姿態(tài)檢測作為機(jī)器人研究領(lǐng)域一個重要的組成部分,對于評定機(jī)器人的運(yùn)動特性和姿態(tài)反饋控制具有非常重要的意義[3]。

    目前,機(jī)器人的姿態(tài)檢測方法按照接觸方式分為間接檢測法和直接檢測法兩種。間接檢測法主要是使用超聲波、激光或者CCD攝像機(jī)進(jìn)行姿態(tài)測量,需要在機(jī)器人上安裝多個外部檢測元件,結(jié)構(gòu)復(fù)雜,成本較高,不適合微小型機(jī)器人。直接檢測法主要是利用陀螺儀、加速度計以及磁強(qiáng)計等慣性敏感元器件和計算機(jī),實(shí)時測量運(yùn)載體的空間位置、姿態(tài)和重力場等參數(shù),慣性測量的好處是響應(yīng)速度快、成本低,隨著MEMS技術(shù)的發(fā)展和對功耗、體積的要求,低成本的微型慣性組件已經(jīng)大量的在機(jī)器人的姿態(tài)測量中使用[4-6]。慣性系統(tǒng)在姿態(tài)解算前必須對其進(jìn)行初始對準(zhǔn),對準(zhǔn)的計算值將作為姿態(tài)解算的初始值,因此它的精度將直接影響姿態(tài)解算的精度。但低成本的慣性器件由于誤差較大,即使在靜基座下也無法準(zhǔn)確辨識地球自轉(zhuǎn)速率等重要信息,需要使用外部信息輔助動態(tài)對準(zhǔn)[7]。

    為解決對機(jī)器人運(yùn)動姿態(tài)進(jìn)行動態(tài)檢測的問題,本文基于北斗/iNEMO慣性模塊,設(shè)計了以ARM為核心的嵌入式組合姿態(tài)解算平臺。在此基礎(chǔ)上,對實(shí)時姿態(tài)解算和初始對準(zhǔn)的算法進(jìn)行了研究。最后通過試驗,驗證了硬件和算法設(shè)計的有效性。

    1 硬件系統(tǒng)總體設(shè)計

    硬件采用ARM為核心的組合姿態(tài)解算平臺,可采集iNEMO慣性模塊輸出量,接收北斗模塊輸出信號,完成姿態(tài)數(shù)據(jù)的解算、融合處理和存儲,系統(tǒng)結(jié)構(gòu)如圖1所示。其中,iNEMO慣性模塊采用的是ST(意法半導(dǎo)體)公司的LSM9DS1芯片,這最一種基于MEMS工藝的9軸運(yùn)動位置傳感器模塊,內(nèi)置一個3軸加速度計、一個3軸陀螺儀和一個3軸磁強(qiáng)計。北斗接收機(jī)采用的和芯星通公司的UM220-III N芯片,是一種雙系統(tǒng)高性能的GNSS 模塊,能夠同時支持 BD2 B1、GPS L1 兩個頻點(diǎn)。主處理器選用的是ST公司的ARM芯片STM32F405,有較強(qiáng)的數(shù)字信號處理能力和多任務(wù)管理能力,有較多的應(yīng)用接口。為了避免處理頻繁的響應(yīng)中斷,協(xié)處理單元由ARM芯片STM32F103和FPGA芯片XC6SLX75共同構(gòu)成。整個系統(tǒng)以無線指令模塊的高電平輸出作為測試的啟動指令。測試開始后,協(xié)處理器ARM從北斗接收機(jī)輸出的NEMA0813信息中提取姿態(tài)解算所需的UTC時間、位置、速度、定位狀態(tài)等信息并暫存入FPGA中;FPGA內(nèi)包含一個雙口RAM存儲器和一個高精度的計數(shù)器,其中計數(shù)器以北斗接收機(jī)有效輸出時的1PPS上升沿脈沖作為清零信號;主處理器ARM根據(jù)FPGA輸出的上升沿信號產(chǎn)生中斷,通過FSMC總線、兩路SPI總線分別讀取FPGA內(nèi)的北斗導(dǎo)航信息、iNEMO慣性模塊的加速度/角速度和地磁信息,利用這些信息完成姿態(tài)數(shù)據(jù)的解算,將結(jié)果存入FLASH芯片中。同時,為防止出現(xiàn)時序的競爭和冒險,并協(xié)調(diào)各個器件的工作,系統(tǒng)以有源晶振的輸出作為唯一的頻率基準(zhǔn)。

    圖1 硬件系統(tǒng)總體框圖

    2 實(shí)時姿態(tài)解算方法

    由于iNEMO慣性模塊的輸出數(shù)據(jù)速率遠(yuǎn)高于北斗接收機(jī)數(shù)據(jù)更新率,因此機(jī)器人姿態(tài)的實(shí)時解算主要針對iNEMO慣性模塊(陀螺儀、加速度計、磁強(qiáng)計)的組合系統(tǒng)。其中,陀螺儀、加速度計、磁強(qiáng)計都能檢測到載體的姿態(tài)角,但由于器件自身的局限性和外界環(huán)境的干擾,單一慣性器件無法實(shí)現(xiàn)對載體姿態(tài)角準(zhǔn)確而完整的檢測,因此必須將這幾種慣性傳感器的測量數(shù)據(jù)進(jìn)行融合。本文實(shí)時姿態(tài)解算的總體思想是,由陀螺儀單獨(dú)定姿、加速度計和磁強(qiáng)計組合定姿分別解算出姿態(tài)角,采用一維卡爾曼濾波對兩類輸出的三個方位的姿態(tài)角數(shù)據(jù)分別進(jìn)行互補(bǔ)融合,進(jìn)而得到更精確的姿態(tài)角,如圖2所示。

    圖2 實(shí)時姿態(tài)解算原理

    2.1 iNEMO陀螺儀解算歐拉角

    (1)

    式中,ψ∈[-π,π]為方位角, θ∈[-π/2,π/2]為俯仰角,γ∈[-π,π]為橫滾角,這三個角稱為載體的歐拉角,用于表示載體的運(yùn)動姿態(tài),此處采用四元數(shù)法求解歐拉角。

    設(shè)載體坐標(biāo)系相對于導(dǎo)航坐標(biāo)系的轉(zhuǎn)動四元數(shù)為:

    (2)

    iNEMO陀螺儀x、y、z軸的角速度輸出為(ωx、ωy、ωz),采樣周期為T,采用一階龍格-庫塔法更新四元數(shù):

    (3)

    (4)

    若用四元數(shù)Q表示載體坐標(biāo)系相對于導(dǎo)航坐標(biāo)系的變換,即:

    (5)

    (6)

    2.2iNEMO加速度計、磁強(qiáng)計解算歐拉角

    iNEMO加速度計x、y、z軸的線性加速度輸出為(ax、ay、az),則可估算俯仰角θ2、橫滾角γ2:

    (7)

    (8)

    iNEMO磁強(qiáng)計x、y、z軸的磁場強(qiáng)度輸出為(mx、my、mz),水平狀態(tài)下可以估算出方位角的主值:

    (9)

    機(jī)器人在運(yùn)動過程中由于顛簸或者上下坡度而不會總處于水平位置上,此處采用幾何旋轉(zhuǎn)法,通過姿態(tài)矩陣建立載體坐標(biāo)系下iNEMO慣性模塊測得的重力矢量和磁場強(qiáng)度矢量的對應(yīng)方程,用加速計估算出的俯仰角θ、橫滾角γ的值對磁強(qiáng)計輸出的進(jìn)行傾角補(bǔ)償:

    (10)

    (11)

    由于地球的地理北向與磁北方向之間存在一個磁偏角,經(jīng)過磁偏角修正,就可以得到載體真實(shí)的方位角ψ2。

    iNEMO加速度計和磁強(qiáng)度計的估算值都通過泰勒級數(shù)的展開完成解算。

    2.3 卡爾曼濾波的姿態(tài)數(shù)據(jù)融合

    由于陀螺儀存在溫度漂移,積分運(yùn)算的累積誤差會迅速放大,因此陀螺儀解算出的歐拉角必須進(jìn)一步融合加速度計和磁強(qiáng)計的數(shù)據(jù)。以方位角為例, 陀螺儀解算出ψ1、加速度和磁強(qiáng)計組合解算出ψ2,對應(yīng)的方差δ1和δ2。為減小計算量,采用一維卡爾曼(Kalman)濾波器,方位角的最佳估計值ψ為[8]:

    (12)

    其方差為:

    (13)

    3 北斗輔助iNEMO慣性模塊精對準(zhǔn)

    機(jī)器人姿態(tài)解算時,首先要對系統(tǒng)進(jìn)行初始化,完成測量平臺的導(dǎo)航數(shù)據(jù)初始對準(zhǔn)工作,消除慣性模塊的初始誤差,并為姿態(tài)解算賦予初始值,為解算的迭代過程打下基礎(chǔ)。初始對準(zhǔn)即在最短的時間內(nèi)以一定的精度確定從載體坐標(biāo)系到導(dǎo)航坐標(biāo)系的初始姿態(tài)變換矩陣[13]。

    對準(zhǔn)分兩步進(jìn)行,粗對準(zhǔn)和精對準(zhǔn)。在粗對準(zhǔn)階段,將iNEMO加速度計、磁強(qiáng)計輸出估算的歐拉角代入方程(1),即可得到粗對準(zhǔn)階段的初始姿態(tài)矩陣。粗對準(zhǔn)可以縮短對準(zhǔn)的時間,但對準(zhǔn)的精度較低,需要精對準(zhǔn)進(jìn)一步提高對準(zhǔn)精度。在精對準(zhǔn)階段,本文采用組合對準(zhǔn)的方式,應(yīng)用拓展卡爾曼濾波方法,從慣導(dǎo)參數(shù)相對于北斗導(dǎo)航系統(tǒng)提供的參數(shù)的偏差中估計出慣性模塊的誤差量,并進(jìn)行校正,用實(shí)時更新的姿態(tài)矩陣替代粗對準(zhǔn)階段得到的姿態(tài)矩陣。

    為有效地提高計算精度和效率,此處直接采用北斗系統(tǒng)參照的地球坐標(biāo)系(e系)實(shí)現(xiàn)對準(zhǔn)。將北斗的導(dǎo)航解直接作為卡爾曼濾波器的輸入,這是一種松組合的數(shù)據(jù)融合方法,因而其數(shù)學(xué)模型并不需要包含北斗狀態(tài)方程。由三軸位置誤差δre、速度誤差δve、姿態(tài)誤差φ,慣性模塊三軸加速度計偏置誤差▽b、陀螺儀漂移εb共同組成15維的系統(tǒng)狀態(tài)向量:

    (16)

    設(shè)理想的地球坐標(biāo)系為e系(簡稱地球系),計算的地球坐標(biāo)系為c系(簡稱計算系)。地球系的位置誤差由于不涉及坐標(biāo)系統(tǒng)轉(zhuǎn)換,其微分形式即是速度誤差,比較簡單。但慣性元件存在測量誤差會導(dǎo)致兩個坐標(biāo)系之間存在一個轉(zhuǎn)動誤差,即平臺歐拉誤差角,需要予以考慮。短時間內(nèi),三軸陀螺儀漂移、加速度計誤差可認(rèn)為為常值,不隨時間變化,其微分為0。在地球坐標(biāo)系下,任意失準(zhǔn)角情況的系統(tǒng)狀態(tài)誤差微分形式為[7]:

    (17)

    則系統(tǒng)的狀態(tài)方程可寫為:

    (18)

    式中,F(xiàn)(15×15)為系統(tǒng)的狀態(tài)轉(zhuǎn)移矩陣,G為系統(tǒng)的動態(tài)噪聲矩陣,W為系統(tǒng)的過程白噪聲矢量。

    在慣性模塊與北斗的松組合模式中,一般取速度和位置為觀測信息。但在靜止與勻速運(yùn)動狀態(tài)下,系統(tǒng)的方位角誤差不可觀測。因此,本系統(tǒng)在選取北斗和慣性模塊中輸出的位置和速度之差的基礎(chǔ)上,將慣性模塊中陀螺儀與磁強(qiáng)計輸出的方位角之差作為觀測量,構(gòu)造量測方程:

    (19)

    式中,rI和vI為慣性模塊輸出的載體3維速度和位置信息,rB、vB為北斗接收機(jī)輸出的載體3維速度和位置信息,ψ1、ψ2分別為慣性模塊中陀螺儀與磁強(qiáng)計輸出的方位角。H為量測轉(zhuǎn)移矩陣,V為量測白噪聲矩陣。

    將狀態(tài)方程(18)和量測方程(19)離散化, 即可利用拓展卡爾曼濾波,通過時間更新和狀態(tài)更新,反饋修正姿態(tài)信息,完成姿態(tài)數(shù)據(jù)的初始對準(zhǔn)。

    4 試驗測試及結(jié)果分析

    4.1 基于Kalman濾波的姿態(tài)信息融合數(shù)據(jù)分析

    使用微系統(tǒng)中的iNEMO慣性模塊進(jìn)行實(shí)時姿態(tài)測量實(shí)驗,陀螺儀和加速度計采樣頻率為100 Hz,磁強(qiáng)計的采用頻率為10 Hz,將兩種姿態(tài)測量結(jié)果進(jìn)行Kalman濾波信息融合,得到的姿態(tài)信息如圖3所示。其中,上圖為陀螺儀輸出的歐拉角,中圖為加速度計和磁強(qiáng)計組合輸出的歐拉角,下圖為運(yùn)用經(jīng)過Kalman濾波后輸出的歐拉角,γ、θ、ψ分別表示系統(tǒng)的橫滾角、俯仰角和方位角。由圖中可以看出,經(jīng)過Kalman濾波信息融合,對兩種方式測得的姿態(tài)誤差進(jìn)行了補(bǔ)償,經(jīng)濾波得到的結(jié)果較為平滑且與系統(tǒng)真實(shí)的姿態(tài)變化相吻合。

    圖3 實(shí)時姿態(tài)解算Kalman濾波結(jié)果

    4.2 系統(tǒng)初始對準(zhǔn)精度分析

    初始對準(zhǔn)在靜基座條件下進(jìn)行。實(shí)驗中,慣組模塊陀螺儀和加速度計采樣頻率為100 Hz、磁強(qiáng)計的采用頻率為10 Hz,北斗接收機(jī)的數(shù)據(jù)更新率為1 Hz。采用慣組輸出進(jìn)行粗對準(zhǔn),得到姿態(tài)角誤差如圖4所示。北斗輔助iNEMO慣性模塊進(jìn)行精對準(zhǔn),得到姿態(tài)角誤差如圖5所示。表1列出了兩種對準(zhǔn)方式下的姿態(tài)角誤差的均值和標(biāo)準(zhǔn)差。從測試結(jié)果可以看出,粗對準(zhǔn)中俯仰角誤差和方位角誤差較大,誤差均值在0.7°左右,精對準(zhǔn)后,輸出的姿態(tài)角誤差的均值和標(biāo)準(zhǔn)差均變小,得到的

    圖4 粗對準(zhǔn)時靜基座的姿態(tài)角誤差

    圖5 精對準(zhǔn)時靜基座的姿態(tài)角誤差

    方位角和俯仰角的誤差控制在±0.1°以內(nèi),橫滾角誤差在±0.15°以內(nèi),表明基于上述模型的初始對準(zhǔn)很好的滿足了姿態(tài)解算精度的要求。

    表1 靜態(tài)測試時的姿態(tài)角誤差結(jié)果

    5 結(jié)論

    本文設(shè)計了一種移動機(jī)器人的姿態(tài)測量系統(tǒng),可以完成對機(jī)器人運(yùn)動姿態(tài)的實(shí)時采集和存儲。研究了對iNEMO慣性模塊測量數(shù)據(jù)進(jìn)行融合的方法,可以較精確地實(shí)時解算出姿態(tài)數(shù)據(jù)。利用北斗輸出的導(dǎo)航信息輔助iNEMO慣性模塊實(shí)現(xiàn)精對準(zhǔn),能夠有效的改善系統(tǒng)的性能。軟硬件的測試結(jié)果驗證了設(shè)計方案的可行性和正確性,該方法方便快捷、功耗低、試驗范圍不受限制,并且用來支持試驗的地面設(shè)備少,節(jié)約經(jīng)費(fèi),設(shè)備可重復(fù)使用,能夠滿足機(jī)器人姿態(tài)解算所要求的高精度、穩(wěn)定性,具有很強(qiáng)的通用性和可拓展性。

    [1] 徐國保,尹怡欣, 周美娟.智能移動機(jī)器人技術(shù)現(xiàn)狀及展望[J]. 機(jī)器人技術(shù)與應(yīng)用,2007(2):29-34.

    [2] Hua M D, Ducard G, Hamel T, et al. Introduction to Nonlinear Attitude Estimation for Aerial Robotic Systems[J]. Aerospace Lab, 2014: AL08-04.

    [3] 任文強(qiáng).多面體機(jī)器人姿態(tài)檢測技術(shù)研究[D].北京:北京交通大學(xué),2016.

    [4] 秦 勇,臧希喆,王曉宇等.基于MEMS 慣性傳感器的機(jī)器人姿態(tài)檢測系統(tǒng)的研究[J].傳感技術(shù)學(xué)報,2007,20(2):298-301.

    [5] Renaudin V, Combettes C. Magnetic, Acceleration Fields and Gyroscope Quaternion (MAGYQ)-Based Attitude Estimation with Smartphone Sensors for Indoor Pedestrian Navigation[J]. Sensors, 2014, 14(12): 22864-22890.

    [6] Di L, Fromm T, Chen Y Q. A Data Fusion System for Attitude Estimation of Low-cost Miniature UAVs[J]. Journal of Intelligent & Robotic Systems, 2012,65(1-4): 621-635.

    [7] 張秋昭, 張書畢, 王 堅,等.露天礦卡車低成本GPS/INS組合導(dǎo)航系統(tǒng)動態(tài)對準(zhǔn)模型[J].煤炭學(xué)報,2013,38(8):1362-1367.

    [8] 姜 蘭.超長航時無人機(jī)持久組合導(dǎo)航系統(tǒng)設(shè)計[D].南京:南京理工大學(xué),2014.

    [9] 楊亞非,譚久彬,鄧正隆.慣導(dǎo)系統(tǒng)初始對準(zhǔn)技術(shù)綜述[J]. 中國慣性技術(shù)學(xué)報,2002,10(2):68-72.

    Method of Robot Attitude Calculation Based on Beidou/ iNEMO Inertial Module

    Chen Kai1, Wang Guanglong2,Wang Cheng2,Bo Hongbo1

    (1.No.66008 Unit of PLA, Tianjin 300250, China; 2.Laboratory of Nanotechnology and Microsystems,Ordnance Engineering College, Shijiazhuang 050003, China)

    This paper researches a combined robot posture measurement system and attitude algorithm based on Beidou/iNEMO. Here designed an embedded combination attitude calculation platform cored by ARM, to calculate posture of robot in real-time under high dynamic condition, by Quaternion and Kalman filter technology fusing measured data of iNEMO inertial measurement unit. Meanwhile, using the Beidou receiver outputted 1pps rising edge pulse as synchronous time of subsequent data fusion, as well as its outputted navigation data aid iNEMO inertial measurement unit, fine alignment can be realized. The feasibility and correctness of the designed test has been verified, it provides robot attitude calculation with high precision, high stability, small volume, low power consumption.

    attitude calculation; inertial measurement unit(IMU); robot; BeiDou system; initial alignment

    2016-04-07;

    2016-06-21。

    國家自然科學(xué)青年基金(61501493)。

    陳 凱(1983-),男,重慶榮昌人,碩士研究生,工程師,主要從事測試計量技術(shù)及儀器的應(yīng)用研究。

    1671-4598(2017)02-0130-04

    10.16526/j.cnki.11-4762/tp.2017.02.036

    TP212.9

    A

    猜你喜歡
    磁強(qiáng)計對準(zhǔn)陀螺儀
    磁強(qiáng)計陣列測量一致性校正
    基于EMD的MEMS陀螺儀隨機(jī)漂移分析方法
    基于矢量磁強(qiáng)計的磁場梯度張量儀誤差校正方法
    組合導(dǎo)航中磁強(qiáng)計干擾估計與補(bǔ)償方法
    基于LabVIEW的微型磁通門磁強(qiáng)計測試系統(tǒng)搭建
    對準(zhǔn)提升組織力的聚焦點(diǎn)——陜西以組織振興引領(lǐng)鄉(xiāng)村振興
    我國著名陀螺儀專家——林士諤
    一種改進(jìn)的速度加姿態(tài)匹配快速傳遞對準(zhǔn)算法
    微機(jī)械陀螺儀概述和發(fā)展
    MEMS三軸陀螺儀中不匹配干擾抑制方法
    1024香蕉在线观看| 亚洲欧美色中文字幕在线| 国产野战对白在线观看| 女人久久www免费人成看片| 永久网站在线| 日韩中文字幕视频在线看片| 美女xxoo啪啪120秒动态图| 国产精品嫩草影院av在线观看| 久久精品国产亚洲av天美| 国产不卡av网站在线观看| 纵有疾风起免费观看全集完整版| 一级片'在线观看视频| 一区二区三区四区激情视频| 欧美国产精品va在线观看不卡| 成年av动漫网址| 91精品三级在线观看| 这个男人来自地球电影免费观看 | 成人18禁高潮啪啪吃奶动态图| 亚洲欧美清纯卡通| 亚洲av在线观看美女高潮| 91精品三级在线观看| 大香蕉久久成人网| 中文字幕最新亚洲高清| 日本爱情动作片www.在线观看| 久久久国产一区二区| 中文字幕人妻丝袜制服| 国产成人精品在线电影| 免费女性裸体啪啪无遮挡网站| 天堂中文最新版在线下载| 国产成人精品久久二区二区91 | 免费日韩欧美在线观看| 亚洲欧美清纯卡通| 超碰成人久久| 少妇熟女欧美另类| 黄色 视频免费看| 国产乱人偷精品视频| 日韩一区二区视频免费看| 亚洲国产日韩一区二区| 免费观看在线日韩| 只有这里有精品99| 一级毛片 在线播放| 久久久精品94久久精品| av国产精品久久久久影院| 我要看黄色一级片免费的| 曰老女人黄片| h视频一区二区三区| 国产女主播在线喷水免费视频网站| 国产一区亚洲一区在线观看| 丝袜脚勾引网站| 日韩视频在线欧美| 亚洲精品日韩在线中文字幕| 国产视频首页在线观看| 久久精品夜色国产| 另类亚洲欧美激情| 成人黄色视频免费在线看| 美女xxoo啪啪120秒动态图| 黄网站色视频无遮挡免费观看| 少妇被粗大的猛进出69影院| 国产日韩欧美亚洲二区| 老汉色av国产亚洲站长工具| 日产精品乱码卡一卡2卡三| 制服人妻中文乱码| 免费黄色在线免费观看| 十八禁网站网址无遮挡| 国产熟女午夜一区二区三区| 精品少妇久久久久久888优播| 一级黄片播放器| 亚洲人成网站在线观看播放| 母亲3免费完整高清在线观看 | 亚洲精品一区蜜桃| 久久精品国产亚洲av高清一级| 在线观看三级黄色| 老司机影院成人| 91成人精品电影| 亚洲第一区二区三区不卡| 各种免费的搞黄视频| 又黄又粗又硬又大视频| 免费看不卡的av| 2022亚洲国产成人精品| 国产精品女同一区二区软件| 自线自在国产av| 90打野战视频偷拍视频| 可以免费在线观看a视频的电影网站 | 欧美激情极品国产一区二区三区| 一级毛片 在线播放| 亚洲四区av| 日韩视频在线欧美| 91aial.com中文字幕在线观看| 欧美日韩视频高清一区二区三区二| 久久精品久久久久久久性| 久久婷婷青草| 一区二区日韩欧美中文字幕| 成年动漫av网址| 久久精品国产亚洲av涩爱| 亚洲,欧美,日韩| 久久av网站| 国产野战对白在线观看| 亚洲第一av免费看| 母亲3免费完整高清在线观看 | 精品亚洲成a人片在线观看| 肉色欧美久久久久久久蜜桃| 亚洲欧洲日产国产| 狠狠婷婷综合久久久久久88av| 观看av在线不卡| 久久国产亚洲av麻豆专区| 亚洲图色成人| 一级片免费观看大全| 91久久精品国产一区二区三区| 亚洲三区欧美一区| 久久久久国产网址| 亚洲av成人精品一二三区| 国产av国产精品国产| 一级片'在线观看视频| 一个人免费看片子| 久久精品国产a三级三级三级| 欧美 日韩 精品 国产| 国产爽快片一区二区三区| 考比视频在线观看| 午夜福利在线免费观看网站| 只有这里有精品99| av网站在线播放免费| 中文字幕av电影在线播放| 久久国内精品自在自线图片| 高清欧美精品videossex| 丝袜脚勾引网站| 国产深夜福利视频在线观看| 国产精品不卡视频一区二区| 大话2 男鬼变身卡| 中文欧美无线码| 最新的欧美精品一区二区| 精品一区二区三区四区五区乱码 | 大片免费播放器 马上看| 一边摸一边做爽爽视频免费| 人人妻人人爽人人添夜夜欢视频| 免费不卡的大黄色大毛片视频在线观看| 69精品国产乱码久久久| 国产男人的电影天堂91| 国产 精品1| 少妇被粗大的猛进出69影院| 日本色播在线视频| 咕卡用的链子| 国产xxxxx性猛交| 精品人妻偷拍中文字幕| 成人漫画全彩无遮挡| 欧美 日韩 精品 国产| 欧美最新免费一区二区三区| 国产亚洲最大av| 国产精品久久久久久精品古装| 日韩,欧美,国产一区二区三区| 国产麻豆69| 一本大道久久a久久精品| 精品国产一区二区三区四区第35| 在线观看免费日韩欧美大片| 最近手机中文字幕大全| 亚洲成av片中文字幕在线观看 | 亚洲精品在线美女| av.在线天堂| www.熟女人妻精品国产| 晚上一个人看的免费电影| 美女脱内裤让男人舔精品视频| 国产亚洲欧美精品永久| 日韩一区二区三区影片| 九草在线视频观看| 精品国产一区二区久久| 美女国产视频在线观看| 夫妻性生交免费视频一级片| 亚洲成人一二三区av| 一边摸一边做爽爽视频免费| 日日啪夜夜爽| 亚洲视频免费观看视频| av在线老鸭窝| av线在线观看网站| 女人高潮潮喷娇喘18禁视频| 精品亚洲成国产av| 日产精品乱码卡一卡2卡三| 国产成人a∨麻豆精品| 少妇精品久久久久久久| 亚洲熟女精品中文字幕| 欧美日韩精品网址| 人人妻人人澡人人爽人人夜夜| 午夜激情久久久久久久| 伊人久久大香线蕉亚洲五| 最近最新中文字幕大全免费视频 | 99香蕉大伊视频| 最近手机中文字幕大全| 国产欧美亚洲国产| 久久久久视频综合| 国产精品蜜桃在线观看| 久久久久久久久久久久大奶| 亚洲国产欧美网| 久热这里只有精品99| 日本色播在线视频| 日韩 亚洲 欧美在线| 国产欧美亚洲国产| 中文字幕最新亚洲高清| 十八禁高潮呻吟视频| 男女边摸边吃奶| 男男h啪啪无遮挡| 日韩欧美一区视频在线观看| 亚洲av男天堂| 蜜桃国产av成人99| 老鸭窝网址在线观看| 亚洲国产精品一区三区| 女人高潮潮喷娇喘18禁视频| 69精品国产乱码久久久| 久久精品夜色国产| 亚洲国产精品一区二区三区在线| 久久午夜福利片| 亚洲国产精品国产精品| 亚洲熟女精品中文字幕| 欧美精品亚洲一区二区| 欧美日韩一区二区视频在线观看视频在线| 色视频在线一区二区三区| 天天躁狠狠躁夜夜躁狠狠躁| 婷婷色综合大香蕉| 秋霞在线观看毛片| 天堂中文最新版在线下载| 综合色丁香网| 精品人妻在线不人妻| 国产欧美日韩综合在线一区二区| 人妻系列 视频| 黑丝袜美女国产一区| 黑人猛操日本美女一级片| 999久久久国产精品视频| 你懂的网址亚洲精品在线观看| 成年人免费黄色播放视频| 一边摸一边做爽爽视频免费| 青春草视频在线免费观看| 这个男人来自地球电影免费观看 | 性少妇av在线| 国产成人a∨麻豆精品| 麻豆乱淫一区二区| 最近手机中文字幕大全| 婷婷成人精品国产| 午夜免费观看性视频| 国产女主播在线喷水免费视频网站| 欧美精品国产亚洲| 五月天丁香电影| h视频一区二区三区| 视频在线观看一区二区三区| 青春草亚洲视频在线观看| 建设人人有责人人尽责人人享有的| 欧美黄色片欧美黄色片| 亚洲熟女精品中文字幕| 亚洲色图综合在线观看| 久久久久久久大尺度免费视频| 久久久久精品人妻al黑| 成人国产麻豆网| 亚洲国产欧美网| 亚洲第一青青草原| 亚洲欧洲日产国产| 日本91视频免费播放| 日本午夜av视频| 丝袜脚勾引网站| 99精国产麻豆久久婷婷| 中文字幕av电影在线播放| 午夜福利一区二区在线看| 久久这里只有精品19| 九色亚洲精品在线播放| 在线看a的网站| 国产精品久久久久久av不卡| 久久久久国产网址| 人人妻人人添人人爽欧美一区卜| 国产精品国产av在线观看| av片东京热男人的天堂| 纯流量卡能插随身wifi吗| 亚洲欧美精品自产自拍| 一级毛片我不卡| 欧美日本中文国产一区发布| 人人妻人人添人人爽欧美一区卜| 80岁老熟妇乱子伦牲交| 中文字幕色久视频| 老女人水多毛片| 99热国产这里只有精品6| 少妇熟女欧美另类| 亚洲精品视频女| av电影中文网址| 大香蕉久久网| 中文字幕av电影在线播放| 免费看不卡的av| 777米奇影视久久| 高清在线视频一区二区三区| 亚洲欧洲国产日韩| 国产精品蜜桃在线观看| 亚洲色图 男人天堂 中文字幕| 777久久人妻少妇嫩草av网站| 国产高清不卡午夜福利| 男男h啪啪无遮挡| 亚洲,欧美,日韩| 欧美国产精品一级二级三级| 精品一品国产午夜福利视频| 啦啦啦在线观看免费高清www| 精品人妻一区二区三区麻豆| 国产成人91sexporn| 超碰97精品在线观看| 伊人亚洲综合成人网| 咕卡用的链子| 大陆偷拍与自拍| 一区二区三区四区激情视频| 亚洲在久久综合| 亚洲精品国产色婷婷电影| 亚洲国产看品久久| 中文精品一卡2卡3卡4更新| 日韩制服丝袜自拍偷拍| 日韩伦理黄色片| 精品国产超薄肉色丝袜足j| 久久久久久久久久久久大奶| 老司机亚洲免费影院| 亚洲成色77777| 亚洲精品av麻豆狂野| 亚洲久久久国产精品| 亚洲欧洲精品一区二区精品久久久 | 久久久久精品性色| 亚洲国产精品一区二区三区在线| 高清欧美精品videossex| 国产男女超爽视频在线观看| 精品卡一卡二卡四卡免费| 国产一区二区三区综合在线观看| 午夜免费鲁丝| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 黄色 视频免费看| 只有这里有精品99| 精品国产超薄肉色丝袜足j| 这个男人来自地球电影免费观看 | 国产精品国产av在线观看| 久久久国产一区二区| 日本91视频免费播放| 老司机影院成人| 在线天堂最新版资源| 丝袜美足系列| 老汉色av国产亚洲站长工具| 日本爱情动作片www.在线观看| 免费在线观看完整版高清| 激情五月婷婷亚洲| 国产不卡av网站在线观看| 欧美国产精品一级二级三级| 国产成人午夜福利电影在线观看| 99热国产这里只有精品6| 免费观看无遮挡的男女| 久久狼人影院| 欧美日韩精品成人综合77777| 精品卡一卡二卡四卡免费| 久久久国产欧美日韩av| 中文字幕人妻丝袜一区二区 | 爱豆传媒免费全集在线观看| 三级国产精品片| 精品酒店卫生间| 日韩中文字幕视频在线看片| 女的被弄到高潮叫床怎么办| 久久女婷五月综合色啪小说| 中国三级夫妇交换| 中文字幕制服av| 性高湖久久久久久久久免费观看| 午夜激情久久久久久久| 18禁动态无遮挡网站| 大话2 男鬼变身卡| 人人妻人人添人人爽欧美一区卜| 精品人妻一区二区三区麻豆| 只有这里有精品99| 老司机亚洲免费影院| videosex国产| 午夜久久久在线观看| 精品亚洲乱码少妇综合久久| 男的添女的下面高潮视频| 一边亲一边摸免费视频| www.自偷自拍.com| 美女大奶头黄色视频| 免费观看无遮挡的男女| 久久国产精品男人的天堂亚洲| 久久久亚洲精品成人影院| 亚洲色图综合在线观看| 99re6热这里在线精品视频| 亚洲精品视频女| 午夜老司机福利剧场| 精品视频人人做人人爽| 亚洲人成电影观看| 九色亚洲精品在线播放| 大陆偷拍与自拍| av一本久久久久| 欧美成人午夜免费资源| 久久免费观看电影| 在线观看免费高清a一片| 乱人伦中国视频| 国产精品无大码| 亚洲精品国产av成人精品| 蜜桃在线观看..| 亚洲伊人久久精品综合| 婷婷成人精品国产| 在线观看免费日韩欧美大片| 女性被躁到高潮视频| 免费高清在线观看日韩| 精品少妇久久久久久888优播| 99久久人妻综合| 免费高清在线观看视频在线观看| 人妻人人澡人人爽人人| 欧美激情 高清一区二区三区| 桃花免费在线播放| 国产爽快片一区二区三区| 91aial.com中文字幕在线观看| 久久韩国三级中文字幕| 18+在线观看网站| 少妇人妻久久综合中文| 妹子高潮喷水视频| 欧美最新免费一区二区三区| 久久热在线av| 最新的欧美精品一区二区| 看免费av毛片| 午夜福利一区二区在线看| 汤姆久久久久久久影院中文字幕| 久久久久精品性色| 国产一区二区激情短视频 | 欧美黄色片欧美黄色片| 精品久久蜜臀av无| 久久久久久人人人人人| 男女下面插进去视频免费观看| 卡戴珊不雅视频在线播放| 亚洲国产精品999| 国产视频首页在线观看| 精品国产一区二区三区四区第35| 91国产中文字幕| 国产 精品1| 叶爱在线成人免费视频播放| 久久精品久久精品一区二区三区| 欧美成人精品欧美一级黄| 日韩一本色道免费dvd| 精品一区在线观看国产| 一区二区三区激情视频| 精品国产一区二区三区久久久樱花| 美女xxoo啪啪120秒动态图| 男女无遮挡免费网站观看| 2022亚洲国产成人精品| 国产免费一区二区三区四区乱码| 考比视频在线观看| 国产精品国产三级专区第一集| 在线 av 中文字幕| 国产午夜精品一二区理论片| 午夜日本视频在线| 亚洲国产精品一区三区| 另类亚洲欧美激情| av免费观看日本| 欧美精品亚洲一区二区| 国产精品av久久久久免费| 一本—道久久a久久精品蜜桃钙片| 亚洲图色成人| 久久免费观看电影| 成人漫画全彩无遮挡| 国产亚洲av片在线观看秒播厂| freevideosex欧美| 国产成人aa在线观看| 成人毛片60女人毛片免费| 成年av动漫网址| 亚洲av中文av极速乱| 亚洲四区av| 天堂中文最新版在线下载| 亚洲精品av麻豆狂野| 国产午夜精品一二区理论片| 国产亚洲欧美精品永久| 水蜜桃什么品种好| 人人妻人人爽人人添夜夜欢视频| 成人漫画全彩无遮挡| 欧美变态另类bdsm刘玥| 成人国产麻豆网| 三上悠亚av全集在线观看| 香蕉国产在线看| 高清不卡的av网站| 午夜免费观看性视频| 色哟哟·www| 久久久久久人妻| 满18在线观看网站| 国产免费一区二区三区四区乱码| 美女高潮到喷水免费观看| 成人国产av品久久久| 国产精品人妻久久久影院| 另类精品久久| 精品国产乱码久久久久久男人| videossex国产| 秋霞伦理黄片| 日韩av免费高清视频| 亚洲欧美色中文字幕在线| 一本—道久久a久久精品蜜桃钙片| 国产亚洲av片在线观看秒播厂| 精品一区二区免费观看| 亚洲三区欧美一区| 日韩一本色道免费dvd| 99久国产av精品国产电影| 边亲边吃奶的免费视频| a 毛片基地| 丝袜美腿诱惑在线| videossex国产| 亚洲成色77777| 丝袜美足系列| 欧美+日韩+精品| 欧美国产精品一级二级三级| 日韩欧美一区视频在线观看| 亚洲av福利一区| 久久婷婷青草| 国产精品三级大全| 在线观看人妻少妇| 一区在线观看完整版| 曰老女人黄片| 成人亚洲精品一区在线观看| 久久久精品94久久精品| 国产在线一区二区三区精| 汤姆久久久久久久影院中文字幕| 侵犯人妻中文字幕一二三四区| 亚洲av电影在线进入| 伦理电影免费视频| 深夜精品福利| 欧美日韩视频高清一区二区三区二| 亚洲第一青青草原| 啦啦啦在线观看免费高清www| 男人添女人高潮全过程视频| 国产亚洲欧美精品永久| 乱人伦中国视频| 欧美日韩一级在线毛片| 婷婷成人精品国产| 国产精品成人在线| 啦啦啦啦在线视频资源| 男人操女人黄网站| 国产欧美日韩一区二区三区在线| 免费在线观看黄色视频的| 纯流量卡能插随身wifi吗| 免费播放大片免费观看视频在线观看| 亚洲精品国产av蜜桃| 精品一区二区三区四区五区乱码 | 国产高清国产精品国产三级| 国产精品一区二区在线观看99| 亚洲欧美一区二区三区国产| 欧美人与性动交α欧美精品济南到 | 亚洲欧美成人精品一区二区| 国产免费一区二区三区四区乱码| 久久久久久久久久久久大奶| 自线自在国产av| www.精华液| 久久精品人人爽人人爽视色| 少妇 在线观看| 国产麻豆69| 久久精品aⅴ一区二区三区四区 | 18禁国产床啪视频网站| 99九九在线精品视频| 一级,二级,三级黄色视频| 99久久人妻综合| a级毛片在线看网站| 久久精品久久久久久久性| 一本—道久久a久久精品蜜桃钙片| 日韩在线高清观看一区二区三区| 性色avwww在线观看| 亚洲精品,欧美精品| 丁香六月天网| 我的亚洲天堂| 久久精品国产亚洲av高清一级| 丝袜在线中文字幕| 777米奇影视久久| 久久久久国产网址| 97人妻天天添夜夜摸| 日本wwww免费看| 国产午夜精品一二区理论片| 日韩欧美精品免费久久| 亚洲人成电影观看| 少妇精品久久久久久久| 免费观看性生交大片5| 国产av码专区亚洲av| 超碰97精品在线观看| 午夜激情久久久久久久| 看十八女毛片水多多多| 精品少妇一区二区三区视频日本电影 | 久久久久久久久久久久大奶| 国产深夜福利视频在线观看| 成人18禁高潮啪啪吃奶动态图| 免费女性裸体啪啪无遮挡网站| 黄色 视频免费看| 日本猛色少妇xxxxx猛交久久| 久久女婷五月综合色啪小说| 青春草国产在线视频| 色吧在线观看| 欧美日本中文国产一区发布| 日韩 亚洲 欧美在线| av片东京热男人的天堂| 亚洲在久久综合| www日本在线高清视频| 男女啪啪激烈高潮av片| 国产精品久久久久久精品古装| 老鸭窝网址在线观看| 日韩精品有码人妻一区| 自线自在国产av| 色吧在线观看| 巨乳人妻的诱惑在线观看| 色哟哟·www| 午夜福利,免费看| 91精品伊人久久大香线蕉| 亚洲欧美精品自产自拍| 精品少妇一区二区三区视频日本电影 | 久久久久网色| 亚洲色图 男人天堂 中文字幕| 在线观看免费视频网站a站| 色播在线永久视频| 国产精品无大码| 男女无遮挡免费网站观看| 99九九在线精品视频| 波多野结衣一区麻豆| 久久久精品免费免费高清| 亚洲国产最新在线播放| 国产精品偷伦视频观看了| 美女中出高潮动态图| 亚洲国产欧美网| 国产 一区精品| 在线观看免费视频网站a站| 性色avwww在线观看| 欧美日韩精品成人综合77777| 18+在线观看网站| 五月天丁香电影| 亚洲综合色惰| 在线天堂中文资源库| 日韩一区二区三区影片| 大片电影免费在线观看免费|