999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

基于ZMP的多姿態多功能下肢輔助裝置設計及穩定性分析

2025-03-04 00:00:00周曉晨孫田雪張保旭賈丙琪趙彥峻
機械傳動 2025年2期
關鍵詞:結構設計

摘要:【目的】為解決下肢運動功能障礙患者出行困難的問題并促進下肢功能恢復,設計了一款多姿態多功能下肢輔助裝置,提供代步、輔助起立和下肢康復訓練功能。【方法】通過Ansys 軟件進行靜力學仿真,評估裝置在使用時的安全性;建立下肢外骨骼運動學及動力學模型,得到各關節的空間位姿和理論力矩,為后續樣機制作提供理論支持;構建人體-輔助裝置系統“零點力矩”(Zero-Moment Point, ZMP)位置模型,計算理論ZMP 點,并將虛擬樣機導入Adams 軟件中仿真,得到支撐足位置變化;通過比較ZMP 點與支撐足位置驗證裝置的穩定性并搭建實物樣機。【結果】對比仿真與試驗數據證實了樣機設計的合理性,為后續研究提供了依據。

關鍵詞:下肢外骨骼;結構設計;ZMP;仿真分析

中圖分類號:TH789 DOI:10. 16578/j. issn. 1004. 2539. 2025. 02. 019

0 引言

隨著醫療水平的提高,人口壽命不斷延長,我國已進入老齡化社會,面臨嚴峻的人口老齡化挑戰[1]。隨著年齡的增長,老年人肌肉流失,下肢運動功能逐漸減弱,行走能力變弱甚至喪失。腦中風疾病和交通事故也會造成患者下肢運動功能損傷[2-3],嚴重影響日常生活。現代康復醫學理論及實踐證明:使用專用康復訓練器械對關節進行持續緩慢的被動訓練,可幫助患者修復受損運動神經,恢復行走功能[4]。下肢輔助裝置在幫助患者恢復肢體功能方面具有先天的優勢,適合開發為助老助殘設備,從而能夠長時間、重復性地幫助患者進行康復訓練,解決傳統康復治療存在的一些弊端 [5-7]。

目前,國內外的下肢輔助裝置可分為移動式和固定式兩種。移動式下肢輔助裝置有ReWalk[8]、HAL[9]、艾動[10]和ExoMotus[11]等,固定式下肢輔助裝置有Lokomat[12]、Flexbot[13]等,這兩種裝置都存在明顯的缺點。患者使用移動式下肢輔助裝置進行康復訓練時,通常需要手杖進行輔助支撐以保持平衡,對患者上肢要求較高;固定式下肢輔助裝置雖然具有較多的主動自由度,但使用時患者軀干被束縛懸掛,活動受限[14],只能在固定范圍內重復枯燥地完成訓練動作。

針對上述兩種裝置存在的問題,本文提出一款輪椅與下肢康復訓練協同的下肢輔助裝置,具有為下肢運動功能障礙患者提供代步、輔助起立和康復訓練的功能。該裝置中的輪椅擴大了患者的活動范圍,可實現快速行走的功能;輔助起立機構既能幫助患者站立,又能在康復行走訓練時起支撐保護作用。為保證患者步行訓練時的穩定性,使用“零點力矩”(Zero-Moment Point, ZMP)理論對裝置穩定性進行了分析[15],建立下肢輔助裝置虛擬樣機模型,進行了康復行走訓練仿真。結果表明,ZMP點落在支撐足范圍內,驗證了下肢輔助裝置的穩定性。隨后,加工實物樣機、搭建樣機試驗平臺,進行了輔助起立和康復訓練試驗,驗證了下肢輔助裝置的合理性。

1 多姿態多功能下肢輔助裝置設計

下肢輔助裝置的結構設計是保證行走穩定性的基礎及關鍵。設計的機械結構如圖1所示,其主要由電動輪椅、輔助起立機構和下肢外骨骼3部分組成。

1. 1 輔助起立機構結構設計

為解決下肢運動功能障礙患者起立困難的問題,設計了輔助起立機構,既能輔助患者進行坐姿和站姿轉換,又能提供一定的支撐力,以保證康復訓練過程中的穩定性。

輔助起立機構的結構如圖2所示,由輔助起立電動機、導軌滑塊、支撐桿、扶手、包膠軸承、背板等構成。扶手的前端和后端分別通過包膠軸承和鉸鏈與支撐桿相連;扶手、支撐桿、輔助起立電動機和導軌滑塊構成剪叉機構,將電動機的水平運動轉化為扶手的垂直運動。

