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

    基于專家PID的帶臂四旋翼無人機控制方法

    2022-08-24 06:31:00陳寶周祖鵬衛(wèi)歡呂延釗睢志成
    計算機應用 2022年8期
    關鍵詞:角速度旋翼擾動

    陳寶,周祖鵬,衛(wèi)歡,呂延釗,睢志成

    (桂林電子科技大學機電工程學院,廣西桂林 541004)

    0 引言

    近年來,隨著科學技術的發(fā)展,導航和控制技術不斷成熟,飛行機器人不斷應用于航空拍攝、災難現場勘查、環(huán)境建模[1]及抓捕罪犯。一種具有作業(yè)能力的帶臂無人機(Unmanned Aerial Vehicle,UAV)成為了新的研究熱點[2],帶機械臂的四旋翼無人機主要由四旋翼無人機和機械臂組成,并具有與環(huán)境交互的能力,通過與環(huán)境的“接觸式”應用,帶臂無人機具有廣泛的應用前景,如危險目標的抓取、飛行目標物的捕獲及抓取物品運送到指定地點等。

    面向潛在的應用,國內外研究機構設計出多種飛行機械臂系統(tǒng),德雷塞爾大學實驗室將四旋翼無人機配備機械臂并對機械臂進行動力學建模,以及搭建了測試平臺[3]。文獻[4]設計了一種應用于林冠采樣的航空機械手,用于解決在森林冠層復雜環(huán)境中的樣本采集,并借助回路中的運動捕捉以補償慣性傳感器的不足,使無人機和機械臂的移動抓取更加精準;文獻[5]提出了一種解決空中操縱器碰撞問題的方法,該方法能夠在空中操縱器碰撞期間和碰撞后保持穩(wěn)定;文獻[6]介紹了一種大工作空間并聯航空機械手的設計,該機械手能夠在側面以及空中機器人下方進行物理交互,以實現空中機械臂的無奇異工作空間和高靈巧度;文獻[7]設計了一種帶機械臂無人直升機控制系統(tǒng),由無人直升機和機械臂組成,并使機械臂拔出固定在地面上的直桿;文獻[8]介紹了四旋翼無人機和單自由度機械臂組成的飛行機械臂控制系統(tǒng),并使機械臂能夠抓取靜止目標;文獻[9]利用一種單自由度抓手或繩索實現目標物體的空中運輸;文獻[10]介紹了機械臂直升機系統(tǒng)的建模和分析,并設計了一種線性二次規(guī)劃控制器(Linear Quadratic Regulator,LQR),同時對系統(tǒng)閉環(huán)性能進行仿真實驗;文獻[11]針對帶臂飛行器飛行過程中的動態(tài)穩(wěn)定性,提出了一種組合控制法,將增益調度跟模型參考自適應控制(Model Reference Adaptive Control,MRAC)有機融合;文獻[12]設計了一種帶臂無人機的控制器,實現其穩(wěn)定飛行,該控制器采用基于反饋線性化和雙積分器相結合進行控制;文獻[13]設計了自適應魯棒阻抗控制器;文獻[14]為飛行機械臂系統(tǒng)設計了模型預測控制器,面向末端執(zhí)行器的軌跡進行跟蹤控制;文獻[15]建立了漂浮在空間中剛體之間接觸的運動方程,對于非合作目標的抓取,引入了阻抗控制,當機械臂與被動目標碰撞時,通過阻抗控制匹配的方法進行有效的調節(jié)。

    本文針對帶二自由度的機械臂無人機控制系統(tǒng)進行設計,利用拉格朗日-歐拉方程對二自由度機械臂的運動學和動力學特性進行全狀態(tài)建模,并通過設計專家PID 控制器對帶臂無人機進行穩(wěn)定性控制,且與常規(guī)的PID 控制器進行對比。通過Matlab/Simulink 對帶臂無人機進行仿真,以驗證帶臂無人機在位置控制和姿態(tài)控制及二自由度機械臂在受干擾下的運動穩(wěn)定情況,其實物如圖1 所示。

    圖1 帶臂四旋翼無人機系統(tǒng)Fig.1 System of quadrotor UAV with manipulator

    1 帶臂無人機系統(tǒng)模型

    1.1 帶臂無人機運動學模型

    帶臂無人機控制系統(tǒng)主要是機械臂與無人機以鉸接方式進行連接的。通過拉格朗日-歐拉方程分析系統(tǒng)的內部運動關系,從機械臂末端的位姿及受力方式向前計算得到機械臂在無人機上的受力情況,進而分析無人機的位姿,簡易模型[16]如圖2 所示。

    由圖2 模型坐標系可推導該系統(tǒng)的運動模型。其中:Σb表示機載坐標系;Σe表示機械臂末端坐標系;ΣI表示大地坐標系;Pb=[x y z]T∈R3和Θb=[φ θ ψ]T∈R3表示無人機在慣性框架中的位置和歐拉角,φ是繞Xb軸的滾轉角(roll),θ是繞Yb軸的俯仰角(pitch),ψ是繞Zb軸的偏航角(yaw);Q=[q1q2]T表示機械臂關節(jié)角度。所有的廣義矢量的坐標變量可表示為L=[PTΘTQT]T∈R8×1。

    圖2 帶臂無人機的坐標系模型Fig.2 Coordinate system model of UAV with manipulator

    本文的運動學是旋翼平臺在慣性坐標系I 下的位姿變換以及機械臂的坐標系E 在慣性系統(tǒng)下位置與姿態(tài)方面的表達。而坐標系B 與大地坐標系I 之間的旋轉變換矩陣為Rb。

    1.2 帶臂無人機動力學模型

    帶臂無人機通過四個升力Fi(i=1,2,3,4)提升無人機并抓取作業(yè),機械臂的交互使無人機的位姿發(fā)生改變,即無人機與機械臂形成整個控制系統(tǒng)。對帶臂無人機動力學建模是根據機械臂角度和無人機位置及偏航角得出帶臂無人機的力矩,并通過動力學的逆問題得出帶臂無人機的實際位姿及角度。

    帶臂無人機動力學模型是通過拉格朗日-歐拉動力學方程分析系統(tǒng)內部的關系與運動情況。主要通過拉格朗日-歐拉方程[16-17]進行計算,拉格朗日-歐拉方程表達式為:

    其中:κ=[x z θ q1q2]T;u主要由無人機在X和Z方向的力矩uf和無人機力矩uλ及機械臂力矩uτ組成;D是慣性矩陣,D∈R5×5;C是離心力和哥氏力構成的矩陣,C∈R5×1;G是重力項構成的矩陣,G∈R5×1。D、C、G表示無人機的X和Z方向的位置和俯仰角及機械臂通過拉格朗日-歐拉方程編寫的動力學結合矩陣。Y方向由無人機的總推力和總質量得到,偏航角及滾轉角取決于力矩和轉動慣量。

    而帶臂四旋翼無人機在運動過程中四個旋翼會產生升力和旋轉力矩,關系可表示為:

    為了實現對帶臂無人機的閉環(huán)控制需對帶臂四旋翼無人機動力學進行正逆運算,正動力學主要是為了計算出無人機位姿[x z θ]及二自由度機械臂的兩個轉角qi(i=1,2)。則方程為:

    2 專家PID控制器設計

    專家PID 控制主要利用受控對象和控制規(guī)律的知識進行控制,對被控制的對象不需要精確的模型,使用專家經驗對系統(tǒng)模型進行設計控制。對于系統(tǒng)模型的控制專家PID具有靈活性、適應性和魯棒性,可根據系統(tǒng)的工作狀態(tài)及誤差情況靈活地選擇相應的控制規(guī)律去控制,并根據專家知識和經驗,能動性地調整控制器的參數,適應對象特性及環(huán)境的變化,通過利用專家規(guī)則,控制系統(tǒng)模型在非線性、大偏差下可靠地工作,使帶臂旋翼無人機在專家PID 控制律下能夠可靠地使旋翼無人機穩(wěn)定飛行及平穩(wěn)地抓取物品。

    專家PID 主要由五個控制律組成[18],通過工作狀態(tài)及誤差去選擇相應的控制律去穩(wěn)定數據,并對參數進行調節(jié)達到對控制系統(tǒng)穩(wěn)定的作用。控制結構如圖3 所示。

    圖3 帶臂無人機的控制系統(tǒng)結構Fig.3 Structure of control system of UAV with manipulator

    專家PID 控制算法如下:

    對期望值和實際值進行誤差分析,令e(k)表示當前的誤差值,e(k-1)、e(k-2)分別表示前一個和前兩個采樣的時刻誤差,表達式分別為:

    設計專家控制規(guī)律,根據誤差及其變化進行定性分析,控制規(guī)律如下。

    1)當|e(k)|>M1時:誤差絕對值很大,應按控制器的輸出最大(或最?。┬问捷敵?,以達到調整誤差的作用,使誤差的最大速度絕對值減小,可以達到快速使系統(tǒng)控制穩(wěn)定的作用。其控制器輸出為:

    其中:c表示系數。

    2)當e(k)Δe(k) >0 時:

    ①如|e(k)|≥M2,應按控制器較大的輸出進行控制,以使絕對值的誤差往減小的方向變化,達到控制的目的。其控制輸出為:

    ②如|e(k)|<M2,表示誤差往大的方向變化但誤差絕對值變化不明顯,可使用一般控制。其控制器輸出為:

    3)當e(k)Δe(k) <0,Δe(k)Δe(k-1) >0 且e(k)=0 時,表示誤差值不是很大,或處于平衡狀態(tài),其控制器輸出可保持不變,輸出為:

    其中:u(k)表示在上一次輸出過程中輸出值很小,可利用上次的輸出值作為這一次的輸出。

    4)當e(k)Δe(k) <0 且Δe(k)Δe(k-1) <0 時:

    ①如|e(k)|≥M2時,表示誤差較大處于極限狀態(tài),應使控制器輸出較大,達到控制的目的。其控制器輸出為:

    ②如|e(k)|<M2時,表示誤差比較小,可使控制器輸出較小,其控制器輸出為:

    5)當|e(k)|≤ε時,表示誤差的絕對值很小,應加入積分作為控制輸出,以減小穩(wěn)態(tài)誤差。控制器輸出為:

    其中:a、b是系數;em(k)為誤差e的第k個絕對值;u(k)為第k次控制器的輸出;u(k-1)為第k-1 次控制器的輸出;k1為增益放大系數,k1>1;k2為抑制系數,0 <k2<1;M1、M2為誤差界限,M1>M2>0;k為周期序數;ε為任意小的正實數。通過對帶臂無人機進行控制器設計,專家PID 能夠很好地對無人機進行穩(wěn)定控制,具有很好的魯棒性。

    3 軌跡規(guī)劃

    機械臂的軌跡規(guī)劃主要以點到點的運動方式進行,即規(guī)劃一條從初始位置q(t0)到最終位置q(tf)的軌跡,機械臂的軌跡有利于抓取物品的流暢性及規(guī)律性。本文采用五次多項式進行軌跡規(guī)劃,在t0初始時刻和tf最終時刻的角度、角速度、角加速度滿足如下:

    五次多項式為:

    其中:t=tf-t0;ai(i=0,1,2,3,4,5)為系數。

    根據初始位置和最終位置的角度、角速度和角加速度可得方程為:

    通過方程組(21)~(22)可得ai(i=0,1,2,3,4,5)系數,即機械臂關節(jié)期望的角度、角速度和角加速度可為:

    根據式(23)可對帶臂無人機的機械臂進行軌跡規(guī)劃。

    4 仿真實驗

    為了驗證在專家PID 控制下帶臂四旋翼無人機的穩(wěn)定性和優(yōu)越性,使用Matlab/Simulink 對帶臂無人機在仿真環(huán)境下進行實驗,帶臂無人機在擾動的情況下仍然處于穩(wěn)定,以驗證算法的可行性。

    經實際測量可知,帶臂四旋翼無人機系統(tǒng)的各項參數和指標[19]如表1 所示。

    表1 帶臂四旋翼無人機系統(tǒng)仿真參數Tab.1 Simulation parameters of quadrotor UAV system with manipulator

    表1 中:Ct和Cq分別是推力系數和扭矩系數;mv是無人機的質量;mi(i=1,2)分別是基座和各關節(jié)的質量;d是轉子到無人機中心的距離;H是無人機中心到地面的距離;Li(i=1,2)分別表示機械臂的臂長;Ixx、Iyy、Izz分別是無人機的轉動慣量。

    對帶機械臂無人機的位置和姿態(tài)進行仿真,以驗證算法的可行性和穩(wěn)定性,在仿真之前設置專家PID 控制器的數值,采樣時間ts=0.001,誤差界限M1=0.8,M2=0.000 5,增益放大系數K1=20,抑制系數K2=0.6,任意小的正實數ε=0.000 1,系數a=0.5,b=0.01,c=0.2 等參數,則帶機械臂無人機在設計參數的專家PID 控制器下的仿真性能如圖4~10 所示。

    從圖4 表示的帶臂無人機在專家PID 和常規(guī)PID 控制器下位置對比曲線可以看出:X方向專家PID 的收斂速度很快在t=5 s 時就已經達到穩(wěn)定而且超調量很小,而常規(guī)PID 在t=22 s 時才達到穩(wěn)定,超調量要更大一些;Y方向專家PID 在t=4.5 s 達到穩(wěn)定,常規(guī)PID 在t=30.5 s 達到穩(wěn)定;Z方向專家PID 在t=5.0 s 達到穩(wěn)定,常規(guī)PID 在t=19.0 s 達到穩(wěn)定。相較于常規(guī)PID,專家PID 收斂更快、穩(wěn)定性更好。圖5 表示帶臂旋翼無人機在偏航角ψ方面的控制,帶臂無人機以π/6 的角度偏航,專家PID 控制器在t=0.9 s 時就達到穩(wěn)定,而常規(guī)PID 在t=5.7 s 才達到穩(wěn)定。圖6 表示帶臂無人機二自由度機械臂關節(jié)角度對比曲線,專家PID 對機械臂的控制比常規(guī)PID 控制更好,機械臂關節(jié)收斂性更好、超調量更小。

    圖4 位置對比曲線Fig.4 Curve comparison of position

    圖5 偏航角ψ對比曲線Fig.5 Curve comparison of yaw angle ψ

    圖6 機械臂關節(jié)角度對比曲線Fig.6 Curve comparison of joint angles of manipulator

    為了驗證專家PID 對擾動的抑制效果,使用階躍信號和正弦波信號分別模擬平均風和陣風的擾動。平均風表示風速不隨時間而變化的風,陣風表示風速隨時間而變化的風,其風速有突發(fā)的特性。用于描述在兩種風的干擾下,帶臂無人機的穩(wěn)定性及收斂性。以偏航角ψ為例,在t=2 s 至t=6 s 的階躍信號為w(t)=1(2 <t<6);在t=4 s 至t=24 s的正弦波信號為w(t)=0.6sin[π(t-t1)/2] (4 <t<24),其中:t1為起始時間,t為干擾時間。設偏航角ψd=0 rad 表示懸停效果,實驗結果如圖7 所示。

    圖7(a)為階躍信號擾動下兩種控制器估計的結果,可以看出,在t=2 s 有突發(fā)擾動,而在擾動開始時專家PID 雖然有一點超調,但影響較小,且很快收斂到真實值,并保持穩(wěn)定;當t=6 s 擾動消失時專家PID 快速恢復到穩(wěn)定值。對于常規(guī)PID,當出現擾動時超調比較大持續(xù)時間長,擾動消失時不能快速收斂到穩(wěn)定值。圖7(b)為正弦波信號擾動下兩種控制器估計對比曲線,可以看出,t=4 s 時有擾動,而在擾動時專家PID 對真實擾動的跟隨性及準確性都比較好,偏差性較??;在t=24 s 擾動消失時雖然有一點超調但影響比較小,而常規(guī)PID 相對于真實擾動有一定的偏差,且擾動消失時超調比較大。

    圖7 外力干擾對比曲線Fig.7 Curve comparison of external interference

    由圖8 所示的帶臂無人機在兩個控制器控制下軌跡曲線對比可知:專家PID 控制的帶臂無人機開始起飛后就快速地到達了期望的目標軌跡,離期望的目標軌跡偏差比較小,最終能到達目標位置;而常規(guī)PID 控制的帶臂無人機開始起飛后需要飛一段時間才慢慢地接近期望目標軌跡,且離期望的目標軌跡有很大的偏差,不能到達最終的目標位置。即專家PID 優(yōu)于常規(guī)PID,收斂性更好。

    圖8 帶臂無人機軌跡曲線對比Fig.8 Comparison of trajectory curves of UAV with manipulator

    而在實際的應用中,可以通過分層控制分別控制無人機和機械臂,上層控制無人機,下層控制機械臂,無人機帶動機械臂到達指定位置,通過機械臂的運動執(zhí)行任務,即機械臂在五次多項式軌跡規(guī)劃下穩(wěn)定到達指定位置。設關節(jié)1 和關節(jié)2 的初始角度q1和q2分別為0,目標位置角度分別為π/6和π/4。機械臂關節(jié)角度如圖9 所示。

    圖9 機械臂角度軌跡曲線對比Fig.9 Curve comparison of manipulator angle trajectories

    由圖9 機械臂角度軌跡曲線知,帶臂無人機的機械臂在五次多項式規(guī)劃下到達指定位置,且在專家PID 控制下機械臂的軌跡跟隨比常規(guī)PID 更好,跟蹤誤差較小。由圖10 機械臂角速度可知,機械臂第1 個關節(jié)的初始角速度為0,然后速度慢慢增大到最高點,并在第5 秒時達到了最大值為0.098 rad/s,最后角速度逐漸減小到0;機械臂第2 個關節(jié)的初始角速度為0,然后慢慢增大到最高點,并在第5 秒時達到了最大值為0.147 rad/s,最后角速度逐漸減小到0??梢钥闯鰴C械臂在專家PID 的控制下,關節(jié)角速度平穩(wěn)地從0 到最大值再逐漸減小到0,顯示了機械臂能夠快速平穩(wěn)地到達指定位置并執(zhí)行任務。

    圖10 機械臂角速度Fig.10 Manipulator angular velocity

    5 結語

    本文對二自由度的帶臂無人機系統(tǒng)全狀態(tài)的運動學和動力學模型進行探討,并在此基礎上對帶臂無人機控制系統(tǒng)進行設計及算法研究,分析帶臂無人機的運動狀態(tài)穩(wěn)定性,并提出了利用專家PID 控制器智能整定帶臂無人機的穩(wěn)定參數,與常規(guī)PID 進行對比,專家PID 收斂速度更快、穩(wěn)定性更好、超調量更小,大幅提高了帶臂無人機的位姿跟蹤效果,使帶臂無人機在抓取物品時能更精準操作。仿真結果表明帶臂旋翼無人機在專家PID 的控制下有更好的抗擾性和魯棒性,改善性能更好,具有很好的優(yōu)越性和全局穩(wěn)定性,有很好的實用價值。

    猜你喜歡
    角速度旋翼擾動
    Bernoulli泛函上典則酉對合的擾動
    改進型自抗擾四旋翼無人機控制系統(tǒng)設計與實現
    大載重長航時油動多旋翼無人機
    (h)性質及其擾動
    基于STM32的四旋翼飛行器的設計
    電子制作(2019年9期)2019-05-30 09:41:48
    圓周運動角速度測量方法賞析
    小噪聲擾動的二維擴散的極大似然估計
    四旋翼無人機動態(tài)面控制
    半捷聯雷達導引頭視線角速度提取
    用于光伏MPPT中的模糊控制占空比擾動法
    電源技術(2015年11期)2015-08-22 08:50:38
    十分钟在线观看高清视频www| 2021少妇久久久久久久久久久| 天堂俺去俺来也www色官网| 国产色爽女视频免费观看| 久久av网站| 满18在线观看网站| 国产成人91sexporn| 日本av手机在线免费观看| 国产精品久久久久久精品电影小说| 97在线人人人人妻| 国产av精品麻豆| 少妇人妻 视频| 9色porny在线观看| 97精品久久久久久久久久精品| 国产又爽黄色视频| 国产精品一区二区在线观看99| 蜜臀久久99精品久久宅男| av免费在线看不卡| 久久国内精品自在自线图片| 国产一区有黄有色的免费视频| 成人亚洲精品一区在线观看| 欧美日韩精品成人综合77777| 观看美女的网站| 亚洲精品成人av观看孕妇| 秋霞伦理黄片| 午夜日本视频在线| 色哟哟·www| av不卡在线播放| 国产乱人偷精品视频| 国产成人精品福利久久| 国产麻豆69| 亚洲精品av麻豆狂野| 久久久精品区二区三区| 久久精品国产亚洲av涩爱| 日产精品乱码卡一卡2卡三| 午夜久久久在线观看| 在现免费观看毛片| 久久精品aⅴ一区二区三区四区 | 日韩av在线免费看完整版不卡| 亚洲国产欧美在线一区| 丝袜脚勾引网站| 国产又爽黄色视频| 亚洲精品aⅴ在线观看| 精品久久蜜臀av无| 婷婷色综合www| 亚洲图色成人| 中文字幕人妻丝袜制服| 午夜激情久久久久久久| 最近的中文字幕免费完整| 免费人成在线观看视频色| 天堂中文最新版在线下载| 美女国产视频在线观看| 一本大道久久a久久精品| 国产日韩欧美亚洲二区| 一区二区三区四区激情视频| 蜜桃国产av成人99| 大香蕉久久网| 欧美日韩精品成人综合77777| av国产精品久久久久影院| 国产成人欧美| 黄色一级大片看看| 搡女人真爽免费视频火全软件| 中文字幕av电影在线播放| 亚洲精品456在线播放app| 亚洲国产精品999| 一区二区三区精品91| 亚洲欧洲精品一区二区精品久久久 | 免费看光身美女| 成人黄色视频免费在线看| 人人妻人人添人人爽欧美一区卜| 999精品在线视频| 久久久久久久久久成人| 亚洲av成人精品一二三区| 丝袜美足系列| 亚洲国产毛片av蜜桃av| 你懂的网址亚洲精品在线观看| 我的女老师完整版在线观看| 国产精品久久久久久精品电影小说| 久久精品熟女亚洲av麻豆精品| 国产免费一级a男人的天堂| 亚洲内射少妇av| 亚洲欧美中文字幕日韩二区| 国产不卡av网站在线观看| 亚洲精品日韩在线中文字幕| av国产精品久久久久影院| 天天躁夜夜躁狠狠躁躁| 久久女婷五月综合色啪小说| 极品人妻少妇av视频| 精品少妇内射三级| 亚洲第一av免费看| 交换朋友夫妻互换小说| www.熟女人妻精品国产 | 国产欧美日韩综合在线一区二区| 国语对白做爰xxxⅹ性视频网站| 亚洲少妇的诱惑av| 国产免费一级a男人的天堂| √禁漫天堂资源中文www| 少妇 在线观看| 欧美老熟妇乱子伦牲交| 草草在线视频免费看| 嫩草影院入口| 久久精品熟女亚洲av麻豆精品| 久久久久久久久久久久大奶| 国产综合精华液| 亚洲精品中文字幕在线视频| 一级毛片黄色毛片免费观看视频| √禁漫天堂资源中文www| 韩国av在线不卡| 亚洲久久久国产精品| 观看av在线不卡| 亚洲精品av麻豆狂野| 久久精品国产a三级三级三级| 黑丝袜美女国产一区| 国产成人免费观看mmmm| 国产精品99久久99久久久不卡 | 日本欧美视频一区| 蜜臀久久99精品久久宅男| 久久久久久久精品精品| 日韩制服丝袜自拍偷拍| 成人毛片a级毛片在线播放| 一级毛片电影观看| 亚洲精品国产av蜜桃| 久久女婷五月综合色啪小说| 亚洲内射少妇av| 麻豆精品久久久久久蜜桃| 免费看光身美女| 极品人妻少妇av视频| 国产亚洲精品久久久com| 亚洲三级黄色毛片| 国产一区二区三区av在线| 中文字幕av电影在线播放| 男女午夜视频在线观看 | 欧美丝袜亚洲另类| 最近中文字幕2019免费版| 国产男人的电影天堂91| 水蜜桃什么品种好| 国产免费现黄频在线看| xxx大片免费视频| 狂野欧美激情性bbbbbb| 精品酒店卫生间| 最近最新中文字幕免费大全7| 咕卡用的链子| 午夜福利乱码中文字幕| 婷婷色综合大香蕉| 亚洲婷婷狠狠爱综合网| 日韩一区二区视频免费看| 免费人妻精品一区二区三区视频| 亚洲 欧美一区二区三区| videos熟女内射| 男女下面插进去视频免费观看 | 久久久久久久久久人人人人人人| 欧美国产精品一级二级三级| 久久久久久久大尺度免费视频| 国产av码专区亚洲av| 久久久国产一区二区| 国产日韩欧美视频二区| 日韩,欧美,国产一区二区三区| 亚洲国产av影院在线观看| 精品一区二区三区视频在线| 欧美激情极品国产一区二区三区 | 宅男免费午夜| 日韩中字成人| 欧美xxxx性猛交bbbb| 99国产精品免费福利视频| 一区二区av电影网| 在线免费观看不下载黄p国产| 欧美人与善性xxx| 一区二区av电影网| 久久精品夜色国产| 夫妻性生交免费视频一级片| 久久影院123| 精品99又大又爽又粗少妇毛片| 日本黄大片高清| 国产日韩一区二区三区精品不卡| 伊人久久国产一区二区| 亚洲精品久久久久久婷婷小说| 中文字幕av电影在线播放| 亚洲国产毛片av蜜桃av| 欧美日韩精品成人综合77777| 亚洲欧美清纯卡通| 99国产精品免费福利视频| 女人精品久久久久毛片| 久久99热6这里只有精品| 国产白丝娇喘喷水9色精品| 一级片'在线观看视频| 久久精品久久久久久久性| 免费av不卡在线播放| 国产成人精品婷婷| 中文字幕亚洲精品专区| 母亲3免费完整高清在线观看 | 春色校园在线视频观看| 免费大片黄手机在线观看| 久久久久人妻精品一区果冻| 亚洲第一av免费看| 少妇人妻 视频| 国产精品99久久99久久久不卡 | 熟女av电影| 国产成人91sexporn| 免费黄色在线免费观看| 久久人妻熟女aⅴ| 青青草视频在线视频观看| 久久久a久久爽久久v久久| 免费高清在线观看视频在线观看| 成人无遮挡网站| 考比视频在线观看| 日本vs欧美在线观看视频| 狂野欧美激情性bbbbbb| 久久久国产精品麻豆| 色哟哟·www| 久久久久久人妻| 中文天堂在线官网| 亚洲欧美清纯卡通| 看十八女毛片水多多多| 国产综合精华液| 国产一区有黄有色的免费视频| 中文欧美无线码| 国语对白做爰xxxⅹ性视频网站| 国产精品不卡视频一区二区| 亚洲国产欧美在线一区| 高清在线视频一区二区三区| 日韩中字成人| 欧美精品亚洲一区二区| 肉色欧美久久久久久久蜜桃| 亚洲 欧美一区二区三区| 男女无遮挡免费网站观看| av福利片在线| 日韩一区二区三区影片| 亚洲av欧美aⅴ国产| 国产免费视频播放在线视频| a级片在线免费高清观看视频| 国产成人免费观看mmmm| 中文字幕免费在线视频6| 日韩人妻精品一区2区三区| 热re99久久国产66热| 日韩一区二区视频免费看| 亚洲美女视频黄频| 亚洲av综合色区一区| 永久网站在线| 美女内射精品一级片tv| 亚洲精品美女久久av网站| 久久人人爽人人爽人人片va| 丝瓜视频免费看黄片| 男女边摸边吃奶| 丝袜美足系列| 岛国毛片在线播放| 丰满迷人的少妇在线观看| 新久久久久国产一级毛片| 黄片无遮挡物在线观看| 少妇熟女欧美另类| av网站免费在线观看视频| 嫩草影院入口| 国产伦理片在线播放av一区| 欧美国产精品一级二级三级| 巨乳人妻的诱惑在线观看| 国产极品粉嫩免费观看在线| 久久久精品区二区三区| 午夜日本视频在线| 我要看黄色一级片免费的| 久久精品国产亚洲av涩爱| 国产激情久久老熟女| 老熟女久久久| 免费观看a级毛片全部| a级毛片黄视频| 人人妻人人澡人人爽人人夜夜| 男女高潮啪啪啪动态图| 久久狼人影院| 99久久中文字幕三级久久日本| 国产成人欧美| 欧美精品一区二区免费开放| 午夜激情av网站| 美女脱内裤让男人舔精品视频| 亚洲内射少妇av| 秋霞在线观看毛片| 亚洲第一区二区三区不卡| 日本爱情动作片www.在线观看| 考比视频在线观看| 国产精品三级大全| 日韩电影二区| 啦啦啦在线观看免费高清www| 久久精品国产亚洲av天美| 亚洲美女黄色视频免费看| 下体分泌物呈黄色| 永久网站在线| 午夜91福利影院| 日本黄色日本黄色录像| 黄片播放在线免费| 赤兔流量卡办理| 精品久久蜜臀av无| 国产精品麻豆人妻色哟哟久久| 99香蕉大伊视频| 侵犯人妻中文字幕一二三四区| 熟女电影av网| 久久午夜福利片| 欧美激情国产日韩精品一区| 人人澡人人妻人| h视频一区二区三区| 一本色道久久久久久精品综合| 久久久久久久国产电影| 日本色播在线视频| 亚洲国产欧美日韩在线播放| 波多野结衣一区麻豆| 亚洲五月色婷婷综合| 午夜福利在线观看免费完整高清在| 久久国产亚洲av麻豆专区| 日韩成人伦理影院| 亚洲欧洲精品一区二区精品久久久 | 成人影院久久| 国产av一区二区精品久久| 国产欧美另类精品又又久久亚洲欧美| 精品亚洲成a人片在线观看| 又大又黄又爽视频免费| 久久久久久久大尺度免费视频| 在线观看免费日韩欧美大片| 性色av一级| 国产永久视频网站| 国产成人aa在线观看| 高清黄色对白视频在线免费看| 美女内射精品一级片tv| 久久久久久久久久人人人人人人| 咕卡用的链子| 熟女电影av网| 我要看黄色一级片免费的| 国产男女超爽视频在线观看| 嫩草影院入口| 大陆偷拍与自拍| 黄色毛片三级朝国网站| 少妇高潮的动态图| 丝袜美足系列| 美女脱内裤让男人舔精品视频| 热re99久久国产66热| 欧美激情极品国产一区二区三区 | 中文字幕另类日韩欧美亚洲嫩草| 性色avwww在线观看| 人人妻人人澡人人爽人人夜夜| 亚洲欧美日韩卡通动漫| 免费大片黄手机在线观看| 黄色怎么调成土黄色| 亚洲国产精品999| 国产亚洲最大av| 久久狼人影院| 欧美国产精品一级二级三级| 免费看不卡的av| 黑人猛操日本美女一级片| 精品国产一区二区三区四区第35| 中文精品一卡2卡3卡4更新| 国产高清三级在线| 黑人高潮一二区| 精品人妻一区二区三区麻豆| 一本—道久久a久久精品蜜桃钙片| 乱人伦中国视频| 日韩伦理黄色片| 纯流量卡能插随身wifi吗| 亚洲精品,欧美精品| 人妻系列 视频| 18在线观看网站| 在线免费观看不下载黄p国产| 一级,二级,三级黄色视频| 国产老妇伦熟女老妇高清| 永久网站在线| 国产精品一区www在线观看| 纵有疾风起免费观看全集完整版| 一二三四在线观看免费中文在 | 免费日韩欧美在线观看| 最近最新中文字幕大全免费视频 | 成人毛片60女人毛片免费| 咕卡用的链子| 性高湖久久久久久久久免费观看| 毛片一级片免费看久久久久| 欧美最新免费一区二区三区| 最近手机中文字幕大全| 亚洲欧美精品自产自拍| 欧美日韩av久久| 亚洲精品aⅴ在线观看| 欧美日韩av久久| a级毛片在线看网站| 国内精品宾馆在线| 免费女性裸体啪啪无遮挡网站| 一级毛片我不卡| 肉色欧美久久久久久久蜜桃| 下体分泌物呈黄色| 欧美日韩成人在线一区二区| 黄色 视频免费看| 91aial.com中文字幕在线观看| 婷婷色av中文字幕| 久久久久精品久久久久真实原创| 高清毛片免费看| 亚洲国产av影院在线观看| 日日啪夜夜爽| 精品国产一区二区久久| 丝袜喷水一区| 少妇的逼好多水| 两个人免费观看高清视频| 一级a做视频免费观看| 成人午夜精彩视频在线观看| 美女xxoo啪啪120秒动态图| 亚洲,欧美精品.| 欧美性感艳星| www.av在线官网国产| 亚洲精品av麻豆狂野| 成人亚洲精品一区在线观看| 欧美人与性动交α欧美软件 | 韩国高清视频一区二区三区| 久久精品国产自在天天线| 丝瓜视频免费看黄片| 亚洲av男天堂| 亚洲第一av免费看| av一本久久久久| 亚洲性久久影院| 亚洲,一卡二卡三卡| 精品酒店卫生间| videossex国产| 亚洲第一区二区三区不卡| 亚洲国产精品成人久久小说| 一区二区三区精品91| 国产精品久久久久久久久免| 亚洲国产精品专区欧美| 久久狼人影院| 美女xxoo啪啪120秒动态图| 国产在线视频一区二区| 美女内射精品一级片tv| 男女免费视频国产| 国产乱来视频区| 黄片无遮挡物在线观看| 久久精品国产亚洲av天美| 中国三级夫妇交换| 两个人免费观看高清视频| 成人漫画全彩无遮挡| 成年人免费黄色播放视频| 久久国产亚洲av麻豆专区| 99国产综合亚洲精品| 老司机亚洲免费影院| 交换朋友夫妻互换小说| 性色av一级| 久久久久精品久久久久真实原创| 男女午夜视频在线观看 | 精品少妇黑人巨大在线播放| 毛片一级片免费看久久久久| 日本黄大片高清| 免费黄色在线免费观看| 99香蕉大伊视频| 亚洲少妇的诱惑av| 国产亚洲午夜精品一区二区久久| 欧美激情国产日韩精品一区| 99九九在线精品视频| 高清黄色对白视频在线免费看| 最近的中文字幕免费完整| 午夜av观看不卡| 熟妇人妻不卡中文字幕| 精品久久蜜臀av无| 中国美白少妇内射xxxbb| 日日摸夜夜添夜夜爱| 亚洲成人av在线免费| av不卡在线播放| 视频在线观看一区二区三区| 国产一级毛片在线| 少妇人妻久久综合中文| 久久狼人影院| 日韩成人av中文字幕在线观看| videosex国产| 精品一区二区三区视频在线| 欧美精品一区二区免费开放| 97在线人人人人妻| 一级毛片电影观看| 少妇高潮的动态图| 久久久久网色| 国产精品偷伦视频观看了| 亚洲美女搞黄在线观看| 高清av免费在线| 人人妻人人爽人人添夜夜欢视频| 日韩av免费高清视频| 久久这里有精品视频免费| 亚洲情色 制服丝袜| 久久精品人人爽人人爽视色| 在线 av 中文字幕| 制服诱惑二区| 涩涩av久久男人的天堂| av一本久久久久| 欧美激情极品国产一区二区三区 | 国产极品天堂在线| 亚洲精品国产av蜜桃| 亚洲精品色激情综合| 亚洲欧美成人精品一区二区| 99久久综合免费| 国产免费一级a男人的天堂| 激情视频va一区二区三区| 精品一区二区三区四区五区乱码 | 免费不卡的大黄色大毛片视频在线观看| 秋霞伦理黄片| av国产精品久久久久影院| 你懂的网址亚洲精品在线观看| 侵犯人妻中文字幕一二三四区| 少妇被粗大猛烈的视频| 精品国产一区二区三区四区第35| 色视频在线一区二区三区| 插逼视频在线观看| 麻豆精品久久久久久蜜桃| 精品国产一区二区久久| 成人午夜精彩视频在线观看| 少妇猛男粗大的猛烈进出视频| 寂寞人妻少妇视频99o| 欧美丝袜亚洲另类| 一本—道久久a久久精品蜜桃钙片| 国产精品人妻久久久久久| 热re99久久国产66热| 交换朋友夫妻互换小说| 纯流量卡能插随身wifi吗| 自拍欧美九色日韩亚洲蝌蚪91| 欧美老熟妇乱子伦牲交| 性色avwww在线观看| 婷婷色麻豆天堂久久| 天天躁夜夜躁狠狠躁躁| 亚洲欧洲日产国产| 免费大片黄手机在线观看| 五月天丁香电影| 91aial.com中文字幕在线观看| 国产亚洲精品第一综合不卡 | 看十八女毛片水多多多| 九色亚洲精品在线播放| 69精品国产乱码久久久| 欧美精品高潮呻吟av久久| 久久综合国产亚洲精品| 免费女性裸体啪啪无遮挡网站| 韩国高清视频一区二区三区| 18禁裸乳无遮挡动漫免费视频| 国产伦理片在线播放av一区| 日韩免费高清中文字幕av| 大香蕉久久成人网| av在线老鸭窝| 色视频在线一区二区三区| 亚洲第一区二区三区不卡| 另类精品久久| 精品午夜福利在线看| 久久久久精品人妻al黑| 精品一区二区三区视频在线| 久久精品久久久久久噜噜老黄| 久久精品熟女亚洲av麻豆精品| 国产成人av激情在线播放| 岛国毛片在线播放| 大片免费播放器 马上看| 成人亚洲欧美一区二区av| 99九九在线精品视频| 全区人妻精品视频| 夜夜骑夜夜射夜夜干| 欧美xxⅹ黑人| 国内精品宾馆在线| 欧美精品一区二区大全| 久久久久精品性色| 久久久久久久久久久久大奶| 国产精品久久久久成人av| 国产高清不卡午夜福利| 久久午夜福利片| av有码第一页| 精品少妇内射三级| 久久鲁丝午夜福利片| 亚洲精品av麻豆狂野| 国产一区二区激情短视频 | 成人毛片a级毛片在线播放| 春色校园在线视频观看| 亚洲国产精品一区二区三区在线| 国产亚洲精品久久久com| 国产在线免费精品| 久久精品国产综合久久久 | 丝袜脚勾引网站| 婷婷色综合大香蕉| 午夜视频国产福利| 老司机亚洲免费影院| 99香蕉大伊视频| 久久久精品免费免费高清| 男男h啪啪无遮挡| 免费高清在线观看视频在线观看| 亚洲国产毛片av蜜桃av| 亚洲,一卡二卡三卡| 日韩av不卡免费在线播放| 新久久久久国产一级毛片| 在线观看美女被高潮喷水网站| 高清黄色对白视频在线免费看| 18禁观看日本| 男女国产视频网站| 日本免费在线观看一区| 女人被躁到高潮嗷嗷叫费观| 欧美另类一区| 在线天堂最新版资源| 婷婷色av中文字幕| 亚洲一级一片aⅴ在线观看| 亚洲精品美女久久av网站| 九色亚洲精品在线播放| 我的女老师完整版在线观看| 欧美日本中文国产一区发布| 狂野欧美激情性bbbbbb| 成人国语在线视频| 岛国毛片在线播放| 一边摸一边做爽爽视频免费| 一本—道久久a久久精品蜜桃钙片| 久久鲁丝午夜福利片| 最近的中文字幕免费完整| 中国三级夫妇交换| 欧美激情极品国产一区二区三区 | 搡女人真爽免费视频火全软件| 日韩成人伦理影院| 美女国产高潮福利片在线看| 精品一区二区免费观看| 色5月婷婷丁香| a级毛片黄视频| 色婷婷av一区二区三区视频| 人妻一区二区av| 久久久久精品性色| 桃花免费在线播放| 精品人妻熟女毛片av久久网站| 日韩 亚洲 欧美在线| 免费不卡的大黄色大毛片视频在线观看| 黄色配什么色好看| 人妻系列 视频| 国产亚洲欧美精品永久| 亚洲av电影在线进入| 久久精品久久久久久久性|