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

    Dynamic stability of quadruped robot walking on slope with trot gait①

    2016-12-06 02:39:52LeiJingtao雷靜桃RenMingming
    High Technology Letters 2016年1期

    Lei Jingtao(雷靜桃), Ren Mingming

    (School of Mechatronic Engineering and Automation, Shanghai University, Shanghai 200072, P.R.China)

    ?

    Dynamic stability of quadruped robot walking on slope with trot gait①

    Lei Jingtao(雷靜桃)②, Ren Mingming

    (School of Mechatronic Engineering and Automation, Shanghai University, Shanghai 200072, P.R.China)

    The dynamic stability of a quadruped robot trotting on slope was analyzed. Compared with crawl gait, trot gait can improve walking speed of quadruped robots. When a quadruped robot trots, each leg is in the alternate state of swing phase or supporting phase, and two legs in the diagonal line are in the same phase. The feet in the supporting phase form a supporting region on the ground. When a quadruped robot walks on slope, the vertical distance from zero moment point (ZMP) to the supporting diagonal line is defined as ZMP offset distance. Whether this distance is less than the maximum offset distance or not, the stability of robot trotting on slope can be judged. The foot trajectory was planned with the sinusoidal function.Based on the kinematic analysis, the ZMP offset distance of quadruped robot under different slope angles, step length and step height was calculated, then the reasonable slope angle, step length and step height for quadruped robot trotting on slope to keep dynamic stability can be determined. On the other hand, the posture angle of quadruped robot should be controlled within the desired range. Computer simulations were executed to verify the theoretical analysis.The study will provide reference for determining reasonable step parameters of the quadruped robot.

    quadruped robot, slope, trot gait, stability, zero moment point (ZMP)

    0 Introduction

    Mobile robots are generally divided into legged, wheeled, tracked and mixed. The quadruped robots have been popularly studied by many research groups[1]because of its advantages over wheeled or tracked locomotion, such as obstacle striding, high terrain adaptability, environmental protection and energy saving. Of course the legged robot has some disadvantages such as more complex structure, which keeps it from being used.

    When walking in an unstructured environment, compared to wheeled or tracked robots, legged robot has some advantages. Four-legged robot has the performance of more stable than biped robot, simple structure than the six-legged robot.

    The locomotion style of animals changes with different environments, which differ from the uneven terrain such as gravel road and mountainous terrain, specialterrain such as stairs and sand beaches, extreme terrain such as snow or swamp.

    It is difficult for a quadruped robot to walk as an animal does. Quadruped robots can not achieve the mobility and dynamic stability similar to four-legged creatures, since a quadruped robot is a multi DOF system. During walking, the robot has the characteristics of changing configuration. On the other hand, the terrain of the unstructured environment is complex, such as, slope, obstacle and so on. The ground reaction forces, sliding and friction and other factors influence stability. It is difficult to achieve dynamic stability for quadruped robots walking in the unstructured environment.

    The stability of a walking robot is divided into static and dynamic according to different gait style. The trot gait belongs to the dynamic one. When a quadruped robot walks with crawl gait, it is easy to keep stable, but the walking speed is lower. As developed recently, a robot with trot gait, can improve effectively the walking speed. It is necessary to analyze its dynamic stability in unstructured environment.

    Some researchers have analyzed the dynamic stability of the quadruped robot with trot gait. The ZMP stability criterion is generally used to analyze the stability of walking robot[2,3]. Ref.[4] proposed a control method for quadruped robot with trot gait pass through the plain and slope transition. Ref.[5] presented that a quadruped robot could perform self-stable running behavior by designing mechanical structure to expand the self-stable region. Ref.[6] presented the centroid trajectory of quadruped robot walking with crawl gait, and discussed on a strategy to maintain the best walking balance. Ref.[7] proposed a kind of walking pattern generator with the selection of every foot placement such as the COG of the robot follows a stable trajectory characterized by a stability margin relative to the current support triangle. Ref.[8] presented a kind of trajectory planning for quadruped robot, which causes the robot to have versatile motion in different directions. The ZMP criterion method is used to analyze dynamic stability. The minimum and maximum of lateral movement of the body is specified to maintain location of the ZMP in the supporting region. Ref.[9] analyzed the stability of the quadruped robot by simulation, and presented that the robot should have proper leg configuration to increase stability by decreasing pitching motion of the robot. Ref.[10] discussed omin-directional walking of the quadruped robot walking on slope. A kind of successive gait-transition method is presented, and the stability is discussed by a measure of stability margin. Ref.[11] presented a center of pressure method to evaluate the stability of the quadruped robot with crawl gait.Turning gait is the most general and important factor for omni-directional walking of walking robot. Ref.[12] presented the turning gait planning for six-legged robot, and stability is analyzed based on the normalized energy stability margin method. Ref.[13] proposed the omni-directional static walking method on a slope, and the stability in gait transition processes was analyzed. The stability analysis of quadruped robot with trot gait is mostly focused on motion over a flat terrainor walking along straight line.Study on the dynamic stability of the quadruped robot walking in the unstructured environment is limited.

    In this study, the dynamic stability of the quadruped robot trotting on slope was analyzed. Firstly, the foot trajectory of the quadruped robot is planned with sinusoidal function. When the quadruped robot walks with trot gait, the vertical distance from ZMP to the supporting diagonal line is defined as the ZMP offset distance, which is used as the dynamic stability criterion.When the quadruped robot walks in different slope conditions, the ZMP offset distance is calculated, respectively. On the one hand, the reasonable gait parameters and slope angle for quadruped robot dynamic stable walking on slope with trot gait can be determined according to the ZMP offset distance. On the other hand, the posture angle of quadruped robot should be controlled to ensure the dynamic walking stability.

    1 Trot gait

    Quadruped robot is composed of one body and four legs, as shown in Fig.1. Each leg is designed with three rotational joints, which are side-swing-hip joint, forward-swing- hip joint and knee joint.

    Fig.1 Quadruped robot

    The gait style of the quadruped robot are crawl, trot, pace and jump. Trot gait is faster but is more unstable than crawl gait. When the quadruped robot walks with trot gait, four legs are classified into two groups, two legs on the diagonal line belong to one group. Two legs that belong to the same group are lifting, swing, falling, supporting synchronously. The leg mechanism is alternately in swing phase and supporting phase.Two legs transfer forward, while the other legs support the body weight and push itself against the ground to generate the needed thrust in the desired direction of motion.

    Two parameters can describe the trot gait, which are duty factor βiand leg phase φi(i=1,2,3,4). The gait formula g=[β1, β2, β3, β4; φ1, φ2, φ3, φ4] is generally used to describe lift and fall sequence of each leg and the interval time between adjacent legs.

    If the periodic gait is adopted and the duty factor is 0.5, the time in supporting phase and swing phase of each leg is equal, as shown in Fig.2.

    Then the gait formula is

    (1)

    The ZMP criterion method is used to analyze dynamic stability of bipedal or quadruped robot.

    Fig.2 Different phase

    Definition:

    (1) The supporting polygon: the region formed by the feet in supporting phase, as shown in Fig.3.

    (2) ZMP offset distance d: the vertical distance from ZMP to the supporting diagonal line.

    (3) Maximum offset distance Δ: the vertical distance from the supporting diagonal line to the polygon boundary.

    The supporting region is shown in Fig.3, which is determined by the foot sizee. During the former half gait cycle, the 2th leg and 4th leg are in supporting phase. During the latter half gait cycle, the 1th leg and 3th leg are in supporting phase, the stable offset distance d is from ZMP to supporting diagonal line.

    Fig.3 The supporting polygon and offset distance

    The stability judgment condition is: If d<Δ, ZMP locate inside the supporting polygon, robot can keep stable. If d>Δ, ZMP locate outside the supporting polygon, robot is unstable. So the dynamic stability of the quadruped robot walks on slope with trot gait should meet the condition: ZMP offset distance d less than maximum offset distance Δ.

    On the one hand, the robot walk on slop, the ZMP should be controlled within the supporting region, which is related to the gait parameters. On the other hand, the posture angle of robot should be controlled.

    2 Kinematics

    2.1 Trajectory planning

    The trajectory function should be continuous and smooth. The combined cycloid function can effectively reduce the contact force between feet and environment[14]. The cycloid function, sinusoidal function or linear function is generally adopted as the foot trajectory, and different trajectories will provide different dynamic performance for quadruped robot.

    The trajectory planning was done in Cartesian coordinates. As the side swing movement of the leg mechanism needn’t be considered,so the trajectory function was divided into two directions. The functions along thex-axis andz-axis direction were planned by the sinusoidal function. According to the boundary constraint conditions:

    px|t=0= 0, px|t=Tsw=S, pz|t=0=0,

    pz|t=Tsw=H

    The trajectory function can be determined as

    (2)

    Where S is step length,H is step height, px, py, pzis the foot displacement along coordinate axis direction, respectively.

    The Denavit-Hartenberg coordinate method was adopted to analyze the kinematics of quadruped robot trotting on the slope, the spatial position and pose relationship between the foot and leg joints can be derived. Because of each leg is alternately in swing phase and supporting phase, so the kinematic of leg in different phase is different.

    2.2 Forward kinematics of leg in swing phase

    Quadruped robot walks on slope, the slope angle is set as α. The coordinate systems of leg in swing phase are shown in Fig.4. The body coordinate system is {B}. The coordinate system {0} and {1} are located at joint 1, whose position coordinates with respect to the body coordinate system {B} is (a,b,c). The foot coordinate system is {4}.

    Fig.4 D-H coordinate systems

    According to the coordinate systems, theD-Hparameters of leg in swing phase are shown in Table 1.

    Table 1 D-H parameters

    The kinematic of leg in swing phase can be described by the transformation matrix:

    (3)

    where sk=sinθk, ck=cosθk, s23=sin(θ2+θ3), c23=cos(θ2+θ3)

    2.3 Forward kinematics of leg in supporting phase

    As the side-swing motion of leg mechanism needn’t to be analyzed, so θ1is always zero. The coordinate systems of leg in supporting phase were shown in Fig.5. θ2, θ3, θ4are hip, knee and ankle joint variable, respectively. The ziaxis direction of the coordinate systems oixiyizi(k=0,1,2,3) is from outside to inside.

    Fig.5 D-H coordinate systems

    According to theD-Hcoordinate systems, theD-Hparameters of the leg mechanism in supporting phase are shown in Table 2.

    Table 2 D-H parameters

    The kinematic of leg in supporting phase can be described by the transformation matrix:

    (4)

    2.4 Inverse kinematics

    The joint angles θ1,θ2,θ3during supporting phase and swing phase can be determined by the inverse transformation method.

    (5)

    where, D=s1py-c1pz+cc1-bs1

    The angular velocity can be calculated by the derivative, which will be used as joint driving functionsfor simulation analysis.

    3 ZMP of robot dynamic walking

    According to the dynamic stability criterion: ZMP offset distance less than the maximum offset distance, the stability of robot walking on slope can be analyzed. The equation of the supporting diagonal line should be first determined.

    3.1 Equation of the supporting diagonal line

    In order to determine the ZMP offset distanced, coordinates of the supporting feet respect to the body coordinate system {B} should be calculated firstly, and the supporting diagonal line equation can be derived. Then the ZMP coordinates in the body coordinate system are derived. The vertical distance for ZMP to the supporting diagonal line is the offset distanced.

    During the former half gait cycle and the latter half gait, the supporting diagonal line is different. The equation of the supporting diagonal line is

    (6)

    3.2 ZMP coordinates

    The quadruped robot walks on the slope with trot gait, the ZMP coordinates should be firstly calculated according to the definition. ZMP is the point that the moment about this point of all force acting on the robot is equal to zero. These forces include the gravity force, external force and inertia force.

    The moment about the ZMP of all force acting on the robot is zero.

    MZMP=0

    (7)

    Take the 4th leg mechanism for example, the force analysis is shown in Fig.6. Fb, F41, F42, F43are the forces acting on the center of mass of body, link1, link2 and link3 of the 4th leg, respectively. rb, r41, r42, r43are the radius vectors from ZMP to the center of mass of body, link1, link2 and link3 of the 4th leg, respectively.

    Fig.6 Force analysis

    Then the moment of all forces acting on the robot about ZMP is

    =0

    (8)

    Set the ZMP coordinates with respect to the body coordinate system {B} is ZMP (xzmp, yzmp). Then Eq.(8) can be expressed as:

    So the ZMP coordinates can be calculated as:

    3.3 ZMP offset distance

    According to the ZMP coordinates and the equation of the supporting diagonal line, the ZMP offset distance could be calculated. According the vertical distance formula of from point (x, y) to line Ax+By+C=0.

    So the ZMP offset distancedcan be derived as:

    (9)

    4 The posture and posture control

    ZMP is controlled within the supporting region, which couldn’t ensure the robot stable. The posture of robot should be also controlled. The posture angle refers to the angles of body coordinate system respect to reference coordinate system, which is denoted by δ, β and γ, as shown in Fig.6. Then the rotate matrix of robot with respect to the reference coordinate system is

    Rxyz(γ, β,δ)=Rot(zB0,δ)Rot(yB0, β)Rot(xB0,γ)

    (10)

    Quadruped robot walks on slope, and the robot posture should be considered. For designing the controller of walking robot, the trajectory and the posture control should be controlled. On the one hand, the desired trajectory related to joint kinematic parameters should be controlled to ensure the ZMP locating within the supporting region. On the other hand, the gyroscope fixed on the body is used to real-time test the posture angle during the quadruped robot walking, and the posture angle should be controlled to ensure the robot to keep the optimal posture.

    5 The influence factors on stability

    The slope angle and the gait parameters have influence on the dynamic stability of quadruped robot.The relationship between the stability and the slop angle, step length and step height will be analyzed. The parameters of the quadruped robot are shown in Table 3.

    According to the foot sizee, the maximum offset distance is 45mm. When ZMP offset distance is less than 45mm, which shows that the ZMP locates within the supporting polygon, and the robot can keep stable.

    Table 3 The parameters of robot

    According to the structure parameters of the robot, the maximum gait parameters of the robot can be determined. The maximum step lengthSmax is 245mm, and the maximum step heightHmax is 120mm.

    5.1 Step length

    The quadruped robot walks on the slop with trot gait. Set the slop angle to 20°, and the step heightHat 50mm. Step lengthSincreases from 100mm to 140mm, ZMP offset distancedchanges with step lengthSduring one step cycle, as shown in Fig.7.

    The stability of the quadruped robot decreases with the increase of step length. When the step length is 136mm, the ZMP offset distance is 44mm, which is less than the maximum offset distance. The robot can keep stable.

    When the slope angle is 20°, and the step height is 50mm, the step length should be less than 136mm to keep the robot stable.

    Fig.7 d change with the step length

    5.2 Step height

    The quadruped robot walks on the slop with trot gait. Set the slop angle to 20°, and step lengthSat 100mm. When step heightHincreases from 40mm to 120mm, offset distancedchanges with step heightHduring one step cycle, as shown in Fig.8.

    Fig.8 d change with step height

    (1) The step length keeps invariable, the stability of the quadruped robot decreases with the increase of the step height.

    (2) Compared with the step length, the step height has less influence on the stability.

    5.3 Slope angle

    The quadruped robot walks on the slop with trot gait. Set the step height at 50mm, and the step lengthSat 100mm. When the slope angle increases from 15° to 40°, ZMP offset distancedchanges with the slope angle during one step cycle, as shown in Fig.9.

    The quadruped robot walks on the slope with trot gait, the stability of the quadruped robot decreases with the slope angle increase. When the slope angle is 30°, the ZMP offset distance is 44mm, which is less than the maximum offset distance 45mm. The robot can keep stable. Otherwise, the robot is unstable.

    Fig.9 d change with slop angle

    So the quadruped robot walks on the slope with trot gait, the step height is 50mm, and the step length is 100mm, the slope angle should be less than 30°.

    It is concluded that:

    (1) The step length keeps invariable, the stability of the quadruped robot decreases with the increase of the step height. But compared with the step length, step height has less influence on the stability.

    (2) The step height keeps invariable, the stability of the quadruped robot decreases with the step length. If the slope angle is 20°, and the step height is 50mm, the step length should be less than 136mm to keep the robot stable.

    (3) If the step length is 100mm, and the step height is 50mm, the slope angle should be less than 30° to keep the robot stable.

    6 Simulation

    The simulation of the quadruped robot walking on slope with trot gait was performed in software ADAMs.

    The virtual prototype model was established in software Solid-works, and inputted into the ADAMs. Each joint was added the angular velocity driving function, which is determined by the inverse kinematics. The contact constraint was added between the feet in supporting phase and environment.

    (11)

    where k is stiffness coefficient depending on the material properties,k=1000kN/mpk, δ is penetration depth of the foot into the ground, pkis exponent of the penetration for the normal elastic force, pk=2.2. and pdis exponent of the penetration for the normal damping contact force, pd=0.5.

    6.1 Simulation in ADAMS software

    Simulations were performed with different gait parameters and different slope angles. The simulation sequence chart of quadruped robot walking on slop is shown in Fig.10.

    Fig.10 The sequence chart

    6.2 Simulation results

    Set the slop angle at 10°, step length at 100mm and step height at 60mm, gait cycle at 10s.The simulation results of foot trajectory and COG were obtained, as shown in Fig.11.

    (a) Foot displacement of 1th and 3th leg along z direction

    (b) Foot displacement of 2th and 4th leg along z direction

    (c) COG displacement

    7 Conclusions

    In this study, dynamic stability of the quadruped robot walking on slope with trot gait was analyzed. The sinusoidal function was planned as the foot trajectory function. The vertical distance from ZMP to the supporting diagonal line was defined as the ZMP offset distanced. The vertical distance between the supporting diagonal line and the boundary of the supporting polygon was defined as the maximum offset distance Δ. Whetherdis less than Δ or not, it was used as the dynamic stability criterion. The slope angle and the gait parameters influence the stability. When the slope angle keeps invariable, the stability of the quadruped robot decreases with the step length and step height increase. The step length influences the stability more than that of step height. When the gait parameters keep invariable, the stability decreases with the slope angle increase. Finally, simulations were performed to verify the theoretical analysis.

    The experiment of robot walking on slope will be conducted in the future. The optimal posture angle will be determined by the experiment, and the control algorithm for dynamic stable walking on slop will be evaluated.

    [ 1] Arikawa K, Hirose S. Development of quadruped walking robot TITAN-VIII. In: Proceedings of the IEEE/RSJ International Conference on Intelligent Robots and Systems, Osaka , Japan,1996. 208-214

    [ 2] Mrinal K, Jonas B, Peter P, et al. Fast, robust quadruped locomotion over challenging terrain. In: Proceedings of the IEEE International Conference on Robotics and Automation Anchorage Convention District, Anchorage, USA, 2010. 2665-2670

    [ 3] Zhang S, Gao J, Duan X, et al. Trot pattern generation for quadruped robot based on the ZMP stability margin. In: Proceedings of the International Conference on Complex Medical Engineering, Beijing, China, 2013. 608-613

    [ 4] Lee J H, Park J H. Control for quadruped robots in trotting on horizontal and slanted surfaces. In: Proceedings of the IEEE 9th Asian Control Conference (ASCC), Istanbul, Turkey, 2013. 1-6

    [ 5] Chatzakos P, Papadopoulos E. Dynamically running quadrupeds self-stable region expansion by mechanical design. In: Proceedings of the IEEE International Conference on Robotics and Automation, Kobe, Japan. 2009. 2749-2754

    [ 6] Kim B H. Centroid-based analysis of quadruped-robot walking balance. In: Proceedings of the International Conference on Advanced Robotics, Munich, Germany, 2009. 1-6

    [ 7] Pongas D, Mistry M, Schaal S. A robust quadruped walking gait for traversing rough terrain. In: Proceedings of the IEEE International Conference on Robotics and Automation, Roma, Italy, 2007. 1474-1479

    [ 8] Reza Y, Vahid J M, Reza O. Dynamically stable trajectory planning for a quadruped robot. In: Proceedings of the 20th Iranian Conference on Electrical Engineering, Tehran, Iran, 2012. 845-850

    [ 9] Meek S, Kim J, Anderson M. Stability of a trotting quadruped robot with passive, underactuated legs. In: Proceedings of the IEEE International Conference on Robotics and Automation (ICRA), Pasadena, USA, 2008. 347-351

    [10] Zhang L, Ma S G, Honda Y, et al. Omni-directional walking of a quadruped robot on a slope. In: Proceedings of the IEEE International Conference on Robotics and Biomimetics, Shenyang, China, 2004. 268-273

    [11] Wang P, Huang B, Sun L. Stability judging method for quadruped bionic robot.JournalofharbinInstituteoftechnology,2008, 4(7):1063-1066

    [12] Shibendu S R, Dilip K P. Effects of turning gait parameters on energy consumption and stability of a six-legged walking robot.RoboticsandAutonomousSystems, 2012, 60(1): 72-82

    [13] Zhang L, Wang L L, Ma S G. Research of a quadruped robot walking on a slope based on 4-legsupporting period. In: Proceedings of the IEEE International Conference on Robotics and Biomimetics, Tianjin, China, 2010. 1346-1351

    [14] Lei J T, Wang F, Yu H Y, et al. Energy efficiency analysis of quadruped robot with trot gait and combined cycloid foot trajectory.ChineseJournalofMechanicalEngineering, 2014, 27(1):138-145

    Lei Jingtao, born in 1970. She received her Ph.D degrees from Beihang University in 2007. She also received her B.S. and M.S. degrees from Henan University of Science and Technology in 1991 and 1996 respectively. Her research interests include the robot mechanisms and robot modular technology.

    10.3772/j.issn.1006-6748.2016.01.001

    ① Supported by the National Natural Science Foundation of China (No. 51375289), Shanghai Municipal National Natural Science Foundation of China (No.13ZR1415500) and Innovation Fund of Shanghai Education Commission (No.13YZ020).

    ② To whom correspondence should be addressed. E-mail: jtlei2000@163.comReceived on Jan. 22, 2015

    成人三级做爰电影| 免费高清在线观看视频在线观看| 国产精品99久久99久久久不卡| 脱女人内裤的视频| 亚洲avbb在线观看| 国产av又大| 久久狼人影院| 少妇猛男粗大的猛烈进出视频| 亚洲激情五月婷婷啪啪| 国产亚洲欧美在线一区二区| 国产亚洲一区二区精品| 日韩中文字幕视频在线看片| 国内毛片毛片毛片毛片毛片| 日本五十路高清| 少妇粗大呻吟视频| 黄色视频在线播放观看不卡| 天天躁日日躁夜夜躁夜夜| 99久久99久久久精品蜜桃| 亚洲国产日韩一区二区| 成年美女黄网站色视频大全免费| 国产精品香港三级国产av潘金莲| 国产精品一二三区在线看| 在线观看www视频免费| 亚洲一区中文字幕在线| 国产精品久久久久久人妻精品电影 | 日韩 亚洲 欧美在线| 五月天丁香电影| 国产成人精品久久二区二区免费| 天堂8中文在线网| 欧美日韩av久久| 男人爽女人下面视频在线观看| 夫妻午夜视频| 夫妻午夜视频| 色视频在线一区二区三区| 亚洲中文日韩欧美视频| 人妻人人澡人人爽人人| 欧美精品亚洲一区二区| 欧美+亚洲+日韩+国产| 九色亚洲精品在线播放| 亚洲五月色婷婷综合| 久久精品国产综合久久久| 制服人妻中文乱码| 亚洲av日韩在线播放| 国产精品九九99| 天天添夜夜摸| 免费观看av网站的网址| 久久99一区二区三区| 欧美在线一区亚洲| 男女免费视频国产| 男人操女人黄网站| 欧美午夜高清在线| www.999成人在线观看| 国产在线一区二区三区精| 制服诱惑二区| 最近最新免费中文字幕在线| av片东京热男人的天堂| 青青草视频在线视频观看| 国产在视频线精品| 每晚都被弄得嗷嗷叫到高潮| 日韩制服丝袜自拍偷拍| 男女无遮挡免费网站观看| 久久女婷五月综合色啪小说| 麻豆国产av国片精品| 国产精品亚洲av一区麻豆| 国产片内射在线| 自拍欧美九色日韩亚洲蝌蚪91| 欧美变态另类bdsm刘玥| 一区福利在线观看| 一级黄色大片毛片| 午夜精品国产一区二区电影| 欧美av亚洲av综合av国产av| www.自偷自拍.com| 91成人精品电影| 日韩精品免费视频一区二区三区| 建设人人有责人人尽责人人享有的| 成人黄色视频免费在线看| 天堂中文最新版在线下载| 亚洲成人手机| 久久精品亚洲av国产电影网| 国产成人精品在线电影| 精品少妇内射三级| 精品卡一卡二卡四卡免费| 免费在线观看完整版高清| 高清欧美精品videossex| 亚洲自偷自拍图片 自拍| 精品少妇内射三级| 国产精品国产三级国产专区5o| 乱人伦中国视频| 久久中文看片网| 日本黄色日本黄色录像| 啦啦啦 在线观看视频| 国产精品九九99| 日本精品一区二区三区蜜桃| 久久青草综合色| 黑人猛操日本美女一级片| 国产色视频综合| 色婷婷久久久亚洲欧美| 免费久久久久久久精品成人欧美视频| 最新在线观看一区二区三区| 欧美黄色淫秽网站| 久久影院123| 男女国产视频网站| xxxhd国产人妻xxx| 亚洲专区字幕在线| 国产一级毛片在线| 国产xxxxx性猛交| 黄色视频不卡| 中文字幕另类日韩欧美亚洲嫩草| 日韩大码丰满熟妇| 一本色道久久久久久精品综合| 国产成人av激情在线播放| 又黄又粗又硬又大视频| 欧美日韩黄片免| 日韩欧美一区视频在线观看| 一个人免费在线观看的高清视频 | 亚洲激情五月婷婷啪啪| 各种免费的搞黄视频| 啦啦啦视频在线资源免费观看| 国产精品1区2区在线观看. | 在线观看人妻少妇| 成人手机av| 免费高清在线观看视频在线观看| 18禁裸乳无遮挡动漫免费视频| 99久久综合免费| 老熟妇乱子伦视频在线观看 | a级片在线免费高清观看视频| av有码第一页| 亚洲av日韩精品久久久久久密| 后天国语完整版免费观看| 日韩大码丰满熟妇| 亚洲一区二区三区欧美精品| 最近最新中文字幕大全免费视频| 免费女性裸体啪啪无遮挡网站| 欧美在线一区亚洲| 久久久精品国产亚洲av高清涩受| 大码成人一级视频| 亚洲精品国产色婷婷电影| 无限看片的www在线观看| 天堂俺去俺来也www色官网| 日本五十路高清| 男女免费视频国产| 伦理电影免费视频| 老汉色∧v一级毛片| 日韩熟女老妇一区二区性免费视频| 亚洲免费av在线视频| 国产精品熟女久久久久浪| 午夜免费观看性视频| 青草久久国产| 久久国产精品人妻蜜桃| 亚洲国产成人一精品久久久| 国产av精品麻豆| 大片电影免费在线观看免费| 婷婷成人精品国产| av有码第一页| 各种免费的搞黄视频| 亚洲欧美精品综合一区二区三区| 欧美老熟妇乱子伦牲交| 别揉我奶头~嗯~啊~动态视频 | cao死你这个sao货| av天堂久久9| 999久久久国产精品视频| 免费女性裸体啪啪无遮挡网站| 欧美国产精品一级二级三级| 999精品在线视频| 国产男女超爽视频在线观看| 国产精品免费大片| 久久久久久免费高清国产稀缺| 视频区欧美日本亚洲| 亚洲国产欧美在线一区| 日本精品一区二区三区蜜桃| 99久久精品国产亚洲精品| 国产精品久久久久久精品电影小说| 亚洲av美国av| bbb黄色大片| 老司机深夜福利视频在线观看 | 国产1区2区3区精品| 中文字幕制服av| 免费观看av网站的网址| 国产精品一区二区在线不卡| www.自偷自拍.com| 大型av网站在线播放| 久久综合国产亚洲精品| 少妇粗大呻吟视频| 欧美xxⅹ黑人| 91大片在线观看| 少妇粗大呻吟视频| 狠狠狠狠99中文字幕| 精品第一国产精品| 久久精品国产a三级三级三级| 在线看a的网站| 国产成人av激情在线播放| 新久久久久国产一级毛片| 久久青草综合色| 狠狠精品人妻久久久久久综合| 91av网站免费观看| 久久影院123| 日本五十路高清| 亚洲一区二区三区欧美精品| 国产1区2区3区精品| 久久狼人影院| 中文字幕人妻熟女乱码| 色婷婷久久久亚洲欧美| 精品少妇黑人巨大在线播放| 爱豆传媒免费全集在线观看| 天堂俺去俺来也www色官网| 大陆偷拍与自拍| 制服诱惑二区| 超碰成人久久| 亚洲成av片中文字幕在线观看| 欧美黄色淫秽网站| 在线观看www视频免费| 欧美日韩亚洲综合一区二区三区_| √禁漫天堂资源中文www| 欧美日韩亚洲高清精品| 国产在线免费精品| 法律面前人人平等表现在哪些方面 | 国产精品麻豆人妻色哟哟久久| 亚洲精品日韩在线中文字幕| 精品国产国语对白av| 久久毛片免费看一区二区三区| 国产成人一区二区三区免费视频网站| 国产精品久久久久久精品电影小说| 亚洲国产日韩一区二区| 2018国产大陆天天弄谢| 飞空精品影院首页| 久久精品国产亚洲av高清一级| 精品国产乱子伦一区二区三区 | 久久人人爽av亚洲精品天堂| 一个人免费看片子| 欧美日韩视频精品一区| 国产老妇伦熟女老妇高清| 亚洲精品美女久久av网站| 久久久国产一区二区| 久久久久久久大尺度免费视频| 国产成人精品在线电影| 亚洲精品国产色婷婷电影| 午夜免费成人在线视频| 日韩大码丰满熟妇| 无限看片的www在线观看| 高清视频免费观看一区二区| 老司机福利观看| 狂野欧美激情性bbbbbb| 国产男女超爽视频在线观看| 成人18禁高潮啪啪吃奶动态图| 午夜福利视频在线观看免费| 后天国语完整版免费观看| 亚洲一区中文字幕在线| 亚洲av日韩精品久久久久久密| 精品一区在线观看国产| 91成人精品电影| tocl精华| 欧美日韩av久久| 一个人免费在线观看的高清视频 | 美女中出高潮动态图| 在线观看免费午夜福利视频| 欧美日韩av久久| 国产亚洲欧美在线一区二区| 亚洲精品国产色婷婷电影| 18禁观看日本| 日韩视频一区二区在线观看| 欧美另类一区| 精品亚洲成国产av| 欧美日韩一级在线毛片| 一级a爱视频在线免费观看| 国产精品一区二区在线观看99| av天堂久久9| 搡老岳熟女国产| 女警被强在线播放| 国产亚洲一区二区精品| 欧美在线一区亚洲| 欧美乱码精品一区二区三区| 欧美日本中文国产一区发布| 欧美国产精品一级二级三级| 一个人免费看片子| 国产欧美日韩综合在线一区二区| 精品少妇一区二区三区视频日本电影| 色94色欧美一区二区| 精品第一国产精品| 一区福利在线观看| 桃花免费在线播放| 成人av一区二区三区在线看 | 天天躁夜夜躁狠狠躁躁| 亚洲五月婷婷丁香| 精品福利观看| 黄频高清免费视频| 亚洲中文日韩欧美视频| www.999成人在线观看| 欧美另类一区| 国产精品久久久av美女十八| 亚洲欧美成人综合另类久久久| 国产视频一区二区在线看| 亚洲精品中文字幕一二三四区 | videos熟女内射| 99香蕉大伊视频| 亚洲精品第二区| 久久精品国产亚洲av高清一级| 亚洲国产欧美网| 捣出白浆h1v1| av天堂久久9| avwww免费| 亚洲精品乱久久久久久| 18禁观看日本| 亚洲色图 男人天堂 中文字幕| 精品福利观看| 亚洲欧美一区二区三区黑人| 99九九在线精品视频| 19禁男女啪啪无遮挡网站| 国产精品一区二区在线不卡| 精品一品国产午夜福利视频| 亚洲成国产人片在线观看| 999久久久国产精品视频| 日本av免费视频播放| 亚洲精品第二区| 国产成人精品无人区| 天天躁日日躁夜夜躁夜夜| 女性生殖器流出的白浆| 色视频在线一区二区三区| 在线永久观看黄色视频| 久久久久国产一级毛片高清牌| 老熟女久久久| 国产极品粉嫩免费观看在线| 免费在线观看影片大全网站| 国产色视频综合| 国产亚洲av高清不卡| 69av精品久久久久久 | 老司机午夜福利在线观看视频 | 成年动漫av网址| 一本色道久久久久久精品综合| 国产在线一区二区三区精| 欧美人与性动交α欧美精品济南到| 国产日韩欧美视频二区| 国产精品 国内视频| 久久久久网色| 亚洲国产精品一区二区三区在线| 国产成+人综合+亚洲专区| 久久女婷五月综合色啪小说| 国产免费视频播放在线视频| 亚洲av美国av| 国产日韩一区二区三区精品不卡| 蜜桃在线观看..| 多毛熟女@视频| 男人添女人高潮全过程视频| 国产成人欧美在线观看 | 国产高清视频在线播放一区 | 午夜视频精品福利| 男人操女人黄网站| 亚洲精品中文字幕在线视频| 欧美精品一区二区免费开放| 精品国产超薄肉色丝袜足j| 这个男人来自地球电影免费观看| 亚洲国产精品成人久久小说| 免费观看av网站的网址| 一本大道久久a久久精品| 国产免费现黄频在线看| 精品高清国产在线一区| 岛国毛片在线播放| 国产av国产精品国产| 99国产精品免费福利视频| 中国国产av一级| 在线观看免费视频网站a站| 国产亚洲欧美精品永久| 搡老岳熟女国产| av又黄又爽大尺度在线免费看| 亚洲欧美色中文字幕在线| 少妇的丰满在线观看| 午夜福利视频在线观看免费| 亚洲成国产人片在线观看| 久久久国产欧美日韩av| 69精品国产乱码久久久| 国产高清视频在线播放一区 | 国产又爽黄色视频| 亚洲综合色网址| 男人添女人高潮全过程视频| 美女扒开内裤让男人捅视频| 国产主播在线观看一区二区| 国产精品免费大片| 久久久精品免费免费高清| 欧美另类一区| 国产91精品成人一区二区三区 | 另类精品久久| 香蕉国产在线看| 妹子高潮喷水视频| 国产熟女午夜一区二区三区| 国产高清视频在线播放一区 | 男人添女人高潮全过程视频| 免费久久久久久久精品成人欧美视频| 两个人免费观看高清视频| av一本久久久久| 欧美午夜高清在线| 丰满迷人的少妇在线观看| 深夜精品福利| 午夜免费鲁丝| 欧美日本中文国产一区发布| 国产1区2区3区精品| 色94色欧美一区二区| 99九九在线精品视频| 国产有黄有色有爽视频| 国产精品熟女久久久久浪| 大型av网站在线播放| 国产精品一区二区精品视频观看| 日韩大码丰满熟妇| 不卡一级毛片| 最黄视频免费看| 亚洲专区中文字幕在线| 精品一区二区三区四区五区乱码| av在线app专区| 狂野欧美激情性xxxx| 男人添女人高潮全过程视频| 男女高潮啪啪啪动态图| 午夜福利,免费看| 一级毛片女人18水好多| 亚洲中文av在线| 99久久综合免费| 久久人妻福利社区极品人妻图片| 女性被躁到高潮视频| 欧美精品人与动牲交sv欧美| 秋霞在线观看毛片| 两性夫妻黄色片| 纯流量卡能插随身wifi吗| 亚洲国产av影院在线观看| 男女午夜视频在线观看| 狂野欧美激情性xxxx| 日韩电影二区| 青春草亚洲视频在线观看| 十八禁网站网址无遮挡| bbb黄色大片| 亚洲国产欧美网| 十八禁网站网址无遮挡| 69av精品久久久久久 | 丝袜美腿诱惑在线| 啦啦啦中文免费视频观看日本| 蜜桃在线观看..| 国产精品秋霞免费鲁丝片| 人妻人人澡人人爽人人| 美女中出高潮动态图| 午夜影院在线不卡| 叶爱在线成人免费视频播放| 黑人猛操日本美女一级片| 久久人人97超碰香蕉20202| 亚洲欧美精品自产自拍| 国产免费视频播放在线视频| 国产伦理片在线播放av一区| 亚洲伊人久久精品综合| 亚洲va日本ⅴa欧美va伊人久久 | 欧美精品一区二区大全| 无限看片的www在线观看| 十八禁网站网址无遮挡| 精品久久蜜臀av无| 99国产极品粉嫩在线观看| 国产黄色免费在线视频| 又黄又粗又硬又大视频| 亚洲视频免费观看视频| 亚洲国产欧美在线一区| 免费黄频网站在线观看国产| 热99re8久久精品国产| 啦啦啦啦在线视频资源| 天天操日日干夜夜撸| 亚洲精品av麻豆狂野| 亚洲成人国产一区在线观看| 国产三级黄色录像| 欧美在线黄色| 99热网站在线观看| 天天躁夜夜躁狠狠躁躁| 国产一区二区三区av在线| 97精品久久久久久久久久精品| 男女下面插进去视频免费观看| av电影中文网址| 午夜福利在线免费观看网站| 一区二区三区四区激情视频| 久久久水蜜桃国产精品网| 两个人看的免费小视频| 日韩欧美一区视频在线观看| 色综合欧美亚洲国产小说| 国产成人系列免费观看| 一本色道久久久久久精品综合| 亚洲精品一区蜜桃| 国产精品秋霞免费鲁丝片| av线在线观看网站| 视频区图区小说| 久久精品国产a三级三级三级| e午夜精品久久久久久久| 国产深夜福利视频在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 国产高清视频在线播放一区 | 亚洲av成人不卡在线观看播放网 | 欧美黑人欧美精品刺激| 天天躁日日躁夜夜躁夜夜| 老鸭窝网址在线观看| 欧美日韩黄片免| 成人黄色视频免费在线看| 妹子高潮喷水视频| 国产一区二区三区av在线| 老司机靠b影院| 别揉我奶头~嗯~啊~动态视频 | 一本—道久久a久久精品蜜桃钙片| 精品亚洲成国产av| 老司机福利观看| h视频一区二区三区| 国产精品免费大片| 欧美在线黄色| 国产男人的电影天堂91| 免费少妇av软件| 爱豆传媒免费全集在线观看| 91成人精品电影| 十分钟在线观看高清视频www| 涩涩av久久男人的天堂| 日本vs欧美在线观看视频| www.自偷自拍.com| 一区二区av电影网| 乱人伦中国视频| 亚洲伊人久久精品综合| 精品人妻1区二区| 夜夜夜夜夜久久久久| 男女午夜视频在线观看| 欧美黄色淫秽网站| 国产一区有黄有色的免费视频| 老司机午夜福利在线观看视频 | 国产无遮挡羞羞视频在线观看| 嫁个100分男人电影在线观看| 另类精品久久| av福利片在线| 精品人妻一区二区三区麻豆| 成人三级做爰电影| 亚洲精品一二三| 亚洲色图综合在线观看| 一个人免费在线观看的高清视频 | 他把我摸到了高潮在线观看 | 自线自在国产av| 亚洲国产中文字幕在线视频| 女人久久www免费人成看片| 曰老女人黄片| 国产国语露脸激情在线看| 在线观看免费视频网站a站| 一级毛片精品| 欧美精品人与动牲交sv欧美| 好男人电影高清在线观看| 国产亚洲午夜精品一区二区久久| 性高湖久久久久久久久免费观看| av免费在线观看网站| 一区二区av电影网| 国产又色又爽无遮挡免| 王馨瑶露胸无遮挡在线观看| 汤姆久久久久久久影院中文字幕| 国产亚洲精品第一综合不卡| 久久综合国产亚洲精品| 精品久久久久久久毛片微露脸 | 美女中出高潮动态图| 国产人伦9x9x在线观看| 欧美午夜高清在线| 搡老乐熟女国产| 丝袜喷水一区| 多毛熟女@视频| 男女床上黄色一级片免费看| 日日摸夜夜添夜夜添小说| 人人妻人人澡人人爽人人夜夜| 高清欧美精品videossex| 交换朋友夫妻互换小说| 欧美性长视频在线观看| 天堂中文最新版在线下载| 99国产精品99久久久久| 日本精品一区二区三区蜜桃| 99久久综合免费| 99精国产麻豆久久婷婷| 一级毛片精品| 老汉色av国产亚洲站长工具| 亚洲av成人一区二区三| 日本a在线网址| 91精品三级在线观看| www.精华液| 在线精品无人区一区二区三| 热re99久久国产66热| 免费观看人在逋| 777久久人妻少妇嫩草av网站| 黄片大片在线免费观看| 亚洲精品国产av蜜桃| 一本色道久久久久久精品综合| 欧美变态另类bdsm刘玥| 久久青草综合色| 国产精品av久久久久免费| 久久精品亚洲熟妇少妇任你| 18在线观看网站| 亚洲精品国产av蜜桃| 在线av久久热| 在线亚洲精品国产二区图片欧美| 欧美黑人精品巨大| 男女国产视频网站| 美女高潮喷水抽搐中文字幕| 高潮久久久久久久久久久不卡| 十八禁网站免费在线| 国产精品 欧美亚洲| 免费在线观看完整版高清| 多毛熟女@视频| 无限看片的www在线观看| 国产有黄有色有爽视频| 欧美日韩av久久| 亚洲精品日韩在线中文字幕| 韩国高清视频一区二区三区| 嫁个100分男人电影在线观看| 久久久久久久大尺度免费视频| 日本vs欧美在线观看视频| h视频一区二区三区| 久久青草综合色| 免费女性裸体啪啪无遮挡网站| 日韩一区二区三区影片| 美女主播在线视频| 亚洲欧美色中文字幕在线| 国产高清videossex| 国产欧美亚洲国产| 高清欧美精品videossex| 18在线观看网站| 色播在线永久视频| 欧美性长视频在线观看| 他把我摸到了高潮在线观看 | 99热网站在线观看| 在线观看免费视频网站a站|