《電子技術應用》
您所在的位置:首頁 > 嵌入式技術 > 設計應用 > 基于四元數的前臂假肢手部自平衡的設計
基于四元數的前臂假肢手部自平衡的設計
2016年電子技術應用第5期
樊炳輝1,張凱麗1,王傳江1,劉圭圭1,張芳麗2
1.山東科技大學機器人研究中心,山東 青島266590;2.友邦電力設備有限公司,黑龍江 哈爾濱150000
摘要: 針對殘疾人的前臂假肢在端水等運動過程中需要手部虎口一直與大地保持平行,即手部實現姿態自平衡問題,提出采用四元數的姿態解算算法,并采用卡爾曼濾波融合MPU6050傳感器中陀螺儀和加速度計的輸出數據,補償了陀螺儀輸出數據的隨機漂移誤差,得到前臂假肢手部姿態的正確估計。該系統由安卓手機、STM32單片機、MPU6050傳感器、高性能數字舵機組成,通過安卓手機語音識別向STM32單片機發送是否保持手部自平衡指令,由STM32實時采集MPU6050中陀螺儀和加速度計的數據并經過算法處理后,通過串口向舵機發送信號。實驗結果表明,這種姿態估計算法有效地解決了前臂假肢在運動過程中手部自平衡問題,為前肢缺失的殘疾人提供了便利。
中圖分類號: TP242.6
文獻標識碼: A
DOI:10.16157/j.issn.0258-7998.2016.05.022
中文引用格式: 樊炳輝,張凱麗,王傳江,等. 基于四元數的前臂假肢手部自平衡的設計[J].電子技術應用,2016,42(5):78-81.
英文引用格式: Fan Binghui,Zhang Kaili,Wang Chuanjiang,et al. Design of hand self-balance for the forearm prosthetic based on quaternion[J].Application of Electronic Technique,2016,42(5):78-81.
Design of hand self-balance for the forearm prosthetic based on quaternion
Fan Binghui1,Zhang Kaili1,Wang Chuanjiang1,Liu Guigui1,Zhang FangLi2
1.Science and Technology Robot Research Center of Shandong University,Qingdao 266590,China; 2.Youbang Electric Limited Company,Harbin 150000,China
Abstract: The forearm prosthetic hand designed for the disabled need achieve it′s first web keeping in parallel with the earth in the processing of movement like drinking and other actions. That is, the hand realizes the attitude self-balancing problem. Aiming at the estimation of the forearm prosthesis′s attitude angles, an attitude solution algorithm based on quaternion is used. Using Kalman filter to fuse the data from gyroscope and accelerometer in MPU6050, an optimal estimation for the attitude of forearm prosthesis is achieved after the compensation of the angular velocity drift error of gyroscope. The system consists of android phones, STM32 microcontroller、MPU6050、high performance of steering gear. Through the voice recognition, android mobile phone sends the instructions whether the forearm prosthetic hand need to keep the balance to STM32. Actuating signal can be send to steering gear through the serial port after data are collected and processed with algorithm from acceleration sensor and gyroscope in MPU6050 via STM32. The experimental results show that the attitude estimation algorithm is effective and the forearm prosthetic hand can keep balance during motion processes. This facilitates the disabled with the lack of the front arm.
Key words : the forearm prosthesis;attitude estimation;quaternion;Kalman filter

0 引言

    世界上每年因戰爭、自然災害、交通事故等原因引起的殘疾人數多達千萬人以上,其中大部分為肢體殘疾,而目前市場上的大部分假肢都是作為裝飾性用品,真正實現功能代償的假肢數量較少。為滿足殘疾人日常生活的需求,本系統對前臂假肢做了姿態規劃,即控制前臂旋轉和手腕俯仰兩個電機,使得假肢末端手爪的虎口平面與大地水平面始終保持平行,如圖1所示。這樣假肢在抓取和移動目標物體的過程中,可以保證運動平穩,對移動水杯等類似操作尤其重要。在前臂假肢的平衡控制過程中,姿態傳感器將檢測到的前臂假肢的姿態信號送入STM32,經算法運算得到控制量驅動舵機來控制前臂的扭轉和手部的俯仰,從而使前臂假肢的手部恢復到平衡狀態??梢钥闯?,獲取實時的前臂假肢的姿態是保持手部平衡的前提。因此,實現一種實時可靠的前臂假肢姿態估計算法尤為重要,它直接影響其手部自平衡性能。

