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

    1T2R并聯(lián)機構(gòu)拓?fù)浣雕钤O(shè)計與運動性能分析

    2021-06-09 10:24:12沈惠平孟慶梅
    農(nóng)業(yè)機械學(xué)報 2021年5期
    關(guān)鍵詞:支鏈耦合度并聯(lián)

    沈惠平 王 達 李 菊 孟慶梅

    (常州大學(xué)現(xiàn)代機構(gòu)學(xué)研究中心,常州 213016)

    0 引言

    三自由度的一平移兩轉(zhuǎn)動(1T2R)并聯(lián)機構(gòu)具有制造方便、靈活性好、精度易保證等優(yōu)點,目前仍是國內(nèi)外學(xué)者的研究熱點之一。SONG等[1]運用基于共形幾何代數(shù)的型綜合方法,綜合了一類一平移兩轉(zhuǎn)動(1T2R)并聯(lián)機構(gòu);汪滿新等[2]研究了基于UP和UPR-SPR型等效運動的1T2R并聯(lián)機構(gòu)的拓?fù)渚C合方法;房立豐等[3]以一平移兩轉(zhuǎn)動并聯(lián)穩(wěn)定平臺為例,研究了少自由度并聯(lián)穩(wěn)定平臺拓?fù)浣Y(jié)構(gòu)設(shè)計方法;ENRIQUE等[4]對3-CUP的1T2R并聯(lián)機構(gòu)進行了運動學(xué)分析;HUSSEIN等[5]對3-CRC的1T2R并聯(lián)機構(gòu)進行了運動學(xué)分析;SUNG等[6]對具有封閉型正解的1T2R并聯(lián)機構(gòu)進行了運動學(xué)分析;SUN等[7]研究了含有限運動與寄生運動的并聯(lián)機構(gòu)拓?fù)渚C合問題;CHEN等[8]研究了無寄生運動的3-PRRU的1T2R并聯(lián)機器人的動力學(xué)建模和性能;GAN等[9]提出基于可重構(gòu)原理可在一平移兩轉(zhuǎn)動(1T2R)運動與純旋轉(zhuǎn)(3R)兩種工作運動模式之間實現(xiàn)可重構(gòu)性;ALI等[10]利用高速艇的運動模擬器(HSB),提出了一種新型的1T2R并聯(lián)機構(gòu);SAIOA等[11-12]分析了2-PRU+1-PRS的1T2R并聯(lián)機構(gòu)的運動學(xué)、奇異性和動力學(xué)問題;JOSé等[13]對具有實時力/位置控制的1T2R型康復(fù)醫(yī)療并聯(lián)機器人進行了模擬實驗;XIE等[14]提出基于1T2R并聯(lián)機構(gòu)的兩種新型鉸接刀頭,以用作模塊化加工單元;車林仙等[15]采用智能算法求解了一平移兩轉(zhuǎn)動2-PUR+1-PSR并聯(lián)機構(gòu)的尺度參數(shù)優(yōu)化設(shè)計問題;CHONG等[16]研究了3-RCU的1T2R并聯(lián)機構(gòu)的運動學(xué)優(yōu)化問題。LIU等[17]研究了具有解析解的1T2R并聯(lián)機構(gòu)的類型綜合。

    求解并聯(lián)機構(gòu)的位置正解一般采用數(shù)值法和解析法,但較少考慮機構(gòu)拓?fù)涮匦耘c運動學(xué)性能之間的關(guān)系。沈惠平等[18]提出按耦合度κ分類求解并聯(lián)機構(gòu)位置正解的方法,對κ為0的并聯(lián)機構(gòu)直接求解其符號式位置正解,對κ不為0的并聯(lián)機構(gòu)可用數(shù)值法求解其位置正解。文獻[19-23]提出低耦合度(或零耦合度)且具有運動解耦性的少自由度并聯(lián)機構(gòu),并求解了其符號式位置正反解,分析了機構(gòu)運動學(xué)性能。

    本文根據(jù)基于方位特征(Position and orientation characteristic,POC)方程的并聯(lián)機構(gòu)拓?fù)湓O(shè)計理論[24],設(shè)計一種一平移兩轉(zhuǎn)動并聯(lián)機構(gòu),因該機構(gòu)耦合度為1,得不到符號式位置正解,故對其進行降耦優(yōu)化設(shè)計[25],得到一種零耦合度的一平移兩轉(zhuǎn)動并聯(lián)機構(gòu),并對其符號式位置正反解、奇異位形及其工作空間進行計算和分析。

    1 機構(gòu)設(shè)計及拓?fù)浞治?/h2>

    1.1 機構(gòu)設(shè)計

    1.1.1POC集計算

    機構(gòu)POC方程[24]為

    (1)

    (2)

    式中MJi——第i個運動副的POC集

    Mbi——第i條支鏈末端的POC集

    Mpa——機構(gòu)動平臺的POC集

    1.1.2支鏈設(shè)計

    (1)支鏈POC集確定

    支鏈末端滿足Mpa的POC集Mbi可取為

    (2)支鏈數(shù)目確定

    一般情況下,支鏈數(shù)目應(yīng)該等于并聯(lián)機構(gòu)自由度數(shù),故設(shè)計該機構(gòu)支鏈數(shù)目為3。

    (3)支鏈拓?fù)浣Y(jié)構(gòu)綜合

    每條支鏈可有1個驅(qū)動副,為實現(xiàn)動平臺長距離的移動,故設(shè)計每條支鏈的驅(qū)動副為P副,并在平行的導(dǎo)軌上。

    第Ⅰ條支鏈可設(shè)計為方位特征支鏈[25],即其末端構(gòu)件的POC集為動平臺的Mpa,即在滿足第Ⅰ條支鏈能實現(xiàn)一平移兩轉(zhuǎn)動的同時,盡可能使其拓?fù)浣Y(jié)構(gòu)簡單。由式(2)可知,在P副的基礎(chǔ)上串聯(lián)2個軸線垂直的R副,故第Ⅰ條支鏈的拓?fù)浣Y(jié)構(gòu)可設(shè)計為

    SOCⅠ:{P⊥R⊥R}

    SOCi:{P-S-S}(i=Ⅱ,Ⅲ)

    1.1.3機構(gòu)構(gòu)成

    根據(jù)支鏈的幾何布置規(guī)則,設(shè)計的一平移兩轉(zhuǎn)動(1T2R)并聯(lián)機構(gòu)如圖1所示,該機構(gòu)由動平臺1、靜平臺0,以及3條支鏈(Ⅰ、Ⅱ、Ⅲ)構(gòu)成。Ⅰ支鏈為方位特征支鏈P31⊥R32⊥R33,P31軸線與P11軸線平行,R32軸線始終與靜平臺0平面平行;Ⅱ、Ⅲ支鏈均為Pi1-Si2-Si3(i=1,2)無約束支鏈,其移動副P11和P21軸線共線。

    1.2 機構(gòu)拓?fù)浞治?/h3>

    1.2.1機構(gòu)POC集驗證

    (1)支鏈拓?fù)浣Y(jié)構(gòu)分別為

    SOCⅠ:{P31⊥R32⊥R33}

    SOCⅡ:{P11-S12-S13}

    SOCⅢ:{P21-S22-S23}

    (2)選取動平臺1上任一點為基點O′。

    (3)確定支鏈末端構(gòu)件的POC集。

    由式(1)、(2)得支鏈末端構(gòu)件POC集分別為

    (4)確定動平臺的POC集,由式(2)得

    Mpa=MbⅠ∩MbⅡ∩MbⅢ=

    因此,動平臺1具有沿P31軸線一維平移以及繞轉(zhuǎn)動副R32和R33軸線兩維轉(zhuǎn)動的輸出特性。

    1.2.2機構(gòu)自由度

    并聯(lián)機構(gòu)全周DOF公式[24]為

    (3)

    (4)

    v=m-n+1

    式中F——機構(gòu)自由度

    fi——第i個運動副的自由度

    m——運動副數(shù)

    v——獨立回路數(shù)n——構(gòu)件數(shù)

    ξLj——第j個獨立回路的獨立位移方程數(shù)

    Mb(j+1)——前j+1條支鏈末端構(gòu)件的POC集

    由第Ⅰ、Ⅱ支鏈構(gòu)成的第1回路(即子并聯(lián)機構(gòu))的獨立位移方程數(shù)為

    該子并聯(lián)機構(gòu)自由度為

    該子并聯(lián)機構(gòu)的POC集為

    Mpa(Ⅰ-Ⅱ)=MbⅠ∩MbⅡ=

    由上述子并聯(lián)機構(gòu)和支鏈Ⅲ構(gòu)成第2回路的獨立位移方程數(shù)為

    機構(gòu)自由度為

    故當(dāng)取靜平臺0的3個移動副P11、P21、P31為驅(qū)動副時,動平臺1可實現(xiàn)沿P31軸線的一維平移以及繞轉(zhuǎn)動副R32和R33軸線的兩維轉(zhuǎn)動。

    1.2.3機構(gòu)耦合度

    由基于序單開鏈[24](Single-open-chain,SOC)的機構(gòu)組成原理可知,任一機構(gòu)可分解為約束度為正、零、負(fù)的3種有序單開鏈(SOC),第j個SOCj的約束度定義為

    (5)

    式中mj——第j個SOCj的運動副數(shù)

    Ij——第j個SOCj的驅(qū)動副數(shù)

    進一步,一組有序的v個SOC可組成1個獨立回路數(shù)為v的(最小)子運動鏈[24](Sub-kinematics chain,SKC),SKC耦合度為

    (6)

    在1.2.2節(jié),已計算出2個回路的獨立位移方程數(shù),即ξL1=6,ξL2=6。因此,由式(5)可得,它們的約束度分別為

    顯然,由第1、2回路構(gòu)成1個SKC,由式(6)得,該SKC耦合度為

    即該機構(gòu)只包含1個SKC,其耦合度κ=1,表明在位置正解求解時,需要在約束度為正值(Δ1=1)的回路上設(shè)置1個虛擬變量,而在約束度為負(fù)值(Δ2=-1)的回路上,建立1個含該虛擬變量的約束方程,可用一維搜索法求其數(shù)值型位置正解;但數(shù)值型位置正解不利于后續(xù)的誤差分析、軌跡規(guī)劃和運動控制,以及動力學(xué)分析。

    因此,可對該機構(gòu)進行拓?fù)浣雕钤O(shè)計,以使該機構(gòu)在保持基本功能(DOF、POC)不變的情況下,具有符號式位置正解。

    1.3 機構(gòu)降耦設(shè)計及其拓?fù)浞治?/h3>

    根據(jù)機構(gòu)拓?fù)浣雕钤O(shè)計方法[25]中的“轉(zhuǎn)動副/球副重合法”,將圖1所示機構(gòu)中支鏈Ⅱ、Ⅲ中的球副S13、S23合并,但考慮雙重球副制造的不方便性,可將S13移至支鏈Ⅲ上,得到的改進機構(gòu)如圖2所示,即可保持機構(gòu)基本功能(DOF、POC)不變,而使其耦合度從1降到0。

    此時,支鏈Ⅱ、Ⅲ組合成為一復(fù)雜支鏈(子并聯(lián)機構(gòu)),可知,其仍為無約束支鏈;同時,方位特征支鏈Ⅰ不變。因此,根據(jù)1.2節(jié)中所述的拓?fù)浞治霾襟E知,機構(gòu)降耦后的POC集、自由度并未發(fā)生改變,即:

    第1回路{P11-S12-S13-S22-P21},顯然ξL1=6,其約束度為

    式中,-3是指應(yīng)扣除繞S12S13、S13S22及S12S22連線的3個局部轉(zhuǎn)動自由度;但繞S12S22連線的局部轉(zhuǎn)動自由度,應(yīng)計入第2回路的約束度計算中。

    第2回路{RS12S22-S23-R33-R32-P31},顯然,ξL2=6,其約束度為

    可知,根據(jù)子運動鏈(SKC)的劃分原則[24],第1、2回路分別構(gòu)成SKC1、SKC2,它們耦合度均為零,即κ1=κ2=0,因此,該并聯(lián)機構(gòu)耦合度為零,其位置正解易求出。

    2 機構(gòu)位置分析

    2.1 基于拓?fù)涮卣鞯臋C構(gòu)位置正解求解原理

    由基于有序單開鏈(SOC)的機構(gòu)組成原理[24]可知,降耦后的機構(gòu)包含2個SKC,每個SKC可分解為一系列約束度分別為正值、零、負(fù)值的單開鏈回路,因此,該機構(gòu)位置正解的求解,可轉(zhuǎn)換為這兩個SKC內(nèi)所含單開鏈的位置求解。對本機構(gòu)而言,SKC內(nèi)僅有約束度為零的單開鏈,其運動位置具有確定性,因此,其位置正解能獨立求解。

    2.2 坐標(biāo)系建立及參數(shù)標(biāo)注

    機構(gòu)運動學(xué)建模如圖3所示,設(shè)靜平臺0為長方形,其寬lOI為a,長為任意正實數(shù),動平臺1上lGF=b。設(shè)靜平臺0上的點O為靜坐標(biāo)系OXYZ原點,X軸沿OC方向,Y軸沿OI方向,Z軸由右手螺旋法則確定;動平臺1上點O′為動坐標(biāo)系O′X′Y′Z′原點,其X′軸沿O′M方向,Y′軸沿FO′方向,Z′軸由右手螺旋法則確定。

    點O′與點G重合,AB、CD、GH都垂直于靜平臺0平面,△EBD的中線EJ與靜平臺0平面夾角為δ。設(shè)動平臺1繞OX、OY軸正方向轉(zhuǎn)動的角度分別為α、β;設(shè)lAB=lCD=d,lBE=lDE=e,lEF=f,lGH=h,并令lFD/lED=(f+e)/e=k0;設(shè)3個驅(qū)動副輸入分別為lOA=l1、lOC=l2、lIH=l3。

    2.3 位置正解求解

    已知:3個移動副的位置分別為l1、l2、l3,求:動平臺上O′=(x,y,z)及轉(zhuǎn)角α和β。

    2.3.1SKC1位置求解

    第1回路A-B-E-D-C中,各點坐標(biāo)為A=(l1,0,0),B=(l1,0,d),C=(l2,0,0),D=(l2,0,d)。由lEB=lED的幾何關(guān)系,易求得點E坐標(biāo)為

    2.3.2SKC2位置求解

    第2回路F-G-H中,各點的坐標(biāo)為G=O′=(l3,a,h),H=(l3,a,0),由lFD/lED=k0,得

    (7)

    由桿長約束條件lFG=b,整理并解得

    {l3-[k0l1+(2-k0)l2]/2}2

    因動坐標(biāo)系中點F在靜坐標(biāo)系OXYZ中的坐標(biāo)為

    其中

    式中Q——動坐標(biāo)系到靜坐標(biāo)系的變換矩陣

    解得

    (8)

    由式(7)、(8)求出動平臺姿態(tài)角α、β、O′(x,y,z)及動平臺上任一點的位置。

    2.4 位置反解求解

    機構(gòu)位置逆解為:已知動平臺1的轉(zhuǎn)角α、β及O′(x,y,z),求l1、l2、l3。

    由式(8)可得

    (9)

    從而得點E坐標(biāo)為

    由lBE=lDE=e,解得

    (10)

    可知,l1、l2各有2組解,故逆解數(shù)為2×2=4,故該機構(gòu)有4種構(gòu)型。

    2.5 正逆解驗算

    設(shè)該機構(gòu)結(jié)構(gòu)參數(shù)分別為:a=60 mm,b=20 mm,d=5 mm,e=75 mm,f=5 mm,h=45 mm。設(shè)此時3個移動副的位置分別為:l1=7.380 3 mm,l2=116.991 6 mm,l3=60.000 0 mm。

    用Matlab求得機構(gòu)位置正解,如表1所示。

    表1 位置正解

    將表1中序號1的數(shù)值代入式(10),用Matlab解得機構(gòu)的所有反解,如表2所示。

    表2 位置反解

    可見,表2中序號2組的逆解數(shù)值,與正解計算給定的3個輸入值一致。因此,正、反解公式推導(dǎo)正確。

    3 機構(gòu)奇異性分析

    3.1 雅可比矩陣

    采用Jacobian法分析該機構(gòu)的奇異位形。將式(10)兩邊同時對時間t求導(dǎo),得

    (11)

    該機構(gòu)動平臺輸出速度v1和輸入速度v2的關(guān)系為

    Jpv1=Jqv2

    (12)

    f11=-2bsinβcosα(XE-l1)/k0+2bsinαYE/k0-

    2bcosβcosα(ZE-d)/k0

    f12=-2bcosβsinα(XE-l1)/k0+

    2bsinβsinα(ZE-d)/k0

    f13=2(XE-l1)/k0

    f21=-2bsinβcosα(XE-l2)/k0+2bsinαYE/k0-

    2bcosβcosα(ZE-d)/k0

    f22=-2bcosβsinα(XE-l2)/k0+

    2bsinβsinα(ZE-d)/k0

    f23=2(XE-l2)/k0f31=0f32=0f33=1

    u11=2(XE-l1)u12=(k0-1)/k0

    u22=-2(XE-l2)/k0u33=1

    依據(jù)Jp、Jq矩陣是否奇異,將機構(gòu)的奇異位形分為如下3類:①當(dāng)det(Jq)=0時,機構(gòu)發(fā)生輸入奇異。②當(dāng)det(Jp)=0時,機構(gòu)發(fā)生輸出奇異。③當(dāng)det(Jq)=det(Jp)=0時,機構(gòu)發(fā)生綜合奇異。

    3.2 奇異位形分析

    當(dāng)機構(gòu)發(fā)生輸入奇異,意味著每條支鏈靠近驅(qū)動桿的兩根桿處于折疊或完全展開狀態(tài)。此時,動平臺的自由度數(shù)減少,則有det(Jq)=0,方程解為

    l1=l2

    (13)

    根據(jù)式(13)可知,滑塊P11和P21重合時才會發(fā)生輸入奇異,但考慮到構(gòu)件間的干涉,因此,此類奇異位形不存在。

    當(dāng)機構(gòu)發(fā)生輸出奇異,意味著每條支鏈靠近動平臺的桿件處于折疊或完全展開的狀態(tài),此時動平臺自由度數(shù)增多,即使鎖住輸入,動平臺也可能存在自由度輸出。設(shè)

    若det(Jp)=0,則向量e1、e2、e3有如下2種情況:

    (1)存在2個向量線性相關(guān)

    若e1=ke2,滿足GF∥△EBD的條件下,機構(gòu)存在奇異位置,如圖4所示。

    若e1=ke3,滿足GF∥EB的條件下,機構(gòu)存在奇異位置,如圖5所示。

    若e2=ke3,滿足點G、F、D共線的條件下,機構(gòu)存在奇異位置,如圖6所示。

    (2)3個向量線性相關(guān)

    設(shè)e3=k1e1+k2e2(k1k2≠0),滿足sinα=0(即FG⊥GH)的條件下,機構(gòu)存在奇異位置,如圖7所示。

    當(dāng)det(Jq)=det(Jp)=0,即輸入奇異和輸出奇異同時發(fā)生,該機構(gòu)不存在綜合奇異位置。

    4 基于符號式位置正解的工作空間分析

    4.1 工作空間計算

    并聯(lián)機構(gòu)的可達工作空間是指在考慮運動副的轉(zhuǎn)角范圍、桿長不干涉的情況下,末端執(zhí)行器的工作區(qū)域,是衡量并聯(lián)機構(gòu)性能的一個重要指標(biāo)。

    傳統(tǒng)的工作空間計算是基于機構(gòu)位置反解公式,采用離散化空間的三維搜索法,即需要預(yù)先估計設(shè)定一個搜索范圍,通過Matlab軟件搜索所有滿足位置反解約束條件的點,由這些點組成的三維圖即為該機構(gòu)的工作空間。但由于預(yù)設(shè)的搜索范圍難于估計其大小,因此,計算量很大。

    本文降耦優(yōu)化機構(gòu)具有符號式位置正解,因此,直接采用位置正解來計算工作空間,具有思路清晰、計算量少、工作空間邊界計算準(zhǔn)確等優(yōu)點。

    當(dāng)滑塊P31固定時(取l3=60 mm),滑塊P11、P21的活動范圍為:-20 mm≤l1≤60 mm,60 mm≤l2≤140 mm。通過Matlab軟件編程,由式(7)、(8)得到該機構(gòu)兩維轉(zhuǎn)動(2R)的工作空間,如圖8所示。

    2R輸出時,工作空間在XOY、XOZ、YOZ平面上的投影圖如圖9所示。

    當(dāng)3個滑塊都運動并具有不同速度時,滑塊P11、P21、P31活動范圍為:-20 mm≤l1≤60 mm,60 mm≤l2≤140 mm、60 mm≤l3≤120 mm,通過Matlab軟件編程,由式(7)、(8)得到該機構(gòu)一平移兩轉(zhuǎn)動(1T2R)的工作空間,如圖10所示(可看作是圖8所示的2R工作空間在X軸方向疊加而成)。

    1T2R工作空間在XOY、XOZ、YOZ平面上的投影圖如圖11所示。

    4.2 工作空間分析

    采用傳統(tǒng)的、基于式(9)、(10)的方法,同樣計算出了動平臺上點F如圖8、9所示兩維轉(zhuǎn)動工作空間,以及如圖10、11所示的一平移兩轉(zhuǎn)動工作空間,但基于符號式位置正解和位置反解求解2R、1T2R工作空間的計算量不同(由程序耗時,可直觀且真實地知道計算量);基于符號式位置正解、位置反解的程序,求得2R和1T2R工作空間的耗時,分別為7.856、129.307 s和 24.913、410.075 s,即基于符號式位置正解計算法,分別是基于位置反解的計算量的31.533%、31.532%。

    由圖8、10可知,用兩種工作空間計算方法求得該機構(gòu)的工作空間完全一致。由圖10、11可知,該機構(gòu)工作空間關(guān)于z=45 mm處的XOY面部分對稱;當(dāng)z>45 mm時,工作空間隨z增大而減小,當(dāng)z≤45 mm時,工作空間隨z減小而先增后減,其形狀較為規(guī)則,且工作空間大?;诜柺轿恢谜獾墓ぷ骺臻g計算,計算極為方便,計算量約為基于位置反解工作空間計算量的31.53%,且工作空間邊界計算準(zhǔn)確。

    5 結(jié)論

    (1)在保持初始設(shè)計1T2R并聯(lián)機構(gòu)的POC、DOF不變的情況下,通過拓?fù)浣雕钤O(shè)計進行優(yōu)化,使其耦合度從1降為0,進而可以直接推導(dǎo)出該機構(gòu)的符號式位置正解,而不需復(fù)雜的數(shù)學(xué)消元等推導(dǎo)過程,故有利于后續(xù)的尺度綜合、誤差分析、軌跡規(guī)劃、運動控制及動力學(xué)分析。

    (2)基于符號式位置正解的機構(gòu)工作空間計算法具有無需預(yù)估工作空間范圍、計算量少、工作空間邊界計算準(zhǔn)確等優(yōu)點,同時表明,該1T2R并聯(lián)機構(gòu)的工作空間大,形狀較為規(guī)則。

    猜你喜歡
    支鏈耦合度并聯(lián)
    中國北方蒸散-降水耦合度時空變化與水熱因子的關(guān)系
    干旱氣象(2022年5期)2022-11-16 04:40:24
    雙速感應(yīng)電機繞組耦合度研究
    防爆電機(2022年1期)2022-02-16 01:13:58
    識別串、并聯(lián)電路的方法
    遼寧省經(jīng)濟與生態(tài)環(huán)境耦合協(xié)調(diào)性分析
    審批由“串聯(lián)”改“并聯(lián)”好在哪里?
    基于耦合度分析的家禽孵化過程模糊解耦控制系統(tǒng)
    并聯(lián)型APF中SVPWM的零矢量分配
    一種軟開關(guān)的交錯并聯(lián)Buck/Boost雙向DC/DC變換器
    臭氧護理皮支鏈皮瓣200例觀察分析
    卵內(nèi)注射支鏈氨基酸對雞胚胎生長發(fā)育和孵化時間的影響
    飼料博覽(2015年4期)2015-04-05 10:34:14
    极品少妇高潮喷水抽搐| 国产精品亚洲av一区麻豆 | 国产成人精品在线电影| 国产亚洲精品第一综合不卡| 自线自在国产av| 国产片特级美女逼逼视频| 久久国产亚洲av麻豆专区| 一级毛片电影观看| 女人精品久久久久毛片| 一级毛片 在线播放| 中文字幕人妻熟女乱码| a级毛片在线看网站| 国产有黄有色有爽视频| 黄频高清免费视频| 中文字幕另类日韩欧美亚洲嫩草| 黄色一级大片看看| 欧美最新免费一区二区三区| 欧美在线一区亚洲| 中文字幕另类日韩欧美亚洲嫩草| 成人亚洲精品一区在线观看| 五月开心婷婷网| 成年人免费黄色播放视频| 天天躁夜夜躁狠狠躁躁| 性色av一级| 一边摸一边做爽爽视频免费| 午夜福利影视在线免费观看| 青春草亚洲视频在线观看| 丝袜脚勾引网站| 欧美精品人与动牲交sv欧美| 亚洲av日韩精品久久久久久密 | 日韩免费高清中文字幕av| 亚洲四区av| 一本大道久久a久久精品| 女人爽到高潮嗷嗷叫在线视频| 大片免费播放器 马上看| netflix在线观看网站| 青春草亚洲视频在线观看| 国产精品女同一区二区软件| 久久久精品94久久精品| 欧美国产精品一级二级三级| av片东京热男人的天堂| 中文字幕色久视频| 亚洲国产欧美在线一区| 色婷婷av一区二区三区视频| 啦啦啦在线观看免费高清www| 成人三级做爰电影| 国产男女内射视频| 久久ye,这里只有精品| 人人妻人人添人人爽欧美一区卜| 日日摸夜夜添夜夜爱| 久久综合国产亚洲精品| 丰满迷人的少妇在线观看| av国产精品久久久久影院| 王馨瑶露胸无遮挡在线观看| 婷婷色综合www| 五月天丁香电影| 免费黄色在线免费观看| 老司机靠b影院| 亚洲精品国产区一区二| 女的被弄到高潮叫床怎么办| 18禁观看日本| 国产亚洲精品第一综合不卡| 欧美97在线视频| 国产不卡av网站在线观看| 人人妻人人爽人人添夜夜欢视频| 亚洲久久久国产精品| 国产成人精品在线电影| 香蕉丝袜av| 伊人亚洲综合成人网| 色精品久久人妻99蜜桃| 久久精品亚洲av国产电影网| av在线观看视频网站免费| 老鸭窝网址在线观看| 中文字幕人妻丝袜制服| 永久免费av网站大全| 成年美女黄网站色视频大全免费| 国产熟女欧美一区二区| 国产成人精品久久二区二区91 | 18禁观看日本| 搡老岳熟女国产| 国产又色又爽无遮挡免| 亚洲av日韩精品久久久久久密 | 欧美日韩精品网址| 国语对白做爰xxxⅹ性视频网站| 9热在线视频观看99| 丝袜在线中文字幕| 91aial.com中文字幕在线观看| 香蕉国产在线看| 最新在线观看一区二区三区 | 午夜福利在线免费观看网站| 电影成人av| 国产精品一国产av| 亚洲国产av影院在线观看| 欧美中文综合在线视频| 精品第一国产精品| 国产人伦9x9x在线观看| 久久久国产一区二区| 制服人妻中文乱码| 久久久精品94久久精品| 成人三级做爰电影| 大码成人一级视频| 亚洲美女视频黄频| 十分钟在线观看高清视频www| 成人三级做爰电影| 久久精品久久精品一区二区三区| 天天躁日日躁夜夜躁夜夜| 五月天丁香电影| 超碰97精品在线观看| 十八禁高潮呻吟视频| 日日啪夜夜爽| 新久久久久国产一级毛片| 亚洲精品一二三| 天天添夜夜摸| 一级,二级,三级黄色视频| 久久国产精品大桥未久av| 久久久精品免费免费高清| 日韩一区二区视频免费看| 看免费成人av毛片| 国产福利在线免费观看视频| 男人操女人黄网站| 大码成人一级视频| 久久99热这里只频精品6学生| 久久国产亚洲av麻豆专区| 97精品久久久久久久久久精品| 亚洲三区欧美一区| 极品少妇高潮喷水抽搐| 久久天堂一区二区三区四区| 国产精品.久久久| 国产精品久久久久久人妻精品电影 | 最近最新中文字幕免费大全7| 黑人欧美特级aaaaaa片| 我的亚洲天堂| 丝袜美足系列| 少妇 在线观看| 国产在线一区二区三区精| 久久久欧美国产精品| 大陆偷拍与自拍| 亚洲av国产av综合av卡| av片东京热男人的天堂| 亚洲精品国产色婷婷电影| 咕卡用的链子| 青春草国产在线视频| 国产野战对白在线观看| 久久亚洲国产成人精品v| 老司机影院毛片| 男女免费视频国产| 男人操女人黄网站| 一本色道久久久久久精品综合| 成人漫画全彩无遮挡| 日韩欧美精品免费久久| 久久精品亚洲熟妇少妇任你| 免费黄网站久久成人精品| 久久久欧美国产精品| 国产精品.久久久| 纵有疾风起免费观看全集完整版| 制服诱惑二区| 日本一区二区免费在线视频| 美女扒开内裤让男人捅视频| 国产成人av激情在线播放| 亚洲综合精品二区| 国产精品99久久99久久久不卡 | 国产av精品麻豆| 婷婷色麻豆天堂久久| 午夜av观看不卡| 国产极品粉嫩免费观看在线| 午夜久久久在线观看| 丰满迷人的少妇在线观看| 咕卡用的链子| 国产亚洲av高清不卡| 免费黄频网站在线观看国产| 国产av一区二区精品久久| 成人免费观看视频高清| 日韩人妻精品一区2区三区| 又黄又粗又硬又大视频| 久久久久网色| 亚洲熟女毛片儿| 国产色婷婷99| 人人妻人人澡人人爽人人夜夜| 狠狠婷婷综合久久久久久88av| 日韩一本色道免费dvd| 最黄视频免费看| 大片电影免费在线观看免费| 国产亚洲最大av| 久久精品久久精品一区二区三区| 人人妻人人添人人爽欧美一区卜| 亚洲综合精品二区| 韩国高清视频一区二区三区| 十八禁高潮呻吟视频| 黄片小视频在线播放| 狠狠婷婷综合久久久久久88av| 欧美日韩成人在线一区二区| 十分钟在线观看高清视频www| 中文字幕av电影在线播放| 亚洲精品国产av成人精品| 啦啦啦啦在线视频资源| 成年动漫av网址| 人体艺术视频欧美日本| 国产精品成人在线| 777米奇影视久久| 中文字幕高清在线视频| 男人爽女人下面视频在线观看| 999久久久国产精品视频| 日韩精品免费视频一区二区三区| 丝袜美腿诱惑在线| 大香蕉久久成人网| 精品午夜福利在线看| 久久97久久精品| 高清欧美精品videossex| 天堂俺去俺来也www色官网| 三上悠亚av全集在线观看| 男人舔女人的私密视频| 国产成人一区二区在线| 亚洲av中文av极速乱| 777米奇影视久久| 国产99久久九九免费精品| 五月天丁香电影| 伊人久久国产一区二区| 熟女少妇亚洲综合色aaa.| 日本wwww免费看| 黑人猛操日本美女一级片| 亚洲四区av| 老熟女久久久| 五月开心婷婷网| 久久99精品国语久久久| 高清av免费在线| 9191精品国产免费久久| 飞空精品影院首页| 久久久欧美国产精品| 中文字幕另类日韩欧美亚洲嫩草| 一区在线观看完整版| 男女高潮啪啪啪动态图| 一边摸一边抽搐一进一出视频| 精品亚洲成国产av| 99re6热这里在线精品视频| 亚洲av中文av极速乱| 日韩av免费高清视频| 国语对白做爰xxxⅹ性视频网站| 久久精品久久久久久久性| 免费看不卡的av| 男的添女的下面高潮视频| 亚洲精品乱久久久久久| 少妇猛男粗大的猛烈进出视频| 欧美亚洲日本最大视频资源| 热99久久久久精品小说推荐| av.在线天堂| 久久久精品国产亚洲av高清涩受| 日本av手机在线免费观看| 亚洲精品国产区一区二| 人妻一区二区av| 久久免费观看电影| 国产激情久久老熟女| 超碰成人久久| 国产亚洲av片在线观看秒播厂| 不卡视频在线观看欧美| 一级,二级,三级黄色视频| 日韩av免费高清视频| 成年人午夜在线观看视频| 欧美av亚洲av综合av国产av | 欧美精品av麻豆av| 大片电影免费在线观看免费| 狂野欧美激情性xxxx| 久久韩国三级中文字幕| 亚洲国产成人一精品久久久| 国产精品二区激情视频| 美女大奶头黄色视频| 亚洲一级一片aⅴ在线观看| 欧美日韩国产mv在线观看视频| 18禁观看日本| 亚洲欧美成人精品一区二区| 亚洲精品,欧美精品| 中文字幕最新亚洲高清| www日本在线高清视频| 国产深夜福利视频在线观看| 亚洲久久久国产精品| 精品免费久久久久久久清纯 | 熟女av电影| 亚洲国产欧美网| 欧美日本中文国产一区发布| 色网站视频免费| 一级毛片 在线播放| 一级毛片我不卡| 高清在线视频一区二区三区| 99国产综合亚洲精品| 日韩av在线免费看完整版不卡| 日韩伦理黄色片| 少妇猛男粗大的猛烈进出视频| 国产乱人偷精品视频| 黑丝袜美女国产一区| 操美女的视频在线观看| av在线播放精品| 国产成人91sexporn| 91精品伊人久久大香线蕉| 九九爱精品视频在线观看| 精品酒店卫生间| 亚洲国产av新网站| 国产免费一区二区三区四区乱码| 亚洲国产精品国产精品| 午夜福利免费观看在线| 国产熟女午夜一区二区三区| 亚洲精品av麻豆狂野| 女人被躁到高潮嗷嗷叫费观| 日韩视频在线欧美| 亚洲欧美精品自产自拍| 一本—道久久a久久精品蜜桃钙片| av卡一久久| 国产一级毛片在线| 免费黄频网站在线观看国产| 十八禁网站网址无遮挡| 久久久久久久久久久免费av| 又大又爽又粗| 久久久久人妻精品一区果冻| 黄色 视频免费看| 黄色一级大片看看| 中文字幕人妻熟女乱码| 国产成人a∨麻豆精品| 精品人妻一区二区三区麻豆| 狠狠婷婷综合久久久久久88av| 国产97色在线日韩免费| av电影中文网址| 午夜影院在线不卡| 精品福利永久在线观看| 男女边吃奶边做爰视频| 侵犯人妻中文字幕一二三四区| 国产伦理片在线播放av一区| 亚洲色图 男人天堂 中文字幕| 精品午夜福利在线看| 欧美日韩成人在线一区二区| 熟女少妇亚洲综合色aaa.| 午夜激情av网站| 国产又色又爽无遮挡免| 美女福利国产在线| 男人爽女人下面视频在线观看| 91成人精品电影| 人人妻人人澡人人爽人人夜夜| 色婷婷av一区二区三区视频| 夜夜骑夜夜射夜夜干| 亚洲精品一二三| 欧美日韩视频精品一区| 亚洲欧洲精品一区二区精品久久久 | 欧美中文综合在线视频| 国产精品一区二区精品视频观看| 国产av一区二区精品久久| 免费人妻精品一区二区三区视频| 成人漫画全彩无遮挡| 18禁观看日本| 免费女性裸体啪啪无遮挡网站| 少妇的丰满在线观看| 少妇 在线观看| 亚洲成av片中文字幕在线观看| 欧美成人精品欧美一级黄| 日韩,欧美,国产一区二区三区| 69精品国产乱码久久久| 王馨瑶露胸无遮挡在线观看| tube8黄色片| 国产av码专区亚洲av| 日韩中文字幕视频在线看片| 国产精品嫩草影院av在线观看| 欧美97在线视频| 国产精品成人在线| xxxhd国产人妻xxx| 黑人巨大精品欧美一区二区蜜桃| 性高湖久久久久久久久免费观看| 成人国产av品久久久| 日本wwww免费看| 免费在线观看视频国产中文字幕亚洲 | 免费观看av网站的网址| 中文字幕色久视频| 97精品久久久久久久久久精品| 亚洲欧美成人精品一区二区| 操出白浆在线播放| 一本大道久久a久久精品| 在线观看免费视频网站a站| 久久精品亚洲av国产电影网| 国产成人免费无遮挡视频| 欧美亚洲日本最大视频资源| 看十八女毛片水多多多| 国产精品秋霞免费鲁丝片| 免费日韩欧美在线观看| h视频一区二区三区| 亚洲精品自拍成人| 高清不卡的av网站| 91成人精品电影| 男的添女的下面高潮视频| 国产高清国产精品国产三级| 欧美成人精品欧美一级黄| 叶爱在线成人免费视频播放| 国产成人欧美在线观看 | 国产欧美亚洲国产| 91aial.com中文字幕在线观看| 久久韩国三级中文字幕| 久久狼人影院| 国产 一区精品| 亚洲欧美一区二区三区久久| 日韩 欧美 亚洲 中文字幕| 免费高清在线观看视频在线观看| 在线免费观看不下载黄p国产| 国产成人一区二区在线| 午夜老司机福利片| 欧美日韩精品网址| 国产在线视频一区二区| 少妇的丰满在线观看| 亚洲熟女毛片儿| 久久久国产精品麻豆| 一本色道久久久久久精品综合| 国产精品免费视频内射| 午夜免费男女啪啪视频观看| 亚洲精品,欧美精品| 一级毛片我不卡| 国产淫语在线视频| 只有这里有精品99| 99九九在线精品视频| 日本91视频免费播放| 大香蕉久久成人网| 女人爽到高潮嗷嗷叫在线视频| 久久婷婷青草| 少妇人妻 视频| 99久久99久久久精品蜜桃| 日本欧美视频一区| 亚洲熟女精品中文字幕| a级毛片黄视频| 美女高潮到喷水免费观看| 精品久久蜜臀av无| 热re99久久精品国产66热6| 欧美日韩一区二区视频在线观看视频在线| 国产视频首页在线观看| 熟妇人妻不卡中文字幕| 男女免费视频国产| 啦啦啦中文免费视频观看日本| 考比视频在线观看| 观看av在线不卡| 夜夜骑夜夜射夜夜干| 欧美激情 高清一区二区三区| e午夜精品久久久久久久| 久久久精品区二区三区| 欧美日韩亚洲综合一区二区三区_| 9色porny在线观看| 亚洲av在线观看美女高潮| 电影成人av| 我的亚洲天堂| 日韩欧美一区视频在线观看| 亚洲av成人精品一二三区| 精品国产超薄肉色丝袜足j| 51午夜福利影视在线观看| 久久影院123| 波多野结衣av一区二区av| 亚洲美女黄色视频免费看| 久久久久网色| 看十八女毛片水多多多| 亚洲欧洲精品一区二区精品久久久 | 丝袜脚勾引网站| 亚洲av日韩精品久久久久久密 | 一区在线观看完整版| av电影中文网址| 亚洲中文av在线| av网站在线播放免费| 在线亚洲精品国产二区图片欧美| 日韩 欧美 亚洲 中文字幕| 中文天堂在线官网| 高清不卡的av网站| 国产成人午夜福利电影在线观看| 黄色 视频免费看| 日韩中文字幕欧美一区二区 | 99国产综合亚洲精品| 国产日韩一区二区三区精品不卡| 99久久人妻综合| 又大又爽又粗| 欧美 亚洲 国产 日韩一| 久久精品人人爽人人爽视色| 免费av中文字幕在线| av在线app专区| 最黄视频免费看| 在现免费观看毛片| 国产一区二区激情短视频 | 亚洲国产精品999| www.精华液| 精品亚洲成国产av| 操出白浆在线播放| 久久精品久久久久久久性| a级毛片在线看网站| 亚洲自偷自拍图片 自拍| 欧美激情极品国产一区二区三区| 久久久久国产一级毛片高清牌| 男人舔女人的私密视频| 欧美精品一区二区免费开放| 狂野欧美激情性bbbbbb| 亚洲国产精品一区二区三区在线| 亚洲,一卡二卡三卡| 亚洲精品一二三| 免费久久久久久久精品成人欧美视频| 国产日韩欧美亚洲二区| 国产一区二区三区av在线| 国产精品一区二区精品视频观看| 狠狠婷婷综合久久久久久88av| 黄色毛片三级朝国网站| 一本一本久久a久久精品综合妖精| 黑人猛操日本美女一级片| 国产精品久久久av美女十八| 国产成人精品无人区| 精品卡一卡二卡四卡免费| 欧美黄色片欧美黄色片| 一级片'在线观看视频| 无遮挡黄片免费观看| 精品一区二区免费观看| 欧美亚洲日本最大视频资源| 亚洲一区二区三区欧美精品| 欧美人与性动交α欧美软件| 久久女婷五月综合色啪小说| 久久综合国产亚洲精品| 另类亚洲欧美激情| 大片电影免费在线观看免费| 亚洲欧美色中文字幕在线| 免费人妻精品一区二区三区视频| 不卡av一区二区三区| 午夜福利免费观看在线| 久久热在线av| 午夜福利一区二区在线看| 各种免费的搞黄视频| 建设人人有责人人尽责人人享有的| 国产精品一区二区在线观看99| 精品国产国语对白av| 伦理电影大哥的女人| 亚洲精品国产av蜜桃| 免费观看人在逋| 天天躁日日躁夜夜躁夜夜| 波野结衣二区三区在线| 午夜福利视频精品| 亚洲成av片中文字幕在线观看| 欧美精品av麻豆av| 考比视频在线观看| 成人三级做爰电影| 欧美在线黄色| 少妇的丰满在线观看| av在线老鸭窝| 国产精品秋霞免费鲁丝片| 国产成人欧美| 免费人妻精品一区二区三区视频| 久久久久精品性色| 在线看a的网站| 在线 av 中文字幕| 满18在线观看网站| 精品久久久精品久久久| 国产高清不卡午夜福利| 午夜日韩欧美国产| 久久av网站| 中国国产av一级| 少妇被粗大的猛进出69影院| av网站在线播放免费| 国产免费又黄又爽又色| 久久久久视频综合| 久久久国产欧美日韩av| 少妇人妻 视频| 中文字幕制服av| 亚洲欧洲国产日韩| 国产欧美日韩一区二区三区在线| 99香蕉大伊视频| 欧美日韩亚洲高清精品| 街头女战士在线观看网站| 一级黄片播放器| 汤姆久久久久久久影院中文字幕| 韩国av在线不卡| 国产精品偷伦视频观看了| 亚洲欧美激情在线| 成人毛片60女人毛片免费| 亚洲一区中文字幕在线| 亚洲国产精品一区二区三区在线| 97人妻天天添夜夜摸| 色精品久久人妻99蜜桃| 一级a爱视频在线免费观看| 中文天堂在线官网| 亚洲一区中文字幕在线| 午夜福利网站1000一区二区三区| 日韩 亚洲 欧美在线| 国产精品一区二区在线观看99| 最近手机中文字幕大全| 国产av码专区亚洲av| 欧美国产精品一级二级三级| 亚洲视频免费观看视频| 久久ye,这里只有精品| a级片在线免费高清观看视频| 丰满迷人的少妇在线观看| 久久精品国产亚洲av涩爱| 看免费av毛片| 啦啦啦 在线观看视频| 国产免费视频播放在线视频| 国产精品av久久久久免费| 99久久99久久久精品蜜桃| 香蕉丝袜av| www日本在线高清视频| 电影成人av| 韩国高清视频一区二区三区| 欧美日韩视频精品一区| 国产精品久久久久久久久免| 精品国产乱码久久久久久小说| 亚洲精品一区蜜桃| 中国三级夫妇交换| 免费少妇av软件| 啦啦啦在线观看免费高清www| 午夜福利乱码中文字幕| 国产精品女同一区二区软件| 黄色视频不卡| 王馨瑶露胸无遮挡在线观看| 在线天堂最新版资源| 午夜激情av网站| 久久精品亚洲av国产电影网| 18禁动态无遮挡网站| 精品一区二区三区四区五区乱码 | 国语对白做爰xxxⅹ性视频网站| 一级毛片黄色毛片免费观看视频| 国产精品熟女久久久久浪| 电影成人av| 一区二区三区激情视频| 制服诱惑二区| 麻豆乱淫一区二区| 街头女战士在线观看网站| 韩国精品一区二区三区|