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

    基于遺傳算法的七自由度機(jī)器人軌跡規(guī)劃

    2016-12-05 09:59:57馬丹妮李傳江張自強(qiáng)
    實(shí)驗(yàn)室研究與探索 2016年9期
    關(guān)鍵詞:適應(yīng)度遺傳算法軌跡

    馬丹妮, 李傳江, 張自強(qiáng)

    (1. 廣東省電信規(guī)劃設(shè)計(jì)院有限公司, 廣東 廣州 528000; 2. 上海師范大學(xué) 信息機(jī)電學(xué)院, 上海 201418)

    ?

    基于遺傳算法的七自由度機(jī)器人軌跡規(guī)劃

    馬丹妮1, 李傳江2, 張自強(qiáng)2

    (1. 廣東省電信規(guī)劃設(shè)計(jì)院有限公司, 廣東 廣州 528000; 2. 上海師范大學(xué) 信息機(jī)電學(xué)院, 上海 201418)

    針對(duì)七自由度串聯(lián)機(jī)器人Robai Cyton Gamma 300軌跡規(guī)劃問題,采用改進(jìn)遺傳算法規(guī)劃?rùn)C(jī)器人各關(guān)節(jié)的運(yùn)動(dòng)軌跡。利用D-H表示法建立起機(jī)器人末端執(zhí)行器的位姿與參考坐標(biāo)系之間的齊次變換矩陣,采用遺傳算法優(yōu)化BP神經(jīng)網(wǎng)絡(luò)求解機(jī)器人的運(yùn)動(dòng)學(xué)求逆解。利用5次B樣條曲線在關(guān)節(jié)空間構(gòu)造機(jī)器人各關(guān)節(jié)隨時(shí)間變化的運(yùn)動(dòng)軌跡。在滿足運(yùn)動(dòng)學(xué)約束條件下,對(duì)傳統(tǒng)遺傳算法在編碼方式、遺傳算子、交叉概率和變異概率等方面進(jìn)行改進(jìn),對(duì)機(jī)器人各關(guān)節(jié)運(yùn)動(dòng)軌跡進(jìn)行時(shí)間最優(yōu)規(guī)劃。運(yùn)用Matlab對(duì)研究進(jìn)行了仿真實(shí)驗(yàn)。結(jié)果表明,經(jīng)改進(jìn)遺傳算法優(yōu)化后的機(jī)器人運(yùn)動(dòng)軌跡時(shí)間明顯縮短,各關(guān)節(jié)的角速度、加速度和加加速度曲線連續(xù)無突變,從而驗(yàn)證了該方法的有效性。

    遺傳算法; Robai Cyton Gamma 300機(jī)器人; Bezier曲線; 軌跡優(yōu)化

    0 引 言

    機(jī)器人軌跡規(guī)劃是指機(jī)器人運(yùn)動(dòng)過程中的運(yùn)動(dòng)軌跡,即運(yùn)動(dòng)時(shí)各關(guān)節(jié)的位移、速度、加速度隨時(shí)間變化的曲線[1-2]。機(jī)器人的軌跡規(guī)劃是機(jī)器人控制的基礎(chǔ),近年來是機(jī)器人領(lǐng)域的研究熱點(diǎn)。

    軌跡規(guī)劃的性能指標(biāo)有很多種,主要分為:時(shí)間最優(yōu)軌跡規(guī)劃[3]、能量最優(yōu)軌跡規(guī)劃[4]、沖擊最優(yōu)軌跡規(guī)劃[5]以及混合最優(yōu)軌跡規(guī)劃[6]。其中,時(shí)間最優(yōu)軌跡規(guī)劃對(duì)提高機(jī)器人工作效率有重要意義,也一直是機(jī)器人軌跡規(guī)劃研究的重點(diǎn)。國(guó)內(nèi)外針對(duì)機(jī)器人的時(shí)間最優(yōu)軌跡規(guī)劃提出了許多不同的規(guī)劃方法,主要有二次規(guī)劃法[7]、動(dòng)態(tài)目標(biāo)規(guī)劃法[8]、PID控制法[9]、迭代法[10]等。

    本文以七自由度串聯(lián)機(jī)器人Robai Cyton Gamma 300為研究對(duì)象,將時(shí)間最短作為軌跡優(yōu)化目標(biāo),采用5次B樣條曲線進(jìn)行機(jī)器人各關(guān)節(jié)軌跡規(guī)劃,在滿足運(yùn)動(dòng)學(xué)參數(shù)的約束下,運(yùn)用改進(jìn)的遺傳算法實(shí)現(xiàn)對(duì)運(yùn)動(dòng)軌跡的時(shí)間優(yōu)化。實(shí)驗(yàn)結(jié)果表明,機(jī)器人運(yùn)動(dòng)總時(shí)間明顯減少,經(jīng)優(yōu)化后機(jī)器人各關(guān)節(jié)的運(yùn)動(dòng)軌跡連續(xù)、無突變,從而驗(yàn)證了該算法的有效性和實(shí)用性。

    1 Robai Cyton Gamma 300機(jī)器人的運(yùn)動(dòng)學(xué)分析

    Robai Cyton Gamma 300是由美國(guó)賓夕法尼亞州的諾百公司生產(chǎn)的一種創(chuàng)新的七自由度機(jī)器人,與人的手臂結(jié)構(gòu)相似,由7個(gè)轉(zhuǎn)動(dòng)關(guān)節(jié)(直流電機(jī))串聯(lián)而成。采用USB或RS485進(jìn)行通信,具體的技術(shù)參數(shù)如下:結(jié)構(gòu)形式為開環(huán)鏈?zhǔn)浇Y(jié)構(gòu),負(fù)載能力300 g,重復(fù)定位精度±0.5 mm,最大展開半徑48 cm,自身高度53.4 cm,自身質(zhì)量1.2 kg,驅(qū)動(dòng)方式為直流電機(jī)驅(qū)動(dòng),控制方式為PTP/CP,操作方式為示教操作,供電電源為12 V/2 A直流電源供電,工作環(huán)境溫度20~35 ℃。位移范圍:θ1=300°,θ2=210°,θ3=300°,θ4=210°,θ5=210°,θ6=210°,θ7=300°。利用機(jī)器人自帶的開發(fā)軟件Cyton Viewer采用D-H法建立機(jī)器人的連桿坐標(biāo)系如圖1所示,其D-H參數(shù)如表1所示。

    圖1 Cyton Viewer中建立的機(jī)器人D-H坐標(biāo)系

    由圖1知,機(jī)器人在幾何結(jié)構(gòu)上不滿足pieper準(zhǔn)則[11-12],無法獲得解析解,運(yùn)動(dòng)學(xué)方程的逆解是高度非線性方程組的求解問題。因此采用遺傳算法優(yōu)化BP神經(jīng)網(wǎng)絡(luò)的初始權(quán)值和閾值,使優(yōu)化后的BP神經(jīng)網(wǎng)絡(luò)能更好地求解機(jī)器人的運(yùn)動(dòng)學(xué)逆解[13]:① 將神經(jīng)網(wǎng)絡(luò)的輸入層、隱含層、輸出層中全部權(quán)值和閾值作為種群個(gè)體;② 把BP神經(jīng)網(wǎng)絡(luò)的預(yù)測(cè)輸出和期望輸出之間的誤差絕對(duì)值之和作為個(gè)體適應(yīng)度值;③ 選擇、交叉、變異等遺傳操作。

    表1 Robai Cyton Gamma 300機(jī)器人的D-H參數(shù)

    2 5次B樣條曲線構(gòu)造運(yùn)動(dòng)軌跡

    (1)

    由deBoor算法[14]可以推導(dǎo)出:

    (2)

    (3)

    式中:j=i-k+l,…,i。

    (4)

    (5)

    3 改進(jìn)GA時(shí)間最優(yōu)軌跡規(guī)劃

    3.1 運(yùn)動(dòng)學(xué)約束條件

    在關(guān)節(jié)空間中,考慮各關(guān)節(jié)所能達(dá)到的最大速度、最大加速度以及最大加加速度的限制,因此定義各關(guān)節(jié)速度、加速度、加加速度的約束條件:

    (6)

    式中,i=0,1,…,n。

    3.2 時(shí)間最優(yōu)軌跡規(guī)劃目標(biāo)的建立

    本文選擇機(jī)器人運(yùn)動(dòng)時(shí)間最短作為5次B樣條曲線軌跡優(yōu)化的目標(biāo)函數(shù):

    (7)

    式中:T為機(jī)器人沿B樣條曲線軌跡運(yùn)動(dòng)的總時(shí)間;hi為機(jī)器人沿每一段軌跡運(yùn)動(dòng)的時(shí)間;m為型值點(diǎn)的個(gè)數(shù)。

    3.3 改進(jìn)GA算法的軌跡優(yōu)化

    為了克服傳統(tǒng)遺傳算法的缺陷與不足,本文對(duì)其在編碼方式、遺傳算子、交叉概率和變異概率等方面進(jìn)行改進(jìn):

    (1) 采用實(shí)數(shù)編碼。直接采用解空間的形式進(jìn)行編碼,能明確地表達(dá)出問題本來的意義,易于引入特定領(lǐng)域的信息,而且能大大縮短串長(zhǎng),無需頻繁地編碼解碼,減少了遺傳算法的計(jì)算量,提高了算法的效率。

    (2) 采用排序選擇法。將每一代種群中各個(gè)個(gè)體的適應(yīng)度值進(jìn)行排序,在排序的基礎(chǔ)上進(jìn)行選擇操作,進(jìn)而更方便、高效地選擇出適應(yīng)度高的個(gè)體參與下一步的交叉、變異等遺傳操作。

    (3) 采用自適應(yīng)的交叉概率和變異概率。對(duì)于適應(yīng)度高于平均適應(yīng)度的個(gè)體,相對(duì)應(yīng)的交叉概率和變異概率降低;而對(duì)于適應(yīng)度低于平均適應(yīng)度的個(gè)體,相對(duì)應(yīng)的交叉概率和變異概率提高,這樣避免了算法陷入局部最優(yōu),從而提高了全局尋優(yōu)能力。具體的調(diào)整公式為[15]:

    (8)

    (9)

    式中:Pc1=0.9,Pc2=0.6,Pm1=0.1,Pm2=0.01;fmax為種群中最大適應(yīng)度值;favg為種群的平均適應(yīng)度值;f′為相互交叉兩個(gè)個(gè)體中較大的適應(yīng)度值;f為變異個(gè)體的適應(yīng)度值。

    (10)

    式中,ti為機(jī)器人運(yùn)動(dòng)到該節(jié)點(diǎn)處的時(shí)刻。

    實(shí)現(xiàn)步驟:

    (1) 編碼。 采用實(shí)數(shù)編碼方式將機(jī)器人沿每一段曲線運(yùn)動(dòng)所需時(shí)間編碼成遺傳算法所需要的染色體;

    (2) 種群初始化。在每段曲線運(yùn)動(dòng)時(shí)間hi的取值范圍內(nèi),隨機(jī)產(chǎn)生一定數(shù)量的個(gè)體;

    (3) 適應(yīng)度函數(shù)。

    (11)

    式中:hmax為給定范圍內(nèi)hi的最大值。該適應(yīng)度函數(shù)能保證不滿足約束條件的個(gè)體被賦予最小的適應(yīng)度值。

    (4) 選擇。對(duì)每代種群中的個(gè)體按適應(yīng)度值由大到小排序,適應(yīng)度值高的個(gè)體優(yōu)先進(jìn)行交叉、變異等遺傳操作。

    (12)

    式中:k為排序操作后個(gè)體適應(yīng)度值在序列中所處的位數(shù)。

    (5) 交叉。選用線性交叉算子。假設(shè)Xk、Yk為兩個(gè)父代個(gè)體,經(jīng)過算術(shù)交叉操作,產(chǎn)生了新的子代個(gè)體Xk+1、Yk+1:

    (13)

    (14)

    (15)

    4 仿真實(shí)驗(yàn)驗(yàn)證

    4.1 GA優(yōu)化BP神經(jīng)網(wǎng)絡(luò)求運(yùn)動(dòng)學(xué)逆解

    用向量θ=(θ1,θ2,θ3,θ4,θ5,θ6,θ7)表示Cyton Gamma 300機(jī)器人完成任意動(dòng)作的各關(guān)節(jié)變量θ,用向量X=(nx,ny,nz,ox,oy,oz,ax,ay,az,px,py,pz)表示機(jī)器人末端執(zhí)行器的位姿。在其運(yùn)動(dòng)空間內(nèi)任取400組末端執(zhí)行器位姿為輸入變量input,和與之對(duì)應(yīng)的各關(guān)節(jié)變量為輸出變量output,其中300組數(shù)據(jù)作為BP網(wǎng)絡(luò)的訓(xùn)練樣本,100組數(shù)據(jù)為預(yù)測(cè)樣本。訓(xùn)練網(wǎng)絡(luò)的輸入變量P為12×300維的矩陣,輸出變量T為7×300維的矩陣。BP網(wǎng)絡(luò)的輸入層神經(jīng)元個(gè)數(shù)為12,輸出層神經(jīng)元個(gè)數(shù)為7。根據(jù)隱層設(shè)計(jì)的經(jīng)驗(yàn)公式得到隱層神經(jīng)元個(gè)數(shù)為12。BP神經(jīng)網(wǎng)絡(luò)和GA優(yōu)化的BP神經(jīng)網(wǎng)絡(luò)的訓(xùn)練結(jié)果分別如圖3、圖4所示。

    圖3 BP網(wǎng)絡(luò)訓(xùn)練結(jié)果

    通過分析和訓(xùn)練結(jié)果可知,經(jīng)過GA優(yōu)化后的BP神經(jīng)網(wǎng)絡(luò)能得到更好的網(wǎng)絡(luò)初始權(quán)值和閾值,訓(xùn)練次數(shù)由221次減少到112次,收斂速度明顯提高,預(yù)測(cè)誤差由0.086 7減小到0.065 969,網(wǎng)絡(luò)預(yù)測(cè)誤差也減小了。實(shí)驗(yàn)證明,經(jīng)優(yōu)化的BP神經(jīng)網(wǎng)絡(luò)能夠更準(zhǔn)確、快速地求解機(jī)器人的各關(guān)節(jié)變量。

    圖4 經(jīng)遺傳算法優(yōu)化后訓(xùn)練結(jié)果

    4.2 改進(jìn)GA優(yōu)化機(jī)器人運(yùn)動(dòng)軌跡

    由表2中第6次優(yōu)化結(jié)果可得出,采用改進(jìn)遺傳算法,在滿足機(jī)器人運(yùn)動(dòng)學(xué)約束條件下,對(duì)5次B樣條各關(guān)節(jié)運(yùn)動(dòng)軌跡進(jìn)行時(shí)間優(yōu)化,總時(shí)間由初始狀態(tài)的18 s縮短到14.003 9 s,縮短了3.996 1 s。經(jīng)改進(jìn)遺傳算法優(yōu)化后的時(shí)間序列生成5次B樣條運(yùn)動(dòng)軌跡見圖6(圖中:s代表關(guān)節(jié)角位移,mm;v代表關(guān)節(jié)角速度,mm/s; a代表關(guān)節(jié)角加速度,mm/s2; J代表關(guān)節(jié)角加加速度,mm/s3),優(yōu)化后的各關(guān)節(jié)的位移、速度和加速度曲線光滑平穩(wěn),加加速度曲線連續(xù)無突變。

    5 結(jié) 語

    本文以七自由度串聯(lián)Robai Cyton Gamma 300機(jī)器人為研究對(duì)象,采用5次B樣條曲線構(gòu)造運(yùn)動(dòng)軌跡。通過對(duì)傳統(tǒng)遺傳算法在編碼方式、選擇算子、交叉概率和變異概率等方面作出改進(jìn),以運(yùn)動(dòng)時(shí)間最短為優(yōu)化目標(biāo),在Matlab中進(jìn)行軌跡優(yōu)化仿真實(shí)驗(yàn),結(jié)果表明機(jī)器人的運(yùn)動(dòng)總時(shí)間由初始狀態(tài)的20 s縮短到14.003 9 s,達(dá)到了運(yùn)動(dòng)軌跡時(shí)間優(yōu)化的目標(biāo)。經(jīng)優(yōu)化后的機(jī)器人各關(guān)節(jié)的角速度、加速度和加加速度變化曲線連續(xù)無突變,從而驗(yàn)證了該方法的可靠性和有效性。

    [1] Saeed B Niku. 機(jī)器人學(xué)導(dǎo)論——分析.控制及應(yīng)用[M]. 2版. 北京:電子工業(yè), 2013:136-241.

    [2] 張 濤. 機(jī)器人引論[M]. 北京:機(jī)械工業(yè)出版社, 2010:303-512.

    [3] Lin C C, Chen M Y. Intelligent mobile mechanical arm design[C]//System Science and Engineering (ICSSE), 2010 International Conference on. IEEE, 2010: 136-141

    [4] Quigley M, Asbeck A, Ng A. A low-cost compliant 7-dof robotic manipulator[C]//Robotics and Automation (ICRA), 2011 IEEE International Conference on. IEEE, 2011: 6051-6058.

    [5] Chettibi T, Lemoine P. Generation of Point to Point Trajectories for Robotic Manipulators Under Electro-Mechanical Constraints[J]. International Review of Mechanical Engineering, IREME, 2007(2):131-143.

    [6] Sakai S, Xu C, Ming A,etal. Motion planning for a high-speed manipulator with mechanical joint stops based on target dynamics and PCH system[C]// Intelligent Robots and Systems, 2009. IROS 2009. IEEE/RSJ International Conference on. IEEE, 2009:279-284.

    [7] Pchelkin S S, Shiriaev A S, Robertsson A,etal. Integrated time-optimal trajectory planning and control design for industrial robot manipulator[C]//Intelligent Robots and Systems (IROS), 2013 IEEE/RSJ International Conference on. IEEE,2013:2521 - 2526.

    [8] Ghasemi M H, Kashiri N, Dardel M. Time-optimal trajectory planning of robot manipulators in point-to-point motion using an indirect method[J]. Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science, 2012, 226:473-484.

    [9] Muller P A, Boucherit R, Liu S. Smooth and time-optimal trajectory planning for robot manipulators[C]//American Control Conference on. IEEE, 2012:5466 - 5471.

    [10] 劉松國(guó). 六自由度串聯(lián)機(jī)器人運(yùn)動(dòng)優(yōu)化與軌跡跟蹤控 制研究 [D]. 杭州: 浙江大學(xué), 2009.

    [11] 余 陽, 林 明, 林永才. 基于混合遺傳算法的工業(yè)機(jī)器人最優(yōu)軌跡規(guī)劃[J]. 計(jì)算機(jī)工程與設(shè)計(jì), 2012, 33(4):1574-1580.

    [12] 錢東海, 馬文羅, 汪建偉,等. 多約束條件下的機(jī)器人時(shí)間最優(yōu)軌跡規(guī)劃[J]. 制造業(yè)自動(dòng)化, 2011, 33(11):1-5.

    [13] Lee J D, Heshmati A. Solution of Inverse Kinematics for SCARA Manipulator Using Adaptive Neuro-Fuzzy Network[J]. International Journal on Soft Computing, 2011, 2(4):501-504.

    [14] Pérez Bailón W, Barrera Cardiel E, Juárez Campos I,etal.Mechanical Energy Optimization in Trajectory Planning for Six DOF Robot Manipulators Based on Eighth-Degree Polynomial Functions and a Genetic Algorithm[C]//Electrical Engineering,Computing Science and Automatic Control Conference on. IEEE,2010:466 - 451.

    [15] Frederik Debrouwere, Wannes Van Loock, Goele Pipeleers,etal.Time-Optimal Path Following for Robots with Trajectory Jerk Constraints using Sequential Convex Programming[C]//2013 IEEE International Conference on Robotics and Automation (ICRA) Karlsruhe,Germany,May 6-10, 2013:1916-1921.

    [16] Liu H, Lai X, Wu W. Time-optimal and jerk-continuous trajectory planning for robot manipulators with kinematic constraints[J]. Robotics and Computer-Integrated Manufacturing, 2013, 29(2):309-317.

    Seven Degree of Freedom Robot Trajectory Planning Scheme Based on Improved Genetic Algorithms

    MADan-ni1,LIChuan-jiang2,ZHANGZi-qiang2

    (1. Guangdong Planning and Designing Institute of Telecommunications CO.LTD, Guangzhou 528000, China; 2. College of Information, Mechanical and Electrical Engineering, Shanghai Normal University, Shanghai 201418, China)

    To solve trajectory planning problem of seven degree of freedom serial robot called Robai Cyton Gamma 300, an improved genetic algorithm was applied to plan robot trajectory of each joint. First, D-H notation was used to establish homogeneous transformation matrix between the actuator's posture at the end of the robot and the reference coordinate system. And genetic algorithm was used to optimize the BP neural network to solve the kinematics inverse solution of robot. Then, five B spline curve was applied in the joint space of each joint’s motion trajectory. Finally, under the condition of kinematic constraint, traditional genetic algorithm was improved on mode of coding, genetic operators, crossover probability and mutation probability, etc, to realize the time optimal planning. By conducting the simulation experiment on Matlab, the results showed the trajectory time of robot was optimized obviously by improved genetic algorithm. Angle velocity and the acceleration and jerk of each joint were continuous, it proved the effectiveness of the proposed method.

    genetic algorithms(GA); Robai Cyton Gamma 300 robot; Bezier(B) spline curve; trajectory optimization

    2015-12-02

    上海市科委基金資助項(xiàng)目(11510502400)

    馬丹妮(1990-),女,河南新縣人,碩士,主要研究方向?yàn)榍度胧较到y(tǒng)與自動(dòng)化裝置。

    Tel.:021-57122955; E-mail:nini502327489@126.com

    張自強(qiáng)(1958-),男,河南洛陽市人,教授,碩士生導(dǎo)師,主要研究方向?yàn)樽詣?dòng)化控制。

    Tel.:021-57122955; E-mail:zhzq@shnu.edu.cn

    TP 273

    A

    1006-7167(2016)09-0033-05

    猜你喜歡
    適應(yīng)度遺傳算法軌跡
    改進(jìn)的自適應(yīng)復(fù)制、交叉和突變遺傳算法
    軌跡
    軌跡
    軌跡
    基于自適應(yīng)遺傳算法的CSAMT一維反演
    一種基于遺傳算法的聚類分析方法在DNA序列比較中的應(yīng)用
    進(jìn)化的軌跡(一)——進(jìn)化,無盡的適應(yīng)
    基于遺傳算法和LS-SVM的財(cái)務(wù)危機(jī)預(yù)測(cè)
    基于空調(diào)導(dǎo)風(fēng)板成型工藝的Kriging模型適應(yīng)度研究
    基于改進(jìn)的遺傳算法的模糊聚類算法
    亚洲综合色惰| 欧美xxⅹ黑人| 日韩欧美国产在线观看| 免费观看的影片在线观看| 日韩av免费高清视频| 日韩一本色道免费dvd| 国产精品一区二区三区四区免费观看| 九九久久精品国产亚洲av麻豆| 欧美三级亚洲精品| 亚洲精华国产精华液的使用体验| 亚洲精品中文字幕在线视频 | 欧美潮喷喷水| 免费看日本二区| 一本久久精品| 91久久精品国产一区二区三区| 日韩亚洲欧美综合| 国产一区二区亚洲精品在线观看| 日本色播在线视频| 寂寞人妻少妇视频99o| 久久久久久久久久久丰满| 成人av在线播放网站| 嫩草影院新地址| 精品人妻偷拍中文字幕| 亚洲人成网站在线观看播放| 成人二区视频| 亚洲熟妇中文字幕五十中出| 成人美女网站在线观看视频| 久久人人爽人人片av| 欧美成人a在线观看| 欧美另类一区| 国产精品av视频在线免费观看| 日本-黄色视频高清免费观看| av在线亚洲专区| 十八禁国产超污无遮挡网站| 国产精品国产三级国产专区5o| av国产久精品久网站免费入址| 国产精品日韩av在线免费观看| 又大又黄又爽视频免费| 国产精品综合久久久久久久免费| 黄色配什么色好看| 可以在线观看毛片的网站| 中文字幕亚洲精品专区| 久99久视频精品免费| 亚洲国产精品成人综合色| 久久人人爽人人爽人人片va| 26uuu在线亚洲综合色| 一区二区三区高清视频在线| 看十八女毛片水多多多| 亚洲av中文字字幕乱码综合| 18禁动态无遮挡网站| 久久国内精品自在自线图片| 搡老乐熟女国产| 国产 亚洲一区二区三区 | 身体一侧抽搐| 人妻制服诱惑在线中文字幕| 80岁老熟妇乱子伦牲交| 久久久久久久久久久丰满| 麻豆精品久久久久久蜜桃| 日本一二三区视频观看| 一夜夜www| 色哟哟·www| 最近最新中文字幕免费大全7| 国产真实伦视频高清在线观看| 永久免费av网站大全| 成人亚洲精品一区在线观看 | 久久精品人妻少妇| 人体艺术视频欧美日本| 高清视频免费观看一区二区 | 久久久久久伊人网av| 国产午夜精品久久久久久一区二区三区| 看十八女毛片水多多多| 97超碰精品成人国产| 久久久久久久亚洲中文字幕| 蜜桃久久精品国产亚洲av| 国产色爽女视频免费观看| 成人亚洲精品av一区二区| 精品一区二区三区视频在线| 好男人视频免费观看在线| 国产精品女同一区二区软件| 亚洲熟妇中文字幕五十中出| 婷婷六月久久综合丁香| 只有这里有精品99| 一区二区三区高清视频在线| 九九在线视频观看精品| 黄色一级大片看看| 中文资源天堂在线| 国产淫语在线视频| 精品人妻偷拍中文字幕| 色网站视频免费| 亚洲欧美成人精品一区二区| 青春草视频在线免费观看| 日日啪夜夜爽| 国产v大片淫在线免费观看| 非洲黑人性xxxx精品又粗又长| 一本一本综合久久| 美女内射精品一级片tv| av在线播放精品| 一级毛片久久久久久久久女| 久久精品综合一区二区三区| 亚洲av免费在线观看| 美女国产视频在线观看| 少妇被粗大猛烈的视频| 日本与韩国留学比较| 久久久亚洲精品成人影院| 天堂中文最新版在线下载 | 如何舔出高潮| 色视频www国产| 亚洲精品成人av观看孕妇| 久久久久网色| 欧美日韩国产mv在线观看视频 | 直男gayav资源| 中文字幕人妻熟人妻熟丝袜美| 国产亚洲精品av在线| 亚洲欧美一区二区三区黑人 | 一区二区三区乱码不卡18| 毛片一级片免费看久久久久| 欧美日韩视频高清一区二区三区二| 一级毛片电影观看| 男女边吃奶边做爰视频| 日韩av免费高清视频| 亚洲自偷自拍三级| 久久久久久伊人网av| 亚洲怡红院男人天堂| 日韩精品有码人妻一区| 五月伊人婷婷丁香| 深夜a级毛片| 久久99蜜桃精品久久| 亚洲激情五月婷婷啪啪| 在线观看av片永久免费下载| 身体一侧抽搐| 三级国产精品欧美在线观看| 亚洲av日韩在线播放| 大话2 男鬼变身卡| 欧美日韩一区二区视频在线观看视频在线 | 美女脱内裤让男人舔精品视频| 女人久久www免费人成看片| 伊人久久精品亚洲午夜| 天天一区二区日本电影三级| 国产伦一二天堂av在线观看| 久久亚洲国产成人精品v| 看黄色毛片网站| 欧美变态另类bdsm刘玥| 天堂av国产一区二区熟女人妻| 国产伦精品一区二区三区四那| 黄色欧美视频在线观看| av黄色大香蕉| 老女人水多毛片| 日日啪夜夜爽| 精品人妻一区二区三区麻豆| 精品久久久久久久久久久久久| av黄色大香蕉| 一级爰片在线观看| 中文字幕久久专区| 欧美潮喷喷水| 国产又色又爽无遮挡免| 在线播放无遮挡| 欧美+日韩+精品| 一级毛片aaaaaa免费看小| 天堂网av新在线| 美女被艹到高潮喷水动态| 成人一区二区视频在线观看| 啦啦啦韩国在线观看视频| 啦啦啦韩国在线观看视频| kizo精华| 久久热精品热| 亚洲熟妇中文字幕五十中出| 亚洲av免费高清在线观看| 内地一区二区视频在线| 女人被狂操c到高潮| 日韩av在线大香蕉| 嫩草影院入口| 日韩av在线大香蕉| 一本久久精品| 成人av在线播放网站| 日韩一区二区三区影片| 成人二区视频| av线在线观看网站| 亚洲av在线观看美女高潮| 日韩一区二区视频免费看| 天堂av国产一区二区熟女人妻| 亚洲精品一二三| 国产一区二区三区综合在线观看 | 免费看不卡的av| 久久久精品免费免费高清| 欧美97在线视频| 小蜜桃在线观看免费完整版高清| 国产伦在线观看视频一区| 国产毛片a区久久久久| 精品国产露脸久久av麻豆 | 婷婷六月久久综合丁香| 日韩欧美三级三区| 色哟哟·www| 天美传媒精品一区二区| 国产av国产精品国产| 午夜免费激情av| 久久精品夜色国产| 超碰97精品在线观看| 久久精品国产鲁丝片午夜精品| 天堂网av新在线| 国产伦理片在线播放av一区| 欧美不卡视频在线免费观看| 亚洲成色77777| 精品一区二区三区人妻视频| 汤姆久久久久久久影院中文字幕 | 国产精品精品国产色婷婷| 国产黄色小视频在线观看| 亚洲国产最新在线播放| 在线 av 中文字幕| 国产精品国产三级国产专区5o| 国模一区二区三区四区视频| 久久久久久久午夜电影| 26uuu在线亚洲综合色| 1000部很黄的大片| 少妇熟女aⅴ在线视频| 真实男女啪啪啪动态图| 日韩一区二区三区影片| 欧美成人一区二区免费高清观看| 久久久久久久午夜电影| 久久精品国产自在天天线| 国产成人精品久久久久久| 亚洲国产日韩欧美精品在线观看| 欧美高清性xxxxhd video| 国产精品久久久久久精品电影小说 | 国产乱人视频| 女人久久www免费人成看片| 亚洲三级黄色毛片| 国产一级毛片在线| 亚洲精品乱码久久久v下载方式| 国产成人精品一,二区| 国产在线男女| 国产一区二区三区综合在线观看 | 国产女主播在线喷水免费视频网站 | 丰满人妻一区二区三区视频av| 一级毛片电影观看| 一个人观看的视频www高清免费观看| 久久国内精品自在自线图片| 99热全是精品| 欧美成人午夜免费资源| 青春草视频在线免费观看| 内射极品少妇av片p| 九九在线视频观看精品| 日韩一本色道免费dvd| 国产不卡一卡二| xxx大片免费视频| 超碰av人人做人人爽久久| 国模一区二区三区四区视频| 午夜激情欧美在线| 高清毛片免费看| 国产欧美日韩精品一区二区| 日韩,欧美,国产一区二区三区| 高清午夜精品一区二区三区| 久久99精品国语久久久| 99热网站在线观看| 国产成年人精品一区二区| 卡戴珊不雅视频在线播放| 日韩欧美精品免费久久| 国产精品美女特级片免费视频播放器| 久久精品国产亚洲网站| 别揉我奶头 嗯啊视频| 波多野结衣巨乳人妻| 免费av不卡在线播放| 男女啪啪激烈高潮av片| 日韩制服骚丝袜av| 亚洲欧洲日产国产| 性插视频无遮挡在线免费观看| 国产人妻一区二区三区在| 亚洲久久久久久中文字幕| 亚洲av中文字字幕乱码综合| 国产成人福利小说| 国产精品女同一区二区软件| 大香蕉97超碰在线| 欧美成人精品欧美一级黄| 大又大粗又爽又黄少妇毛片口| 国内精品一区二区在线观看| av网站免费在线观看视频 | 少妇裸体淫交视频免费看高清| 国产亚洲精品久久久com| 人妻系列 视频| 久久久久久久亚洲中文字幕| 九草在线视频观看| 在线观看av片永久免费下载| 国产精品久久久久久精品电影| 嫩草影院精品99| 男插女下体视频免费在线播放| 51国产日韩欧美| 日本欧美国产在线视频| 好男人视频免费观看在线| 91在线精品国自产拍蜜月| 男女国产视频网站| 国产亚洲91精品色在线| 我的女老师完整版在线观看| 男女下面进入的视频免费午夜| 国产真实伦视频高清在线观看| av网站免费在线观看视频 | av.在线天堂| 欧美日韩综合久久久久久| 免费av不卡在线播放| 国产白丝娇喘喷水9色精品| 国产一区二区三区综合在线观看 | 中文字幕久久专区| 色网站视频免费| 欧美日韩精品成人综合77777| 亚洲最大成人av| 精品欧美国产一区二区三| 高清在线视频一区二区三区| 日本熟妇午夜| 亚洲真实伦在线观看| 麻豆国产97在线/欧美| 国产在视频线精品| 丝瓜视频免费看黄片| 不卡视频在线观看欧美| 日韩精品青青久久久久久| 又粗又硬又长又爽又黄的视频| 日韩欧美精品免费久久| 99热这里只有精品一区| 99re6热这里在线精品视频| 午夜福利在线在线| 国产成人福利小说| 男人和女人高潮做爰伦理| 伊人久久国产一区二区| 伦理电影大哥的女人| 深夜a级毛片| 成人漫画全彩无遮挡| 国产黄色视频一区二区在线观看| 黄色欧美视频在线观看| 日本免费在线观看一区| 国产精品av视频在线免费观看| 亚洲欧美日韩卡通动漫| 亚洲成人一二三区av| 日韩在线高清观看一区二区三区| 一区二区三区高清视频在线| 亚洲av福利一区| 免费观看的影片在线观看| 精品人妻一区二区三区麻豆| 久久久久久久久久久免费av| 精品久久久久久久久亚洲| 免费看美女性在线毛片视频| 秋霞在线观看毛片| 少妇的逼水好多| 久久精品人妻少妇| 久久久久久国产a免费观看| 老司机影院成人| 国产大屁股一区二区在线视频| 在线a可以看的网站| 菩萨蛮人人尽说江南好唐韦庄| 国产成人a∨麻豆精品| 中文字幕免费在线视频6| 黄色配什么色好看| 国产不卡一卡二| 女人十人毛片免费观看3o分钟| 国产精品伦人一区二区| 久久久久久久久久成人| 日韩欧美精品v在线| 超碰av人人做人人爽久久| 欧美3d第一页| 十八禁网站网址无遮挡 | 国产午夜精品久久久久久一区二区三区| 久久久精品欧美日韩精品| 麻豆乱淫一区二区| 国产成人午夜福利电影在线观看| 日本黄色片子视频| 久久精品国产亚洲网站| 在线播放无遮挡| 国产三级在线视频| 久久人人爽人人爽人人片va| 精品久久久久久久久久久久久| 精品一区二区三卡| 精品国内亚洲2022精品成人| 成人毛片a级毛片在线播放| 97超视频在线观看视频| 久久久久九九精品影院| 国产大屁股一区二区在线视频| 一级毛片aaaaaa免费看小| 亚洲av日韩在线播放| 久久国产乱子免费精品| 色综合站精品国产| 一级黄片播放器| 伦精品一区二区三区| 日本色播在线视频| 91久久精品国产一区二区成人| 中文字幕制服av| 国产精品1区2区在线观看.| 亚洲av男天堂| 国产精品嫩草影院av在线观看| 成人欧美大片| 2022亚洲国产成人精品| 校园人妻丝袜中文字幕| 18禁动态无遮挡网站| 黄色日韩在线| 99热网站在线观看| 亚洲av免费高清在线观看| 午夜免费观看性视频| 神马国产精品三级电影在线观看| 深爱激情五月婷婷| 中文字幕制服av| 精品不卡国产一区二区三区| 亚洲自偷自拍三级| 国产黄色视频一区二区在线观看| 尾随美女入室| 又大又黄又爽视频免费| 丝袜美腿在线中文| 久久精品综合一区二区三区| 一二三四中文在线观看免费高清| 国产久久久一区二区三区| av.在线天堂| 国国产精品蜜臀av免费| 久久99热6这里只有精品| 岛国毛片在线播放| 高清日韩中文字幕在线| 午夜爱爱视频在线播放| 色5月婷婷丁香| 亚洲欧美成人综合另类久久久| 午夜免费激情av| 日韩电影二区| 国产精品福利在线免费观看| 一个人观看的视频www高清免费观看| 中文在线观看免费www的网站| 亚洲欧美精品自产自拍| 久久精品综合一区二区三区| 国产单亲对白刺激| 三级经典国产精品| 亚洲av福利一区| 天堂俺去俺来也www色官网 | 搞女人的毛片| 久久久久久伊人网av| 能在线免费观看的黄片| 久久精品国产亚洲av涩爱| 男女下面进入的视频免费午夜| 久久国内精品自在自线图片| 人人妻人人看人人澡| eeuss影院久久| 国产亚洲av嫩草精品影院| 中文字幕av在线有码专区| 美女大奶头视频| 18禁裸乳无遮挡免费网站照片| a级一级毛片免费在线观看| 哪个播放器可以免费观看大片| 精品久久久精品久久久| a级一级毛片免费在线观看| 亚洲精品成人久久久久久| 久久综合国产亚洲精品| 久久国内精品自在自线图片| 中文字幕av在线有码专区| 日韩伦理黄色片| 亚洲最大成人手机在线| 亚洲性久久影院| 一级毛片 在线播放| 少妇熟女欧美另类| 岛国毛片在线播放| 国产精品一区二区三区四区久久| 成年免费大片在线观看| 三级男女做爰猛烈吃奶摸视频| 欧美另类一区| 亚洲不卡免费看| 91久久精品电影网| 午夜免费激情av| 性插视频无遮挡在线免费观看| 欧美潮喷喷水| 男人舔女人下体高潮全视频| 六月丁香七月| 午夜免费男女啪啪视频观看| 久久这里只有精品中国| 最后的刺客免费高清国语| 国产淫语在线视频| 亚洲熟女精品中文字幕| 亚洲av免费高清在线观看| 欧美激情国产日韩精品一区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产一区二区三区av在线| 欧美日韩国产mv在线观看视频 | 日韩强制内射视频| 中文字幕免费在线视频6| 国产精品综合久久久久久久免费| av在线播放精品| 欧美97在线视频| 日韩伦理黄色片| 久久久久精品性色| 成人高潮视频无遮挡免费网站| 国产av不卡久久| 又爽又黄无遮挡网站| 亚洲国产高清在线一区二区三| 少妇裸体淫交视频免费看高清| 免费黄频网站在线观看国产| 亚洲乱码一区二区免费版| 春色校园在线视频观看| 亚洲性久久影院| 亚洲国产欧美在线一区| 国产精品蜜桃在线观看| 我的老师免费观看完整版| 丝袜美腿在线中文| 免费看不卡的av| 99热网站在线观看| 最近手机中文字幕大全| 一级爰片在线观看| av播播在线观看一区| 亚洲四区av| 日韩,欧美,国产一区二区三区| 97人妻精品一区二区三区麻豆| 精品久久久久久久末码| av专区在线播放| 国产乱人视频| 女人十人毛片免费观看3o分钟| 成人午夜精彩视频在线观看| 免费av毛片视频| 免费看av在线观看网站| 亚洲精品aⅴ在线观看| 午夜福利在线观看吧| 熟妇人妻久久中文字幕3abv| 一边亲一边摸免费视频| 精品99又大又爽又粗少妇毛片| 国产黄色免费在线视频| 国产亚洲精品久久久com| 色综合色国产| 欧美一区二区亚洲| 国产探花在线观看一区二区| 日韩制服骚丝袜av| 男人狂女人下面高潮的视频| 中文字幕av成人在线电影| 最近手机中文字幕大全| 少妇高潮的动态图| 久久综合国产亚洲精品| 国产伦精品一区二区三区四那| 免费电影在线观看免费观看| 欧美成人精品欧美一级黄| 男人舔奶头视频| 亚洲欧美清纯卡通| 久久久久久久久久久免费av| 深爱激情五月婷婷| 80岁老熟妇乱子伦牲交| 亚洲自偷自拍三级| 高清在线视频一区二区三区| 精品久久久久久久久亚洲| 麻豆精品久久久久久蜜桃| 三级国产精品欧美在线观看| av国产久精品久网站免费入址| 不卡视频在线观看欧美| 中国国产av一级| 国产免费又黄又爽又色| 日日撸夜夜添| 亚洲人成网站在线播| 亚洲欧美精品专区久久| 亚洲婷婷狠狠爱综合网| 国产成人freesex在线| 欧美性猛交╳xxx乱大交人| 久久久精品免费免费高清| 亚洲av福利一区| 在线观看人妻少妇| 国产一区二区在线观看日韩| 超碰97精品在线观看| 国产极品天堂在线| 亚洲成人久久爱视频| 亚洲av男天堂| 久久精品国产亚洲av天美| 美女大奶头视频| 国产亚洲一区二区精品| 久久热精品热| 亚洲国产色片| 亚洲av日韩在线播放| 国产亚洲av嫩草精品影院| 国产老妇女一区| 久久午夜福利片| 亚洲色图av天堂| 女人十人毛片免费观看3o分钟| 伦精品一区二区三区| 18禁在线播放成人免费| 成人午夜精彩视频在线观看| 床上黄色一级片| 久久精品熟女亚洲av麻豆精品 | 毛片一级片免费看久久久久| 黄色配什么色好看| 一区二区三区四区激情视频| 国产成人一区二区在线| 99九九线精品视频在线观看视频| 久久久久久久久中文| a级毛片免费高清观看在线播放| 美女被艹到高潮喷水动态| 亚洲内射少妇av| 可以在线观看毛片的网站| 高清av免费在线| 久久精品熟女亚洲av麻豆精品 | 最后的刺客免费高清国语| av天堂中文字幕网| 能在线免费看毛片的网站| 亚洲成人av在线免费| 欧美zozozo另类| 国产精品一区二区三区四区免费观看| av在线播放精品| 大香蕉久久网| 国产av国产精品国产| 日韩,欧美,国产一区二区三区| 最近中文字幕高清免费大全6| 天天一区二区日本电影三级| 午夜福利视频精品| 国产探花极品一区二区| 毛片一级片免费看久久久久| 人人妻人人澡人人爽人人夜夜 | 亚洲精品中文字幕在线视频 | 久久精品国产亚洲av涩爱| 久久99热这里只频精品6学生| 久久鲁丝午夜福利片| 丝袜美腿在线中文| 黄色日韩在线| 一区二区三区免费毛片| 国产精品久久久久久久电影| 禁无遮挡网站| 国产精品综合久久久久久久免费| 欧美xxxx性猛交bbbb| 免费看a级黄色片| 青春草国产在线视频| 深夜a级毛片| 超碰97精品在线观看| 国产不卡一卡二| 久久久久久久久久久丰满| 国产高清不卡午夜福利| 能在线免费看毛片的网站|