ck5-t1.gif

    姿態解算的算法有多種,常用的算法有歐拉角法、方向余弦法與四元數法,由于歐拉角法求解時,方程中存在奇點,方向余弦法計算量大,因此現在常用的方法為四元數法,其計算量小,無奇點[1-2]。姿態估計算法必須考慮到傳感器容易受溫度及噪聲的影響而產生不同程度的漂移,它會影響姿態估計的準確性。針對這一問題,本文采用卡爾曼濾波的方法來融合陀螺儀和加速度計輸出的數據,修正了陀螺儀輸出量中的隨機漂移分量[3]。

1 系統硬件設計

    隨著安卓智能手機的發展與普及,為了讓前臂假肢的控制更加方便和標準化,提出了以安卓手機離線識別語音及應用藍牙通信技術實現安卓手機無線控制假肢運動的方案。本文介紹的前臂假肢可實現手部平衡,除此之外,在不需要手部保持平衡的環境下,安卓手機通過語音識別可控制前臂假肢關閉手部自平衡功能和實現手部上抬、下落、內轉、外轉和張合等動作。該系統由安卓手機、藍牙模塊、MPU6050、STM32微控制器、MX-28高性能舵機組成,系統的控制流程如圖2所示。

ck5-t2.gif

    其中,選用MPU6050做為前臂假肢的姿態檢測系統,MPU6050在單個產品封裝內集成了3軸陀螺儀和3軸加速度計。由于前臂假肢不需獲知偏航角,所以本文采用了MPU6050來獲取前臂假肢的俯仰角和橫滾角。整個控制過程由MPU6050測量姿態變化,并通過前臂扭轉以及手腕俯仰來根據姿態變化實時進行調整,使其達到自平衡狀態。ck5-t3.gif

    MPU6050安裝在假肢的腕部,其安裝位置如圖3所示。將MPU6050安裝在控制前臂旋轉關節的上部并固定在連桿上,左右位置可以變化,但要確保其所在的平面與手爪虎口所在平面平行。安裝完成后,前臂繞MPU6050的X軸旋轉,手腕繞MPU6050的Y軸俯仰。

2 四元數姿態解算算法

    四元數的數學概念是1843年由哈密頓首先提出的。它是由一個實部單位和3個虛部單位i、j、k組成,其形式為:q=q0+q1i+q2j+q3k。四元數可以表示旋轉,其標量部分表示旋轉角度的大小,矢量部分表示轉軸的方向[4-5],應用四元數時必須進行歸一化,也就是令其模值為1。

    前臂假肢的定位由固定的內部的三維坐標確定,一個坐標系到另一個坐標系的變換,可以通過繞不同坐標軸的3次連續轉動來實現。順時針繞橫滾軸,即繞載體坐標系的X軸旋轉,得到橫滾角φ;順時針繞俯仰軸,即繞載體坐標系的Y軸旋轉,得到俯仰角θ;順時針繞偏航軸,即繞載體坐標系的?爪軸旋轉,得到偏航角γ[6]。

    四個參數的初始值是由歐拉角的初始值經過以下計算得到的[6-9]

    ck5-gs1-3.gif

    利用四元數微分方程可以實現四元數的更新,其中,ω是陀螺儀測量的角速度。

    應用四階龍格—庫塔算法解微分方程,得到4個參數的更新值。具體過程如下[12-14]

    設:

    ck5-gs4.gif

    則:

    ck5-gs5.gif

其中,T為更新周期,q(t)為初始四元數,q(t+T)為更新后的四元數,K1是時間段T開始時的斜率,K2是時間段中點的斜率,K3也是中點的斜率,K4是時間段終點的斜率。

    計算得到這4個參數之后,可由姿態矩陣T更新歐拉角[15]:

    ck5-gs6-8.gif

    得到姿態角,即可得到前臂假肢的姿態估計。

    前臂假肢在平衡的情況下,歐拉角的理想輸出值應該是0°,但由于陀螺儀有隨機漂移誤差,歐拉角會一直積累[16]。實驗結果圖4所示。

ck5-t4.gif

    針對這一問題,提出了應用卡爾曼濾波的數據融合算法來解決。卡爾曼濾波可以將陀螺儀的輸出數據和加速度計的輸出數據融合起來。陀螺儀用于測量角速度,動態跟蹤性能好,但陀螺儀受溫度、不穩定力矩等因素的影響,會產生隨機漂移誤差;加速度計可用于測量加速度,通過測得的加速度值與重力加速度的比值得到傾角值,加速度計靜態性能好,但是動態跟蹤性能差。可以看出,陀螺儀和加速度計優缺點互補,結合起來會有好的效果[17-18]。

