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

    差分粒子群算法在PUMA機(jī)器人逆運(yùn)動(dòng)學(xué)求解中的應(yīng)用

    2020-06-19 08:45:58汪院林袁銳波袁安華
    軟件導(dǎo)刊 2020年4期

    汪院林 袁銳波 袁安華

    摘要:傳統(tǒng)逆運(yùn)動(dòng)學(xué)求解主要從逆運(yùn)動(dòng)學(xué)方程出發(fā),基于一定的數(shù)學(xué)理論推導(dǎo),不能完全實(shí)現(xiàn)計(jì)算機(jī)程序化,且精度與計(jì)算效率較低:為改善這一缺陷,基于機(jī)器人正向運(yùn)動(dòng)學(xué)方程,借助MATLAB工具,使用蒙特卡洛法仿真分析出PUMA560機(jī)器人的工作空間,任取一點(diǎn)末端執(zhí)行器位姿作為逆運(yùn)動(dòng)學(xué)求解的已知位姿矩陣T,結(jié)合差分粒子群仿生智能算法作為逆運(yùn)動(dòng)學(xué)求解的主要理論算法。將計(jì)算出的旋轉(zhuǎn)關(guān)節(jié)變量θ1-θ6代入正運(yùn)動(dòng)學(xué)方程,得出末端位姿矩陣于;通過計(jì)算分析T與于相關(guān)角度誤差,兩矩陣所對應(yīng)的位置向量與姿態(tài)向量誤差精度為0.001數(shù)量級,完全滿足目前機(jī)器人定位要求?;诓罘至W尤豪碚摰臋C(jī)器人逆運(yùn)動(dòng)學(xué)求解方法計(jì)算收斂速度更快,能高度實(shí)現(xiàn)計(jì)算機(jī)程序化,誤差精度高,提高計(jì)算效率。

    關(guān)鍵詞:正運(yùn)動(dòng)學(xué);差分粒子群算法;收斂速度;位姿矩陣;逆運(yùn)動(dòng)學(xué)解

    DOI: 10. 11907/rjdk.191724

    開放科學(xué)(資源服務(wù))標(biāo)識碼(OSID):

    中圖分類號:TP319

    文獻(xiàn)標(biāo)識碼:A

    文章編號:1672-7800(2020)004-0203-05

    0 引言

    近年來,機(jī)器人技術(shù)發(fā)展迅速,在工業(yè)等領(lǐng)域占據(jù)重要地位,機(jī)器人運(yùn)動(dòng)學(xué)研究也日益增多。機(jī)器人運(yùn)動(dòng)學(xué)主要由正運(yùn)動(dòng)學(xué)與逆運(yùn)動(dòng)學(xué)兩部分組成:機(jī)器人正運(yùn)動(dòng)學(xué)主要以基坐標(biāo)為基礎(chǔ),根據(jù)各關(guān)節(jié)旋轉(zhuǎn)角度確定末端執(zhí)行器的位姿;逆運(yùn)動(dòng)學(xué)與其相反,以末端關(guān)節(jié)為出發(fā)點(diǎn)最終回到基坐標(biāo)。正運(yùn)動(dòng)學(xué)由多個(gè)關(guān)節(jié)的位姿矩陣進(jìn)行連乘運(yùn)算,最終確定末端執(zhí)行器唯一的位姿矩陣,而逆運(yùn)動(dòng)學(xué)末端位姿可以對應(yīng)機(jī)器人系統(tǒng)關(guān)節(jié)的多個(gè)解[1]。機(jī)器人逆運(yùn)動(dòng)學(xué)問題隨著運(yùn)動(dòng)關(guān)節(jié)的增多而更加復(fù)雜,它實(shí)質(zhì)上是解一個(gè)非線性方程問題。實(shí)際應(yīng)用中,機(jī)器人所要執(zhí)行的任務(wù)通常用其基坐標(biāo)系描述[2]。因此,機(jī)器人逆運(yùn)動(dòng)學(xué)求解在控制領(lǐng)域顯得尤為重要。

    目前,機(jī)器人逆運(yùn)動(dòng)方程求解方法主要有解析法、幾何法、迭代法、幾何一解析法[3]等。Paul等[4]對PUMA機(jī)器人提出了當(dāng)d6=0時(shí)逆運(yùn)動(dòng)的解析解;Lee等[5]提出應(yīng)用幾何學(xué)方法求解各關(guān)節(jié)變量;王奇志等[6]提出當(dāng)d6≠0時(shí)的一種解析解法;戴齊等[2]對PUMA機(jī)器人逆運(yùn)動(dòng)解利用幾何一解析法進(jìn)行研究?;趥鹘y(tǒng)理論方法雖然模型建立簡單,但計(jì)算效率低下,難于析出機(jī)器人逆運(yùn)動(dòng)學(xué)的所有解,且計(jì)算機(jī)程序的可移植性較差;吳振宇等[7]利用遺傳粒子群算法在解決非線性問題時(shí)的優(yōu)勢,對PUMA機(jī)器人的逆運(yùn)動(dòng)解提出遺傳粒子群算法理論,但該方法需要大量的樣本數(shù)據(jù)進(jìn)行訓(xùn)練,且對目標(biāo)函數(shù)進(jìn)行最優(yōu)搜索時(shí)收斂速度較慢,會(huì)陷入局部最優(yōu)狀態(tài),存在一定不足。

    本文提出一種新的逆運(yùn)動(dòng)求解方法:基于差分進(jìn)化算法( Differential Evolution Algorithm.DE),該方法結(jié)構(gòu)簡單,收斂速度快,魯棒性較強(qiáng),是與粒子群優(yōu)化算法( ParticleSwarm Optimization,PSO)同時(shí)考慮位置與速度兩個(gè)參數(shù)相結(jié)合的一種混合智能算法[8]。從正向運(yùn)動(dòng)學(xué)方程出發(fā),在各關(guān)節(jié)角變量范圍內(nèi)取出相應(yīng)角度值,建立PUMA機(jī)器人工作空間范圍模型。通過圖像方式驗(yàn)證三維空間尺寸是否與機(jī)器人結(jié)構(gòu)尺寸一致,圖形各投影平面中有沒有明顯的空洞現(xiàn)象?;谒心┒藞?zhí)行器位姿點(diǎn),取其中某一點(diǎn)位姿,將差分粒子群混合算法作為PUMA機(jī)器人逆運(yùn)動(dòng)求解主要理論,彌補(bǔ)了傳統(tǒng)逆運(yùn)動(dòng)求解方法的不足,提高了計(jì)算效率和精度。

    1 正運(yùn)動(dòng)學(xué)模型

    1.1 D-H法坐標(biāo)系及幾何參數(shù)

    PUMA型機(jī)器人有6個(gè)自由度,且6個(gè)關(guān)節(jié)均為旋轉(zhuǎn)類型,可采用標(biāo)準(zhǔn)D-H法建立桿件坐標(biāo)系描述各關(guān)節(jié)之間的特征參數(shù)和運(yùn)動(dòng)關(guān)系。該矩陣為4x4的齊次變換矩陣,通過該矩陣可得出與PUMA560機(jī)器人幾何尺寸相關(guān)的運(yùn)動(dòng)學(xué)方程。標(biāo)準(zhǔn)D-H法得到的相鄰桿件之間的齊次變換矩陣如下:

    各連桿幾何參數(shù)如表1所示。

    1.2 運(yùn)動(dòng)學(xué)方程

    PUMA560機(jī)器人由多個(gè)旋轉(zhuǎn)關(guān)節(jié)組成,其中包含旋轉(zhuǎn)基座、大小臂、手腕等幾部分[9]。前3個(gè)關(guān)節(jié)用于對末端執(zhí)行器進(jìn)行定位,后3個(gè)關(guān)節(jié)主要用于保證末端執(zhí)行器的姿態(tài)。PUMA機(jī)器人運(yùn)動(dòng)學(xué)正解方程如下:

    若給出關(guān)節(jié)旋轉(zhuǎn)角及各連桿幾何參數(shù),通過式(2)可求出PUMA機(jī)器人末端執(zhí)行器位姿;若僅已知末端關(guān)節(jié)位姿,則需要確定中間各關(guān)節(jié)旋轉(zhuǎn)角,這需要進(jìn)行大量計(jì)算方可求出。

    2 PUMA機(jī)器人工作空間仿真

    2.1 蒙特卡洛法工作空間分析

    目前,確定機(jī)器人工作空間的方法主要有數(shù)值法、解析法及圖解法。圖解法方法較簡單,但受關(guān)節(jié)數(shù)目限制,對多自由度機(jī)器人進(jìn)行工作空間分析時(shí)需采用分組處理[10];解析法主要采用微分幾何學(xué)原理,通過繪制函數(shù)解析式的多次包絡(luò)曲線確定工作空間,表示復(fù)雜,多用于三自由度以下機(jī)器人的工作空間分析[11];數(shù)值法多建立在極值理論及相關(guān)優(yōu)化理論之上,通常計(jì)算量龐大,用于圖形邊界為曲面且凹陷時(shí),計(jì)算可靠性低[12]。本文提出一種建立于數(shù)值法基礎(chǔ)上且引進(jìn)隨機(jī)分布方法的蒙特卡洛法,它主要依托隨機(jī)抽樣的數(shù)學(xué)理論將工程問題轉(zhuǎn)化為數(shù)學(xué)模型,應(yīng)用范圍較廣。在工程中,主要用于一些隨機(jī)無規(guī)律性的物理問題,而在機(jī)器人領(lǐng)域多用來進(jìn)行誤差分析。然而,機(jī)器人工作空間的仿真計(jì)算也能通過蒙特卡洛法理論完成,可利用計(jì)算機(jī)輔助計(jì)算[13]。

    蒙特卡洛法求解機(jī)器人工作空間,主要是在機(jī)器人各關(guān)節(jié)的旋轉(zhuǎn)變量范圍內(nèi)按一定步長隨機(jī)取值,然后把各旋轉(zhuǎn)關(guān)節(jié)的取值代入方程式(2)中,最終計(jì)算出自由端執(zhí)行器的空間坐標(biāo)值,而末端關(guān)節(jié)所有隨機(jī)值的集合則構(gòu)成機(jī)器人工作空間。將末端執(zhí)行器得到的坐標(biāo)點(diǎn)輸出到圖形設(shè)備上,可得出機(jī)器人工作空間的三維仿真圖形。

    蒙特卡洛法主要步驟如下[14]:

    (5)設(shè)定所需要輸出的點(diǎn)集數(shù),通過循環(huán)方式重復(fù)步驟(3)、步驟(4),直到滿足條件為止。

    上述算法可借助MATLAB軟件程序化完成,最終將末端位置空間點(diǎn)通過三維繪圖函數(shù)plot3輸出。利用蒙特卡洛法進(jìn)行機(jī)器人工作空間仿真步驟如圖2所示。

    2.2 PUMA560機(jī)器人工作空間仿真分析

    根據(jù)機(jī)器人運(yùn)動(dòng)學(xué)正解方程及蒙特卡洛法理論,將該算法編寫成MATLAB程序進(jìn)行仿真分析。PUMA560機(jī)器人各關(guān)節(jié)角度取值點(diǎn)為50 000個(gè),最終輸出50 000個(gè)點(diǎn)作為三維工作空間,仿真結(jié)果如圖3-圖6所示。

    從以上三維工作空間仿真圖及各平面投影圖分析可得出:仿真分析出的工作空間均未超出各關(guān)節(jié)變量旋轉(zhuǎn)范圍,且各圖形尺寸范圍與PUMA560機(jī)器人結(jié)構(gòu)尺寸相一致。從圖中看出各點(diǎn)集均勻,未出現(xiàn)空洞等現(xiàn)象。因此,該仿真結(jié)果與實(shí)際相符。

    3 逆運(yùn)動(dòng)學(xué)求解

    粒子群優(yōu)化算法(PSO)是來源于觀察鳥群覓食而誕生的一種仿生算法,廣泛應(yīng)用于非線性問題、多極值和不可微等場合的優(yōu)化,需要設(shè)置的參數(shù)較少[16]。PSO算法在搜索空間范圍內(nèi)尋找最優(yōu)解,主要與粒子適應(yīng)度和群體兩個(gè)參數(shù)有關(guān),而粒子群的每一個(gè)體可表示某優(yōu)化問題的一個(gè)概率解。每個(gè)粒子主要由速度和位置兩大屬性決定。粒子適應(yīng)度是將粒子位置代入目標(biāo)函數(shù)得到的數(shù)值,粒子適應(yīng)度是判斷粒子優(yōu)劣的標(biāo)準(zhǔn)。在對問題進(jìn)行搜索迭代時(shí),粒子主要依據(jù)兩個(gè)參數(shù)作出變換:①粒子受自身經(jīng)驗(yàn)影響而得到的最優(yōu)解,即個(gè)體歷史極值pBest;②粒子受全局影響而尋到的最優(yōu)解,即全局歷史極值gBest。用以下公式對粒子位置及速度進(jìn)行更新。

    式中,v(t+1)、x(t+1)分別為下一時(shí)刻粒子的速度與位置;C1、c2為學(xué)習(xí)因子,代表個(gè)體經(jīng)驗(yàn)與群體經(jīng)驗(yàn)對粒子速度影響所占權(quán)重;rand()為隨機(jī)數(shù),范圍為[0,1];ω為慣性權(quán)重,是與時(shí)間有關(guān)的動(dòng)態(tài)調(diào)整參數(shù)。ω取值對算法影響至關(guān)重要,ω較大能提高全局搜索能力,但收斂速度很慢,較小則易陷入局部搜索。本文對ω的取值隨迭代次數(shù)從0.9線性遞減為0.4[17],變化公式如下:

    差分進(jìn)化算法(DE)是一種并行求解的進(jìn)化算法,與遺傳算法類似,用相同的算術(shù)操作(如:變異、交叉和選擇算子)對初始群體或父代群體進(jìn)行更新。DE首先進(jìn)行變異操作,在搜索范圍內(nèi)任意選取一個(gè)體作為基準(zhǔn)父代個(gè)體,并在父代個(gè)體外隨機(jī)選取兩個(gè)向量進(jìn)行差值計(jì)算,加權(quán)后與父代個(gè)體求和;然后經(jīng)交叉算子在父代個(gè)體與變異個(gè)體間進(jìn)行選擇,產(chǎn)生實(shí)驗(yàn)個(gè)體;最終結(jié)合目標(biāo)函數(shù)對父代個(gè)體和實(shí)驗(yàn)個(gè)體作選擇,生成新的子代個(gè)體[18]。

    其中,F(xiàn)(x)為對應(yīng)問題的目標(biāo)函數(shù)值。

    本文根據(jù)PUMA560機(jī)器人逆運(yùn)動(dòng)學(xué)方程建立數(shù)學(xué)模型(目標(biāo)函數(shù)):機(jī)器人末端執(zhí)行器位姿矩陣為T,利用差分粒子群算法優(yōu)化得出第t次搜索時(shí)各關(guān)節(jié)變量所對應(yīng)的矩陣T,數(shù)學(xué)模型可表達(dá)為兩矩陣T與T之間的空間距離,使其為極小佰,計(jì)算如下:

    本文結(jié)合DE與PSO各自的優(yōu)點(diǎn),提出一種差分粒子群算法(DEPSO),用于解決PUMA560機(jī)器人運(yùn)動(dòng)學(xué)問題。該算法以PSO流程為主,通過DE算法增強(qiáng)了粒子群體的多樣性,可避免PSO算法陷入局部最優(yōu)解。DEPSO算法步驟如下:①在粒子群體范圍內(nèi),初始化粒子位置與速度,設(shè)計(jì)數(shù)器t=0;②計(jì)算粒子適應(yīng)度,用于DE算法選擇粒子,計(jì)數(shù)器t=t+1;③粒子根據(jù)式(8)進(jìn)行DE變異操作并作選擇;④更新局部最優(yōu)解pBest與全局最優(yōu)解,根據(jù)式(5).式(6)計(jì)算粒子位置與速度;⑤粒子根據(jù)式(9)進(jìn)行DE交叉計(jì)算;⑥計(jì)算粒子適應(yīng)度;⑦選擇新一代粒子;⑧更新局部最優(yōu)解和全局最優(yōu)解;⑨若計(jì)數(shù)器大于給定循環(huán)次數(shù),則退出計(jì)算,否則,調(diào)至步驟②。

    DEPSO算法基本流程如圖7所示。

    4 實(shí)驗(yàn)結(jié)果與分析

    PUMA560機(jī)器人各桿件參數(shù)如表1所示。根據(jù)蒙特卡洛法機(jī)器人正運(yùn)動(dòng)學(xué)求出所有末端關(guān)節(jié)位姿矩陣,任取其中一點(diǎn)作為逆運(yùn)動(dòng)學(xué)求解的末端位姿矩陣T,如下所示:

    根據(jù)矩陣T與于,可從位置與姿態(tài)矩陣分析求解精度。實(shí)驗(yàn)?zāi)┒岁P(guān)節(jié)矩陣T對應(yīng)的位置坐標(biāo)為[ 600.227,622.914,164.221];關(guān)節(jié)變量θ1-θ6所對應(yīng)矩陣T的位置坐標(biāo)為[600.2283,622.9124,164.2217],計(jì)算得出T與于之間的誤差精度為0.001mm數(shù)量級。

    姿態(tài)矩陣所對應(yīng)角度為T和于各自的反余弦值[20]。對T與于進(jìn)行分析可得出:PUMA560機(jī)器人指定末端位姿矩陣法向向量各自的方向角[106.5334°,76.3895°,21.6705°],所求關(guān)節(jié)角對應(yīng)位姿矩陣的法向向量方向角為[106.5313°,76.3855°,21.6714°];機(jī)器人指定末端位姿矩陣滑動(dòng)向量方向角為[61.0559°,143.2778°,69.4417°],所求關(guān)節(jié)角對應(yīng)的位姿矩陣滑動(dòng)向量方向角為[61.0585°,143.2788°,69.4397°];PUMA560機(jī)器人指定末端位姿矩陣所對應(yīng)的接近向量方向角各自為[34.1543°,56.6547°,83.4417°],所求關(guān)節(jié)角對應(yīng)位姿矩陣的接近向量方向角各自為[34.1507°,56.6577°,83.4445°】。通過以上計(jì)算可看出,利用差分粒子群算法計(jì)算PUMA560機(jī)器人得出的關(guān)節(jié)變量所對應(yīng)的末端姿態(tài)矩陣與指定末端姿態(tài)矩陣的誤差精度在0.001°數(shù)量級,由此可知,將本文的差分粒子群算法用于求解PUMA560機(jī)器人逆運(yùn)動(dòng)學(xué)求解完全滿足定位精度要求。

    5 結(jié)語

    本文從正運(yùn)動(dòng)學(xué)方程出發(fā),根據(jù)蒙特卡洛法求解出PUMA560機(jī)器人工作空間的末端關(guān)節(jié)位姿矩陣,任取一點(diǎn)用于逆運(yùn)動(dòng)學(xué)求解過程。結(jié)合差分優(yōu)化算法與粒子群算法優(yōu)點(diǎn),提出一種差分粒子群混合智能算法,將其作為逆運(yùn)動(dòng)學(xué)求解方法。該方法與其它傳統(tǒng)逆運(yùn)動(dòng)學(xué)求解方法相比較,更適合于計(jì)算機(jī)計(jì)算,具有定位精度高、收斂速度快、計(jì)算高效等優(yōu)點(diǎn)。通過實(shí)例驗(yàn)證可知,該方法用于機(jī)器人逆運(yùn)動(dòng)求解問題滿足當(dāng)前機(jī)器人定位誤差要求,是智能算法與逆運(yùn)動(dòng)數(shù)值法相結(jié)合的新理論。

    參考文獻(xiàn):

    [1]JOHN J,CRAIG.機(jī)器人學(xué)導(dǎo)論[M].以,贠超,譯.北京:機(jī)械工業(yè)出版社.2006.

    [2] 戴齊,姚先啟.一種求解PUMA機(jī)械手運(yùn)動(dòng)學(xué)逆問題的分解法[J].西南交通大學(xué)學(xué)報(bào),1989(4):30-35.

    [3] 陳寧,焦恩璋.PUMA機(jī)械手逆運(yùn)動(dòng)方程求解新方法[J].南京林業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版),2003,27(4):23-26.

    [4]PAUL R P.Robot manipulators: mathematics, programming and con-trol: the computer control of robot manipulators[M].Cambridge: MITPress.1981.

    [5] FU K S,CONZALEZ R C,LEE C S G.Robotics: control, sensing,vision and intelligence [J]. Me Graw-Hill. 1987(6):78-82.

    [6] 王奇志,徐心和,尹朝萬.PUMA機(jī)械手逆運(yùn)動(dòng)方程新的推導(dǎo)方法及求解[J].機(jī)器人,1998,20(2):81-87.

    [7]吳振宇,姚明江,馮林,等.基于遺傳粒子群的PUMA機(jī)器人逆運(yùn)動(dòng)求解[J].計(jì)算機(jī)仿真,2011,28( 10):173-176.

    [8]池元成,方杰,蔡國飆.基于差分進(jìn)化和粒子群優(yōu)化算法的混合優(yōu)化算法[J].計(jì)算機(jī)工程與設(shè)計(jì),2009.30( 12):641-647.

    [9] 布魯諾·西西里安諾.機(jī)器人學(xué):建模、規(guī)劃與控制[M].張國良,譯.西安:西安交通大學(xué)出版社,2013.

    [10]蔡蒂,謝存禧,張鐵,等.基于蒙特卡洛法的噴涂機(jī)器人工作空間分析及仿真[J].機(jī)械設(shè)計(jì)與制造,2009(3):161-162.

    [11] KOHLI D. HSU M S.The jacobian analysis of workspaces of mechan-ical manipulators[J]. Mechanism& Machine Theory, 1987, 22(3):265-275.

    [12]趙燕江,張永德,姜金剛.基于Matlah的機(jī)器人工作空間求解方法[J].機(jī)械科學(xué)與技術(shù),2009,28(12):1657-1661.

    [13] RASTEGAR J,PEREL D.Ceneration of manipulator workspaceboundary geometry using the monte carlo method and interactive com-puter graphics[J]. Journal of mechanical design, 1990, 112 (3):452-454.

    [14] 郭明,周國斌.多關(guān)節(jié)機(jī)器人工作空間的分析與評價(jià)方法[J].機(jī)器人,1988,10 (4):7-12.

    [15] 沙漠,鄧子龍.基于MATLAB對于PUMA560機(jī)器人的運(yùn)動(dòng)空間分析研究[J].機(jī)械制造與自動(dòng)化,2016,45(2):156-159.

    [16] KAO Y T,ZAHARA E.A hybrid genetic algorithm and particleswarm optimization for multimodal functions[J]. Applied Soft Com-puting Journal, 2008.8 (2): 849-857.

    [17] 王志,胡小兵,何雪海.一種新的差分與粒子群算法的混合算法[J].計(jì)算機(jī)工程與應(yīng)用,2012,48(6):46-48.

    [18] 欒麗君,譚立靜,牛奔.一種基于粒子群優(yōu)化算法和差分進(jìn)化算法的新型混合全局優(yōu)化算法[J].信息與控制,2007,36(6): 708-714.

    [19]謝宏,禹文科,楊鵬,等,多子群分層差分粒子群算法的逆運(yùn)動(dòng)求解方法[J].電子測量與儀器學(xué)報(bào),2015,29(10):95-99.

    [20] 芮梃.PUMA機(jī)器人逆運(yùn)動(dòng)模擬退火粒子群求解方法[J].計(jì)算機(jī)工程與應(yīng)用,2010,46(3):27-29.

    (責(zé)任編輯:杜能鋼)

    作者簡介:汪院林(1994-),男,昆明理工大學(xué)機(jī)電工程學(xué)院碩士研究生,研究方向?yàn)榱黧w傳動(dòng)與控制;袁銳波(1968-),男,博士,昆明理工大學(xué)機(jī)電工程學(xué)院教授,研究方向?yàn)殡娨罕壤欧刂?、?fù)雜機(jī)電系統(tǒng)集成與控制;袁安華(1991-),男,昆明理工大學(xué)機(jī)電工程學(xué)院碩士研究生,研究方向?yàn)榱黧w傳動(dòng)與控制。本文通訊作者:袁銳波。

    亚洲精品美女久久久久99蜜臀| 国产免费av片在线观看野外av| 色婷婷久久久亚洲欧美| e午夜精品久久久久久久| 琪琪午夜伦伦电影理论片6080| 亚洲电影在线观看av| 丁香六月欧美| 中文字幕精品亚洲无线码一区 | www.www免费av| 国产91精品成人一区二区三区| 久热爱精品视频在线9| 久久国产精品人妻蜜桃| 麻豆成人午夜福利视频| 国产伦人伦偷精品视频| 18禁美女被吸乳视频| 天天添夜夜摸| 日本免费a在线| 一边摸一边抽搐一进一小说| 欧美日韩一级在线毛片| 一级作爱视频免费观看| 成人三级做爰电影| 啦啦啦韩国在线观看视频| 国产激情偷乱视频一区二区| 少妇 在线观看| 日韩免费av在线播放| 日韩中文字幕欧美一区二区| 在线观看66精品国产| 国内少妇人妻偷人精品xxx网站 | 在线观看午夜福利视频| tocl精华| 亚洲一区高清亚洲精品| 欧美久久黑人一区二区| 热re99久久国产66热| 级片在线观看| 成人一区二区视频在线观看| 欧美乱色亚洲激情| 麻豆成人av在线观看| 久久欧美精品欧美久久欧美| av有码第一页| 欧美人与性动交α欧美精品济南到| 久久热在线av| 日韩欧美 国产精品| 久久精品国产亚洲av香蕉五月| 日本精品一区二区三区蜜桃| 哪里可以看免费的av片| 国产三级在线视频| 色在线成人网| 在线免费观看的www视频| 欧美激情极品国产一区二区三区| 久久久久久大精品| 91av网站免费观看| 国产精品久久久久久人妻精品电影| 久9热在线精品视频| 亚洲成国产人片在线观看| 成人亚洲精品av一区二区| 中文在线观看免费www的网站 | 亚洲中文日韩欧美视频| 欧美日韩精品网址| 久久精品人妻少妇| 男女下面进入的视频免费午夜 | 热99re8久久精品国产| 国内精品久久久久久久电影| 老司机福利观看| 亚洲在线自拍视频| 757午夜福利合集在线观看| 身体一侧抽搐| 国产在线观看jvid| 亚洲av熟女| 亚洲欧美激情综合另类| 熟妇人妻久久中文字幕3abv| 国内揄拍国产精品人妻在线 | 在线国产一区二区在线| 97人妻精品一区二区三区麻豆 | 中文字幕另类日韩欧美亚洲嫩草| 亚洲国产精品成人综合色| 后天国语完整版免费观看| 国产v大片淫在线免费观看| 国产在线精品亚洲第一网站| 欧美日本视频| 美女大奶头视频| 搡老妇女老女人老熟妇| 国产v大片淫在线免费观看| 岛国在线观看网站| 日韩欧美三级三区| 99热6这里只有精品| 韩国av一区二区三区四区| 人人妻人人澡欧美一区二区| 国产精品久久久久久人妻精品电影| 亚洲国产中文字幕在线视频| 亚洲精品一卡2卡三卡4卡5卡| 日韩欧美一区视频在线观看| 后天国语完整版免费观看| 窝窝影院91人妻| 天天躁夜夜躁狠狠躁躁| 久久中文字幕人妻熟女| 国产精品 欧美亚洲| 国产视频一区二区在线看| 亚洲男人天堂网一区| 久久午夜亚洲精品久久| 国产亚洲精品av在线| 波多野结衣av一区二区av| 欧美zozozo另类| 午夜成年电影在线免费观看| 少妇粗大呻吟视频| 性色av乱码一区二区三区2| 国产成人欧美在线观看| 黄色成人免费大全| 19禁男女啪啪无遮挡网站| 精品一区二区三区四区五区乱码| 欧美精品啪啪一区二区三区| 岛国视频午夜一区免费看| 久久性视频一级片| 亚洲国产精品久久男人天堂| 精品卡一卡二卡四卡免费| 在线观看66精品国产| 久久人人精品亚洲av| 国产私拍福利视频在线观看| 久久久精品欧美日韩精品| 日韩大码丰满熟妇| 夜夜夜夜夜久久久久| 久久久久久久久免费视频了| 亚洲一区中文字幕在线| 亚洲精品在线观看二区| 天堂√8在线中文| 麻豆成人午夜福利视频| 在线十欧美十亚洲十日本专区| 国产主播在线观看一区二区| 欧美日韩福利视频一区二区| 天天一区二区日本电影三级| 亚洲国产欧洲综合997久久, | 欧美av亚洲av综合av国产av| 在线观看免费日韩欧美大片| 国产亚洲精品综合一区在线观看 | 亚洲欧美一区二区三区黑人| 久久国产精品人妻蜜桃| 免费看a级黄色片| 久久精品国产亚洲av高清一级| 日韩欧美在线二视频| 久久精品91蜜桃| 女生性感内裤真人,穿戴方法视频| 两人在一起打扑克的视频| 少妇熟女aⅴ在线视频| 久久精品国产亚洲av香蕉五月| 成人手机av| 久久久久久久久中文| 国产成人系列免费观看| 国产熟女午夜一区二区三区| 少妇熟女aⅴ在线视频| 午夜免费观看网址| 亚洲免费av在线视频| 白带黄色成豆腐渣| 久久香蕉国产精品| 精品久久久久久成人av| 日韩欧美 国产精品| 免费看美女性在线毛片视频| 精品久久久久久久久久久久久 | 99riav亚洲国产免费| 久久青草综合色| 91成人精品电影| 母亲3免费完整高清在线观看| 男女那种视频在线观看| 国产一区二区三区视频了| av欧美777| 国产成人影院久久av| 国产单亲对白刺激| 精品不卡国产一区二区三区| 黑丝袜美女国产一区| 国产激情久久老熟女| 精品一区二区三区视频在线观看免费| 极品教师在线免费播放| 91老司机精品| 精华霜和精华液先用哪个| 18美女黄网站色大片免费观看| 国产激情欧美一区二区| 色综合亚洲欧美另类图片| 国产精品美女特级片免费视频播放器 | 亚洲精品国产精品久久久不卡| 夜夜夜夜夜久久久久| 亚洲精品一卡2卡三卡4卡5卡| 在线天堂中文资源库| 夜夜夜夜夜久久久久| 久久国产精品影院| 叶爱在线成人免费视频播放| 夜夜躁狠狠躁天天躁| 亚洲美女黄片视频| 中文字幕人妻熟女乱码| 国语自产精品视频在线第100页| 欧美黑人欧美精品刺激| 色综合欧美亚洲国产小说| 精品国产超薄肉色丝袜足j| 成在线人永久免费视频| 看免费av毛片| 国产精品久久久久久精品电影 | 久热这里只有精品99| 精品熟女少妇八av免费久了| 婷婷精品国产亚洲av| 久久精品影院6| 国产成人欧美在线观看| 久久久久久免费高清国产稀缺| 深夜精品福利| 国产av又大| 久久青草综合色| 日本免费一区二区三区高清不卡| 国产一区二区激情短视频| 久久精品国产99精品国产亚洲性色| 色婷婷久久久亚洲欧美| 国产精品久久久久久精品电影 | 日韩精品青青久久久久久| 久久香蕉国产精品| 亚洲男人天堂网一区| 色综合站精品国产| 在线观看日韩欧美| 热re99久久国产66热| 国产一区在线观看成人免费| 亚洲精品粉嫩美女一区| 给我免费播放毛片高清在线观看| 日本在线视频免费播放| 亚洲 欧美 日韩 在线 免费| av欧美777| 久久天躁狠狠躁夜夜2o2o| 欧美一级a爱片免费观看看 | 少妇被粗大的猛进出69影院| 18禁裸乳无遮挡免费网站照片 | 日本免费一区二区三区高清不卡| 成在线人永久免费视频| 99热这里只有精品一区 | 制服丝袜大香蕉在线| 亚洲成a人片在线一区二区| 国产亚洲精品综合一区在线观看 | 久久狼人影院| 一级a爱片免费观看的视频| 免费观看人在逋| 亚洲欧美精品综合久久99| 国产精品亚洲一级av第二区| 久久久久久久午夜电影| 成人三级做爰电影| 精品不卡国产一区二区三区| 国产一区二区在线av高清观看| www国产在线视频色| 精品不卡国产一区二区三区| 波多野结衣高清作品| 欧美色视频一区免费| 国产激情偷乱视频一区二区| 日韩成人在线观看一区二区三区| 亚洲黑人精品在线| 人人澡人人妻人| 亚洲精品中文字幕在线视频| 国产精品精品国产色婷婷| 亚洲人成伊人成综合网2020| 日日摸夜夜添夜夜添小说| 美女国产高潮福利片在线看| 国产一卡二卡三卡精品| 国产色视频综合| 两个人看的免费小视频| 在线观看www视频免费| 婷婷六月久久综合丁香| 国产精品日韩av在线免费观看| 午夜福利在线观看吧| 国产亚洲欧美98| 我的亚洲天堂| av免费在线观看网站| 久久亚洲精品不卡| 国产成人影院久久av| 窝窝影院91人妻| 丝袜在线中文字幕| 亚洲熟妇中文字幕五十中出| 亚洲色图 男人天堂 中文字幕| 美女扒开内裤让男人捅视频| 国产精品一区二区精品视频观看| www日本黄色视频网| 99久久无色码亚洲精品果冻| 丰满人妻熟妇乱又伦精品不卡| 在线观看66精品国产| 亚洲精品色激情综合| 免费在线观看黄色视频的| 免费在线观看视频国产中文字幕亚洲| 国产一区二区三区在线臀色熟女| 亚洲av中文字字幕乱码综合 | 人妻久久中文字幕网| 18禁观看日本| 亚洲一卡2卡3卡4卡5卡精品中文| 精品国产超薄肉色丝袜足j| 人成视频在线观看免费观看| 岛国在线观看网站| 久久精品国产亚洲av高清一级| 色综合婷婷激情| 亚洲av美国av| 亚洲第一电影网av| 久久99热这里只有精品18| 欧美丝袜亚洲另类 | 婷婷精品国产亚洲av| 少妇裸体淫交视频免费看高清 | 757午夜福利合集在线观看| 一个人免费在线观看的高清视频| 国产精品电影一区二区三区| 国产精品二区激情视频| 制服人妻中文乱码| 免费电影在线观看免费观看| 午夜久久久久精精品| 99re在线观看精品视频| 精品午夜福利视频在线观看一区| 国产成人系列免费观看| 久久精品国产99精品国产亚洲性色| 欧美+亚洲+日韩+国产| 精品欧美国产一区二区三| 女同久久另类99精品国产91| 亚洲五月天丁香| 少妇熟女aⅴ在线视频| 亚洲性夜色夜夜综合| 国产激情欧美一区二区| 日韩三级视频一区二区三区| 国产激情欧美一区二区| 欧美激情极品国产一区二区三区| 亚洲精品av麻豆狂野| 亚洲,欧美精品.| www.精华液| 美女免费视频网站| 欧美成人免费av一区二区三区| 女人高潮潮喷娇喘18禁视频| 免费在线观看成人毛片| 亚洲全国av大片| 久久九九热精品免费| 侵犯人妻中文字幕一二三四区| 99国产精品一区二区蜜桃av| www日本在线高清视频| 国产精品av久久久久免费| 真人做人爱边吃奶动态| 91成年电影在线观看| 免费在线观看日本一区| 免费高清视频大片| 成人欧美大片| 亚洲精品粉嫩美女一区| 久久 成人 亚洲| 亚洲国产欧洲综合997久久, | 女生性感内裤真人,穿戴方法视频| 国产亚洲欧美精品永久| 黄色 视频免费看| 此物有八面人人有两片| 在线永久观看黄色视频| 最近最新免费中文字幕在线| 国产1区2区3区精品| av欧美777| 一进一出抽搐动态| 又紧又爽又黄一区二区| avwww免费| 久99久视频精品免费| 亚洲精品一卡2卡三卡4卡5卡| 日本免费a在线| 色尼玛亚洲综合影院| 亚洲全国av大片| 一个人观看的视频www高清免费观看 | 亚洲成人免费电影在线观看| 久久中文字幕人妻熟女| 亚洲欧洲精品一区二区精品久久久| 精品国产超薄肉色丝袜足j| 国产熟女午夜一区二区三区| 国产高清videossex| 久久精品人妻少妇| 成人亚洲精品一区在线观看| 成年免费大片在线观看| 国产精品98久久久久久宅男小说| 国产91精品成人一区二区三区| 久久久精品国产亚洲av高清涩受| 免费观看精品视频网站| 国产乱人伦免费视频| 不卡一级毛片| 国产精品日韩av在线免费观看| 操出白浆在线播放| 黄色成人免费大全| 欧美日韩中文字幕国产精品一区二区三区| 老司机福利观看| 欧美亚洲日本最大视频资源| 免费看美女性在线毛片视频| 亚洲avbb在线观看| 在线观看www视频免费| 欧美日韩乱码在线| 怎么达到女性高潮| 岛国在线观看网站| 亚洲全国av大片| 国产精品久久久久久人妻精品电影| 免费观看精品视频网站| 欧美在线黄色| 国内精品久久久久久久电影| 午夜激情av网站| 国产三级黄色录像| 一卡2卡三卡四卡精品乱码亚洲| 亚洲三区欧美一区| 国产精品爽爽va在线观看网站 | 黄片播放在线免费| 国产真人三级小视频在线观看| 一级作爱视频免费观看| av免费在线观看网站| 国产精品自产拍在线观看55亚洲| 午夜久久久久精精品| 久久人人精品亚洲av| 男人操女人黄网站| 欧美不卡视频在线免费观看 | www日本黄色视频网| 日韩有码中文字幕| 国产亚洲精品综合一区在线观看 | 香蕉丝袜av| 日韩国内少妇激情av| 亚洲aⅴ乱码一区二区在线播放 | 欧美精品啪啪一区二区三区| 亚洲专区字幕在线| 国产三级在线视频| 少妇裸体淫交视频免费看高清 | 母亲3免费完整高清在线观看| 精品国产国语对白av| 男女床上黄色一级片免费看| 国产精品久久电影中文字幕| 嫩草影视91久久| 禁无遮挡网站| 欧美精品啪啪一区二区三区| 黑丝袜美女国产一区| 高潮久久久久久久久久久不卡| 久久婷婷人人爽人人干人人爱| 国产区一区二久久| 每晚都被弄得嗷嗷叫到高潮| 色精品久久人妻99蜜桃| 香蕉久久夜色| 亚洲精品粉嫩美女一区| 最近最新中文字幕大全电影3 | 99热6这里只有精品| 最好的美女福利视频网| 成熟少妇高潮喷水视频| 欧美午夜高清在线| 日日摸夜夜添夜夜添小说| 嫩草影院精品99| 亚洲专区字幕在线| 国产精品久久视频播放| 每晚都被弄得嗷嗷叫到高潮| 免费在线观看成人毛片| av在线播放免费不卡| 中文字幕最新亚洲高清| 长腿黑丝高跟| 满18在线观看网站| 欧美乱色亚洲激情| 男女视频在线观看网站免费 | 国产私拍福利视频在线观看| 可以在线观看毛片的网站| 久久久国产精品麻豆| 国产免费男女视频| 日本 欧美在线| 亚洲中文字幕日韩| 久9热在线精品视频| 日本五十路高清| 中出人妻视频一区二区| 在线永久观看黄色视频| 亚洲 欧美一区二区三区| 欧美性猛交黑人性爽| 日韩国内少妇激情av| 国产精品久久视频播放| 午夜视频精品福利| 久久精品国产99精品国产亚洲性色| 久久亚洲精品不卡| 脱女人内裤的视频| 国产伦人伦偷精品视频| 精品少妇一区二区三区视频日本电影| 亚洲人成电影免费在线| 无遮挡黄片免费观看| 亚洲男人的天堂狠狠| 久久中文看片网| 日日夜夜操网爽| 中文字幕另类日韩欧美亚洲嫩草| 黑人巨大精品欧美一区二区mp4| 欧美成狂野欧美在线观看| 日本免费一区二区三区高清不卡| 免费看日本二区| 国产欧美日韩一区二区三| 窝窝影院91人妻| 制服丝袜大香蕉在线| 一级a爱视频在线免费观看| 久久久久国内视频| 男人舔奶头视频| 99久久国产精品久久久| 天天躁夜夜躁狠狠躁躁| 亚洲中文av在线| www.www免费av| 久久久久免费精品人妻一区二区 | av视频在线观看入口| 中文资源天堂在线| 国产一区二区三区在线臀色熟女| 久久 成人 亚洲| 女生性感内裤真人,穿戴方法视频| 国产激情偷乱视频一区二区| 亚洲国产毛片av蜜桃av| 国产亚洲欧美98| 韩国精品一区二区三区| 黄色片一级片一级黄色片| 999久久久精品免费观看国产| 丝袜人妻中文字幕| 精品一区二区三区av网在线观看| 国产成+人综合+亚洲专区| 极品教师在线免费播放| 色婷婷久久久亚洲欧美| 亚洲黑人精品在线| 日日摸夜夜添夜夜添小说| or卡值多少钱| av电影中文网址| 51午夜福利影视在线观看| 亚洲一区中文字幕在线| 在线观看午夜福利视频| 久久精品91蜜桃| 麻豆成人午夜福利视频| 国内少妇人妻偷人精品xxx网站 | 亚洲精品久久成人aⅴ小说| 色综合欧美亚洲国产小说| 久久亚洲精品不卡| 黄色丝袜av网址大全| 久久精品国产清高在天天线| 伦理电影免费视频| 18禁黄网站禁片免费观看直播| 国产成人欧美在线观看| 美女大奶头视频| 国内精品久久久久久久电影| 亚洲精品国产区一区二| 真人一进一出gif抽搐免费| av免费在线观看网站| 欧美黑人巨大hd| 国产精品一区二区精品视频观看| 欧美黑人巨大hd| 一级黄色大片毛片| 18禁裸乳无遮挡免费网站照片 | www.999成人在线观看| 欧美最黄视频在线播放免费| 亚洲av成人av| 国产激情偷乱视频一区二区| 精品久久久久久久毛片微露脸| 91九色精品人成在线观看| 国产精品精品国产色婷婷| 国产成+人综合+亚洲专区| av视频在线观看入口| 国语自产精品视频在线第100页| 久久久久国内视频| 亚洲精品在线美女| 别揉我奶头~嗯~啊~动态视频| 日本免费一区二区三区高清不卡| 18禁黄网站禁片午夜丰满| 97碰自拍视频| 亚洲成a人片在线一区二区| 午夜福利18| 亚洲精品国产一区二区精华液| 国产精品香港三级国产av潘金莲| 一卡2卡三卡四卡精品乱码亚洲| av在线播放免费不卡| 国产激情久久老熟女| 亚洲美女黄片视频| 在线av久久热| 国产三级在线视频| 久久久国产成人精品二区| 日韩成人在线观看一区二区三区| 女人被狂操c到高潮| 亚洲第一av免费看| svipshipincom国产片| 免费av毛片视频| 男女那种视频在线观看| 老汉色av国产亚洲站长工具| 黄色毛片三级朝国网站| 淫秽高清视频在线观看| 又黄又粗又硬又大视频| 99久久无色码亚洲精品果冻| 国产亚洲精品久久久久久毛片| 制服诱惑二区| 999久久久国产精品视频| 亚洲欧美激情综合另类| 欧美av亚洲av综合av国产av| www国产在线视频色| 国产av一区二区精品久久| 国产又黄又爽又无遮挡在线| 一二三四社区在线视频社区8| 中文字幕高清在线视频| 欧美激情高清一区二区三区| 日本熟妇午夜| a级毛片在线看网站| 欧美日韩黄片免| 可以在线观看的亚洲视频| 亚洲美女黄片视频| 国产成人啪精品午夜网站| 日韩欧美三级三区| av在线天堂中文字幕| 精品国产国语对白av| 黄频高清免费视频| 国产精品98久久久久久宅男小说| 一本久久中文字幕| 悠悠久久av| av片东京热男人的天堂| 巨乳人妻的诱惑在线观看| 久久久久精品国产欧美久久久| 国产av不卡久久| 美国免费a级毛片| 亚洲aⅴ乱码一区二区在线播放 | 亚洲一区中文字幕在线| 十八禁人妻一区二区| 日本免费一区二区三区高清不卡| 成人国语在线视频| 日本 av在线| 啦啦啦观看免费观看视频高清| a级毛片a级免费在线| 又黄又爽又免费观看的视频| 麻豆av在线久日| 午夜视频精品福利| 少妇的丰满在线观看| 免费av毛片视频| 日日摸夜夜添夜夜添小说| 国产黄色小视频在线观看| 又紧又爽又黄一区二区| 操出白浆在线播放| 啦啦啦免费观看视频1| 国产高清videossex| 亚洲狠狠婷婷综合久久图片| 久久香蕉精品热| 色综合欧美亚洲国产小说| 日韩视频一区二区在线观看| 国产午夜福利久久久久久|