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

    六由度并聯(lián)機器人位姿控制分析與仿真

    2014-07-08 02:19:31王栿棟鄒華勇
    科技創(chuàng)新與應(yīng)用 2014年21期
    關(guān)鍵詞:坐標(biāo)仿真位姿

    王栿棟++鄒華勇

    摘 要:利用剛體空間六個自由度運動原理,以Stewart平臺為例,介紹了并聯(lián)機構(gòu)正反解的基本原理。并針對并聯(lián)機構(gòu)更困難的正解問題,利用Matlab設(shè)計一個圖形加數(shù)字坐標(biāo)仿真程序,可以方便地求得并聯(lián)機構(gòu)符合實際客觀條件的位置解,并直觀地觀察各桿件之間的空間關(guān)系。

    關(guān)鍵詞:正解;反解;坐標(biāo);位姿;仿真

    引言

    機器人從20世紀(jì)誕生以來,以迅速成為一門應(yīng)用十分廣泛的技術(shù)。機器人技術(shù)在工程中的應(yīng)用日益廣泛,即工業(yè)機器人。其中研究比較早,目前工程界應(yīng)用比較成熟的是串聯(lián)機器人。它是一種能夠自動定位控制、多自由度、各大大小小的機械臂以串聯(lián)的形式聯(lián)接起來。這種典型的開環(huán)機構(gòu)具有結(jié)構(gòu)簡單,成本低,控制簡單,運動空間大等優(yōu)點。但伴隨著其結(jié)構(gòu)特性的缺點則是承載能力有限、機構(gòu)剛度差,各運動機構(gòu)疊加誤差大,運動慣性大。而在80年代以來隨機計算機技術(shù)的發(fā)展,并聯(lián)機器人的研究成為了新的熱點。并聯(lián)機器人與串聯(lián)機器人正好可以形成互補,它的優(yōu)點極為串聯(lián)機器人的缺點,相應(yīng)它的缺點則是串聯(lián)機器人的優(yōu)點,因此并聯(lián)機器人在飛行模擬器,空間對接器、裝配生產(chǎn)線等需要高精度,高穩(wěn)定性,高速的場合應(yīng)用越來越廣泛。但串聯(lián)機器人與并聯(lián)機器人的基礎(chǔ)理論基本一樣,所以文章以六自由度并聯(lián)機器人最典型的結(jié)構(gòu)——Stewart平臺簡析六自由度機器手入手,在運動學(xué)的位置正反解。文章總結(jié)的內(nèi)容中略去了有關(guān)機器人的數(shù)學(xué)基礎(chǔ)知識和基礎(chǔ)力學(xué)知識,假設(shè)已經(jīng)具備一定的矢量代數(shù),矩陣論,和工學(xué)力學(xué)知識,直接進(jìn)入對機器人本身的運動進(jìn)行了分析。并約定,一般黑色字體變量為矢量,非黑字體變量為標(biāo)量。

    1 六自由度并聯(lián)機器人運動學(xué)

    研究操作手的幾何學(xué)要區(qū)分兩個問題,即在同一系統(tǒng)下的運動學(xué)正問題和逆問題,所謂的正問題就是給定關(guān)節(jié)變量前提下,求機器人執(zhí)行終端的位姿,這個求解過程就叫正解。而逆問題反過來,設(shè)定定機器人執(zhí)行終端的位置,反求各關(guān)節(jié)變量,這個求解過程叫反解。

    見圖1a所示并聯(lián)操作手,關(guān)于這個并聯(lián)機構(gòu)的反解其實是容易的,這里不再敘述,有興趣者可參閱本文列出得參考文獻(xiàn)?,F(xiàn)在主要對其行運動學(xué)正問題分析,但在應(yīng)用到如圖1的并聯(lián)機器人中時,他的正解運動學(xué)卻并不那么簡單。

    在圖1a中,考慮三角形AiSiBi,i=1~3,角標(biāo)i表示第i對腿。當(dāng)6個腿的長度固定,把平臺M去掉,三角形AiSiBi只能繞AiBi軸旋轉(zhuǎn)。這樣,我們可以用一個長度為li的腿代替長度為qia和qib的一對腿,這條腿通過一個繞AiBi的轉(zhuǎn)動關(guān)節(jié)和六角基座平臺β相連接。結(jié)構(gòu)簡化后的結(jié)果如圖2所示,它在運動學(xué)上的等價于圖1a的原來結(jié)構(gòu)。

    圖1 6自由度飛行器模擬器

    a)基本結(jié)構(gòu) b)三角形活動平臺、六角形固定平臺

    然后引入坐標(biāo)系Fi,它的原點設(shè)在第i個腿和基座平臺β連接的點口Oi處,并作如下約定:

    對i=1,2,3,

    Oi是連接旋轉(zhuǎn)關(guān)節(jié)的中心的集合;

    Xi的方向是從Ai指向Bi

    Yi的選擇是由Zi垂直于6邊形固定平臺指向上來確定的,也就是X與Y坐標(biāo)向量的矢量積方向朝上。

    之后來確定從六邊形的中心O到移動平臺的3個頂點S1,S2 和S3的位置矢量。

    我們需要確定li和Oi。參考圖2和圖3,其中ai和bi表示Ai,和Bi相對O為原點基座坐標(biāo)的矢量。對于i=1,2,3,有

    因此ui是由Ai指向Bi的單位矢量。而且,原點Oi的位置由矢量oi確定,oi的表述如下:

    (1)

    進(jìn)而,設(shè)si為在坐標(biāo)系Fi(Oi,Xi,Yi,Zi)中Si的位置矢量,則

    (2)

    圖2 等效簡化機構(gòu)

    圖3 一條腿與一對腿的等效 圖4 坐標(biāo)系F0與坐標(biāo)系Fi的關(guān)系

    現(xiàn)在,將一個坐標(biāo)系空間F0定義為原點 O和X,Y,軸在六邊形固定平面上,其上的 Xi和 Yi的關(guān)系如圖4所示.當(dāng)在坐標(biāo)系F0中表達(dá),si的形式如下

    (2.3)

    這里,[Ri]0是在坐標(biāo)系F0中表達(dá)的從坐標(biāo)系F0到坐標(biāo)系Fi的旋轉(zhuǎn)矩陣,給出如下

    (2.4)

    參考圖(4)有

    (2.5)

    (2.6)

    將式(4)~式(6)代入式(3)得到

    (2.7)

    這里,Oi由式(8)給出。

    因為活動三角平臺三個頂點之間的距離是固定的,所以位置矢量S1,S2,S3必須滿足以下約束

    (8a) (8b) (8c)

    展開后,式(8a-c)取形式

    (9a)

    (9b)

    (9c)

    Di,Ei,F(xiàn)i都是已知數(shù)據(jù),i為1~5, 他們可以通過計算得到,經(jīng)過整理后具體形式如下。

    其中 ,是一個2*2反對稱矩陣,這樣三個方程已經(jīng)

    列好,并對應(yīng)三個未知數(shù)?準(zhǔn)1,?準(zhǔn)2,?準(zhǔn)3。方程數(shù)與未知數(shù)數(shù)量一致,在理論上是可以有解的,但這個數(shù)據(jù)運算也好,方程本身的復(fù)雜度也好,要把它迅速正確求出代數(shù)解是很困難的。如果?準(zhǔn)1,?準(zhǔn)2,?準(zhǔn)3三個角符合實際情況的解(0~?仔)能夠得到,則通過式(7)可得到確定三角運動平臺三個頂點的矢量,從而根據(jù)三點確定一面的原理,來確定并聯(lián)機器人執(zhí)行端(三角活動平臺)的空間位姿。接下去就是如何快速正確地解這方程。

    2 仿真程序

    2.1 相關(guān)仿真參數(shù)簡述

    按照上述原理,編寫計算程序,實現(xiàn)坐標(biāo)顯示和仿真圖示。如圖5a所示,綠色的正六邊形代表固定底盤,其大小由邊長b確定。上方正三角形是活動平臺,其大小由邊長a確定。L1~L6表示六條機器腿的長度,同一顏色為一組。a1~a3表示活動平臺三個頂點,其坐標(biāo)通過計算已顯示在左上角。通過計算還可以識別活動平臺法線朝向,如果是朝上則顯示青色(如圖5a),如果朝下則顯示品紅色(如圖5b),一般多數(shù)Stewart平臺的活動平臺正常使用情況桿多為法線方向一致朝上。而在此處是為了分析活動平臺所有可能的情況而設(shè)計的仿真程序。

    除此之外還可以利用圖形界面的放大和旋轉(zhuǎn)功能來觀察各桿件之間的位置關(guān)系,從而判斷其互相干涉情況。

    2.2 仿真過程羅列

    請看圖6。整個仿真過程,分為仿真計算、仿真圖示兩個大部分。其中仿真計算是最根本,而且執(zhí)行過程相對耗時,仿真圖示是根據(jù)初始輸入?yún)?shù)把計算結(jié)果來顯示出來。

    圖6 仿真計算結(jié)構(gòu)簡圖

    2.3 仿真計算結(jié)構(gòu)

    為了使得整個仿真過程快速度,作者采用了數(shù)值迭代法。因為在理論上采用符號計算(復(fù)數(shù)域求解)是最全最精確的,而事實上在matlab中符號計算對于參數(shù)較多,相對復(fù)雜的方程組,計算相當(dāng)耗時。而且返回的結(jié)果不一定齊全。

    為了使仿真計算結(jié)果沒有遺漏,作者在數(shù)值迭代計算中合理分布了三個方程共9個初值,使得計算結(jié)果沒有遺漏,并且避免三角函數(shù)因周期而產(chǎn)生的重復(fù)解。

    為了使仿真計算結(jié)果符合客觀實際情況,作者采用了數(shù)值篩選,目的之一是去除不符合客觀物理空間復(fù)數(shù)解。目的之二是去除不符合并聯(lián)機構(gòu)實際結(jié)構(gòu)的負(fù)數(shù)解。

    整個計算機構(gòu)請看圖7所示。中間計算解值是指圖(3)中α1~α3三個角度。

    圖3.3仿真計算結(jié)構(gòu)

    2.4 仿真舉例

    當(dāng)輸入?yún)?shù)為a=6,b=9,L1=8,L2=8,L3=8,L4=8,L5=13,L6=15.時候仿真結(jié)果如下,見圖8,結(jié)果顯示,本初始條件系下,并聯(lián)結(jié)構(gòu)有兩個位姿解,其中一個法線朝上,一個法線朝下。相應(yīng)的活動平臺頂點坐標(biāo)也顯示在左上角。

    3 結(jié)束語

    通過對并聯(lián)機構(gòu)運動學(xué)的分析,設(shè)計用來求并聯(lián)機構(gòu)正解位置的仿真程序,這一來不但可以快速正確地求得在不同機構(gòu)參數(shù)下的位置解,而且可以通過圖示清楚直觀地了解各工作臺面的和機器腿在三維空間中的實際姿態(tài)。

    參考文獻(xiàn)

    [1]宋偉剛.機器人學(xué)[J].科學(xué)出版社,2007.

    [2]黃真,孔令富,方躍法.并聯(lián)機器人機構(gòu)學(xué)理論及控制[J].機械工業(yè)出版社,1997.

    [3]唐錫寬,金德聞.機械動力學(xué)[J].高等教育出版社,1983.endprint

    摘 要:利用剛體空間六個自由度運動原理,以Stewart平臺為例,介紹了并聯(lián)機構(gòu)正反解的基本原理。并針對并聯(lián)機構(gòu)更困難的正解問題,利用Matlab設(shè)計一個圖形加數(shù)字坐標(biāo)仿真程序,可以方便地求得并聯(lián)機構(gòu)符合實際客觀條件的位置解,并直觀地觀察各桿件之間的空間關(guān)系。

    關(guān)鍵詞:正解;反解;坐標(biāo);位姿;仿真

    引言

    機器人從20世紀(jì)誕生以來,以迅速成為一門應(yīng)用十分廣泛的技術(shù)。機器人技術(shù)在工程中的應(yīng)用日益廣泛,即工業(yè)機器人。其中研究比較早,目前工程界應(yīng)用比較成熟的是串聯(lián)機器人。它是一種能夠自動定位控制、多自由度、各大大小小的機械臂以串聯(lián)的形式聯(lián)接起來。這種典型的開環(huán)機構(gòu)具有結(jié)構(gòu)簡單,成本低,控制簡單,運動空間大等優(yōu)點。但伴隨著其結(jié)構(gòu)特性的缺點則是承載能力有限、機構(gòu)剛度差,各運動機構(gòu)疊加誤差大,運動慣性大。而在80年代以來隨機計算機技術(shù)的發(fā)展,并聯(lián)機器人的研究成為了新的熱點。并聯(lián)機器人與串聯(lián)機器人正好可以形成互補,它的優(yōu)點極為串聯(lián)機器人的缺點,相應(yīng)它的缺點則是串聯(lián)機器人的優(yōu)點,因此并聯(lián)機器人在飛行模擬器,空間對接器、裝配生產(chǎn)線等需要高精度,高穩(wěn)定性,高速的場合應(yīng)用越來越廣泛。但串聯(lián)機器人與并聯(lián)機器人的基礎(chǔ)理論基本一樣,所以文章以六自由度并聯(lián)機器人最典型的結(jié)構(gòu)——Stewart平臺簡析六自由度機器手入手,在運動學(xué)的位置正反解。文章總結(jié)的內(nèi)容中略去了有關(guān)機器人的數(shù)學(xué)基礎(chǔ)知識和基礎(chǔ)力學(xué)知識,假設(shè)已經(jīng)具備一定的矢量代數(shù),矩陣論,和工學(xué)力學(xué)知識,直接進(jìn)入對機器人本身的運動進(jìn)行了分析。并約定,一般黑色字體變量為矢量,非黑字體變量為標(biāo)量。

    1 六自由度并聯(lián)機器人運動學(xué)

    研究操作手的幾何學(xué)要區(qū)分兩個問題,即在同一系統(tǒng)下的運動學(xué)正問題和逆問題,所謂的正問題就是給定關(guān)節(jié)變量前提下,求機器人執(zhí)行終端的位姿,這個求解過程就叫正解。而逆問題反過來,設(shè)定定機器人執(zhí)行終端的位置,反求各關(guān)節(jié)變量,這個求解過程叫反解。

    見圖1a所示并聯(lián)操作手,關(guān)于這個并聯(lián)機構(gòu)的反解其實是容易的,這里不再敘述,有興趣者可參閱本文列出得參考文獻(xiàn)?,F(xiàn)在主要對其行運動學(xué)正問題分析,但在應(yīng)用到如圖1的并聯(lián)機器人中時,他的正解運動學(xué)卻并不那么簡單。

    在圖1a中,考慮三角形AiSiBi,i=1~3,角標(biāo)i表示第i對腿。當(dāng)6個腿的長度固定,把平臺M去掉,三角形AiSiBi只能繞AiBi軸旋轉(zhuǎn)。這樣,我們可以用一個長度為li的腿代替長度為qia和qib的一對腿,這條腿通過一個繞AiBi的轉(zhuǎn)動關(guān)節(jié)和六角基座平臺β相連接。結(jié)構(gòu)簡化后的結(jié)果如圖2所示,它在運動學(xué)上的等價于圖1a的原來結(jié)構(gòu)。

    圖1 6自由度飛行器模擬器

    a)基本結(jié)構(gòu) b)三角形活動平臺、六角形固定平臺

    然后引入坐標(biāo)系Fi,它的原點設(shè)在第i個腿和基座平臺β連接的點口Oi處,并作如下約定:

    對i=1,2,3,

    Oi是連接旋轉(zhuǎn)關(guān)節(jié)的中心的集合;

    Xi的方向是從Ai指向Bi

    Yi的選擇是由Zi垂直于6邊形固定平臺指向上來確定的,也就是X與Y坐標(biāo)向量的矢量積方向朝上。

    之后來確定從六邊形的中心O到移動平臺的3個頂點S1,S2 和S3的位置矢量。

    我們需要確定li和Oi。參考圖2和圖3,其中ai和bi表示Ai,和Bi相對O為原點基座坐標(biāo)的矢量。對于i=1,2,3,有

    因此ui是由Ai指向Bi的單位矢量。而且,原點Oi的位置由矢量oi確定,oi的表述如下:

    (1)

    進(jìn)而,設(shè)si為在坐標(biāo)系Fi(Oi,Xi,Yi,Zi)中Si的位置矢量,則

    (2)

    圖2 等效簡化機構(gòu)

    圖3 一條腿與一對腿的等效 圖4 坐標(biāo)系F0與坐標(biāo)系Fi的關(guān)系

    現(xiàn)在,將一個坐標(biāo)系空間F0定義為原點 O和X,Y,軸在六邊形固定平面上,其上的 Xi和 Yi的關(guān)系如圖4所示.當(dāng)在坐標(biāo)系F0中表達(dá),si的形式如下

    (2.3)

    這里,[Ri]0是在坐標(biāo)系F0中表達(dá)的從坐標(biāo)系F0到坐標(biāo)系Fi的旋轉(zhuǎn)矩陣,給出如下

    (2.4)

    參考圖(4)有

    (2.5)

    (2.6)

    將式(4)~式(6)代入式(3)得到

    (2.7)

    這里,Oi由式(8)給出。

    因為活動三角平臺三個頂點之間的距離是固定的,所以位置矢量S1,S2,S3必須滿足以下約束

    (8a) (8b) (8c)

    展開后,式(8a-c)取形式

    (9a)

    (9b)

    (9c)

    Di,Ei,F(xiàn)i都是已知數(shù)據(jù),i為1~5, 他們可以通過計算得到,經(jīng)過整理后具體形式如下。

    其中 ,是一個2*2反對稱矩陣,這樣三個方程已經(jīng)

    列好,并對應(yīng)三個未知數(shù)?準(zhǔn)1,?準(zhǔn)2,?準(zhǔn)3。方程數(shù)與未知數(shù)數(shù)量一致,在理論上是可以有解的,但這個數(shù)據(jù)運算也好,方程本身的復(fù)雜度也好,要把它迅速正確求出代數(shù)解是很困難的。如果?準(zhǔn)1,?準(zhǔn)2,?準(zhǔn)3三個角符合實際情況的解(0~?仔)能夠得到,則通過式(7)可得到確定三角運動平臺三個頂點的矢量,從而根據(jù)三點確定一面的原理,來確定并聯(lián)機器人執(zhí)行端(三角活動平臺)的空間位姿。接下去就是如何快速正確地解這方程。

    2 仿真程序

    2.1 相關(guān)仿真參數(shù)簡述

    按照上述原理,編寫計算程序,實現(xiàn)坐標(biāo)顯示和仿真圖示。如圖5a所示,綠色的正六邊形代表固定底盤,其大小由邊長b確定。上方正三角形是活動平臺,其大小由邊長a確定。L1~L6表示六條機器腿的長度,同一顏色為一組。a1~a3表示活動平臺三個頂點,其坐標(biāo)通過計算已顯示在左上角。通過計算還可以識別活動平臺法線朝向,如果是朝上則顯示青色(如圖5a),如果朝下則顯示品紅色(如圖5b),一般多數(shù)Stewart平臺的活動平臺正常使用情況桿多為法線方向一致朝上。而在此處是為了分析活動平臺所有可能的情況而設(shè)計的仿真程序。

    除此之外還可以利用圖形界面的放大和旋轉(zhuǎn)功能來觀察各桿件之間的位置關(guān)系,從而判斷其互相干涉情況。

    2.2 仿真過程羅列

    請看圖6。整個仿真過程,分為仿真計算、仿真圖示兩個大部分。其中仿真計算是最根本,而且執(zhí)行過程相對耗時,仿真圖示是根據(jù)初始輸入?yún)?shù)把計算結(jié)果來顯示出來。

    圖6 仿真計算結(jié)構(gòu)簡圖

    2.3 仿真計算結(jié)構(gòu)

    為了使得整個仿真過程快速度,作者采用了數(shù)值迭代法。因為在理論上采用符號計算(復(fù)數(shù)域求解)是最全最精確的,而事實上在matlab中符號計算對于參數(shù)較多,相對復(fù)雜的方程組,計算相當(dāng)耗時。而且返回的結(jié)果不一定齊全。

    為了使仿真計算結(jié)果沒有遺漏,作者在數(shù)值迭代計算中合理分布了三個方程共9個初值,使得計算結(jié)果沒有遺漏,并且避免三角函數(shù)因周期而產(chǎn)生的重復(fù)解。

    為了使仿真計算結(jié)果符合客觀實際情況,作者采用了數(shù)值篩選,目的之一是去除不符合客觀物理空間復(fù)數(shù)解。目的之二是去除不符合并聯(lián)機構(gòu)實際結(jié)構(gòu)的負(fù)數(shù)解。

    整個計算機構(gòu)請看圖7所示。中間計算解值是指圖(3)中α1~α3三個角度。

    圖3.3仿真計算結(jié)構(gòu)

    2.4 仿真舉例

    當(dāng)輸入?yún)?shù)為a=6,b=9,L1=8,L2=8,L3=8,L4=8,L5=13,L6=15.時候仿真結(jié)果如下,見圖8,結(jié)果顯示,本初始條件系下,并聯(lián)結(jié)構(gòu)有兩個位姿解,其中一個法線朝上,一個法線朝下。相應(yīng)的活動平臺頂點坐標(biāo)也顯示在左上角。

    3 結(jié)束語

    通過對并聯(lián)機構(gòu)運動學(xué)的分析,設(shè)計用來求并聯(lián)機構(gòu)正解位置的仿真程序,這一來不但可以快速正確地求得在不同機構(gòu)參數(shù)下的位置解,而且可以通過圖示清楚直觀地了解各工作臺面的和機器腿在三維空間中的實際姿態(tài)。

    參考文獻(xiàn)

    [1]宋偉剛.機器人學(xué)[J].科學(xué)出版社,2007.

    [2]黃真,孔令富,方躍法.并聯(lián)機器人機構(gòu)學(xué)理論及控制[J].機械工業(yè)出版社,1997.

    [3]唐錫寬,金德聞.機械動力學(xué)[J].高等教育出版社,1983.endprint

    摘 要:利用剛體空間六個自由度運動原理,以Stewart平臺為例,介紹了并聯(lián)機構(gòu)正反解的基本原理。并針對并聯(lián)機構(gòu)更困難的正解問題,利用Matlab設(shè)計一個圖形加數(shù)字坐標(biāo)仿真程序,可以方便地求得并聯(lián)機構(gòu)符合實際客觀條件的位置解,并直觀地觀察各桿件之間的空間關(guān)系。

    關(guān)鍵詞:正解;反解;坐標(biāo);位姿;仿真

    引言

    機器人從20世紀(jì)誕生以來,以迅速成為一門應(yīng)用十分廣泛的技術(shù)。機器人技術(shù)在工程中的應(yīng)用日益廣泛,即工業(yè)機器人。其中研究比較早,目前工程界應(yīng)用比較成熟的是串聯(lián)機器人。它是一種能夠自動定位控制、多自由度、各大大小小的機械臂以串聯(lián)的形式聯(lián)接起來。這種典型的開環(huán)機構(gòu)具有結(jié)構(gòu)簡單,成本低,控制簡單,運動空間大等優(yōu)點。但伴隨著其結(jié)構(gòu)特性的缺點則是承載能力有限、機構(gòu)剛度差,各運動機構(gòu)疊加誤差大,運動慣性大。而在80年代以來隨機計算機技術(shù)的發(fā)展,并聯(lián)機器人的研究成為了新的熱點。并聯(lián)機器人與串聯(lián)機器人正好可以形成互補,它的優(yōu)點極為串聯(lián)機器人的缺點,相應(yīng)它的缺點則是串聯(lián)機器人的優(yōu)點,因此并聯(lián)機器人在飛行模擬器,空間對接器、裝配生產(chǎn)線等需要高精度,高穩(wěn)定性,高速的場合應(yīng)用越來越廣泛。但串聯(lián)機器人與并聯(lián)機器人的基礎(chǔ)理論基本一樣,所以文章以六自由度并聯(lián)機器人最典型的結(jié)構(gòu)——Stewart平臺簡析六自由度機器手入手,在運動學(xué)的位置正反解。文章總結(jié)的內(nèi)容中略去了有關(guān)機器人的數(shù)學(xué)基礎(chǔ)知識和基礎(chǔ)力學(xué)知識,假設(shè)已經(jīng)具備一定的矢量代數(shù),矩陣論,和工學(xué)力學(xué)知識,直接進(jìn)入對機器人本身的運動進(jìn)行了分析。并約定,一般黑色字體變量為矢量,非黑字體變量為標(biāo)量。

    1 六自由度并聯(lián)機器人運動學(xué)

    研究操作手的幾何學(xué)要區(qū)分兩個問題,即在同一系統(tǒng)下的運動學(xué)正問題和逆問題,所謂的正問題就是給定關(guān)節(jié)變量前提下,求機器人執(zhí)行終端的位姿,這個求解過程就叫正解。而逆問題反過來,設(shè)定定機器人執(zhí)行終端的位置,反求各關(guān)節(jié)變量,這個求解過程叫反解。

    見圖1a所示并聯(lián)操作手,關(guān)于這個并聯(lián)機構(gòu)的反解其實是容易的,這里不再敘述,有興趣者可參閱本文列出得參考文獻(xiàn)?,F(xiàn)在主要對其行運動學(xué)正問題分析,但在應(yīng)用到如圖1的并聯(lián)機器人中時,他的正解運動學(xué)卻并不那么簡單。

    在圖1a中,考慮三角形AiSiBi,i=1~3,角標(biāo)i表示第i對腿。當(dāng)6個腿的長度固定,把平臺M去掉,三角形AiSiBi只能繞AiBi軸旋轉(zhuǎn)。這樣,我們可以用一個長度為li的腿代替長度為qia和qib的一對腿,這條腿通過一個繞AiBi的轉(zhuǎn)動關(guān)節(jié)和六角基座平臺β相連接。結(jié)構(gòu)簡化后的結(jié)果如圖2所示,它在運動學(xué)上的等價于圖1a的原來結(jié)構(gòu)。

    圖1 6自由度飛行器模擬器

    a)基本結(jié)構(gòu) b)三角形活動平臺、六角形固定平臺

    然后引入坐標(biāo)系Fi,它的原點設(shè)在第i個腿和基座平臺β連接的點口Oi處,并作如下約定:

    對i=1,2,3,

    Oi是連接旋轉(zhuǎn)關(guān)節(jié)的中心的集合;

    Xi的方向是從Ai指向Bi

    Yi的選擇是由Zi垂直于6邊形固定平臺指向上來確定的,也就是X與Y坐標(biāo)向量的矢量積方向朝上。

    之后來確定從六邊形的中心O到移動平臺的3個頂點S1,S2 和S3的位置矢量。

    我們需要確定li和Oi。參考圖2和圖3,其中ai和bi表示Ai,和Bi相對O為原點基座坐標(biāo)的矢量。對于i=1,2,3,有

    因此ui是由Ai指向Bi的單位矢量。而且,原點Oi的位置由矢量oi確定,oi的表述如下:

    (1)

    進(jìn)而,設(shè)si為在坐標(biāo)系Fi(Oi,Xi,Yi,Zi)中Si的位置矢量,則

    (2)

    圖2 等效簡化機構(gòu)

    圖3 一條腿與一對腿的等效 圖4 坐標(biāo)系F0與坐標(biāo)系Fi的關(guān)系

    現(xiàn)在,將一個坐標(biāo)系空間F0定義為原點 O和X,Y,軸在六邊形固定平面上,其上的 Xi和 Yi的關(guān)系如圖4所示.當(dāng)在坐標(biāo)系F0中表達(dá),si的形式如下

    (2.3)

    這里,[Ri]0是在坐標(biāo)系F0中表達(dá)的從坐標(biāo)系F0到坐標(biāo)系Fi的旋轉(zhuǎn)矩陣,給出如下

    (2.4)

    參考圖(4)有

    (2.5)

    (2.6)

    將式(4)~式(6)代入式(3)得到

    (2.7)

    這里,Oi由式(8)給出。

    因為活動三角平臺三個頂點之間的距離是固定的,所以位置矢量S1,S2,S3必須滿足以下約束

    (8a) (8b) (8c)

    展開后,式(8a-c)取形式

    (9a)

    (9b)

    (9c)

    Di,Ei,F(xiàn)i都是已知數(shù)據(jù),i為1~5, 他們可以通過計算得到,經(jīng)過整理后具體形式如下。

    其中 ,是一個2*2反對稱矩陣,這樣三個方程已經(jīng)

    列好,并對應(yīng)三個未知數(shù)?準(zhǔn)1,?準(zhǔn)2,?準(zhǔn)3。方程數(shù)與未知數(shù)數(shù)量一致,在理論上是可以有解的,但這個數(shù)據(jù)運算也好,方程本身的復(fù)雜度也好,要把它迅速正確求出代數(shù)解是很困難的。如果?準(zhǔn)1,?準(zhǔn)2,?準(zhǔn)3三個角符合實際情況的解(0~?仔)能夠得到,則通過式(7)可得到確定三角運動平臺三個頂點的矢量,從而根據(jù)三點確定一面的原理,來確定并聯(lián)機器人執(zhí)行端(三角活動平臺)的空間位姿。接下去就是如何快速正確地解這方程。

    2 仿真程序

    2.1 相關(guān)仿真參數(shù)簡述

    按照上述原理,編寫計算程序,實現(xiàn)坐標(biāo)顯示和仿真圖示。如圖5a所示,綠色的正六邊形代表固定底盤,其大小由邊長b確定。上方正三角形是活動平臺,其大小由邊長a確定。L1~L6表示六條機器腿的長度,同一顏色為一組。a1~a3表示活動平臺三個頂點,其坐標(biāo)通過計算已顯示在左上角。通過計算還可以識別活動平臺法線朝向,如果是朝上則顯示青色(如圖5a),如果朝下則顯示品紅色(如圖5b),一般多數(shù)Stewart平臺的活動平臺正常使用情況桿多為法線方向一致朝上。而在此處是為了分析活動平臺所有可能的情況而設(shè)計的仿真程序。

    除此之外還可以利用圖形界面的放大和旋轉(zhuǎn)功能來觀察各桿件之間的位置關(guān)系,從而判斷其互相干涉情況。

    2.2 仿真過程羅列

    請看圖6。整個仿真過程,分為仿真計算、仿真圖示兩個大部分。其中仿真計算是最根本,而且執(zhí)行過程相對耗時,仿真圖示是根據(jù)初始輸入?yún)?shù)把計算結(jié)果來顯示出來。

    圖6 仿真計算結(jié)構(gòu)簡圖

    2.3 仿真計算結(jié)構(gòu)

    為了使得整個仿真過程快速度,作者采用了數(shù)值迭代法。因為在理論上采用符號計算(復(fù)數(shù)域求解)是最全最精確的,而事實上在matlab中符號計算對于參數(shù)較多,相對復(fù)雜的方程組,計算相當(dāng)耗時。而且返回的結(jié)果不一定齊全。

    為了使仿真計算結(jié)果沒有遺漏,作者在數(shù)值迭代計算中合理分布了三個方程共9個初值,使得計算結(jié)果沒有遺漏,并且避免三角函數(shù)因周期而產(chǎn)生的重復(fù)解。

    為了使仿真計算結(jié)果符合客觀實際情況,作者采用了數(shù)值篩選,目的之一是去除不符合客觀物理空間復(fù)數(shù)解。目的之二是去除不符合并聯(lián)機構(gòu)實際結(jié)構(gòu)的負(fù)數(shù)解。

    整個計算機構(gòu)請看圖7所示。中間計算解值是指圖(3)中α1~α3三個角度。

    圖3.3仿真計算結(jié)構(gòu)

    2.4 仿真舉例

    當(dāng)輸入?yún)?shù)為a=6,b=9,L1=8,L2=8,L3=8,L4=8,L5=13,L6=15.時候仿真結(jié)果如下,見圖8,結(jié)果顯示,本初始條件系下,并聯(lián)結(jié)構(gòu)有兩個位姿解,其中一個法線朝上,一個法線朝下。相應(yīng)的活動平臺頂點坐標(biāo)也顯示在左上角。

    3 結(jié)束語

    通過對并聯(lián)機構(gòu)運動學(xué)的分析,設(shè)計用來求并聯(lián)機構(gòu)正解位置的仿真程序,這一來不但可以快速正確地求得在不同機構(gòu)參數(shù)下的位置解,而且可以通過圖示清楚直觀地了解各工作臺面的和機器腿在三維空間中的實際姿態(tài)。

    參考文獻(xiàn)

    [1]宋偉剛.機器人學(xué)[J].科學(xué)出版社,2007.

    [2]黃真,孔令富,方躍法.并聯(lián)機器人機構(gòu)學(xué)理論及控制[J].機械工業(yè)出版社,1997.

    [3]唐錫寬,金德聞.機械動力學(xué)[J].高等教育出版社,1983.endprint

    猜你喜歡
    坐標(biāo)仿真位姿
    Buck開關(guān)變換器的基本參數(shù)設(shè)計及仿真分析
    試析PLC控制下的自動化立體倉庫仿真情況分析
    基于共面直線迭代加權(quán)最小二乘的相機位姿估計
    基于CAD模型的單目六自由度位姿測量
    基于MADYMO的航空座椅約束系統(tǒng)優(yōu)化設(shè)計
    科技視界(2016年18期)2016-11-03 21:44:44
    中國體態(tài)假人模型與FAA Hybrid Ⅲ 型假人模型沖擊差異性分析
    科技視界(2016年18期)2016-11-03 20:31:49
    用單位圓記憶誘導(dǎo)公式
    MAPGIS的制圖技巧
    基于物聯(lián)網(wǎng)技術(shù)的“蔬菜”溯源體系探索
    小型四旋翼飛行器位姿建模及其仿真
    美女主播在线视频| 午夜福利视频精品| 日本色播在线视频| 久久女婷五月综合色啪小说| 熟妇人妻不卡中文字幕| 欧美 亚洲 国产 日韩一| 成人影院久久| 性少妇av在线| 1024视频免费在线观看| 色94色欧美一区二区| 天天躁狠狠躁夜夜躁狠狠躁| 国产无遮挡羞羞视频在线观看| 人妻人人澡人人爽人人| 在线观看国产h片| 日日撸夜夜添| 日韩av在线免费看完整版不卡| 久久精品人人爽人人爽视色| 秋霞伦理黄片| 中文字幕亚洲精品专区| 99久久99久久久精品蜜桃| av又黄又爽大尺度在线免费看| 男人操女人黄网站| 免费观看av网站的网址| 久久精品国产亚洲av涩爱| 在线观看免费视频网站a站| 天天躁夜夜躁狠狠躁躁| 亚洲欧美激情在线| 99久久99久久久精品蜜桃| 亚洲情色 制服丝袜| 伦理电影免费视频| 女人精品久久久久毛片| 久久性视频一级片| 免费观看人在逋| 九草在线视频观看| h视频一区二区三区| av线在线观看网站| 一二三四在线观看免费中文在| 免费日韩欧美在线观看| 久久国产亚洲av麻豆专区| √禁漫天堂资源中文www| 国产一区二区在线观看av| 在线观看免费日韩欧美大片| 亚洲一码二码三码区别大吗| 亚洲精品美女久久av网站| 美女午夜性视频免费| 韩国精品一区二区三区| 性色av一级| 午夜日韩欧美国产| 亚洲天堂av无毛| 久久免费观看电影| 国产午夜精品一二区理论片| 亚洲欧美激情在线| 欧美日韩亚洲国产一区二区在线观看 | 亚洲天堂av无毛| 久久久亚洲精品成人影院| 久久久国产精品麻豆| 99re6热这里在线精品视频| 日韩一区二区视频免费看| 精品亚洲成a人片在线观看| 中文字幕亚洲精品专区| 韩国av在线不卡| 老司机影院成人| 欧美 亚洲 国产 日韩一| 在线观看免费高清a一片| 久久综合国产亚洲精品| 黑人猛操日本美女一级片| 又粗又硬又长又爽又黄的视频| 99久久精品国产亚洲精品| 日本av免费视频播放| 亚洲精品日韩在线中文字幕| 亚洲精品aⅴ在线观看| 在线天堂最新版资源| 人人妻人人澡人人爽人人夜夜| 亚洲国产精品999| 日韩精品免费视频一区二区三区| 国产在线一区二区三区精| 在线观看www视频免费| 下体分泌物呈黄色| 午夜福利影视在线免费观看| 最近中文字幕高清免费大全6| 男女国产视频网站| 亚洲av男天堂| 亚洲国产av影院在线观看| 91aial.com中文字幕在线观看| 电影成人av| 51午夜福利影视在线观看| 99精国产麻豆久久婷婷| 国产男女超爽视频在线观看| 国产成人一区二区在线| 80岁老熟妇乱子伦牲交| 一区二区日韩欧美中文字幕| 美女午夜性视频免费| 天天影视国产精品| 人人妻,人人澡人人爽秒播 | 免费av中文字幕在线| 高清av免费在线| 自拍欧美九色日韩亚洲蝌蚪91| 日韩不卡一区二区三区视频在线| 一区在线观看完整版| 制服诱惑二区| 91aial.com中文字幕在线观看| 国产精品国产三级专区第一集| 中文天堂在线官网| 亚洲第一区二区三区不卡| 精品国产一区二区三区久久久樱花| 日本黄色日本黄色录像| 在线观看免费高清a一片| 欧美中文综合在线视频| 777米奇影视久久| 777米奇影视久久| 高清av免费在线| 国产成人啪精品午夜网站| 2018国产大陆天天弄谢| 亚洲国产中文字幕在线视频| 久久久久国产一级毛片高清牌| 亚洲av男天堂| 成人黄色视频免费在线看| 高清视频免费观看一区二区| 亚洲欧美激情在线| 欧美精品av麻豆av| 亚洲一区二区三区欧美精品| 久久精品熟女亚洲av麻豆精品| 老司机在亚洲福利影院| 免费av中文字幕在线| 欧美精品人与动牲交sv欧美| 国产色婷婷99| 免费日韩欧美在线观看| 777米奇影视久久| 男的添女的下面高潮视频| 国产探花极品一区二区| 美女高潮到喷水免费观看| 中文字幕另类日韩欧美亚洲嫩草| 免费观看av网站的网址| 亚洲国产欧美一区二区综合| 又大又爽又粗| 日本午夜av视频| bbb黄色大片| 18禁裸乳无遮挡动漫免费视频| 大码成人一级视频| 午夜激情久久久久久久| 久久午夜综合久久蜜桃| 亚洲国产欧美网| 高清欧美精品videossex| a 毛片基地| 国产日韩欧美视频二区| 欧美 日韩 精品 国产| 自线自在国产av| 亚洲成人免费av在线播放| 欧美97在线视频| 国产又爽黄色视频| 国产精品香港三级国产av潘金莲 | 精品人妻一区二区三区麻豆| 国产亚洲精品第一综合不卡| 少妇精品久久久久久久| 日韩一区二区视频免费看| a级片在线免费高清观看视频| 亚洲,欧美精品.| 美女高潮到喷水免费观看| 亚洲av电影在线进入| 男女之事视频高清在线观看 | 日本欧美视频一区| 成人午夜精彩视频在线观看| 国产不卡av网站在线观看| 一本一本久久a久久精品综合妖精| 亚洲精品美女久久久久99蜜臀 | 亚洲久久久国产精品| 黑丝袜美女国产一区| 国产成人精品久久二区二区91 | 国产男女内射视频| 欧美另类一区| 日本av手机在线免费观看| 午夜日韩欧美国产| 女性被躁到高潮视频| 97人妻天天添夜夜摸| 亚洲欧美精品综合一区二区三区| 麻豆乱淫一区二区| 亚洲伊人色综图| 国产xxxxx性猛交| 深夜精品福利| 欧美黑人精品巨大| 久久鲁丝午夜福利片| 国产熟女欧美一区二区| 母亲3免费完整高清在线观看| 一级,二级,三级黄色视频| 哪个播放器可以免费观看大片| 久久久精品国产亚洲av高清涩受| av又黄又爽大尺度在线免费看| 观看美女的网站| 91精品国产国语对白视频| 又黄又粗又硬又大视频| 亚洲欧美激情在线| 自线自在国产av| 国产精品香港三级国产av潘金莲 | 亚洲欧美成人精品一区二区| 久久女婷五月综合色啪小说| 中文天堂在线官网| 婷婷色综合大香蕉| 日韩大片免费观看网站| 国产高清国产精品国产三级| 亚洲成人手机| 在现免费观看毛片| 999精品在线视频| videos熟女内射| kizo精华| 色精品久久人妻99蜜桃| 免费在线观看黄色视频的| 最近最新中文字幕免费大全7| av一本久久久久| 卡戴珊不雅视频在线播放| 在现免费观看毛片| 久久精品久久久久久久性| 亚洲av成人不卡在线观看播放网 | 国产97色在线日韩免费| 欧美久久黑人一区二区| 精品久久久精品久久久| 老司机在亚洲福利影院| 亚洲av在线观看美女高潮| 免费看不卡的av| 女性生殖器流出的白浆| av在线观看视频网站免费| 色婷婷av一区二区三区视频| 国产视频首页在线观看| 午夜日韩欧美国产| 九色亚洲精品在线播放| 久久久久久久大尺度免费视频| 日本欧美视频一区| 久久久久精品人妻al黑| 欧美xxⅹ黑人| 亚洲天堂av无毛| 中文精品一卡2卡3卡4更新| 午夜免费鲁丝| 老司机亚洲免费影院| 999精品在线视频| 国产淫语在线视频| 伦理电影大哥的女人| 日韩av不卡免费在线播放| 看非洲黑人一级黄片| 国产亚洲午夜精品一区二区久久| 国产1区2区3区精品| 欧美日韩精品网址| 精品亚洲乱码少妇综合久久| 亚洲国产精品999| 在线看a的网站| 日韩精品有码人妻一区| 国产精品 欧美亚洲| 精品国产国语对白av| 国产精品一二三区在线看| 在线观看一区二区三区激情| 精品少妇黑人巨大在线播放| 亚洲一码二码三码区别大吗| 婷婷色综合大香蕉| 午夜久久久在线观看| 久久精品国产综合久久久| 亚洲精品在线美女| 日韩视频在线欧美| 黄色一级大片看看| 男女边摸边吃奶| 中文字幕色久视频| 国产99久久九九免费精品| 9热在线视频观看99| 日日啪夜夜爽| 在线精品无人区一区二区三| 国产成人精品无人区| 国产麻豆69| 国产精品女同一区二区软件| 又黄又粗又硬又大视频| 纯流量卡能插随身wifi吗| 久久天堂一区二区三区四区| 欧美日韩一区二区视频在线观看视频在线| 欧美精品一区二区免费开放| 一边亲一边摸免费视频| av国产久精品久网站免费入址| 久久久精品区二区三区| 久热这里只有精品99| 精品国产乱码久久久久久男人| 看非洲黑人一级黄片| 久久午夜综合久久蜜桃| 老司机影院成人| 国产乱人偷精品视频| 高清视频免费观看一区二区| 午夜福利网站1000一区二区三区| 精品一区二区三区四区五区乱码 | 在线观看免费日韩欧美大片| 美女扒开内裤让男人捅视频| 亚洲综合精品二区| 国产欧美日韩一区二区三区在线| 国产av一区二区精品久久| 男女之事视频高清在线观看 | 啦啦啦 在线观看视频| 最新的欧美精品一区二区| 一级片'在线观看视频| 嫩草影院入口| 久久久久久久大尺度免费视频| 人妻一区二区av| 伊人亚洲综合成人网| 中文字幕色久视频| av国产久精品久网站免费入址| 亚洲av综合色区一区| 热99久久久久精品小说推荐| 人成视频在线观看免费观看| 爱豆传媒免费全集在线观看| 午夜91福利影院| 国产欧美日韩一区二区三区在线| 在线精品无人区一区二区三| 18在线观看网站| 男女高潮啪啪啪动态图| www日本在线高清视频| 男女午夜视频在线观看| 成人免费观看视频高清| 亚洲美女黄色视频免费看| 欧美日韩视频高清一区二区三区二| 成年人午夜在线观看视频| 男女午夜视频在线观看| 亚洲av电影在线进入| 久久精品久久久久久噜噜老黄| 精品少妇久久久久久888优播| 国产精品国产三级国产专区5o| 秋霞伦理黄片| av在线老鸭窝| 国产不卡av网站在线观看| 1024视频免费在线观看| 日本猛色少妇xxxxx猛交久久| 婷婷色综合www| 亚洲国产精品国产精品| 久久女婷五月综合色啪小说| 少妇人妻 视频| 久久精品亚洲熟妇少妇任你| 久久av网站| 十分钟在线观看高清视频www| 欧美乱码精品一区二区三区| 久久国产精品男人的天堂亚洲| 伦理电影大哥的女人| 日本午夜av视频| 黄片播放在线免费| 日日爽夜夜爽网站| 另类精品久久| 成年av动漫网址| 亚洲精品美女久久av网站| 汤姆久久久久久久影院中文字幕| 国产日韩一区二区三区精品不卡| 建设人人有责人人尽责人人享有的| 亚洲成人国产一区在线观看 | 老鸭窝网址在线观看| 美女扒开内裤让男人捅视频| 国产片特级美女逼逼视频| 日韩精品免费视频一区二区三区| 久久99一区二区三区| 一边亲一边摸免费视频| 国产片特级美女逼逼视频| avwww免费| 男女下面插进去视频免费观看| 国产男人的电影天堂91| www日本在线高清视频| 国产色婷婷99| 在线看a的网站| 久久久国产一区二区| 777久久人妻少妇嫩草av网站| 巨乳人妻的诱惑在线观看| 你懂的网址亚洲精品在线观看| 性少妇av在线| 美国免费a级毛片| 日本欧美视频一区| 久久久久久久久久久久大奶| 多毛熟女@视频| 久久人人爽人人片av| 午夜免费男女啪啪视频观看| 亚洲精品自拍成人| 黑人猛操日本美女一级片| 亚洲五月色婷婷综合| 午夜福利,免费看| 亚洲天堂av无毛| 高清不卡的av网站| 国产亚洲精品第一综合不卡| 国产免费福利视频在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 久久久久精品性色| av视频免费观看在线观看| 国产精品国产三级专区第一集| 国产精品女同一区二区软件| 婷婷色av中文字幕| 黄色一级大片看看| 亚洲国产av影院在线观看| 亚洲精品国产av成人精品| 午夜影院在线不卡| 亚洲激情五月婷婷啪啪| 欧美日韩视频高清一区二区三区二| 亚洲国产成人一精品久久久| 免费高清在线观看日韩| 国产一区亚洲一区在线观看| 午夜激情av网站| 在线精品无人区一区二区三| e午夜精品久久久久久久| 免费日韩欧美在线观看| 亚洲精品久久成人aⅴ小说| 亚洲精品久久久久久婷婷小说| 亚洲精品久久午夜乱码| 丁香六月天网| 欧美在线一区亚洲| 国产精品一区二区在线不卡| 熟妇人妻不卡中文字幕| 亚洲第一区二区三区不卡| 国产片内射在线| 欧美亚洲日本最大视频资源| 90打野战视频偷拍视频| 国产亚洲欧美精品永久| 一边摸一边做爽爽视频免费| 国产在线视频一区二区| 国产成人91sexporn| 黄色视频在线播放观看不卡| 亚洲人成网站在线观看播放| 如日韩欧美国产精品一区二区三区| 丁香六月欧美| 丝袜喷水一区| 满18在线观看网站| 9色porny在线观看| 国产乱人偷精品视频| 天天操日日干夜夜撸| 亚洲精品久久午夜乱码| 精品亚洲成a人片在线观看| 国产精品久久久久久人妻精品电影 | 久久久国产精品麻豆| 国产精品人妻久久久影院| 少妇人妻精品综合一区二区| 日韩,欧美,国产一区二区三区| 久久久久视频综合| 久久久久国产精品人妻一区二区| 捣出白浆h1v1| 免费不卡黄色视频| 国产日韩欧美视频二区| 久久国产精品大桥未久av| 大码成人一级视频| 麻豆精品久久久久久蜜桃| 大陆偷拍与自拍| 亚洲欧洲精品一区二区精品久久久 | 青春草国产在线视频| 亚洲av电影在线进入| 亚洲国产中文字幕在线视频| 国产成人欧美| 久久人妻熟女aⅴ| 久久久久久久久免费视频了| 丁香六月天网| 汤姆久久久久久久影院中文字幕| 大片免费播放器 马上看| 秋霞伦理黄片| 日韩伦理黄色片| 中文乱码字字幕精品一区二区三区| 欧美激情 高清一区二区三区| 成人国产麻豆网| 中文乱码字字幕精品一区二区三区| 国产97色在线日韩免费| 一区二区三区四区激情视频| 午夜福利视频在线观看免费| 日韩av在线免费看完整版不卡| 啦啦啦视频在线资源免费观看| 久久久久久久大尺度免费视频| 18禁国产床啪视频网站| 免费在线观看黄色视频的| 欧美 亚洲 国产 日韩一| av视频免费观看在线观看| 免费观看人在逋| 美女扒开内裤让男人捅视频| 国产精品一区二区精品视频观看| 国产精品久久久av美女十八| 国产精品偷伦视频观看了| 欧美人与性动交α欧美软件| 亚洲成av片中文字幕在线观看| 日本vs欧美在线观看视频| 韩国av在线不卡| 色网站视频免费| 9色porny在线观看| 欧美精品高潮呻吟av久久| 一区在线观看完整版| 国产片特级美女逼逼视频| 国产一区二区三区综合在线观看| xxx大片免费视频| 三上悠亚av全集在线观看| 欧美精品人与动牲交sv欧美| netflix在线观看网站| 1024香蕉在线观看| 成人国产av品久久久| 久久影院123| 在线观看人妻少妇| 亚洲天堂av无毛| 老司机在亚洲福利影院| 久久女婷五月综合色啪小说| 成人18禁高潮啪啪吃奶动态图| 国产一区有黄有色的免费视频| 两性夫妻黄色片| 五月天丁香电影| 性高湖久久久久久久久免费观看| 欧美少妇被猛烈插入视频| 精品人妻一区二区三区麻豆| 91国产中文字幕| 亚洲精品日本国产第一区| 国产精品.久久久| av天堂久久9| 国产精品久久久av美女十八| 欧美人与性动交α欧美精品济南到| 在线 av 中文字幕| 成人免费观看视频高清| 制服人妻中文乱码| 巨乳人妻的诱惑在线观看| 日韩伦理黄色片| 久久久精品区二区三区| 亚洲一码二码三码区别大吗| 男男h啪啪无遮挡| 人人妻人人爽人人添夜夜欢视频| 大片电影免费在线观看免费| 伊人久久国产一区二区| 日本色播在线视频| 一边摸一边做爽爽视频免费| av.在线天堂| 美女福利国产在线| 亚洲男人天堂网一区| 国产精品欧美亚洲77777| 国产在线一区二区三区精| 哪个播放器可以免费观看大片| 国产精品久久久久久精品电影小说| videosex国产| 色视频在线一区二区三区| 久久毛片免费看一区二区三区| 日韩中文字幕视频在线看片| 国产深夜福利视频在线观看| 别揉我奶头~嗯~啊~动态视频 | 天天躁夜夜躁狠狠久久av| 国产精品久久久久成人av| 在线天堂中文资源库| av天堂久久9| 久久精品aⅴ一区二区三区四区| 极品人妻少妇av视频| 欧美日韩亚洲高清精品| 久久性视频一级片| 少妇 在线观看| 满18在线观看网站| 亚洲四区av| 精品国产乱码久久久久久男人| 久久国产亚洲av麻豆专区| 高清视频免费观看一区二区| 久久午夜综合久久蜜桃| 少妇被粗大的猛进出69影院| 国产男女内射视频| 国产精品久久久av美女十八| 99久久精品国产亚洲精品| 国产深夜福利视频在线观看| 美女主播在线视频| 在线亚洲精品国产二区图片欧美| 久久精品亚洲av国产电影网| 久热爱精品视频在线9| 国产成人精品久久二区二区91 | 中文字幕亚洲精品专区| 91aial.com中文字幕在线观看| 黄网站色视频无遮挡免费观看| a级毛片在线看网站| 欧美国产精品va在线观看不卡| 黑丝袜美女国产一区| 毛片一级片免费看久久久久| 精品人妻熟女毛片av久久网站| 国产人伦9x9x在线观看| 久久精品熟女亚洲av麻豆精品| 精品久久蜜臀av无| 在线观看人妻少妇| 国产成人91sexporn| 丝瓜视频免费看黄片| 日本色播在线视频| 亚洲成人一二三区av| 91国产中文字幕| 亚洲精品国产区一区二| 久久99精品国语久久久| 少妇人妻精品综合一区二区| a级毛片在线看网站| 777久久人妻少妇嫩草av网站| 国产福利在线免费观看视频| 日韩欧美精品免费久久| 亚洲成色77777| netflix在线观看网站| 午夜福利在线免费观看网站| 午夜日韩欧美国产| 国产精品一区二区在线观看99| 久久精品国产综合久久久| 成人亚洲精品一区在线观看| 在线观看免费高清a一片| 性高湖久久久久久久久免费观看| 丁香六月欧美| 亚洲av日韩精品久久久久久密 | 99热国产这里只有精品6| 久久久久久久国产电影| 又大又黄又爽视频免费| 国产午夜精品一二区理论片| 综合色丁香网| 欧美黑人欧美精品刺激| 天堂8中文在线网| 成年人免费黄色播放视频| 一区二区三区四区激情视频| 国产欧美亚洲国产| 黄色毛片三级朝国网站| 女人久久www免费人成看片| 美女福利国产在线| 欧美最新免费一区二区三区| 国产免费一区二区三区四区乱码| 亚洲人成77777在线视频| 女人爽到高潮嗷嗷叫在线视频| 叶爱在线成人免费视频播放| 丰满少妇做爰视频| 免费观看av网站的网址| 99久国产av精品国产电影| 在线看a的网站| 欧美日韩精品网址| 亚洲精品自拍成人| 国产日韩欧美视频二区| 91国产中文字幕| 欧美亚洲日本最大视频资源| 哪个播放器可以免费观看大片| 久久久久久久国产电影| 国产一区二区在线观看av| 国产亚洲精品第一综合不卡|