3 實驗結果分析

3.1 四元數法和卡爾曼濾波實驗結果分析

    將MPU6050靜止平放在桌上,采集其經四元數姿態解算和卡爾曼濾波后輸出的數據,通過MATLAB工具對實測數據進行仿真分析,如圖5所示,俯仰角和橫滾角都在附近,對比圖4發現卡爾曼濾波有效解決了陀螺儀的隨機漂移誤差,可以準確地獲取靜態時前臂假肢的姿態。

ck5-t5.gif

    將MPU6050安裝到前臂假肢上,令前臂假肢繞MPU6050的X軸從平衡位置順時針逐步旋轉到-40°,一段時間后,回到平衡位置,再沿X軸逆時針逐步旋轉到40°,采集這一過程MPU6050經四元數姿態解算和卡爾曼濾波后輸出的數據,通過MATLAB工具對實測數據進行仿真分析,橫滾角發生相應的變化,俯仰角基本保持不變,如圖6(a)所示。令前臂假肢繞MPU6050的Y軸重復以上實驗,根據圖6(b)可以看出,俯仰角發生相應的變化,橫滾角基本保持不變。由于實驗過程中手動旋轉前臂假肢,操作者偶有抖動和難做到勻速轉動,因此,圖6中的波形中偶有小的波動。實驗表明,本文采用的姿態估計可以實現前臂假肢動作時姿態的正確估計,經過STM32運算處理后,實現前臂假肢在運動過程中保持手部平衡。

ck5-t6.gif

3.2 前臂假肢手部自平衡實驗結果分析

    對前臂假肢系統組裝并調試,通過安卓手機軟件向假肢發出保持手部平衡的命令。具體測試步驟:連接好線路,打開電源,觀察指示燈正常亮起,觸擊控制界面的藍牙連接按鈕,觀察指示燈,以確保手機藍牙和藍牙模塊的配對連接成功;否則,檢查系統的連接是否正確。藍牙配對成功后,向安卓手機喊出“平衡”指令,手機通過藍牙將指令傳到單片機,從而控制前臂假肢運動。實驗結果如圖7、圖8所示。

ck5-t7.gif

ck5-t8.gif

    圖7(a)~圖7(e)展示了前臂假肢從下往上轉動過程中,其手部可以實現自平衡,手部虎口一直保持與大地平行。

    圖8中以箭頭指向的地點為參考點,可以看出圖8(a)~圖8(c)為前臂假肢由內往外轉動過程,同樣假肢手部能實現自平衡。

4 結論

    此前臂假肢采用四元數法的姿態解算算法,并利用卡爾曼濾波融合MPU6050中陀螺儀和加速度計的數據,得到其正確的姿態估計。經實驗測試,該系統運行穩定、可靠,所設計的姿態估計算法能夠很好地跟蹤前臂假肢的真實姿態,實現了假肢手部的姿態自平衡控制。此前臂假肢可以輔助殘疾人完成喝水等動作,方便了他們的生活,也減輕了家屬的負擔。

參考文獻

[1] 杜海龍,張榮輝,劉平,等.捷聯慣導系統姿態解算模塊的實現[J].光學精密工程,2008,16(10):1956-1962.

[2] 張榮輝,賈宏光,陳濤,等.基于四元數法的捷聯式慣性導航系統的姿態解算[J].光學精密工程,2008,16(10):1963-1970.

[3] 楊凌霄,李曉陽.基于卡爾曼濾波的兩輪自平衡車姿態解算方法[J].計算機仿真,2014,31(6):406-409.

[4] 葉锃鋒,馮恩信.基于四元數和卡爾曼濾波的兩輪車姿態穩定方法[J].傳感技術學報,2012,25(4):524-528.

[5] 呂印新,肖前貴,胡壽松.基于四元數互補濾波的無人機姿態解算[J].燕山大學學報,2014,38(2):175-180.

[6] 張思奇,趙憶文.便攜式移動機器人姿態解算方法[J].微計算機信息,2007,23(10-2):188-190.

[7] 祝喬,汪旭東,崔家瑞.基于四元數的航天器姿態非線性控制器[J].信息與控制,2012,41(2):170-179.

[8] 賈秋玲,閆建國,王建民.空間多剛體系統姿態的協同控制[J].信息與控制,2005,34(4):385-388.