輔助起立機構要支撐人體上半身的質量,基于Ansys 軟件對其進行靜力學分析。設定人體質量為75 kg,選擇Q235管材作為材料,通過軟件計算得到圖3所示的應力分析圖。最大應力值為169. 8 MPa,低于材料的屈服強度,證明該結構設計符合強度要求,不會發生明顯變形。

1. 2 下肢外骨骼結構設計

為滿足運動功能障礙患者康復訓練的需求,在輪椅的基礎上增加了下肢外骨骼結構,可為患者提供額外動力,輔助患者進行下肢康復訓練。由于外骨骼部分和人體下肢直接相連,在設計時需充分考慮穿戴舒適性和人機協同的靈活性。結合生物力學特點,將復雜的下肢運動模型抽象為3個運動關節:髖關節、膝關節和踝關節[16]。

由于人體下肢運動主要發生在矢狀面,本文設計外骨骼時為簡化結構和外骨骼控制系統,將髖關節和膝關節設計為主動關節,保留矢狀面的屈曲/伸展自由度,由推桿電動機進行驅動。踝關節設計為被動關節,在小腿構件末端安裝桿端關節軸承,保留完整自由度。為兼顧坐姿和站姿,將髖關節和膝關節驅動電動機分別固定在輪椅框架和大腿桿上,髖關節運動角度屈/伸變化范圍為-100°~5°,膝關節運動角度屈/伸變化范圍為-80°~0°,符合關節自由度需求。

由《中國成年人人體尺寸》[17]可知,大腿和小腿長度分別占身高的25%和23. 8%。本文設計外骨骼時在大腿和小腿處增加彈簧尺寸快調機構,使其適用于身高在1. 5~1. 9 m 的患者。下肢外骨骼結構如圖4所示。

2 下肢輔助裝置運動/動力學及穩定性建模

2. 1 輔助起立機構運動學建模

患者起立動作由輔助起立機構、髖關節、膝關節和踝關節共同作用完成。為保證輔助起立過程與人體自由狀態下起立過程相匹配,需對輔助起立過程進行運動學建模。假設人體關于矢狀面完全對稱,選擇踝關節作為基準點,將人體同外骨骼簡化為四連桿系統,建立起立運動學模型,如圖5所示。圖中,lub、lt 和ls 分別表示人體軀干長度、大腿長度和小腿長度。

將踝關節、膝關節和髖關節分別記作M0、M1 和M2,如圖5所示,則膝關節和髖關節的坐標分別為

2. 2 下肢外骨骼運動/動力學建模

患者下肢和外骨骼通過綁帶相連接,由綁帶帶動下肢運動。為保證外骨骼運動規律與人體步態相近,需對下肢外骨骼進行運動學建模。由于下肢外骨骼左右對稱,本文以左側腿為研究對象,分析各關節在矢狀面的運動。以腰部中心點為基坐標系原點,根據D-H法建立模型坐標系,如圖6所示。各連桿的D-H參數如表1所示。

基于D-H法,下肢外骨骼各連桿變化矩陣為

式中,s、c分別為sin與cos的縮寫。

患者進行康復訓練時,要運動到預期位置,下肢外骨骼各關節必須輸出相應的驅動力矩,需對下肢外骨骼進行動力學建模。本文采用拉格朗日法進行建模,忽略內部作用力,簡化建模過程。人體下肢運動時,左、右腿動作基本相同,僅在時間上相差半個周期。因此,本文僅對左側腿進行建模。以髖關節旋轉中心為坐標原點,如圖7所示,踝關節為被動關節,不考慮其關節轉矩變化。

圖7中,l1、l2 分別為大腿和小腿的長度;la、lb分別為大腿和小腿質心到相應關節中心的距離;θa、θb 分別為髖關節和膝關節的旋轉角度;m1、m2 分別為大腿和小腿的質量;I1、I2 分別為繞質心旋轉的轉動慣量。

大腿的動能K1 和勢能P1 分別為

K1 =1/2I1θ?2a(7)

P1 = m1 glacos θa (8)

小腿質心的位置坐標(x2,y2)為

2. 3 下肢輔助裝置穩定性建模

目前,機器人行走穩定性的判定依據主要有零點力矩(Zero-Moment Point, ZMP)和重心地面投影點(Center of Gravity, COG)兩種[18]。當機器人步速低于10 s/步時,處于靜態步行狀態,機器人慣性可忽略不計,機器人重心在地面上的投影始終位于支撐板多邊形內,機器人保持穩定。當步速提高后,機器人處于動態步行狀態,慣性力增強,COG依據不再適用。因此, 使用ZMP 來判斷機器人系統的穩定性[19]。

