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

    Comparative Analysis of Characteristics of the Coupled and Decoupled Parallel Mechanisms

    2010-03-01 01:48:46ZENGDaxingHOUYuleiLUWenjuanandHUANGZhen

    ZENG Daxing, HOU Yulei, LU Wenjuan, and HUANG Zhen

    College of Mechanical Engineering, Yanshan University, Qinhuangdao 066004, China

    1 Introduction

    Just as mentioned on philosophy, things are always the unity of opposites, and this is also confirmed in the research process of mechanism. When the serial mechanism is forging ahead in its development, the parallel mechanism (PM), in which the fixed base and the moving platform are connected by two or more limbs, becomes the research focus of the international academic and engineering field gradually.

    Compared with the serial mechanism, the parallel mechanism has some outstanding properties in certain aspects, such as a large payload to weight ratio,considerable stiffness, low inertia, high dynamic performances, and so on. The advantages of the PM are closely related to the coupling between the limbs. While it is no other than coupling that makes the PM possess some disadvantages, including smaller workspace, lower dexterity, more complex control command, and multiple singularities. The difficulties in the structural design, static and dynamic analysis and development of control system of the PM are also caused by the existence of the coupling,which restrict the applications of the PM.

    Then the decoupled parallel mechanism (DPM), in which there exists a one-to-one correspondence relationship between the input and output variables, attracts many scholars’ interests[1–3]. In this paper, the kinematic pairs will be denoted with following symbols: P for prismatic pair, R for revolute pair, U for universal pair, and C for cylindrical pair. HUANG and LI[4]invented a fully-isotropic 3-CPR three-degrees-of-freedom (3-DOF) decoupled translational PM. CARRICATO, et al, discussed the general problem of the topological synthesis and classification of translational PMs, and investigated both their constraint and direct singularities[5], and presented a pointing parallel mechanism U-PUR-PRRU with fully decoupled degrees of freedom[6].KONG and GOSSELIN proposed a 3-CRR 3-DOF translational parallel manipulator with linear actuators[7]and three classes of input-output decoupled parallel manipulators with 2 to 4 DOFs using a geometric approach[8]. KIM and TSAI[9]introduced a 3-DOF Cartesian parallel manipulator and performed its optimal design. LI, et al[10], presented a 3-DOF translational parallel manipulator named R-CUBE which has decoupled motion in three axes. GOGU[11]presented singularity-free fully-isotropic parallel wrists with 2-DOF and proposed a method for structural synthesis based on the theory of linear transformations.

    Generally speaking, the DPM not only keeps the principal advantages of the PM, but also avoids the main drawbacks of the general coupled PM. To reveal the separated characteristics of the coupled and decoupled PMs,taking the 3-RPUR PM and the 3-CPR PM as example respectively, this paper focuses on their comparative analysis.

    The organization of this paper is as follows: Following the introduction, the structure composition of the 3-RPUR PM and the 3-CPR PM are described and their modes of motion are analyzed in section 2. Displacement analysis,including the forward and inverse displacement solution, is performed in section 3. Section 4 concerns performance analysis, including workspace, dexterity, velocity, payload capability and stiffness etc. The paper is concluded in section 5, summarizing the present work.

    2 Mechanism Features Analysis

    As shown in Fig. 1, the 3-RPUR PM is composed of an upper platform, a lower platform and three uniform limbs connecting the two platforms with R, P, U and R in turn.The connection points are located on a circumcircle symmetrically 120° apart. The axis of the R pair connected with the lower platform is parallel with one of the two axes of the U pair, the motion direction of the P pair is perpendicular to these two parallel axes, and the other axis of the U pair is parallel with the axis of the R pair connected with the upper platform. Taking the upper and lower platform as the moving platform and the fixed base respectively, we set up the moving coordinate system o-xyz and the stationary coordinate system O-XYZ with their origins located at the geometrical center of the upper and lower platform respectively. In addition, the P pairs of three limbs are selected as the input pairs. The circumcircles of the upper and lower platform are denoted as r and R,respectively. The distance between Ai(i=1, 2, 3) and Bi(i=1, 2, 3) is denoted as li(i=1, 2, 3), and the length of the pole BiCi(i=1, 2, 3) is supposed as L.

    Fig. 1. Schematic diagram of 3-RPUR 3-DOF PM

    The coordinate system oi-xiyizi(i=1, 2, 3) is built with its origin located on the center of the R pair of the ith limb connected with the fixed base, in which yiaxis is consistent with the axis of this R pair, and ziaxis is perpendicular to the fixed base plane. Selecting the limb A1B1C1as example,we can obtain the screw system of the limb in o1-x1y1z1as follows:

    where a, b, c, d, t, l, e, f, g, h, p and q are all finite non-zero real number.

    The reciprocal screw of the screw system in Eq. (1) is

    Similarly, the reciprocal screwsandof the other two limbs A2B2C2and A3B3C3can be solved, and they are also perpendicular to the plane formed by the two rotational axes of U pair of the corresponding limb respectively. Obviously,andare noncoplanar which accounts for that they are linear independent. Therefore the three rotations of the moving platform are constrained completely, that is, the mechanism only has three translations.

    The mobility can also be obtained by the modified Kutzbach-Grübler criterion[12]. What should be noted is that the 3-RPUR PM has not common constraint, passive DOF,and redundant constraints, and then its DOF is

    where M is the DOF of the PM; d the order of the mechanism, and d= 6?λ ; λ the common constraint;n denotes the number of links; g the number of kinematic pairs; fithe freedom of the ith pair; ν the number of redundant constraints; ? the passive degree of freedom.

    Fig. 2 shows the schematic diagram of the 3-CPR PM with three identical limbs. Each limb connects the moving platform and the fixed base with the R pair and the C pair,respectively. The axis of the R pair is parallel with that of the C pair, and the motion direction of the P pair, which locates between the R pair and the C pair, is perpendicular to these two parallel axes. In addition, the axes of the C pairs of the three limbs keep perpendicular one another.The stationary coordinate system O-XYZ is set up with its origin located at the intersection of the three C pairs, and the X, Y and Z axis are collinear with the axis of the three C pairs respectively. The moving coordinate system o-xyz is set up with its origins located at the geometrical center of the triangle B1B2B3, and the x, y and z axes are parallel with that of the stationary coordinate system respectively. The translation of the three C pairs are selected as the input, and the input quantity are denoted as L11, L21and L31respectively. The vertical distance from the origin of the coordinate system o-xyz to the axis of the R pair connecting with the moving platform is denoted as m, and the length of the link AiBi(i=1, 2, 3) is denoted as Li(i=1, 2, 3).

    Fig. 2. Schematic diagram of 3-CPR 3-DOF PM

    Taking the limb A1B1as example, we can obtain its screw system as follows:

    where a, b, c and d are all finite non-zero real number.

    Then the reciprocal screws of the screw system of the limb A1B1can be solved as

    Similarly, the reciprocal screwsandof the limb A2B2which around X axis and Z axis respectively, andandof the limb A3B3which around X axis and Y axis respectively, can also be solved.

    Then it is obvious that these six reciprocal screws constrain all the rotation of the moving platform and there are three redundant constraints, so this PM only has three translations.

    By the modified Kutzbach-Grübler criterion, the DOF of the 3-CPR PM can be solved as

    3 Displacement Analysis

    3.1 Inverse displacement solution

    For the 3-RPUR PM, the coordinates of points Aiand Ci(i=1, 2, 3) with respect to O-XYZ can be expressed as follows:

    On one hand, the trajectory of point Bi(i=1, 2, 3) locates on a circle of which takes Ci(i=1, 2, 3) as the centre and L as radius, on the other hand, Bialso locates on a circle of which takes Ai(i=1, 2, 3) as the centre and li(i=1, 2, 3) as radius. Then the input parameter1l, l2and3l can be obtained as

    where

    When derivatives are taken with respect to time for Eqs.(7), (8) and (9), the relationship between the input and output velocity can be expressed with matrix form as

    where

    If matrix Jiis not singular, the Jacobian matrix of the 3-RPUR PM can be solved as

    For the 3-CPR PM, the relationship between the input and output can be directly obtained from Fig. 2, that is

    Taking derivatives with respect to time for Eq. (12), the Jacobian matrix of the 3-CPR PM can be obtained as

    3.2 Forward displacement solution

    For the 3-RPUR PM, the forward displacement solution needs to solve the simultaneous Eqs. (7), (8) and (9) with knownNevertheless, it is comparative difficult to obtain the closed solution by using current methods.

    While for the 3-CPR PM, if the driving inputis determined, the location coordinates of the moving platform are

    From the above derivation, it can be seen that, for the 3-RPUR PM, the input-output equations are highly coupled and the expression of Jacobian matrix is complex, which will result in the difficulty in its inverse and forward displacement solution, and bring out some disadvantages for the performance analysis. While the Jacobian matrix of the 3-CPR PM is a constant matrix, and the inverse and forward displacement solution of this mechanism is so simple to solve, which is convenient for further research and practical application.

    4 Perfromance Analysis

    4.1 Normalization of structural parameters

    The theory of physical model of the solution space[13–15]provides a tool for the research of the relationship between the performance evaluation indices and the structural parameters. By using the performance indices atlases, we can determine the domain of structural parameters with good values of the indices conveniently.

    The three limbs of the 3-RPUR PM possess the same structure and distribute symmetrically, then without consideration of the twist angle between the axis of the pair and the moving platform or the fixed base, this mechanism has three variable parameters r, R and L. For the simplification consideration, we normalize the structural parameters with linear dimensions as follows:

    Eq. (15) should subject to

    In the following analysis, the initial values of structural parameters of the 3-RPUR PM are supposed as r′=0.6 ,

    For the similar consideration, the dimensionless input parameters of 3-RPUR PM can be expressed as

    For the 3-CPR PM, there is only one variable parameters m, which does not influence the output and the Jacobian matrix of the PM without consideration of the twist angles between the axes of pairs.

    For the comparability between the 3-RPUR and the 3-CPR PM, the same variable range of the input is set as follows:

    4.2 Workspace analysis

    As one of the important kinematics indices, the size of the workspace determines the range of motion of the mechanism. For a PM with good application value in practical engineering, the shape of its workspace should be regular, and the mechanism should have better kinematics performances and do not exist singular configuration in the whole workspace.

    According to Eqs. (7), (8), (9), (15) and (16), a mushroom workspace of the 3-RPUR PM is obtained and shown in Fig. 3. The workspace distributes symmetrically 120° apart and the thickness of the center is about 0.6.

    Fig. 3. Workspace of the 3-RPUR PM

    On the basis of the theory of physical model of the solution space, the workspace of the 3-CPR PM is described as shown in Fig. 4. This workspace is just like a cube of which each side is 1, and there does not exist singular point in the whole workspace.

    By comparing Fig. 3 and Fig. 4, it can be seen that the available workspace of the 3-CPR PM is bigger than that of the 3-RPUR PM, so the application fields of the former should be wider than the latter from this point.

    Fig. 4. Workspace of the 3-CPR PM

    4.3 Dexterity analysis

    Generally speaking, the condition number of the Jacobian matrix is taken as the performance evaluation index for the dexterity analysis[16], and the bigger the condition number is, the worse the dexterity is. The condition number of the Jacobian matrix can be defined as[17]

    where σmax( J ) and σmin( J )are the maximum and minimum singular values of Jacobian matrix, respectively.

    According to Eqs. (10), (11), and (19), the dexterity atlases of the 3-RPUR PM are shown in Fig. 5.

    Fig. 5. Dexterity atlases of 3-PRUR PM

    In fact, the condition number defined in Eq. (19) only describes the dexterity of the mechanism in a certain pose,to give an overall evaluation on the dexterity of a certain workspace, the global conditioning index, proposed by GOSSELIN, et al[18], is defined as follows:

    whereJη denotes the global conditioning index, andW denotes the workspace of the mechanism.

    The global conditioning index is essentially the average of the reciprocal of the condition number of the Jacobian matrix in a certain workspace. The biggerJη is, the higher the dexterity of the mechanism is.

    According to Eq. (20), the global conditioning index of the 3-RPUR PM can be calculated out, which is equal to 0.046 414 and is far smaller than 1. This shows that the dexterity of this mechanism is very poor.

    While for the 3-CPR PM, according to Eq. (13), the condition number of the Jacobian matrix is equivalent to 1,which indicates that the global conditioning index is equal to 1 and keeps invariable in the whole workspace. So it can be said that the 3-CPR PM possesses isotropic kinematic performances and its dexterity is very good.

    4.4 Velocity analysis

    Denoting the input and output velocity vector of the 3-RPUR PM asand V respectively, we have

    According to Eqs. (10), (11) and (22), the velocity index atlases of the 3-RPUR PM are plotted and shown in Fig. 6.

    Being similar to the condition number discussed above,the velocity extremum only describes the maximum or minimum of the velocity of the mechanism on a certain point. For the purpose of evaluating the velocity performance in a certain workspace, the global velocity indices[19]are defined as follows:

    where ?Vmaxand ?Vminare the maximum and minimum global velocity indices, respectively; W denotes the workspace of the mechanism.

    Fig. 6. Atlases of the velocity index of the 3-RPUR PM

    Then the maximum and minimum global velocity indices of the 3-RPUR PM in the whole workspace can be solved,which are equal to 68.393 and 0.543 7, respectively.

    Owing to the Jacobian matrix of the 3-CPR PM is a unit matrix, the maximum and minimum of the output velocity vector are both equal to 1, and keep invariable in the workspace. Therefore the maximum and minimum global velocity indices of the 3-CPR PM are also equivalent to 1.

    From the above analysis, it can be seen that the minimum global velocity index of the 3-CPR PM is bigger than that of the 3-RPUR PM, while the condition of the maximum global velocity index is just opposite to that of the minimum one. So the coupled 3-RPUR PM is more suitable for the application which requires high velocity in some directions and low velocity in the other directions,and the decoupled 3-CPR PM is more suitable for the situation requiring the uniform velocity in all directions.

    4.5 Payload capability analysis

    The relationship between the generalized force vector f and the external force vector F applied to the end actuator of the mechanism can be expressed as

    where G is the force Jacobian matrix, and

    According to Eqs. (10), (11) and (25), the payload capability index atlases of the 3-RPUR PM are plotted and shown in Fig. 7.

    The maximum or minimum obtained from Eq. (25) only denotes the payload capability of the mechanism when the moving platform locates on a certain position in the workspace, then to describe the payload capability in the whole workspace, the global payload capability indices are defined as follows[20]:

    where ?Fmaxand ?Fminare the maximum and minimum global payload capability indices, respectively; W denotes the workspace of the mechanism.

    Then the maximum and minimum global payload capability indices of the 3-RPUR PM in the whole workspace can be solved, which are equal to 1.865 7 and 0.091 338, respectively. The maximum and minimum global payload capability indices of the 3-CPR PM are both equivalent to 1.

    From the above analysis, it can be seen that the minimum global payload capability index of the 3-CPR PM is far bigger than that of the 3-RPUR PM, while the maximum global payload capability index of the 3-RPUR PM is a little bigger than that of the 3-CPR PM. So the coupled 3-RPUR PM can be used for the machine which requires large payload capability in some directions and small payload capability in the other directions, and the decoupled 3-CPR PM can be used for the situation requiring the uniform payload capability in all directions,and should possess wider application fields.

    Fig. 7. Atlases of the payload capability index of the 3-RPUR PM

    4.6 Stiffness analysis

    The relationship between the deformation D of the end actuator and the general force vector F applied to the mechanism can be expressed as

    where C is the flexibility matrix, and

    According to Eqs. (10), (11) and (28), the stiffness index atlases of the 3-RPUR PM are plotted and shown in Fig. 8.

    Fig. 8. Atlases of the stiffness index of 3-PRUR PM

    The maximum or minimum deformation obtained from Eq. (28) only reflects the deformation situations of the mechanism in a certain position, then to evaluate the stiffness performance of the mechanism in the whole workspace, the global stiffness indices are defined as follows[21]:

    where ?Dmaxand ?Dminare the maximum and minimum global stiffness indices, respectively; W denotes the workspace of the mechanism.

    Then the maximum and minimum global stiffness indices of the 3-RPUR PM in the whole workspace can be solved, which are equal to 1 228.5 and 0.300 38,respectively. The maximum and minimum global stiffness indices of the 3-CPR PM are both equivalent to 1.

    From the above analysis, it can be seen that the maximum global stiffness index of the 3-RPUR PM is far bigger than that of the 3-CPR PM, while the minimum global stiffness index of the 3-RPUR PM is smaller than that of the 3-CPR PM. This indicates that the coupled 3-RPUR PM is more suitable for the application requiring large stiffness in some directions.

    5 Conclusions

    (1) Taking the 3-RPUR PM and the 3-CPR PM as example, the comparative analysis of characteristics of the coupled and decoupled PMs, including the kinematics,workspace, dexterity, velocity, payload capability and stiffness, has been performed.

    (2) Compared with the coupled PM, the Jacobian matrix of the DPM is a constant matrix, which is independent of the pose of the mechanism. The condition number of the Jacobian matrix is equal to 1, and the mechanism possesses excellent force and motion transmission capabilities.

    (3) Both the forward and the inverse displacement solutions of the DPM are straightforward and almost do not need to calculate. The shape of the workspace of the DPM is regular, and it is free from singularities in the whole workspace.

    (4) The velocity, payload capability and stiffness indices of the coupled PM are better than that of the DPM in some directions. While the performance indices of the DPM,including dexterity, velocity, payload capability and stiffness, keep invariable in the workspace.

    (5) The higher the decoupling extent of the PM is, the easier the kinematics and dynamics analysis are, which is more favorable for the control and the improvement of the movement precision.

    (6) The coupled and decoupled PMs have their own advantages and disadvantages, and which one is more suitable for a certain application should depend on the specific situation.

    [1] PATARINSKI S P, UCHIYAMA M. Position/orientation decoupled parallel manipulator[C]//Proceedings of International Conference on Advanced Robotics, Tokyo, Japan, 1993: 153–158.

    [2] INNOCENTI C, PARENTI-CASTELLI V. Direct kinematics of the 6–4 fully parallel manipulator with position and orientation uncoupled[C]//Proceedings of the European Robotics and Intelligent Systems, Corfou, 1991: 23–28.

    [3] RYU D, CHO C, KIM M, et al. Design of a 6 DOF haptic master for teleoperation of a mobile manipulator[C]//Proceedings of IEEE International Conference on Robotics and Automation, Taipei,China, 2003: 3 243–3 248.

    [4] HUANG Zhen, LI Qinchuan. 3-DOF translational parallel robot mechanism: China, 01104454.3[P]. 2001-09-05.

    [5] CARRICATO M, PARENTI-CASTELLI V. Singularity-free fully-isotropic translational parallel mechanisms[J]. The International Journal of Robotics Research, 2002, 21(2): 161–174.

    [6] CARRICATO M, PARENTI-CASTELLI V. A novel fully decoupled two-degrees-of-freedom parallel wrist[J]. The International Journal of Robotics Research, 2004, 23(6): 661–667.

    [7] KONG Xianwen, GOSSELIN C M. Kinematics and singularity analysis of a novel type of 3-CRR 3-DOF translational parallel manipulator[J]. The International Journal of Robotics Research,2002, 21(9): 791–798.

    [8] KONG Xianwen, GOSSELIN C M. Type synthesis of input-output decoupled parallel manipulators[J]. Transactions of the Canadian Society for Mechanical Engineering, 2004, 28(2): 185–196.

    [9] KIM H S, TSAI L W. Design optimization of a Cartesian parallel manipulator[C]//Proceedings of the ASME Design Engineering Technical Conference, Montreal, Quebec, Canada, September 29,2002, 5B: 865–872.

    [10] LI Weimin, GAO Feng, ZHANG Jianjun. R-CUBE, a decoupled parallel manipulator only with revolute joints[J]. Mechanism and Machine Theory, 2005, 40(4): 467–473.

    [11] GOGU G. Fully-isotropic over-constrained parallel wrists with two degrees of freedom[C]//Proceedings of the 2005 IEEE International Conference on Robotics and Automation, Barcelona, Spain, 2005:4 025–4 030.

    [12] HUANG Zhen, ZHAO Yongsheng, ZHAO Tieshi. Advanced spatial mechanism[M]. Beijing: Higher Education Press, 2006. (in Chinese)

    [13] DAVIES T H, BAKER J E, THOMPSON A G R. A finite,3-dimensional atlas of 4-bar linkages[J]. Mechanism and Machine Theory, 1979, 14(6): 389–403.

    [14] BARKER C R, WEEKS G A. A physical model of the solution space for four bar mechanisms[C]//Proceedings-OSU Applied Mechanism Conference, St Louis, MO, USA, 1983: 26–33.

    [15] YANG Jihou. The space model and dimensional types of the four-bar mechanisms[J]. Mechanism and Machine Theory, 1987,22(1): 71–76.

    [16] POND G, CARRETERO J A. Formulating Jacobian matrices for the dexterity analysis of parallel manipulators[J]. Mechanism and Machine Theory, 2006, 41(12): 1 505–1 519.

    [17] SALISBURY J K, CRAIG J J. Articulated hands: force control and kinematic issues[J]. The International Journal of Robotics Research,1982, 1(1): 4–17.

    [18] GOSSELIN C M, ANGELES J. A global performance index for the kinematic optimization of robotic manipulators[J]. Journal of Mechanical Design, 1991, 113(3): 220–226.

    [19] GAO Feng, LIU Xinjun, GRUVER W A. Performance evaluation of two-degree-of-freedom planar parallel robots[J]. Mechanism and Machine Theory, 1998, 33(6): 661–668.

    [20] LIU Xinjun, WANG Jinsong, PRITSCHOW G. Performance atlases and optimum design of planar 5R symmetrical parallel mechanisms[J]. Mechanism and Machine Theory, 2006, 41(2):119–144.

    [21] LIU Xinjun, JIN Zhenlin, GAO Feng. Optimum design of 3-DOF spherical parallel manipulators with respect to the conditioning and stiffness indices[J]. Mechanism and Machine Theory, 2000, 35(9):1 257–1 267.

    亚洲av第一区精品v没综合| 国产精品免费一区二区三区在线| 欧美日韩瑟瑟在线播放| 国产熟女xx| 亚洲一区中文字幕在线| 美女国产高潮福利片在线看| 丝袜美足系列| 国产伦一二天堂av在线观看| 99久久99久久久精品蜜桃| 国产精品98久久久久久宅男小说| 亚洲三区欧美一区| 国产高清视频在线播放一区| 国产一区二区三区视频了| 天天躁夜夜躁狠狠躁躁| 韩国精品一区二区三区| 搡老岳熟女国产| 亚洲视频免费观看视频| 精品乱码久久久久久99久播| 亚洲熟妇中文字幕五十中出 | 国产精品久久久人人做人人爽| 两个人免费观看高清视频| 国产高清激情床上av| 久久人人精品亚洲av| 欧美+亚洲+日韩+国产| 亚洲国产精品sss在线观看 | 老熟妇仑乱视频hdxx| 伊人久久大香线蕉亚洲五| 国产深夜福利视频在线观看| 亚洲av成人不卡在线观看播放网| 国产深夜福利视频在线观看| av中文乱码字幕在线| 天天影视国产精品| 国产免费av片在线观看野外av| 亚洲成人国产一区在线观看| 啦啦啦 在线观看视频| 99精国产麻豆久久婷婷| 国产精品久久视频播放| 久久这里只有精品19| 欧洲精品卡2卡3卡4卡5卡区| 黑人欧美特级aaaaaa片| 国产成人免费无遮挡视频| 超碰成人久久| 色哟哟哟哟哟哟| 国产精品1区2区在线观看.| 色老头精品视频在线观看| 国产视频一区二区在线看| 国产欧美日韩一区二区三区在线| av中文乱码字幕在线| 男女午夜视频在线观看| 精品久久久久久,| 亚洲激情在线av| 国产三级黄色录像| 日本a在线网址| 国产精品久久视频播放| 午夜老司机福利片| 精品日产1卡2卡| 很黄的视频免费| 久久久国产成人精品二区 | 午夜福利欧美成人| 男女床上黄色一级片免费看| 成人免费观看视频高清| 国产精品一区二区三区四区久久 | 高清在线国产一区| 精品卡一卡二卡四卡免费| 久热爱精品视频在线9| 伊人久久大香线蕉亚洲五| 波多野结衣av一区二区av| 另类亚洲欧美激情| 美女高潮到喷水免费观看| 搡老乐熟女国产| x7x7x7水蜜桃| 女人被狂操c到高潮| 欧美久久黑人一区二区| 色综合婷婷激情| a级毛片黄视频| 老司机靠b影院| 日韩精品免费视频一区二区三区| 中文字幕人妻丝袜一区二区| 国产成人精品无人区| 女人精品久久久久毛片| 日日摸夜夜添夜夜添小说| 成年版毛片免费区| 亚洲精品国产区一区二| 国产一区在线观看成人免费| 亚洲,欧美精品.| 国产精品自产拍在线观看55亚洲| a级毛片在线看网站| 国产成人影院久久av| 很黄的视频免费| 91大片在线观看| 老司机深夜福利视频在线观看| 999久久久国产精品视频| 18禁国产床啪视频网站| 欧美日韩乱码在线| 久久久久国内视频| 久久精品国产99精品国产亚洲性色 | 精品福利永久在线观看| 国产成人精品无人区| 中文欧美无线码| 亚洲精品国产色婷婷电影| 香蕉久久夜色| 久久久久久免费高清国产稀缺| 午夜视频精品福利| 搡老岳熟女国产| 久久精品影院6| 久久亚洲真实| www国产在线视频色| 麻豆一二三区av精品| 国产亚洲精品久久久久久毛片| 国产精品影院久久| 岛国视频午夜一区免费看| 99精国产麻豆久久婷婷| 操出白浆在线播放| 国产片内射在线| 又黄又粗又硬又大视频| 男人舔女人的私密视频| 亚洲熟妇熟女久久| 亚洲九九香蕉| 欧美精品一区二区免费开放| 国产单亲对白刺激| 精品久久久久久久久久免费视频 | 国产主播在线观看一区二区| 成在线人永久免费视频| 亚洲人成电影观看| 亚洲成a人片在线一区二区| 免费观看精品视频网站| 精品国产亚洲在线| 久久精品国产亚洲av香蕉五月| av在线天堂中文字幕 | 午夜福利影视在线免费观看| 午夜精品国产一区二区电影| 日韩欧美一区视频在线观看| 90打野战视频偷拍视频| 欧美乱妇无乱码| 久久热在线av| 热re99久久精品国产66热6| 一区二区日韩欧美中文字幕| 性色av乱码一区二区三区2| 成人国语在线视频| 国产高清videossex| 一本综合久久免费| 亚洲专区中文字幕在线| 亚洲av美国av| 99riav亚洲国产免费| 丰满人妻熟妇乱又伦精品不卡| 老司机在亚洲福利影院| 日韩 欧美 亚洲 中文字幕| 日本精品一区二区三区蜜桃| 日韩av在线大香蕉| 在线观看免费午夜福利视频| 男女高潮啪啪啪动态图| 午夜福利影视在线免费观看| 大型黄色视频在线免费观看| 精品一品国产午夜福利视频| 男女之事视频高清在线观看| 一区在线观看完整版| 男女高潮啪啪啪动态图| 欧美性长视频在线观看| 久久久久亚洲av毛片大全| 美女高潮喷水抽搐中文字幕| cao死你这个sao货| av网站在线播放免费| 成人永久免费在线观看视频| 国产精品影院久久| 国产成人精品在线电影| 性少妇av在线| 亚洲国产欧美网| 亚洲av成人不卡在线观看播放网| 亚洲av日韩精品久久久久久密| 99精品久久久久人妻精品| 国产成人免费无遮挡视频| 悠悠久久av| 久久久水蜜桃国产精品网| 久9热在线精品视频| 50天的宝宝边吃奶边哭怎么回事| 国产精品久久久久成人av| 亚洲 国产 在线| 69av精品久久久久久| 亚洲精品国产区一区二| 亚洲第一青青草原| 日日干狠狠操夜夜爽| www国产在线视频色| 亚洲专区字幕在线| 人人妻,人人澡人人爽秒播| 90打野战视频偷拍视频| 两人在一起打扑克的视频| 国产成人免费无遮挡视频| 99久久久亚洲精品蜜臀av| 一二三四社区在线视频社区8| 制服诱惑二区| 欧美大码av| 亚洲欧美日韩另类电影网站| 国产一卡二卡三卡精品| 中文字幕另类日韩欧美亚洲嫩草| bbb黄色大片| 91在线观看av| 午夜精品国产一区二区电影| 女生性感内裤真人,穿戴方法视频| 母亲3免费完整高清在线观看| 精品福利永久在线观看| 午夜免费激情av| 搡老乐熟女国产| 久久久久久久久中文| 国产成+人综合+亚洲专区| 波多野结衣av一区二区av| 大型av网站在线播放| www.精华液| 精品少妇一区二区三区视频日本电影| av欧美777| 国产黄色免费在线视频| 在线av久久热| 欧美精品亚洲一区二区| 日本 av在线| 国产成人影院久久av| 亚洲欧美日韩高清在线视频| 亚洲精品一卡2卡三卡4卡5卡| 亚洲片人在线观看| 麻豆av在线久日| 国产精品永久免费网站| 亚洲熟女毛片儿| 男女床上黄色一级片免费看| 国产精品美女特级片免费视频播放器 | 日韩三级视频一区二区三区| 亚洲欧美激情综合另类| 超碰成人久久| 国产不卡一卡二| 交换朋友夫妻互换小说| 成人三级做爰电影| 亚洲专区国产一区二区| 亚洲精品久久成人aⅴ小说| 久久精品成人免费网站| 久久午夜综合久久蜜桃| 老汉色∧v一级毛片| 88av欧美| 欧美乱妇无乱码| 免费一级毛片在线播放高清视频 | 99re在线观看精品视频| 亚洲国产欧美一区二区综合| 精品久久久久久成人av| 女同久久另类99精品国产91| 亚洲成人国产一区在线观看| 青草久久国产| 日韩大尺度精品在线看网址 | 韩国精品一区二区三区| 大陆偷拍与自拍| 久久久国产成人免费| 国产99久久九九免费精品| 亚洲伊人色综图| а√天堂www在线а√下载| 久久久国产一区二区| 少妇 在线观看| 亚洲精品美女久久av网站| 欧美日韩一级在线毛片| 国产1区2区3区精品| 日本免费一区二区三区高清不卡 | 精品福利永久在线观看| 亚洲男人的天堂狠狠| 久久国产精品影院| 国产成年人精品一区二区 | 成人国产一区最新在线观看| 精品国产乱码久久久久久男人| 午夜久久久在线观看| 变态另类成人亚洲欧美熟女 | 亚洲精品美女久久久久99蜜臀| 欧美人与性动交α欧美精品济南到| 老司机福利观看| 精品国内亚洲2022精品成人| 麻豆久久精品国产亚洲av | 熟女少妇亚洲综合色aaa.| 日日摸夜夜添夜夜添小说| 午夜视频精品福利| 精品一品国产午夜福利视频| а√天堂www在线а√下载| 黑人猛操日本美女一级片| 亚洲欧美日韩另类电影网站| 国产一区二区在线av高清观看| 美女 人体艺术 gogo| www日本在线高清视频| 欧美黑人精品巨大| 两性夫妻黄色片| 午夜91福利影院| 免费少妇av软件| 亚洲国产欧美一区二区综合| 国产精品一区二区免费欧美| 天天躁夜夜躁狠狠躁躁| 久久人妻熟女aⅴ| 美女扒开内裤让男人捅视频| 身体一侧抽搐| 国产精品自产拍在线观看55亚洲| 成人三级做爰电影| 欧美成人午夜精品| 国产成人精品在线电影| 成熟少妇高潮喷水视频| 亚洲成人免费电影在线观看| 美女高潮到喷水免费观看| 久久人妻av系列| 黄网站色视频无遮挡免费观看| 亚洲精品美女久久av网站| 日本撒尿小便嘘嘘汇集6| 岛国在线观看网站| 国产精品久久久人人做人人爽| 99精品久久久久人妻精品| 热99国产精品久久久久久7| 精品乱码久久久久久99久播| 精品一区二区三区四区五区乱码| 日韩国内少妇激情av| 宅男免费午夜| 美女 人体艺术 gogo| 亚洲国产精品合色在线| 在线永久观看黄色视频| 欧美日韩福利视频一区二区| 男男h啪啪无遮挡| 好男人电影高清在线观看| 男女之事视频高清在线观看| 精品第一国产精品| 精品人妻1区二区| 视频在线观看一区二区三区| av天堂久久9| 国产无遮挡羞羞视频在线观看| 男女之事视频高清在线观看| 国产精品一区二区免费欧美| 人人妻,人人澡人人爽秒播| 18禁裸乳无遮挡免费网站照片 | 亚洲久久久国产精品| 亚洲免费av在线视频| 欧美国产精品va在线观看不卡| 精品人妻在线不人妻| 日本a在线网址| 最近最新中文字幕大全电影3 | 国产av又大| 中文字幕精品免费在线观看视频| 另类亚洲欧美激情| 日本vs欧美在线观看视频| 最新美女视频免费是黄的| 国产男靠女视频免费网站| 久久久久国内视频| 国产熟女xx| 日韩欧美三级三区| 日韩一卡2卡3卡4卡2021年| 国产亚洲av高清不卡| 色综合站精品国产| 老司机福利观看| 亚洲一区二区三区色噜噜 | 精品久久蜜臀av无| 色哟哟哟哟哟哟| 男女下面插进去视频免费观看| 国产精品一区二区免费欧美| 欧美日本亚洲视频在线播放| 国产成人影院久久av| 好看av亚洲va欧美ⅴa在| 99热只有精品国产| av片东京热男人的天堂| 欧美日韩av久久| 女性生殖器流出的白浆| 国产亚洲欧美精品永久| 欧美日韩黄片免| 夜夜爽天天搞| 久久国产亚洲av麻豆专区| 午夜a级毛片| 欧美日韩视频精品一区| 一本综合久久免费| av中文乱码字幕在线| 国产精品成人在线| 中文字幕精品免费在线观看视频| 十八禁网站免费在线| 亚洲免费av在线视频| 亚洲五月色婷婷综合| 成人免费观看视频高清| www国产在线视频色| 12—13女人毛片做爰片一| av超薄肉色丝袜交足视频| 黑人操中国人逼视频| 黄色女人牲交| 国产麻豆69| 一边摸一边做爽爽视频免费| 中文字幕人妻丝袜制服| 中文字幕最新亚洲高清| 亚洲美女黄片视频| 麻豆av在线久日| 老司机深夜福利视频在线观看| 国产一区二区三区综合在线观看| www.精华液| 午夜免费成人在线视频| 午夜精品在线福利| 国产乱人伦免费视频| 天天添夜夜摸| 免费观看精品视频网站| 另类亚洲欧美激情| 日本撒尿小便嘘嘘汇集6| 亚洲人成网站在线播放欧美日韩| av在线天堂中文字幕 | 亚洲精品国产区一区二| 日日干狠狠操夜夜爽| 黑人操中国人逼视频| 国产不卡一卡二| 亚洲一码二码三码区别大吗| 可以在线观看毛片的网站| 国产高清激情床上av| 操美女的视频在线观看| 亚洲成人免费电影在线观看| 亚洲欧美一区二区三区黑人| 脱女人内裤的视频| 18禁国产床啪视频网站| 亚洲成a人片在线一区二区| 纯流量卡能插随身wifi吗| 母亲3免费完整高清在线观看| 亚洲专区字幕在线| 女性生殖器流出的白浆| 久久伊人香网站| 999精品在线视频| 亚洲五月婷婷丁香| 99国产精品一区二区三区| 亚洲中文av在线| 丰满迷人的少妇在线观看| 亚洲片人在线观看| 首页视频小说图片口味搜索| 国产成人欧美| 亚洲欧美一区二区三区黑人| 麻豆一二三区av精品| 深夜精品福利| 国产av一区在线观看免费| 日本 av在线| 成人三级做爰电影| 国产精品影院久久| 欧美精品一区二区免费开放| 亚洲精品一卡2卡三卡4卡5卡| 国产伦人伦偷精品视频| 色播在线永久视频| 免费在线观看完整版高清| 精品久久久久久成人av| 涩涩av久久男人的天堂| 久99久视频精品免费| 热re99久久精品国产66热6| a级毛片在线看网站| 无限看片的www在线观看| 午夜福利免费观看在线| 久久久久久久久免费视频了| 久久亚洲精品不卡| 桃色一区二区三区在线观看| 天堂俺去俺来也www色官网| 亚洲欧美一区二区三区黑人| 看黄色毛片网站| 99国产精品免费福利视频| 欧美大码av| 日韩欧美一区二区三区在线观看| 欧美中文综合在线视频| 亚洲伊人色综图| 国产成人av教育| 国产高清视频在线播放一区| 男人的好看免费观看在线视频 | 欧美激情 高清一区二区三区| 人妻久久中文字幕网| 亚洲成人精品中文字幕电影 | 无限看片的www在线观看| 90打野战视频偷拍视频| 国产1区2区3区精品| 欧美日本中文国产一区发布| 久久久国产成人精品二区 | 97超级碰碰碰精品色视频在线观看| 老汉色∧v一级毛片| 久久精品亚洲av国产电影网| 亚洲精品一卡2卡三卡4卡5卡| а√天堂www在线а√下载| 看片在线看免费视频| 亚洲av日韩精品久久久久久密| 99精品在免费线老司机午夜| 在线观看免费视频网站a站| 国产精品 国内视频| 在线十欧美十亚洲十日本专区| 亚洲人成77777在线视频| 50天的宝宝边吃奶边哭怎么回事| 亚洲激情在线av| 亚洲精品国产色婷婷电影| 两个人免费观看高清视频| 亚洲avbb在线观看| 一级a爱片免费观看的视频| 亚洲色图 男人天堂 中文字幕| 不卡一级毛片| 视频在线观看一区二区三区| av电影中文网址| 亚洲性夜色夜夜综合| 亚洲中文日韩欧美视频| 久久久久久大精品| 91精品国产国语对白视频| av有码第一页| 欧美最黄视频在线播放免费 | 欧美日韩视频精品一区| 夜夜看夜夜爽夜夜摸 | 久久久久久久久免费视频了| 精品一品国产午夜福利视频| 久久草成人影院| 亚洲国产毛片av蜜桃av| aaaaa片日本免费| 欧美+亚洲+日韩+国产| 中文字幕另类日韩欧美亚洲嫩草| 久久精品国产综合久久久| 免费高清在线观看日韩| 中文字幕高清在线视频| 美女高潮到喷水免费观看| 真人做人爱边吃奶动态| 国产黄a三级三级三级人| 怎么达到女性高潮| 欧美激情 高清一区二区三区| 91大片在线观看| 国产欧美日韩一区二区三| 国产av一区在线观看免费| 他把我摸到了高潮在线观看| 天堂影院成人在线观看| 他把我摸到了高潮在线观看| 国产aⅴ精品一区二区三区波| 国产欧美日韩一区二区三| 法律面前人人平等表现在哪些方面| 亚洲片人在线观看| 黑人猛操日本美女一级片| av视频免费观看在线观看| 久久伊人香网站| 久久精品国产亚洲av香蕉五月| 免费av中文字幕在线| 午夜福利一区二区在线看| 色尼玛亚洲综合影院| 国产一卡二卡三卡精品| 一级片'在线观看视频| 欧美在线黄色| 欧美另类亚洲清纯唯美| av视频免费观看在线观看| 91大片在线观看| 欧美激情 高清一区二区三区| 国产日韩一区二区三区精品不卡| 日韩欧美三级三区| 午夜免费鲁丝| 国产午夜精品久久久久久| 中文字幕av电影在线播放| 亚洲人成77777在线视频| 一本大道久久a久久精品| 大型av网站在线播放| 亚洲av五月六月丁香网| 十八禁人妻一区二区| 妹子高潮喷水视频| 国产av又大| 91字幕亚洲| 桃色一区二区三区在线观看| 波多野结衣av一区二区av| 欧美成人免费av一区二区三区| 亚洲国产精品sss在线观看 | 亚洲专区国产一区二区| www.www免费av| 久久婷婷成人综合色麻豆| 视频区图区小说| 新久久久久国产一级毛片| 中文字幕色久视频| 男女下面进入的视频免费午夜 | av欧美777| 国产精品一区二区在线不卡| 欧美黑人欧美精品刺激| 窝窝影院91人妻| 国产av精品麻豆| 激情在线观看视频在线高清| 午夜福利在线观看吧| 久热爱精品视频在线9| 香蕉国产在线看| 99久久久亚洲精品蜜臀av| 很黄的视频免费| 国产真人三级小视频在线观看| 多毛熟女@视频| 亚洲男人天堂网一区| 精品高清国产在线一区| 性色av乱码一区二区三区2| 国产精品98久久久久久宅男小说| 亚洲五月婷婷丁香| 大型黄色视频在线免费观看| 在线观看免费日韩欧美大片| 亚洲av日韩精品久久久久久密| 搡老熟女国产l中国老女人| a级毛片在线看网站| 男人舔女人的私密视频| 精品久久蜜臀av无| 在线永久观看黄色视频| 91老司机精品| 女生性感内裤真人,穿戴方法视频| 美女福利国产在线| 大型av网站在线播放| 超碰97精品在线观看| 国产亚洲欧美98| 巨乳人妻的诱惑在线观看| 亚洲欧洲精品一区二区精品久久久| 国产男靠女视频免费网站| 久久精品人人爽人人爽视色| 好看av亚洲va欧美ⅴa在| 最近最新免费中文字幕在线| 99久久久亚洲精品蜜臀av| 国产成人欧美在线观看| 日本 av在线| 亚洲精品久久午夜乱码| 在线观看免费高清a一片| 国产免费av片在线观看野外av| 精品第一国产精品| 国产成人精品在线电影| 久久久久国产精品人妻aⅴ院| 欧美日韩亚洲综合一区二区三区_| 亚洲人成电影免费在线| 免费看a级黄色片| 午夜成年电影在线免费观看| 日本精品一区二区三区蜜桃| 精品卡一卡二卡四卡免费| 夜夜爽天天搞| 免费在线观看影片大全网站| 操美女的视频在线观看| 国产又色又爽无遮挡免费看| 人人澡人人妻人| 亚洲性夜色夜夜综合| av免费在线观看网站| 久久性视频一级片| 神马国产精品三级电影在线观看 | 久久久久亚洲av毛片大全| 国产精品成人在线|