[9] 趙國文.機器人運動學中的四元數矩陣和基變換矩陣[J].機器人,1993,15(2):59-64.

[10] 張天光,王秀萍,王麗霞,等,譯.捷聯慣性導航技術(第二版)[M].北京:國防工業出版社,2007:29-33.

[11] 秦永元.慣性導航[M].北京:科學出版社,2006:287-305.

[12] 崔璐璐.基于MEMS器件的姿態測量系統研究與實現[D].大連:大連理工大學,2009.

[13] 王勵揚,霍昆明,何文濤,等.四階龍格庫塔算法在捷聯慣性導航中的應用[J].計算機仿真,2014,31(11):56-59.

[14] 張慧春,吳簡彤,何昆鵬,等.四階龍格-庫塔法在捷聯慣導系統姿態解算中的應用[J].應用科技,2005,32(6):37-39.

[15] 鄧正隆.慣性技術[M].哈爾濱:哈爾濱工業大學出版社,2006:14-164.

[16] 金光明,張國良,陳林鵬.MEMS陀螺儀靜態漂移模型與濾波方法研究[J].傳感器與微系統,2007,26(11):48-50.

[17] 馮智勇,曾瀚,張力,等.基于陀螺儀及加速度計信號融合的姿態角度測量[J].西南師范大學學報(自然科學版),2011,36(4):137-141.

[18] 劉星.多維MEMS慣性傳感器的姿態解算算法研究[D].哈爾濱:哈爾濱工程大學,2013.