ZMP作為判斷機器人穩定行走的依據,由Vuko?bratovic于1972年提出[20]。根據機器人所受合力計算機器人實際ZMP點(地面作用力的力矩水平分量為0的作用點),如圖8所示,ZMP點在人體行走的支撐區域內時,行走過程穩定。

當下肢輔助裝置輔助患者行走時,至少有1個足底與地面保持接觸。為了保證步行的穩定性,患者和輔助裝置組成的系統的ZMP點應在支撐多邊形內。利用ZMP理論分析系統動態穩定性時,需考慮人的參與對系統的影響。因此,ZMP在支撐多邊形內的位置須有一定余量。

假設下肢輔助裝置與患者動作一致,將兩者視為整體建立模型,系統支撐力、行走依靠的反作用力和足底與地面間的摩擦力如圖9所示。O 點為坐標原點;G 為重心;外骨骼足底和地面接觸面以四邊形表示;A點相對于坐標系原點的矢量記作r;將足底受的作用力分解為作用于A點的力矢量F和力矩矢量TA。

力矢量F 相對于原點的力矩T 為

T = r × F + TA (18)

系統總動量S 與地面支撐力的關系為

S?= mg + F (19)

式中, m 為系統質量; g 為重力加速度, g =[ 0 0 -g ]T。

系統總角動量L 與地面作用力矩的關系為

L? = G × mg + T (20)

式中,G = [ x y z ]T。

將式(18)和式(19)代入式(20), 得力矩矢量TA為

TA = L? - G × mg + (S?- Mg ) × r (21)

令式(21)中TA = 0, 得ZMP的位置為

將式(23)和式(24)代入式(22),得簡化狀態下的ZMP位置為

式中,Gc = [ xc yc zc ]T;x?c,z?c 分別為外骨骼運動時x 軸和z 軸方向的加速度。

3 ZMP虛擬樣機仿真及分析

建立下肢輔助裝置的虛擬樣機,進行康復行走訓練仿真,以判斷裝置在康復訓練時的穩定性。將該裝置的SolidWorks 模型導入到Adams 軟件中,并在軟件中定義各個零件的材料、構件之間的運動副以及與地面的接觸力[21-22],在髖關節、膝關節和踝關節處添加驅動函數。本文使用OpenSim軟件中以每秒1步行走所獲得的各關節運動角度數據作為關節驅動數據。如圖11所示,利用Cubic Fitting Method 對關節驅動添加Spline 樣條函數CUBSPL(time, 0,spline_,0)*1d,模型如圖12所示。

仿真步數設為3 000,得到的下肢輔助裝置運動過程如圖13所示。

由式(25)可知,若將步行輔助機構系統簡化為一個質點,簡化狀態下ZMP位置需使用重心的坐標。采用Adams軟件的后處理模塊導出系統的重心位置和重心加速度,如圖14所示。分析可知,x 軸方向是人體前進方向,在x 軸方向上,重心位置變化為持續上升曲線,符合人體前進時重心位置的變化;z 軸方向是人體站立方向,在z 軸方向上,重心位置變化為周期性變化曲線,變化范圍為1 180~1 240 mm,位于腰部附近,同時,z 軸上重心加速度趨近于0,僅在0. 5 s和1 s處出現突變,這是由于左右腿支撐交換導致的。

重心在y 軸方向上的位置變化如圖15所示。由于在y 軸方向上有輪椅作為支撐,因此,在y 軸方向上重心位置變化非常小,可忽略不計。在計算過程中,可通過比較左右足底的壓力確定支撐腿,以判斷支撐足位置。在模型中,支撐足腳跟到腳尖的x 坐標軸范圍為-100~150 mm。將式(25)計算出的ZMP位置與行走過程中的支撐足位置進行比較,結果如圖16所示。圖中,縱坐標表示ZMP點在x 軸方向上的位置變化;虛線表示支撐足隨時間的變化。由圖可知,ZMP點始終不超出支撐足的范圍,可以判定使用下肢輔助裝置行走時狀態穩定。

4 實物樣機試驗

