• <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)的遺傳算法的模糊聚類算法
    五月玫瑰六月丁香| 黄色怎么调成土黄色| 久久精品久久精品一区二区三区| 51国产日韩欧美| 中文在线观看免费www的网站| 国产欧美亚洲国产| 精品国产乱码久久久久久小说| av在线天堂中文字幕| 婷婷色综合www| 亚洲久久久久久中文字幕| 国产毛片在线视频| 亚洲欧美精品专区久久| 搡老乐熟女国产| 国产精品99久久99久久久不卡 | 国国产精品蜜臀av免费| 日韩 亚洲 欧美在线| 老师上课跳d突然被开到最大视频| 男女无遮挡免费网站观看| 日本欧美国产在线视频| 91午夜精品亚洲一区二区三区| 看十八女毛片水多多多| 可以在线观看毛片的网站| 亚洲精品,欧美精品| av在线蜜桃| 久久久久久久久久人人人人人人| 精品酒店卫生间| 亚洲av中文av极速乱| 又大又黄又爽视频免费| 国产成人福利小说| 婷婷色综合www| 老司机影院成人| 男人爽女人下面视频在线观看| 在线免费十八禁| 国产精品一二三区在线看| 蜜桃亚洲精品一区二区三区| 69人妻影院| 2021少妇久久久久久久久久久| 蜜桃亚洲精品一区二区三区| 波多野结衣巨乳人妻| 国产精品一区二区在线观看99| 国产精品一区二区三区四区免费观看| 婷婷色麻豆天堂久久| 亚洲第一区二区三区不卡| 午夜福利高清视频| 日韩不卡一区二区三区视频在线| 国产一区二区在线观看日韩| 日韩av在线免费看完整版不卡| 久久人人爽av亚洲精品天堂 | 亚洲国产欧美人成| 久久99热这里只频精品6学生| 2021少妇久久久久久久久久久| 久久精品久久久久久久性| 又爽又黄a免费视频| 97超碰精品成人国产| .国产精品久久| 91久久精品电影网| 久久久午夜欧美精品| 夫妻午夜视频| 国产精品一区www在线观看| 亚洲三级黄色毛片| 2022亚洲国产成人精品| av在线播放精品| 国产精品伦人一区二区| 久久女婷五月综合色啪小说 | 天堂俺去俺来也www色官网| 亚洲精品一二三| 久久久午夜欧美精品| 国产精品精品国产色婷婷| 亚洲精品成人av观看孕妇| 亚洲av中文字字幕乱码综合| 自拍偷自拍亚洲精品老妇| 欧美性猛交╳xxx乱大交人| 人妻系列 视频| 久久久精品免费免费高清| 久久久久久九九精品二区国产| 丰满人妻一区二区三区视频av| 久久久成人免费电影| 三级国产精品片| 最近的中文字幕免费完整| 色婷婷久久久亚洲欧美| 成年免费大片在线观看| 直男gayav资源| 国产欧美日韩精品一区二区| 欧美97在线视频| 国产日韩欧美在线精品| tube8黄色片| 午夜免费观看性视频| 黄色配什么色好看| 亚洲成色77777| 午夜福利高清视频| 亚洲第一区二区三区不卡| 日本黄大片高清| 亚洲真实伦在线观看| 国产精品一及| 欧美国产精品一级二级三级 | 成人综合一区亚洲| 一级av片app| 丝瓜视频免费看黄片| 内地一区二区视频在线| 欧美成人精品欧美一级黄| 亚洲在线观看片| 在线播放无遮挡| 九九在线视频观看精品| 国产乱人视频| 日韩欧美精品v在线| 97人妻精品一区二区三区麻豆| 亚洲精品乱久久久久久| 18禁在线播放成人免费| 午夜老司机福利剧场| 亚洲欧美中文字幕日韩二区| 亚洲精品乱码久久久v下载方式| 国产一区有黄有色的免费视频| 搞女人的毛片| 欧美日韩一区二区视频在线观看视频在线 | 午夜激情福利司机影院| 中文字幕免费在线视频6| 久久精品国产a三级三级三级| 国产精品人妻久久久久久| 国产成人精品久久久久久| 国产一区二区在线观看日韩| 777米奇影视久久| 少妇猛男粗大的猛烈进出视频 | 亚洲精品久久久久久婷婷小说| 精品视频人人做人人爽| 在线播放无遮挡| 天堂中文最新版在线下载 | 成人漫画全彩无遮挡| 在线 av 中文字幕| 国产在线一区二区三区精| 国产成人一区二区在线| 人妻夜夜爽99麻豆av| 一区二区av电影网| 国产精品女同一区二区软件| 成年人午夜在线观看视频| 直男gayav资源| 人妻制服诱惑在线中文字幕| 亚洲一区二区三区欧美精品 | 日本免费在线观看一区| 青青草视频在线视频观看| 免费黄网站久久成人精品| 全区人妻精品视频| 国产精品女同一区二区软件| 久久久a久久爽久久v久久| 日韩一本色道免费dvd| 超碰av人人做人人爽久久| 高清毛片免费看| 免费观看在线日韩| 亚洲av一区综合| 欧美成人午夜免费资源| 亚洲精品视频女| 久久6这里有精品| av黄色大香蕉| 免费高清在线观看视频在线观看| 国语对白做爰xxxⅹ性视频网站| 大又大粗又爽又黄少妇毛片口| 亚洲精品久久久久久婷婷小说| 国精品久久久久久国模美| 丝袜喷水一区| av国产免费在线观看| 国模一区二区三区四区视频| 男人添女人高潮全过程视频| 国产一区二区三区av在线| 97热精品久久久久久| av专区在线播放| 日本色播在线视频| 亚洲国产精品成人久久小说| 大香蕉久久网| 综合色av麻豆| 高清毛片免费看| 两个人的视频大全免费| 一个人观看的视频www高清免费观看| 午夜精品一区二区三区免费看| 蜜臀久久99精品久久宅男| 性插视频无遮挡在线免费观看| 国产日韩欧美在线精品| 在线观看一区二区三区| 久久久久精品性色| 亚洲国产日韩一区二区| 精品久久国产蜜桃| 精品亚洲乱码少妇综合久久| 久久99蜜桃精品久久| 大片免费播放器 马上看| 久久这里有精品视频免费| 欧美日韩视频精品一区| 久久国内精品自在自线图片| 伊人久久精品亚洲午夜| 国产男女内射视频| 别揉我奶头 嗯啊视频| 伦理电影大哥的女人| av国产免费在线观看| 夜夜骑夜夜射夜夜干| 热99久久久久精品小说推荐| 色吧在线观看| 我要看黄色一级片免费的| 天天躁狠狠躁夜夜躁狠狠躁| 18在线观看网站| av网站在线播放免费| 在线免费观看不下载黄p国产| 九草在线视频观看| 国产视频首页在线观看| 天天影视国产精品| 国产国语露脸激情在线看| 国产老妇伦熟女老妇高清| 久久精品国产亚洲av涩爱| 最新在线观看一区二区三区 | 九九爱精品视频在线观看| 国产女主播在线喷水免费视频网站| 午夜免费鲁丝| 80岁老熟妇乱子伦牲交| 七月丁香在线播放| 国产精品蜜桃在线观看| 男女边摸边吃奶| 男男h啪啪无遮挡| 久久久精品免费免费高清| 少妇被粗大猛烈的视频| 亚洲成人手机| 在线观看人妻少妇| 99精国产麻豆久久婷婷| 久久这里只有精品19| 亚洲精品久久午夜乱码| 建设人人有责人人尽责人人享有的| 熟女av电影| 街头女战士在线观看网站| 亚洲国产欧美在线一区| 美女扒开内裤让男人捅视频| 国产日韩一区二区三区精品不卡| 亚洲色图综合在线观看| 国产精品熟女久久久久浪| 在线观看人妻少妇| 成年美女黄网站色视频大全免费| 日韩伦理黄色片| 女人精品久久久久毛片| 国产激情久久老熟女| 制服人妻中文乱码| 久久久久久久精品精品| 精品少妇内射三级| 一区二区三区激情视频| 亚洲精品国产av蜜桃| 国产精品久久久久成人av| 欧美另类一区| 精品国产一区二区久久| 母亲3免费完整高清在线观看| 深夜精品福利| 国产一区有黄有色的免费视频| 久久久久精品人妻al黑| www.自偷自拍.com| 欧美日韩av久久| 亚洲精品久久午夜乱码| 中文精品一卡2卡3卡4更新| 中文欧美无线码| 亚洲人成77777在线视频| 国产精品香港三级国产av潘金莲 | 一边亲一边摸免费视频| 成人手机av| 免费女性裸体啪啪无遮挡网站| 亚洲精华国产精华液的使用体验| 亚洲伊人色综图| 青春草国产在线视频| 99久国产av精品国产电影| 久久99热这里只频精品6学生| 熟女少妇亚洲综合色aaa.| 国产麻豆69| av在线观看视频网站免费| 丝袜美腿诱惑在线| 亚洲在久久综合| 久热这里只有精品99| 精品亚洲乱码少妇综合久久| 亚洲成人一二三区av| 777久久人妻少妇嫩草av网站| 好男人视频免费观看在线| 欧美少妇被猛烈插入视频| 免费在线观看黄色视频的| 亚洲精品久久久久久婷婷小说| 国产成人精品久久久久久| 伦理电影大哥的女人| 午夜福利,免费看| 一区福利在线观看| 99久久人妻综合| 亚洲av综合色区一区| 欧美精品亚洲一区二区| 美国免费a级毛片| 久久久久久久精品精品| 欧美中文综合在线视频| 亚洲欧美成人综合另类久久久| videosex国产| 成年人免费黄色播放视频| 午夜福利,免费看| 赤兔流量卡办理| 亚洲av日韩在线播放| 亚洲精品国产色婷婷电影| 蜜桃国产av成人99| 天天躁夜夜躁狠狠躁躁| 亚洲精品中文字幕在线视频| 久热爱精品视频在线9| 五月天丁香电影| 亚洲自偷自拍图片 自拍| av国产久精品久网站免费入址| 啦啦啦视频在线资源免费观看| 我的亚洲天堂| 亚洲国产精品一区三区| 国产激情久久老熟女| 又粗又硬又长又爽又黄的视频| 操美女的视频在线观看| 亚洲欧美成人综合另类久久久| 99香蕉大伊视频| 99精品久久久久人妻精品| 国产精品一国产av| 只有这里有精品99| 日韩av在线免费看完整版不卡| 男女下面插进去视频免费观看| 在线看a的网站| 久久99一区二区三区| 最近最新中文字幕大全免费视频 | 人妻 亚洲 视频| 两个人免费观看高清视频| 久热这里只有精品99| 欧美黑人精品巨大| 午夜免费观看性视频| 久久精品aⅴ一区二区三区四区| 久久婷婷青草| 赤兔流量卡办理| 国产成人午夜福利电影在线观看| 午夜免费男女啪啪视频观看| 午夜福利乱码中文字幕| 久久久精品免费免费高清| 成人午夜精彩视频在线观看| avwww免费| 国产精品久久久人人做人人爽| av视频免费观看在线观看| 精品一品国产午夜福利视频| 秋霞伦理黄片| 亚洲精品一二三| 深夜精品福利| 美女高潮到喷水免费观看| 亚洲av电影在线进入| 老司机影院毛片| 精品国产露脸久久av麻豆| 中文欧美无线码| 少妇被粗大猛烈的视频| 精品久久久久久电影网| 亚洲免费av在线视频| 色综合欧美亚洲国产小说| 久久久国产一区二区| 久久精品亚洲熟妇少妇任你| 80岁老熟妇乱子伦牲交| 欧美在线一区亚洲| 日韩中文字幕视频在线看片| 中文字幕高清在线视频| 亚洲av综合色区一区| 免费在线观看完整版高清| 大香蕉久久成人网| av又黄又爽大尺度在线免费看| 激情五月婷婷亚洲| 校园人妻丝袜中文字幕| 精品国产国语对白av| av女优亚洲男人天堂| 国产亚洲一区二区精品| 999久久久国产精品视频| 亚洲国产毛片av蜜桃av| 日韩欧美一区视频在线观看| 亚洲国产精品999| 成人亚洲精品一区在线观看| 精品免费久久久久久久清纯 | 国产精品嫩草影院av在线观看| 少妇精品久久久久久久| 精品国产乱码久久久久久男人| 国产成人免费观看mmmm| 99九九在线精品视频| 老熟女久久久| 夫妻性生交免费视频一级片| 亚洲国产av新网站| 午夜91福利影院| 1024香蕉在线观看| 国产精品久久久久久人妻精品电影 | 国产精品熟女久久久久浪| 久久久精品区二区三区| 人人妻,人人澡人人爽秒播 | 观看美女的网站| 久久久久国产一级毛片高清牌| 亚洲av福利一区| 99香蕉大伊视频| 自拍欧美九色日韩亚洲蝌蚪91| 国产片内射在线| 好男人视频免费观看在线| 最近2019中文字幕mv第一页| 天天躁夜夜躁狠狠躁躁| 国产精品香港三级国产av潘金莲 | 一个人免费看片子| 久久av网站| 一级毛片黄色毛片免费观看视频| 日韩大码丰满熟妇| e午夜精品久久久久久久| 精品一区在线观看国产| 卡戴珊不雅视频在线播放| a级毛片黄视频| 亚洲av日韩精品久久久久久密 | 婷婷色av中文字幕| 搡老岳熟女国产| 国产精品二区激情视频| 亚洲七黄色美女视频| 99香蕉大伊视频| www.精华液| 好男人视频免费观看在线| 一边摸一边做爽爽视频免费| 精品亚洲成a人片在线观看| 欧美激情 高清一区二区三区| 丝袜美腿诱惑在线| 两性夫妻黄色片| 成年女人毛片免费观看观看9 | 在线观看国产h片| 美国免费a级毛片| 999精品在线视频| 日韩欧美一区视频在线观看| 亚洲人成77777在线视频| 只有这里有精品99| 久久久久国产一级毛片高清牌| 国产99久久九九免费精品| 黄色视频在线播放观看不卡| 国产乱人偷精品视频| 尾随美女入室| 高清黄色对白视频在线免费看| 日韩精品有码人妻一区| 亚洲国产精品一区三区| 热99国产精品久久久久久7| 国产一级毛片在线| av视频免费观看在线观看| 一级毛片黄色毛片免费观看视频| 99香蕉大伊视频| 在线观看国产h片| 五月天丁香电影| 精品亚洲乱码少妇综合久久| 99re6热这里在线精品视频| 一级片免费观看大全| 99久久精品国产亚洲精品| 校园人妻丝袜中文字幕| 免费黄色在线免费观看| 欧美在线黄色| 久久国产精品大桥未久av| 亚洲天堂av无毛| 在线观看国产h片| 日韩 亚洲 欧美在线| 国产xxxxx性猛交| xxxhd国产人妻xxx| 最近的中文字幕免费完整| 搡老岳熟女国产| 欧美日韩综合久久久久久| 如日韩欧美国产精品一区二区三区| 美女扒开内裤让男人捅视频| 日韩免费高清中文字幕av| 韩国精品一区二区三区| 看免费成人av毛片| 国产一区二区三区av在线| 国产片内射在线| 久久国产精品男人的天堂亚洲| 男女高潮啪啪啪动态图| 亚洲婷婷狠狠爱综合网| 国产精品久久久久久久久免| 免费在线观看完整版高清| 视频在线观看一区二区三区| 国产男女内射视频| 色婷婷久久久亚洲欧美| 天天躁狠狠躁夜夜躁狠狠躁| 欧美最新免费一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| 精品人妻一区二区三区麻豆| 亚洲美女视频黄频| 精品一区二区三区av网在线观看 | 欧美xxⅹ黑人| 久久这里只有精品19| 国精品久久久久久国模美| 18禁动态无遮挡网站| 欧美激情极品国产一区二区三区| 街头女战士在线观看网站| 久久综合国产亚洲精品| 国产日韩欧美视频二区| 国产不卡av网站在线观看| 亚洲av国产av综合av卡| 欧美国产精品一级二级三级| 丝袜脚勾引网站| 老汉色av国产亚洲站长工具| 欧美人与性动交α欧美精品济南到| 人人妻人人添人人爽欧美一区卜| 黄网站色视频无遮挡免费观看| 午夜福利,免费看| 国产 一区精品| 亚洲,一卡二卡三卡| 久久人人97超碰香蕉20202| 女性生殖器流出的白浆| 成人国产麻豆网| 嫩草影院入口| av在线老鸭窝| 纵有疾风起免费观看全集完整版| 亚洲精品国产一区二区精华液| 亚洲欧美中文字幕日韩二区| 秋霞伦理黄片| 久久99精品国语久久久| 国产精品欧美亚洲77777| 国产精品免费视频内射| 亚洲,欧美精品.| 免费观看a级毛片全部| 日日摸夜夜添夜夜爱| 久久精品国产a三级三级三级| 国产高清不卡午夜福利| 成人国产av品久久久| 日韩大码丰满熟妇| 人人妻人人添人人爽欧美一区卜| 久久久久久人人人人人| 天天添夜夜摸| 又粗又硬又长又爽又黄的视频| 一二三四在线观看免费中文在| 亚洲男人天堂网一区| 69精品国产乱码久久久| 99九九在线精品视频| 美女中出高潮动态图| 超色免费av| 丁香六月欧美| 国产一区有黄有色的免费视频| 一级片免费观看大全| 天美传媒精品一区二区| 国产在线一区二区三区精| 精品免费久久久久久久清纯 | 国产精品免费视频内射| 亚洲精品美女久久久久99蜜臀 | 国产成人免费观看mmmm| 免费黄频网站在线观看国产| 狂野欧美激情性xxxx| 伊人久久国产一区二区| 亚洲av欧美aⅴ国产| 亚洲av电影在线观看一区二区三区| 18禁动态无遮挡网站| 亚洲第一区二区三区不卡| 国产高清国产精品国产三级| 三上悠亚av全集在线观看| 亚洲第一区二区三区不卡| 国产精品一区二区在线不卡| 成年av动漫网址| 国产成人精品久久二区二区91 | netflix在线观看网站| 日韩一区二区视频免费看| 免费观看a级毛片全部| 欧美 亚洲 国产 日韩一| 黄片小视频在线播放| 日韩制服骚丝袜av| 丝袜脚勾引网站| 欧美激情极品国产一区二区三区| 国产一区亚洲一区在线观看| 97精品久久久久久久久久精品| 看免费成人av毛片| 色综合欧美亚洲国产小说| 丝袜在线中文字幕| 在线看a的网站| av在线播放精品| 美女中出高潮动态图| 婷婷色综合大香蕉| 欧美国产精品va在线观看不卡| 欧美激情 高清一区二区三区| 99热网站在线观看| 建设人人有责人人尽责人人享有的| 亚洲成人av在线免费| 搡老乐熟女国产| 国产亚洲av片在线观看秒播厂| 婷婷色麻豆天堂久久| 日韩欧美精品免费久久| 国产伦理片在线播放av一区| 日韩大片免费观看网站| 国产日韩欧美视频二区| 欧美日韩一区二区视频在线观看视频在线| 国产高清国产精品国产三级| 国产亚洲av片在线观看秒播厂| 成人手机av| 国产片特级美女逼逼视频| 国产一区二区激情短视频 | 亚洲av综合色区一区| 亚洲一码二码三码区别大吗| 久久精品久久精品一区二区三区| 久久狼人影院| 国产一区二区在线观看av| 黄网站色视频无遮挡免费观看| 国产伦人伦偷精品视频| 黄色毛片三级朝国网站| 国产一卡二卡三卡精品 | 亚洲av电影在线观看一区二区三区| 天天躁狠狠躁夜夜躁狠狠躁| 久久久久国产一级毛片高清牌| 成人午夜精彩视频在线观看| 精品一区在线观看国产| 考比视频在线观看| 久久久精品94久久精品| av视频免费观看在线观看| 国产精品 欧美亚洲| 一区二区日韩欧美中文字幕| 熟女少妇亚洲综合色aaa.| 日韩电影二区| 男女边吃奶边做爰视频| 午夜福利一区二区在线看| 色吧在线观看| av线在线观看网站| 一区在线观看完整版| 亚洲熟女毛片儿| 日韩av不卡免费在线播放| 日本欧美国产在线视频| 色网站视频免费| av卡一久久| 19禁男女啪啪无遮挡网站| 久久热在线av| 久热这里只有精品99| 成人三级做爰电影| 一个人免费看片子| 日本vs欧美在线观看视频| 男人舔女人的私密视频| 亚洲天堂av无毛| 啦啦啦啦在线视频资源| 丝袜人妻中文字幕| 欧美精品亚洲一区二区| 亚洲一码二码三码区别大吗|