此內容為AET網站原創,未經授權禁止轉載。
热re99久久精品国产66热_欧美小视频在线观看_日韩成人激情影院_庆余年2免费日韩剧观看大牛_91久久久久久国产精品_国产原创欧美精品_美女999久久久精品视频_欧美大成色www永久网站婷_国产色婷婷国产综合在线理论片a_国产精品电影在线观看_日韩精品视频在线观看网址_97在线观看免费_性欧美亚洲xxxx乳在线观看_久久精品美女视频网站_777国产偷窥盗摄精品视频_在线日韩第一页
  • <strike id="ygamy"></strike>
  • 
    
      • <del id="ygamy"></del>
        <tfoot id="ygamy"></tfoot>
          <strike id="ygamy"></strike>
          欧美精品三级在线观看| 亚洲激情不卡| 伊人久久大香线蕉av超碰演员| 欧美日韩大片一区二区三区| 欧美日韩高清不卡| 久久人人看视频| av成人福利| 国产真实精品久久二三区| 久久久国际精品| 亚洲午夜精品一区二区三区他趣| 国产精品午夜在线| 国产精品爱久久久久久久| 国产欧美在线看| 国产欧美一区二区三区在线老狼| 欧美日本一区二区视频在线观看| 欧美亚洲综合另类| 国产精品视频网| 久久一区激情| 一区二区激情小说| 欧美激情网站在线观看| 亚洲性人人天天夜夜摸| 国产精品久久久久久一区二区三区| 欧美一级在线视频| 欧美激情视频在线免费观看 欧美视频免费一| 亚洲国产欧美一区二区三区久久| 国产欧美精品一区二区色综合| 亚洲国产一二三| 国产精品一区二区你懂的| 国产精品久久久久久久久久久久久久| 亚洲精品国产精品国自产观看浪潮| 麻豆freexxxx性91精品| 在线观看成人小视频| 欧美日韩国产影片| 久久久久久9| 欧美电影美腿模特1979在线看| 宅男噜噜噜66国产日韩在线观看| 一区在线电影| 国产精品久久久一区麻豆最新章节| 久久久精品国产99久久精品芒果| 91久久夜色精品国产网站| 1204国产成人精品视频| av72成人在线| 欧美日韩国产欧美日美国产精品| 国产精品久久久久久福利一牛影视| 99re这里只有精品6| 另类天堂视频在线观看| 欧美国产综合一区二区| 国产一区二区三区在线观看视频| 午夜精品久久久久久久白皮肤| 伊人精品在线| 国产一区二区电影在线观看| 亚洲国产一区二区三区高清| 狠狠色狠狠色综合日日tαg| 久久久久久香蕉网| 国产精品一二三| 久久婷婷激情| 午夜精品久久久久久久久| 欧美日韩八区| 国产欧美精品一区二区三区介绍| 欧美高清在线观看| 国产亚洲二区| 国产精品中文在线| 亚洲欧美日韩精品一区二区| 亚洲天堂成人在线视频| 亚洲欧美激情精品一区二区| 韩国av一区二区三区| 亚洲精选成人| 亚洲国产三级网| 99精品视频免费全部在线| 免费成人在线观看视频| 最新国产乱人伦偷精品免费网站| 国产日韩在线亚洲字幕中文| 亚洲免费久久| 亚洲七七久久综合桃花剧情介绍| 欧美日韩在线观看一区二区三区| 久久久久一区二区| 欧美激情综合色综合啪啪| 久久夜色精品亚洲噜噜国产mv| 亚洲一区二区三区四区五区黄| 欧美日韩一区三区| 欧美高清视频一二三区| 国产视频久久网| 新狼窝色av性久久久久久| 国产精品v日韩精品| 久久精品国产一区二区三区| 艳妇臀荡乳欲伦亚洲一区| 亚洲美女视频| 欧美一级黄色网| 欧美一级二区| 欧美成人一区二区三区| 国产日韩精品电影| 亚洲自拍偷拍视频| 欧美成人r级一区二区三区| 国产日韩欧美精品在线| 亚洲国产日韩欧美综合久久| 在线亚洲一区| 久久露脸国产精品| 亚洲欧美日韩一区二区三区在线观看| 国产精品一区2区| 国产欧美日韩视频一区二区| 伊人久久av导航| 国产日韩一区二区三区在线播放| 亚洲丰满少妇videoshd| 亚洲成人资源网| 另类酷文…触手系列精品集v1小说| 免费h精品视频在线播放| 国产精品久久久久久久久久妞妞| 宅男噜噜噜66国产日韩在线观看| 欧美成人第一页| 在线免费不卡视频| 久久香蕉国产线看观看网| 久久高清一区| 在线观看一区| 一区二区在线视频播放| 欧美久久久久久久| 欧美性做爰毛片| 国产精品久久久久9999高清| 亚洲精品123区| 国产欧美91| 欧美日韩国内| 亚洲欧美日韩精品久久奇米色影视| 性欧美在线看片a免费观看| 久久激情五月激情| 亚洲国产成人久久综合一区| 午夜一区二区三视频在线观看| 99人久久精品视频最新地址| 亚洲专区一区二区三区| 黄色精品一二区| 先锋影院在线亚洲| 国产亚洲精品综合一区91| 国产日韩欧美视频在线| 欧美大片一区二区| 久久久久九九九| 国产在线日韩| 韩国视频理论视频久久| 国产精品女人网站| 一区二区日韩免费看| 欧美大片一区二区| 裸体一区二区| 亚洲欧美一级二级三级| 亚洲精品久久久久久一区二区| 亚洲欧美一区在线| 在线一区二区视频| 国产精品久久久久9999吃药| 久久久免费精品视频| 亚洲精一区二区三区| 亚洲欧洲一区二区在线播放| 欧美成人精精品一区二区频| 亚洲午夜一二三区视频| 国产日韩欧美另类| 亚洲日本无吗高清不卡| 牛夜精品久久久久久久99黑人| 欧美国产精品日韩| 99热在这里有精品免费| 欧美日韩午夜| 久久国产日本精品| 欧美激情视频在线播放| 久久久久久综合| 国产精品国产馆在线真实露脸| 免费看成人av| 国产精品青草久久久久福利99| 亚洲人成毛片在线播放| 亚洲第一页在线| 亚洲一区在线直播| 亚洲一级片在线看| 国产一区二区三区四区老人| 国产精品高清在线| 欧美日本视频在线| 久久精品欧美| 中文国产亚洲喷潮| 久久成人18免费观看| 亚洲视频欧美在线| 制服诱惑一区二区| 一本色道久久综合精品竹菊| 欧美三级资源在线| 亚洲性人人天天夜夜摸| 亚洲欧美日韩精品久久亚洲区| 国产一区二区0| 午夜精彩视频在线观看不卡| 国产一本一道久久香蕉| 欧美国产在线电影| 欧美激情中文字幕一区二区| 久久久欧美精品sm网站| 欧美日韩国产首页| 亚洲午夜视频在线观看| 欧美日韩高清免费| 美女诱惑一区| 亚洲精品午夜精品| 国产精品国产三级国产aⅴ9色| 国产手机视频一区二区| 欧美激情一区| 国产精品尤物福利片在线观看| 久久精品国产在热久久| 久久成人综合网| 亚洲美女av网站| 国产日韩欧美一区二区| 国产精品日韩欧美大师| 亚洲一区二区日本| 久久亚洲精品欧美| 亚洲大胆av| 国产精品免费看片| 国产一区二区高清视频| 欧美日韩国产精品一区二区亚洲| 一区二区不卡在线视频 午夜欧美不卡'| 欧美日韩大陆在线| 亚洲午夜精品福利| 久久综合狠狠综合久久综青草| 亚洲免费小视频| 亚洲狠狠婷婷| 国产女人18毛片水18精品| 激情综合自拍| 亚洲自拍偷拍色片视频| 国产一区欧美| 欧美日韩一区三区四区| 伊人一区二区三区久久精品| 欧美在线不卡视频| 欧美视频福利| 免费看黄裸体一级大秀欧美| 久久夜精品va视频免费观看| 午夜久久久久久| 亚洲国产精品专区久久| 狠狠狠色丁香婷婷综合久久五月| 亚洲第一伊人| 国产婷婷精品| 久久漫画官网| 国产精品v片在线观看不卡| 欧美www在线| 国产精品www色诱视频| 亚洲久久一区| 免费日韩av片| 久久久久国产精品午夜一区| 欧美色中文字幕| 亚洲激情国产精品| 亚洲男人影院| 亚洲国产精品第一区二区三区| 欧美成人情趣视频| 国产精品免费观看视频| 欧美日韩免费观看一区二区三区| 久久国产福利| 99国产精品国产精品毛片| 欧美在线观看日本一区| 久热精品在线| 伊人精品久久久久7777| 欧美日韩国产123| 欧美一区视频在线| 亚洲国产精品成人综合色在线婷婷| 1769国内精品视频在线播放| 性色av一区二区三区红粉影视| 亚洲国产精品一区二区www| 美女黄毛**国产精品啪啪| 欧美成人免费播放| 国内综合精品午夜久久资源| 亚洲影院色在线观看免费| 国产亚洲福利社区一区| 亚洲欧美自拍偷拍| 亚洲影视中文字幕| 一区二区三区鲁丝不卡| 一色屋精品视频在线观看网站| 先锋影音久久| 欧美激情中文字幕在线| 亚洲福利视频在线| 欧美激情中文字幕乱码免费| 欧美日韩成人一区二区三区| 欧美日韩免费高清一区色橹橹| 一区二区三区四区蜜桃| 女同性一区二区三区人了人一| 久久一区免费| 午夜精品久久久久久99热| 久久久久久久一区二区三区| 亚洲人体一区| 中文在线资源观看视频网站免费不卡| 久久久久一区| 尤物网精品视频| 夜夜爽99久久国产综合精品女不卡| 黄网动漫久久久| 亚洲国产裸拍裸体视频在线观看乱了中文| 免费观看国产成人| 欧美日韩午夜激情| 欧美国产高潮xxxx1819| 亚洲片在线资源| 榴莲视频成人在线观看| 国产一区二区三区直播精品电影| 欧美中在线观看| 国产日韩在线亚洲字幕中文| 亚洲综合电影| 欧美日韩亚洲成人| 亚洲免费精彩视频| 亚洲天堂视频在线观看| 亚洲午夜精品在线| 亚洲三级网站| 亚洲国产精品悠悠久久琪琪| 99www免费人成精品| 欧美视频中文一区二区三区在线观看| 国产精品久久国产精麻豆99网站| 欧美日韩精品国产| 欧美一级在线播放| 久久偷看各类wc女厕嘘嘘偷窃| 国产精品高潮呻吟久久| 亚洲国产日韩欧美在线图片| 欧美三级电影精品| 国产亚洲成精品久久| 欧美日韩一区成人| 久久久亚洲精品一区二区三区| 欧美四级在线| 久久免费午夜影院| 久久久久国色av免费看影院| 国产欧美日韩一区二区三区在线观看| 欧美一区观看| 亚洲高清精品中出| 亚洲视频播放| 欧美三级视频在线播放| 亚洲一区区二区| 欧美视频中文字幕在线| 亚洲高清在线观看一区| 欧美一级久久久| 免费成人黄色片| 亚洲欧美在线免费| 欧美中文字幕在线| 欧美黑人一区二区三区| 国产精品xxxav免费视频| 久久久久成人精品免费播放动漫| 欧美一区亚洲一区| 欧美视频日韩| 亚洲欧美日韩电影| 日韩视频一区| 欧美一区二区视频网站| 久久久99精品免费观看不卡| 亚洲精品国产精品国产自|