為進一步驗證下肢輔助裝置設計的合理性和輔助運動的效果,在理論研究的基礎上搭建了下肢輔助裝置的樣機進行真人試驗,由身高180 cm、質量90 kg的試驗員進行試驗。代步是下肢輔助裝置的重要功能之一,主控制器根據傳感器確認下肢輔助裝置是否處于輪椅狀態,確認后接通電源。經測試,代步模式下輔助裝置可以承載人體進行前進、后退和轉彎等功能,前進速度不超過0. 5 m/s;由于下肢輔助裝置比傳統輪椅的尺寸大,故其最小回轉半徑較大,為1. 3 m。

輔助起立機構是實現從代步功能向康復訓練功能轉換的關鍵,需驗證其是否可以輔助人體完成從坐姿到站姿的轉換。試驗過程如圖17所示。輔助起立關節角度變化如圖18所示。

分析圖18可得,輔助起立時仿真過程和試驗過程的關節角度數據變化基本吻合,產生偏差的原因在于使用輔助起立機構站立時,軀干固定在輪椅靠背上,軀干不需要前移也可保證穩定,因此,髖關節和膝關節角度變化略有區別。

該裝置另一重要功能是輔助康復訓練。將控制系統所需的傳感器安裝在下肢輔助裝置上,調試控制系統及傳感器系統保證其正常運行,進行真人康復輔助訓練試驗,如圖19所示。

使用下肢輔助裝置對試驗員進行10次康復輔助訓練,每次為2 min,將髖關節和膝關節處傳感器采集到的數據進行處理分析,并對比下肢輔助裝置理想軌跡與實際軌跡,如圖20所示。

由圖20可得,實際康復訓練時,髖關節和膝關節角度變化與理想狀態下的趨勢一致,沒有出現劇烈抖動和機構干涉的情況。但實際軌跡與期望軌跡之間仍然存在一定的偏差,關節角度運動范圍有一定的收窄,這是由于推桿電動機無法做到編碼電動機的無級調速,只能按照一定的梯度進行調速且無法做到精確的速度控制,但整體上的誤差在可接受的范圍之內。

5 結論

1) 從功能多樣性角度出發,為下肢運動障礙患者設計了一款多姿態多功能下肢輔助裝置,解決了固定式和移動式下肢輔助裝置的弊端,既能實現坐姿和站姿的轉換,又為患者提供了代步、輔助起立和康復訓練的功能。

2) 對下肢輔助裝置進行靜力學分析,保證其結構符合強度要求。對外骨骼部分進行運動學和動力學建模,得到外骨骼各關節的空間位置及關節驅動力矩。建立人體-輔助裝置系統ZMP位置模型,得到簡化狀態下的ZMP位置計算式,通過使用Adams軟件進行運動仿真,得到該裝置運動時重心位置和加速度的變化,計算得知,ZMP點不超出支撐足的范圍,驗證了下肢輔助裝置的穩定性。

3) 搭建實物樣機對下肢輔助裝置設計及功能進行試驗驗證,結果表明,各項功能均符合設計預期效果,具有較高的穩定性和安全性。后續會進行更多的樣機試驗,進一步驗證下肢輔助裝置的可行性。

參考文獻

[1] 董燕,黃健元. 中國人口性別年齡結構研究成果述評[J]. 鹽城工學院學報(社會科學版),2020,33(3):50-55.

DONG Yan,HUANG Jianyuan. A review of studies on the genderagestructure of China’s population[J]. Journal of Yancheng Instituteof Technology(Social Science Edition),2020,33(3):50-55.

[2] COLEMAN E R,MOUDGAL R,LANG K,et al. Early rehabilitationafter stroke:a narrative review[J]. Current Atherosclerosis Reports,2017,19(12):1-12.

[3] 盧明毅,吳冰潔,許建文,等. 步行訓練在脊髓損傷致下肢運動功能障礙康復中的臨床應用進展[J]. 中國臨床新醫學,2020,13(1):99-102.

LU Mingyi,WU Bingjie,XU Jianwen,et al. Advances in researchon rehabilitation of lower limb dyskinesia caused by spinal cord injuryvia walking training[J]. Chinese Journal of New ClinicalMedicine,2020,13(1):99-102.

[4] 朱旭,林松. 髖關節康復訓練仿真與康復器械個性化參數配置[J]. 機電一體化,2018,24(增刊1):39-45.

ZHU Xu,LIN Song. Simulation of hip rehabilitation training andrehabilitation equipment personalized parameter configuration[J].Mechatronics,2018,24(Suppl. 1):39-45.

[5] 黃高,張偉民,MARCO Ceccarelli,等. 一種新的康復與代步外骨骼機器人研究[J]. 自動化學報,2016,42(12):1933-1942.

HUANG Gao,ZHANG Weimin,MARCO Ceccarelli,et al. Researchof a new rehabilitation and assisting robot[J]. Acta AutomaticaSinica,2016,42(12):1933-1942.

[6] 石男強,劉剛峰,鄭天驕,等. 下肢康復機器人的研究進展與臨床應用[J]. 信息與控制,2021,50(1):43-53.

SHI Nanqiang,LIU Gangfeng,ZHENG Tianjiao,et al. Researchprogress and clinical application of lower limb rehabilitation robot[J]. Information and Control,2021,50(1):43-53.

[7] 賈丙琪,畢文龍,魏笑,等. 一種多功能下肢外骨骼機器人的設計與仿真分析[J]. 機械傳動,2021,45(1):59-64.

JIA Bingqi,BI Wenlong,WEI Xiao,et al. Design and simulationanalysis of multi-purpose lower limb exoskeleton robot[J]. Journalof Mechanical Transmission,2021,45(1):59-64.

[8] AWAD L N,ESQUENAZI A,FRANCISCO G E,et al. The Re‐Walk ReStore? soft robotic exosuit:a multi-site clinical trial of thesafety,reliability,and feasibility of exosuit-augmented post-strokegait rehabilitation[J]. Journal of Neuroengineering and Rehabilitation,2020,17(1):1-80.

[9] SCZESNY-KAISER M,TROST R,AACH M,et al. A randomizedand controlled crossover study investigating the improvement ofwalking and posture functions in chronic stroke patients usingHAL exoskeleton-The HALESTRO study (HAL-exoskeletonSTROke Study)[J]. Frontiers in Neuroscience,2019,13:259.

[10] 李龍飛,朱凌云,茍向鋒. 可穿戴下肢外骨骼康復機器人研究現狀與發展趨勢[J]. 醫療衛生裝備,2019,40(12):89-97.

LI Longfei,ZHU Lingyun,GOU Xiangfeng. Current status and developmenttrend of wearable lower-limb exoskeleton rehabilitationrobot[J]. Chinese Medical Equipment Journal,2019,40(12):89-97.

[11] GE Y F,DAN Y P,WANG A H,et al. Lower limb rehabilitation robotcontrol based on human gait data and plantar reaction force[C]//2020 International Conference on Advanced MechatronicSystems(ICAMechS). IEEE,2020:286-289.

[12] NAM K Y,KIM H J,KWON B S,et al. Robot-assisted gait training(Lokomat) improves walking function and activity in peoplewith spinal cord injury:a systematic review[J]. Journal of Neuroengineeringand Rehabilitation,2017,14(1):1-13.

[13] ZHANG A,HENG I,ZIA F. Acquiring \"working\" skills and experiencethrough hands-on multidisciplinary design projects[C]//ASME International Mechanical Engineering Congress and Exposition.American Society of Mechanical Engineers,2012,45219:85-89.

[14] 史小華,王洪波,孫利,等. 外骨骼型下肢康復機器人結構設計與動力學分析[J]. 機械工程學報,2014,50(3):41-48.

SHI Xiaohua,WANG Hongbo,SUN Li,et al. Design and dynamicanalysis of an exoskeletal lower limbs rehabilitation robot[J]. Journalof Mechanical Engineering,2014,50(3):41-48.

[15] MOKHTARI M,TAGHIZADEH M,MAZARE M. Hybrid adaptiverobust control based on CPG and ZMP for a lower limb exoskeleton[J]. Robotica,2020,39(2):181-199.

[16] 汪步云,王月朋,梁藝,等. 下肢外骨骼助力機器人關節驅動設計及試驗分析[J]. 機械工程學報,2019,55(23):55-66.

WANG Buyun,WANG Yuepeng,LIANG Yi,et al. Design on articularmotion amp; servo driving with experimental analysis for lowerlimb exoskeleton robot[J]. Journal of Mechanical Engineering,2019,55(23):55-66.

[17] 全國人類工效學標準化技術委員會. 中國成年人人體尺寸:GB/T 10000—2023[S]. 北京:中國標準出版社,2023:1.

National Technical Committee for Standardization of Ergonomics.Human dimensions of Chinese adults:GB/T 10000—2023[S].Beijing:Standards Press of China,2023:1.

[18] KIM S,HIROTA K,NOZAKI T,et al. Human motion analysis andits application to walking stabilization with COG and ZMP[J].IEEE Transactions on Industrial Informatics,2018,14(11):5178-5186.

[19] 魏揚帆,周川,郭健,等. 基于CPG的四足機器人坡面穩定行走控制研究[J]. 控制工程,2021,28(6):1055-1060.

WEI Yangfan,ZHOU Chuan,GUO Jian,et al. CPG-based stablewalking control of the quadruped robot on the slope[J]. ControlEngineering of China,2021,28(6):1055-1060.

[20] VUKOBRATOVIC M,BOROVAC B. Zero-moment point—thirtyfiveyears of its life[J]. International Journal of Humanoid Robotics,2004,1(1):157-173.

[21] HE Y,LI N,WANG C,et al. Development of a novel autonomouslower extremity exoskeleton robot for walking assistance[J]. Frontiersof Information Technology amp; Electronic Engineering,2019,20:318-329.

[22] WEI X,BI W L,XU G X,et al. Human-machine resonance analysisof multi-position lower limb rehabilitation exoskeleton[J]. InternationalJournal of Computational and Engineering,2021,1(6):1-4.

基金項目:國家自然科學基金項目(51905319);國家自然科學基金青年基金項目(51505263)

猜你喜歡
結構設計
高層建筑連體結構設計
人防結構設計疏漏的思考
高層建筑結構設計分析探討
高層建筑結構設計分析探討
一種承壓防爆接線盒結構設計
防爆電機(2020年4期)2020-12-14 03:11:14
一種新型蒸汽發生器結構設計
中國核電(2017年1期)2017-05-17 06:10:04
一種多方向導光按鍵結構設計
電子制作(2016年1期)2016-11-07 08:43:05
一種新型衛生防疫車結構設計
專用汽車(2016年5期)2016-03-01 04:14:49
一種新型流動車管所車結構設計
專用汽車(2016年5期)2016-03-01 04:14:48
BIM結構設計應用
主站蜘蛛池模板: 一级在线毛片| 91成人免费观看| 91青青草视频在线观看的| 欧美色丁香| 永久天堂网Av| 999福利激情视频| 国产精品女熟高潮视频| 亚洲永久色| 在线高清亚洲精品二区| 97视频免费在线观看| 久久综合亚洲鲁鲁九月天| 国产内射一区亚洲| 亚洲高清中文字幕在线看不卡| 亚洲人成网站在线播放2019| 亚洲天堂区| 亚洲精品欧美重口| 无码 在线 在线| 喷潮白浆直流在线播放| 欧美.成人.综合在线| 色婷婷视频在线| 亚洲国产成人麻豆精品| 国产福利微拍精品一区二区| 亚洲日韩久久综合中文字幕| 日韩欧美视频第一区在线观看| 91精品伊人久久大香线蕉| 欧美精品亚洲精品日韩专区va| 国产欧美专区在线观看| 在线精品亚洲国产| 精品剧情v国产在线观看| 亚洲女同欧美在线| 亚洲欧美在线精品一区二区| 国产小视频免费观看| 国产精品亚欧美一区二区| a色毛片免费视频| 亚洲第一精品福利| 久久国产热| 国产综合在线观看视频| 男人天堂亚洲天堂| 日韩国产一区二区三区无码| 尤物视频一区| 国模沟沟一区二区三区| 成人韩免费网站| 91久久国产综合精品女同我| 性喷潮久久久久久久久| 永久在线播放| 免费无码一区二区| 国产第一页屁屁影院| 免费在线色| 亚洲一区第一页| 99激情网| 国产区在线看| 国产黄在线观看| jizz国产在线| 日韩国产亚洲一区二区在线观看| 亚洲国产日韩在线观看| 98超碰在线观看| 色综合热无码热国产| a天堂视频| 色偷偷男人的天堂亚洲av| 在线视频精品一区| 91久久偷偷做嫩草影院电| 久久精品人人做人人| av手机版在线播放| 成人免费午夜视频| 亚洲一区二区日韩欧美gif| 成人免费黄色小视频| 欧美中文字幕在线二区| 国产精品毛片在线直播完整版| 最新国产精品第1页| 午夜日b视频| 国产天天射| 老司机aⅴ在线精品导航| 夜夜高潮夜夜爽国产伦精品| 免费看一级毛片波多结衣| 日韩AV手机在线观看蜜芽| 潮喷在线无码白浆| 99久久免费精品特色大片| 亚洲成AV人手机在线观看网站| 亚洲最新地址| 中文无码伦av中文字幕| 亚洲精品无码不卡在线播放| www.99在线观看|