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

    Analysis of reaction torque-based control of a redundant free-floating space robot

    2017-11-17 08:32:02MingheJINChengZHOUYechaoLIUHongLIU
    CHINESE JOURNAL OF AERONAUTICS 2017年5期

    Minghe JIN,Cheng ZHOU,Yechao LIU,Hong LIU

    State Key Laboratory of Robotics and System,Harbin Institute of Technology,Harbin 150001,China

    Analysis of reaction torque-based control of a redundant free-floating space robot

    Minghe JIN,Cheng ZHOU*,Yechao LIU,Hong LIU

    State Key Laboratory of Robotics and System,Harbin Institute of Technology,Harbin 150001,China

    Owing to the dynamics coupling between a free-floating base and a manipulator,the non-stationary base of a space robot will face the issue of base disturbance due to a manipulator’s motion.The reaction torque acted on the satellite base’s centroid is an important index to measure the satellite base’s disturbance.In this paper,a comprehensive analysis of the reaction torque is made,and a novel way to derive the analytical form of the reaction torque is proposed.In addition,the reaction torque null-space is derived,in which the manipulator’s joint motion is dynamically decoupled from the motion of the satellite base,and its novel expression demonstrates the equivalence between the reaction torque null-space and the reaction null-space.Furthermore,the reaction torque acted as an optimization index can be utilized to achieve satellite base disturbance minimization in the generalized Jacobian-based end-effector Cartesian path tracking task.Besides,supposing that the redundant degrees of freedom are abundant to achieve reaction torque-based active control,the reaction torque can be used to realize satellite base attitude control,that is,base attitude adjustment or maintenance.Moreover,because reaction torque-based control is a second-order control scheme,joint torque minimization can be regarded as the optimization task in reaction torque-based active or in-active control.A real-time simulation system of a 7-DOF space robot under Linux/RTAI is developed to verify and test the feasibility and reliability of the proposed ideas.Our extensive empirical results demonstrate that the corresponding analysis about the reaction torque is correct and the proposed methods arefeasible.

    1.Introduction

    A space robot is an important tool for space activities like target capturing,on-orbit servicing,outer space exploration,and so on.It seems to be an extension of human beings’arms,that is to say,our limitations on the field of activities are expanded by the advent of a space robot.Hence,the last 30 years have witnessed an increasing interest towards robotic applications in space.1–5

    Unlike a terrestrial manipulator,a space robot is always freefrom the influence of an external force.The coupling relationship between a satellite base and a mounted manipulator will alter the satellite’s position and posture resulting in difficulty in the design of a space robot’s control systems.Besides,the conventional treatment on this coupling of a free-floating space robot is always from the view of one-order differential equations:the rotational momentum conservation equation and the translational momentum conservation equation.6–8The former is just an integrable expression which can be transformed into a position relationship of each link of a space robot,and hence the conclusion of a free-floating space robot’s invariant centroid is derived.Whereas,the latter translational momentum conservation is just a non-integrable relationship and it is also the non-holonomic constraint of a free-floating space robot. Meanwhile, in practical applications, a satellite’s attitude overwhelms its position,as the attitude is so important when considering solar supplement and information communication.9

    Furthermore,the way to deal with the satellite attitude is always a hot spot in the field of planning and control of a space robot,and to some degree,the satellite attitude adjustment or maintenance will determine the task implementation.10,11Therefore,path planning algorithms and motion control laws are always the footstones of autonomous control of a space robot.In literature,Dubowsky and Torres planned the trajectory of a space manipulator using an enhanced disturbance map(EDM)to minimize the base attitude disturbance.12They took a 2-DOF manipulator for example;however,it is difficult to obtain the EDMs of manipulators with more DOFs.Vafa and Dubowsky used a virtual manipulator model to develop path planning that reduced the base disturbance,which is called the self-correcting path planning algorithm.13In this method,a system is considered as a linear system with the assumption that the movement of joints is small enough.Nakamura and Mukherjee utilized Lyapunov function to achieve the regulations of both the satellite orientation and the manipulator joint angle simultaneously.14However,the stability of this method was mot demonstrated strictly and the planned joint angles were not smooth.Fernands et al.proposed near-optimalnon-holonomic motion planning to achieve attitude control inspired by the fact that a falling cat can change its orientation in midair.15Nenchev et al.originally utilized the notation of reaction null-space(RNS)to achieve attitude control.16,17To sum up,the RNS method is the only method that can achieve both the base attitude regulation and end-effector trajectory planning simultaneously,and it has been verified in the flight experiments of ETS-VII.

    The above proposed algorithms are offline or online motion control of a space robot on the basis of one-order differential motion equations.As for a space robot’s autonomous control,online coordination of the motion between a satellite and a manipulator is crucial.In spite of the fact that the online application of the RNS can be achieved,the end-effector Cartesian path tracking task is limited to position or posture control.The application of the method is quiet a conservative way,and the extended ETS-VII mission is just the path tracking of predefined reactionless paths.In addition,in Cocuzza’s studies,18–20the importance of a reaction torque has been emphasized and the corresponding three algorithms have been proposed,that is,least-squares based reaction torque control,null-space based reaction torque control,and weighted pseudoinverse based reaction torque control aimed at achieving a minimum satellite base attitude disturbance.However,the reaction torque in his studies is the torque acted on the hinge point of the first joint and the satellite base,besides the reaction torque is viewed as the optimization task in his proposed algorithms,and the subsequent experiments are made on a ground-based 3-DOF manipulator on an air-bearing table.

    In addition,the dynamic coupling and mounted manipulator reaction control issues are of vital importance in path planning and autonomous control of a free-floating space robot.

    The main contribution of this paper is to give a comprehensive analysis of the reaction torque acted on a satellite base’s centroid.Firstly,a novel solution to get the analytical expression about the satellite base reaction torque of a free-floating space robot is proposed.Then,when an abundant DOF exists,the reaction torque can be utilized to achieve satellite attitude adjustment or maintenance,that is to say,the reaction torque is acted as a control task to achieve satellite base active control.Besides,when more DOFs are occupied,the reaction torque is viewed as an optimization task to achieve a satellite attitude disturbance minimization task.Moreover,due to the usage of second-order motion equations,the joint driving torque can be included in the optimization task.

    The paper is organized as follows.In Section 2,a kinematics and dynamics model is built,and the inverse kinematics of a space robot is also explained.Then a novel approach to obtain the analytical expression of the reaction torque is developed in Section 3.Section 4 illustrates the reaction torque based inactive control for satellite attitude disturbance minimization.Section 5 shows the reaction torque based active control for satellite attitude adjustment/maintenance with joint torque optimization.In Section 6,a real-time simulator under Linux/RTAI is built and a set of experiments is made to verify the proposed algorithms.The conclusions are summarized in Section 7.

    2.Motion equation of a redundant space robot

    2.1.Dynamic and kinematic model of a space robot

    Considering the linear and angular velocities of the baseis the base linear velocity and wbis the base angular velocity)and the angular velocity of each joint˙θ∈Rn×1,the equation of motion of the space manipulator system is presented as6

    where Hb∈ R6×6is the inertia matrix of the base,Hbm∈ R6×nis the coupled inertia matrix,Hm∈ Rn×nis the inertia matrix of the manipulator,cb∈R6is the velocity-dependent non-linear term of the base,cm∈Rnis the velocity-dependent nonlinear term of the manipulator,and τ∈ Rnis the manipulator joint torque.

    The motion of the free-floating space robot observes the law of momentum conservation,that is,

    where miis the mass of body i,riis the vector from the inertia coordinate system to the centre of the mass,Iiis the inertia matrix of body i,and wiis the angular velocity of rigid body i.

    From the rotational momentum conservation equation,the following expression can be obtained:

    where IS∈ Rl×lis the rotation related inertia matrix of the base.IM∈ Rl×nis the coupling inertia related matrix of the manipulator,in which l=3 stands for the angular momentum solution.L0is the initial momentum.

    Furthermore,assuming that the dimension of the Cartesian space task is m,we can get the following equation from the translational momentum conservation equation:

    From Eq.(5),we have

    Substituting Eq.(7)in Eq.(6)obtains

    2.2.Inverse kinematics of a redundant space robot

    Assuming that the initial of a space robot is still,L0=0 and ˙x0=0 are obtained.On the basis of Eq.(8),the generalized Jacobian based inverse kinematics of a space robot is as follows:13

    where˙ξ is a null-space arbitrary vector.The corresponding second-order form of Eq.(9)is as follows:

    The above Eqs.(9)and(10)are the first-and second-order control laws to achieve generalized Jacobian based endeffector Cartesian path tracking control of a space robot.

    Moreover,from Eqs.(5)and(8),we can get the following equations respectively:

    Combining Eqs.(11)and(12),we can get

    where JEis the extended Jacobian.

    Eq.(13)gives a forward kinematic of a space robot,and its general solution is given as

    which we call the extended Jacobian(EJ)algorithm.

    Furthermore,the acceleration level control law is as follows:

    which is the general solution of Eq.(13).Meanwhile,the second-order forward kinematics is obtained as

    The above Eq.(16)is an important conclusion for the following derivation in Section 5.

    The inverse kinematic based motion control of a space robot can solve satellite base attitude control,zero reaction maneuver,or generalized Jacobian based end-effector Cartesian task control.Besides,the redundant DOF available beyond thatisrequired forthe end-effector Cartesian path tracking is characterized by the null-space of the main task related Jacobian,which can be utilized to assist in the realization of minimizing a position-dependent scalar performance criterion H(θ).Thus the corresponding nullspace vector minimizes˙ξ=k?H(θ)/?θ.Furthermore,the joint limits avoidance and the singularity avoidance index are as follows:21

    joint limits avoidance:

    singularity avoidance:

    3.Reaction torque of a free-floating space robot

    The reaction and joint driving torques are the dynamics quantity of a space robot,and the reaction torque will cause a satellite base disturbance,while the joint torque drives a joint to arrive at the desired position.These torques are shown in Fig.1.

    Fig.1 Schematic diagram of the reaction and joint driving torques of a space robot.

    The satellite base disturbance is directly caused by the reaction force/torque acted on the satellite base’s centroid.The satellite position is changed by the reaction force,while the satellite attitude is altered by the reaction torque.Because the satellite base rotation is a significant issue which leads to communication problems with a ground control centre,the reaction torque is considered prior.The reaction torque can be calculated by the Newton–Euler iterative algorithm.Thus the analytical form of the reaction torque can be derived,but the corresponding derivation will be too sophisticated.The numerical calculation of the reaction torque is very important in the building of the ground verification systems of a space robot,while its analytical form is of vital importance in the development of a space robot’s control algorithm.

    Consider the motion equation of the free-floating base

    By recalling the rotational momentum conservation,from the differential expression of Eq.(5),we can obtain

    and

    From Eq.(21),we can see that the reaction torque acted on the base’s mass centre consists of Mr(θwith respect toand Nr(θ)with respect toθ.

    It is obvious that the above two ways to calculate the reaction torque are equivalent.Besides,the latter is simpler and easier,and this analytical form is obtained by recalling the momentum conservation equation.The derivation Eqs.(19)–(21)also give the relationship between the reaction torque and the rotational momentum conservation of a space robot.

    By recalling the reaction null-space of a free-floating space robot,the rotational momentum conservation equation with a zero attitude disturbance

    yields the following null-space solution:

    In the above Eq.(22),the coupling inertia matrix from the rotational momentum conservation equation comprises a nontrivial kernel,and a set of reactionless joint velocities is as follows:

    As for˙θ∈NRNS,the manipulator joint motion is decoupled from the motion of the satellite base.This reactionless kernel is well known as the RNS.

    Besides,the reactionless control of a space robot can be achieved by the bounded reaction torque,T0=0.Then we can obtain

    which yields the following null-space solution:

    Then,a set of reactionless joint accelerations is as follows:

    This reactionless kernel is well known as the reaction torque null-space(RTNS),in which the manipulator joint motion is dynamically decoupled from the motion of the satellite base.

    Furthermore,the reaction torque can be utilized to achieve satellite attitude control when abundant redundancy exists,besides,joint torque optimization can be achieved in this acceleration level control algorithm.If the end-effector task occupies major DOFs,satellite attitude adjustment/maintenance cannot be achieved,and we aim to obtain satellite disturbance minimization control in this case.

    4.Reaction torque based inactive control for satellite attitude disturbance minimization

    4.1.Null-space based solution for satellite base disturbance minimization

    When the end-effector needs to be sent into the desired position and posture for target capturing,that is to say,major DOFs are occupied to enforce the Cartesian main task,the re will be no abundant DOF to achieve reaction torque based active control.As for the limited redundant DOF,the dynamics task like reaction torque based active control cannot be realized,so the reaction torque based optimization control is adopted in this part to realize satellite minimum disturbance.The redundancy resolution through reaction torque plus driving torque optimization can be concluded into the following quadratic programming(QP)problem:

    where w is the joint torque related weight in the multidynamics-task optimization.

    The cost function is equivalent to the following leastsquares problem that minimizes

    Furthermore,the solution of the above equation can be obtained according to the least-squares solution as

    Thus,the solution of Eq.(28)can be derived by substituting Eq.(30)into Eq.(10)as

    The secondary optimization task is the least-squares problem for reaction torque plus joint driving torque minimization.Assuming that the redundant DOF is limited,the optimization of the joint driving torque works a little in the null-space of a main task.In other words,the optimization effect is up to the dimension of the null-space of a main task.

    4.2.Lagrange multipliers based solution for satellite base disturbance minimization

    Another way to solve the multi-torque optimization under an equality constraint is from the view of Lagrange multipliers.With the utilization of Lagrange multipliers λ,the augmented cost function can be expressed as follows:

    Substituting the lower part of Eq.(27)into the above equation obtains

    The solution is obtained by vanishing the partial derivatives of L(¨θ,λ)with respect to ¨θ and λ as

    where B=MTM ∈ Rn×n.The above Eq.(34)is the Karush-Kuhn-Tucker(KKT)condition of the convex optimization problem Eq.(27),and from ?L/?¨θ to obtain

    Then,substituting Eq.(35)into the expression of ?L/?λ obtains

    Furthermore,the Lagrange multiplier based solution of Eq.(27)with the combination of Eqs.(35)and(36)is

    where is the right pseudoinverse solution of M,andis the weighted pseudoinverse of JG.

    4.3.Equivalence between the proposed solutions

    The above two solutions of the same constrained least squares problem are equivalent,and the necessary and sufficient condition of the equivalence verification is as follows:

    Proof.In order to simplify the derivation,let the matrixand the matrixHence,the proof of Eq.(38)is transformed into the verification that X is the pseudoinverse of A.For the matrix A,the pseudoinverse A+is defined by the following four properties:

    Evaluation of the expression AXA results in

    Evaluation of the expression XAX results in

    Evaluation of the expression (AX)Tresults in

    Evaluation of the expression (XA)Tresults in

    Thus,Eq.(39)is proven.However,if w=0 and M ∈ R3×n,the n M#M= (MTM)-1MTM≠0 and B is not a positive definite,that is to say,the se two ways exist a little miss-distance in this case.□

    In addition,the above proof shows that different ways to achieve the minimization of the same criterion with the same constraints yields the same results.However,the algorithmic singularity exists in the null-space based solution,which will lead to instability in the real applications of the algorithm,and we will give an intuitional explanation in Section 6.

    5.Reaction torque based active control for satellite attitude adjustment/maintenance

    5.1.Reaction torque based active control for satellite attitude adjustment/maintenance

    Generally speaking,the reaction torque is the origin that brings about satellite disturbance,but in this part,the reaction torque will be viewed as an input variable in the control systems of a space robot,and the planned desired reaction torque can realize satellite attitude regulation.A combination of Eq.(16)and the acceleration form of Eq.(6)will lead to a linear system,which is formulated by the joint acceleration dependent expression as

    where wdand wrare the desired and real angular velocities of a satellite base,respectively.In the above equation,the desired reaction torque is up to the desired base angular velocity.

    The general solution of the above reaction torque control and end-effector position control is as follows:

    The difference between Eqs.(15)and(47)is the special solution,while the general solutions belong to the same set of solutions. To a non-redundant manipulator,that is to say,the same set of solutions are specified vectors,and thus the se two are equivalent.Furthermore,the desired reaction torque set to be 0 will guarantee reactionless manipulation as

    5.2.Reaction torque based active control with joint torque optimization

    The second-order reaction torque based control can include an optimization task about dynamics quantity.The local joint torque optimization problem can be expressed as

    There are two ways to solve the above equation.Thefirst way is a null-space based solution,and the cost function is equivalent to the following least-squares problem that minimizes

    Substituting Eq.(47)in Eq.(50),we have

    Furthermore,the solution of the above equation can be obtained according to the least-squares solution as

    Thus,the solution of Eq.(49)can be derived by substituting Eq.(52)into Eq.(47)as

    The above derivation is the null-space based method,and the null-space vector is obtained to achieve the joint torque optimization task.

    Another way to solve the joint torque optimization under an equality constraint is from the view of Lagrange multipliers.With the utilization of Lagrange multipliers,the augmented cost function can be expressed as follows:

    The solution is obtained by vanishing the partial derivatives of)with respect toand λ as

    In order to verify the equivalence between Eqs.(53)and(56),it is sufficient that the following equations exist:

    The proof of the se equations is analogy to the derivation from Eq.(38)to Eq.(45).

    The above ways are the two solutions to solve the leastsquares problem with an equality constraint(LSE).Furthermore,a convenient space robot behavior can be obtained when minimizing the joint torque in terms of a norm weighted by the squared inverse of the inertia matrix as

    Meanwhile,the other torque–level control method for a redundant robot is the dynamically consistent redundancy resolution method to solve the following LSE problem:where τ+and τ-are the upper and lower joint torque bounds,respectively.The Lagrange multipliers based method can give a solution of Eq.(57)and Eq.(58).

    Fig.2 Body fixed frame of a 7-DOF space robot.

    Table 1 Mass properties of the space manipulator.

    In addition,the reaction torque based active control can achieve satellite base attitude adjustment or maintenance,and the dynamics criteria optimization can be introduced.

    6.Algorithm verification in a real-time simulation system

    6.1.Control systems of a redundant space robot

    In order to verify the proposed algorithms,we take a seven-DOF space robot for example.The body fixed frame of the space robot is shown in Fig.2,and the kinematics and dynamics parameters are listed in Table 1.

    To a space robot,ground hardware-based verification methods such as air-bearing table,22airplaneflying or freefalling motion,23and suspension system are costly built.24Recently,a kind of hybrid simulation method following the principles of dynamics simulation and kinematics equivalence has been presented,which is called hardware-in-the -loop simulation systems.25In this part,we build a real-time simulation system of a redundant space robot under Linux/RTAI to explore the capabilities and limitations of the proposed method.This real-time numerical simulator is conceived and realized as an aid to the design of the controller of a new space robot in development.

    Fig.3 Functional architecture of the real time simulator.

    The real-time verification system runs under the Linux operating system with the real time application interface(RTAI)extension,which is a hard real-time extension of Linux.It provides task scheduling and synchronization and inter-task communication services among others.Thefunctional architecture of the simulator is shown in Fig.3.The real-time simulation system consists of 4 modules:

    (1)Forward dynamics simulator of the space robot

    The dynamic model of the space robot is built of fline in the Sim Mechanics,and the MATLAB Real-Time Workshop is utilized to port the virtual prototype model into c codes for online simulation.

    (2)Simulator of the controller

    The controller of the real space robot system consists of a central controller and a joint controller.The central controller is responsible for online visual based autonomous planning,Cartesian path planning,joint space path planning,computed torque control,and so on.The functional diagram of the central controller simulator is shown in Fig.4.Whereas,the simulator of the joint controller is used to simulate the servo system of the space robot,and its functional diagram is shown in Fig.5.

    (3)Task scheduling and synchronization

    The modules mentioned above are made into three realtime tasks and run in three independent threads.Another thread for scheduling the three tasks which is called clock task is created to make the m run periodically,and the synchronous execution of the tasks is achieved by a blocking inter-process/thread communication semaphore.When one executes its period and writes data into the shared memory,the semaphore value is set to one and the next is allowed to run.The shared memory works as a virtual bus.

    (4)Monitor

    The monitor part runs in another Windows PC installed with Labview which can modify control parameters online and save data to of fline analyze fur the r.It is connected with the main simulator by TCP/IP communication.The control parameters could be modified and sent back to the model.

    To sum up,the real-time simulator of a free-floating space robot is built through the analogy with the real control system of the space robot in development.This simulator can be used to verify and develop the path planning algorithm,motion control algorithm,and dynamics based control algorithm,and to respond to offline databases.

    6.2.Experiments results

    The performance of the proposed Cartesian continuous path planning is illustrated by the following simulation examples.The initial joint angles of the space robot are(0°,-90°,0°,0°,0°,0°,0°)and the base attitude is(0°,0°,0°)(Cardan angles).The control cycle is 0.01 s.Furthermore,the following equation is utilized to plan the desired trajectory according to the desired position and posture curves:

    Fig.4 Functional diagram of the central controller simulator.

    Fig.5 Functional diagram of the joint controller simulator.

    where s is the position,sfand s0are the final and initial values,respectively,and tfis the terminal time.

    6.2.1.Example 1:Satellite base attitude disturbance minimization

    In this example,a comparison to analyze the proposed nullspace based solution and the Lagrange multiplier based weighted pseudoinverse solution is made.The main task is the end-effector path tracking task,while the optimization task is the reaction torque plus joint torque minimization.Suppose that the desired end-effector Cartesian position (x,y,z)and posture (α,β,γ)are divided into 4 cases:the first case is(-100,100,50)mm and(2°,2°,2°),and the other three cases are 1.5 times,double,and 2.5 times of the first case separately,that is,(-150,150,75)mm and(3°,3°,3°),(-200,200,100)mm and(4°,4°,4°),and(-250,250,125)mm and(5°,5°,5°).Thesefour cases are noted as 1×,1.5×,2×,and 2.5×,respectively,in Fig.6 and Fig.7,which show the corresponding tracking errors(Δx,Δy,Δz)and (Δα,Δβ,Δγ).

    Fig.6 Tracking errors in Eq.(31)method.

    Fig.7 Tracking errors in Eq.(38)method.

    The results shown in Figs.6 and 7 demonstrate that the tracking errors in Eq.(38)method are smaller than those in Eq.(31)method,that is,the latter is more stable than the former.This is because the latter solution is without algorithmic singularity while the former is easy to divergent with the increasing task dependence.

    Moreover,in order to verify the local optimization algorithm of the reaction torque plus the joint torque,assume that the end-effector Cartesian desired position (x,y,z)and posture (α,β,γ)are(-200,100,100)mm and(5°,5°,0°),respectively,the terminal time is 20 s,and in this process,the position and posture of the end-effector are shown in Fig.8,while the base relative orientation (ψ1,ψ2,ψ3)and the base angular velocity (w1,w2,w3)are shown in Fig.9.Figs.10 and 11 show the reaction torque (Tx,Ty,Tz) and the jointtorque τi,i=1,2,···,7 curves,respectively.In the following,case 1 represents path tracking control without torque optimization and case 2 represents path tracking control with torque optimization.

    Fig.8 Position and posture of the end-effector.

    Fig.9 Relative orientation and angular velocity of the base.

    Fig.10 Reaction torque acted on the free-floating base.

    Fig.11 Joint torque of 7 joints in the base attitude minimization control.

    Fig.13 Attitude of the base and the position of the end-effector.

    Fig.14 Joint torque of 7 joints in base attitude adjustment control.

    Fig.12 Tracking errors of base attitude adjustment and end-effector position tracking tasks.

    From Fig.8,it can be seen that the optimization of the reaction torque plus the joint torque will have a little in fluence on the end-effector path tracking.Besides,from Fig.9,the final base attitude is(0.094,0.123,-0.062)rad if no optimization is made,but with reaction torque based optimization control,the base disturbance is(0.057,0.078,-0.039)rad,and the corresponding base angular velocity is reduced,which is caused by the optimization of the reaction torque shown in Fig.10.Moreover,Fig.11 shows the curves of the joint torque of 7 joints,and it can be recognized that the joint torque is reduced to some degree.From the above results,we can see that the reaction torque is an important index to measure the satellite base attitude disturbance.

    6.2.2.Example 2:Satellite base attitude adjustment

    In this example,a set of experiments is performed to make it clear about the difference between the reaction torque based active control and the rotational momentum based solution for base attitude control of a space robot.Four situations are considered:the first situation is that the desired base attitude (ψ1,ψ2,ψ3)is(1°,1°,0°)and the desired end-effector displacement(x,y,z)is(20,45,30)mm,while the other cases are double,triple,and quadruple of the first case separately.Take the first situation for example,the desired reaction torque is computed by the desired angular velocity and angular acceleration.Fig.12 shows the tracking errors.

    Fig.12 reports the results obtained from the reaction torque based solution which overwhelm those obtained from the momentum based solution.In the four cases,the accuracy of the former solution is higher than that of the latter solution,and the former solution is more stable.This is because the former solution is dynamics-level control considering the centrifugal force,while the latter solution is just a second-order kinematics based control.

    Furthermore,assume that the desired satellite base attitude(ψ1,ψ2,ψ3)is(1°,1°,0°)and the end-effector desired position (x,y,z)is(-40,-80,50)mm.Fig.13 shows the curves of the main tasks,and Fig.14 shows the optimization effect of the joint torque.

    From the above results,the real motion of the space robot agrees with the desired planned trajectory,and it is obvious that the reaction torque based active control can achieve satellite base attitude adjustment effectively.The joint torque optimization will have a little influence on the base attitude control.The multiple dynamics related task is implemented simultaneously,among which the reaction torque is regarded as the active control input value and the joint driving torque is viewed as the optimization task.

    7.Conclusions

    (1)In this paper,a novel way is proposed to derive the analytical form of a reaction torque acted on a satellite base’s centroid.Besides,a reaction torque null-space is obtained on the basis of the analytical form of the reaction torque.Furthermore,an in-depth analysis of the relationship between the RTNS and the RNS is made.

    (2)The reaction torque can be viewed as an optimization index to achieve satellite base attitude disturbance minimization,and this reaction torque based inactive control strategy is of vital importance in the generalized Jacobian based Cartesian task control when most of the DOFs are occupied.

    (3)If redundant DOFs are abundant to achieve reaction torque based active control,the reaction torque can be utilized to achieve satellite attitude control,and the relationship between the reaction torque based base attitude control and rotational momentum based base attitude control is elaborated in this paper.

    (4)A real-time simulator under Linux/RTAI of a freefloating space robot is built to verify the online usage of the proposed method.

    Acknowledgements

    This project was supported in part by the National Program on Key Basic Research Project 973 Program under Grant 2013CB733103 and the Program for New Century Excellent Talents in University under Grand NCET-10-0058.

    1.Flores-Abad A,Ma O,Pham K,Ulrich S.A review of space robotics technologies for on-orbit servicing.Prog Aerosp Sci 2014;68:1–26.

    2.Xu W,Liang B,Xu Y.Survey of modeling,planning,and ground verification of space robotic systems.Acta Astronaut 2011;68(11):1629–49.

    3.Moosavian SAA,Papadopoulos E.Free-flying robots in space:an overview of dynamics modeling,planning and control.Robotica 2007;25(5):537–47.

    4.Wen Z,Wang Y,Kuijper A,Di N,Luo J,Zhang L,et al.On-orbit real-time robust cooperative target identification in complex background.Chin J Aeronaut 2015;28(5):1451–63.

    5.Huang P,Wang D,Meng Z,Zhang F,Guo J.Adaptive postcapture backstepping control for tumbling tethe red space robot-target combination.J Guidance Control Dyn 2016;39(1):150–6.

    6.Dubowsky S,Papadopoulos E.The kinematics,dynamics,and control offree-flying and free-floating space robotic systems.IEEE Trans Robot Autom 1993;9(5):531–43.

    7.Zhang L,Jia Q,Chen G,Sun H.Pre-impact trajectory planning for minimizing base attitude disturbance in space manipulator systemsfora capturetask.Chin J Aeronaut2015;28(4):1199–208.

    8.Fernandes C,Gurvits L,Li ZX.Attitude control of space platform/manipulator system using internal motion.Space robotics:dynamics and control.Springer;1993.p.131–63.

    9.Jin M,Zhou C,Liu Y,Liu H.Cartesian path planning for base attitude adjustment of space robot.In:IEEE international conference on mechatronics and automation;2015 Aug 2–5;Beijing,China;IEEE Press;2015.p.582–7.

    10.Huang P,Wang M,Meng Z,Zhang F,Liu Z.Attitude takeover control for post-capture of target spacecraft using space robot.Aerosp Sci Technol 2016;51:171–80.

    11.Huang P,Wang D,Meng Z,Zhang F,Guo J.Adaptive postcapture backstepping control for tumbling tethe red space robot–target combination.J Guidance Control Dyn 2015;39(1):150–6.

    12.Dubowsky S,Torres M.Path planning for space manipulators to minimize spacecraft attitude disturbances.In:IEEE international conference on robotics and automation 1991.IEEE Press;1991.p.2522–8.

    13.Vafa Z,Dubowsky S.On the dynamics of space manipulators using the virtual manipulator,with applications to path planning.Space Robotics:Dynamics and Control.Springer,1993.p.45–76.

    14.Nakamura Y,Mukherjee R.Nonholonomic motion planning offree- flying space robots via a bi-directional ap-proach.Space Robotics:Dynamics and Control.Springer.1993.p.101–30.

    15.Fernandes C,Gurvits L,Li Z.Near-optimal nonholonomic motion planning for a system of coupled rigid bodies.IEEE Trans.Autom.Control 1994;39(3):450–63.

    16.Nenchev D,Umetani Y,Yoshida K.Analysis of a redundant freeflying spacecraft/manipulator system.IEEE Trans Robot Autom 1992;8(1):1–6.

    17.Nenchev DN,Yoshida K,Vichitkulsawat P,Uchiyama M.Reaction null-space control of flexible structure mounted manipulator systems.IEEE Trans Robot Autom 1999;15(6):1011–23.

    18.Cocuzza S,Pretto I,Debei S.Reaction torque control of redundant space robotic systems for orbital maintenance and simulated microgravity tests.Acta Astronaut 2010;67(3):285–95.

    19.Cocuzza S,Pretto I,Debei S.Least-squares-based reaction control of space manipulators.J Guidance Control Dyn 2012;35(3):976–86.

    20.Cocuzza S,Pretto I,Debei S.Novel reaction control techniques for redundant space manipulators:Theory and simulated microgravity tests.Acta Astronaut 2011;68(11):1712–21.

    21.Liegeois BA.Automatic supervisory control of the configuration and behavior of multibody mechanisms.IEEE Trans Syst Man Cybern 2010;7(12):868–71.

    22.Yoshida Kazuya.Experimental study on the dynamics and control of a space robot with experimental free-floating robot satellite.Adv Robot 1994;9(6):583–602.

    23.Menon C,Aboudan A,Cocuzza S,et al.Free-flying robot tested on parabolic flights:Kinematic control.J Guidance Control Dyn 2005;28(4):623–30.

    24.Fujii H,Uchiyama K,Yoneoka H,et al.Ground-based simulation of space manipulators using test bed with suspension system.Chem Eng Commun 1996;19(5):985–91.

    25.Krenn R,Schaefer B.Limitations of hardware-in-the -loop simulations of space robotics dynamics using industrial robots.European Space Agency-Publications-ESA SP,1999;440:681–6.

    9 September 2016;revised 12 October 2016;accepted 22 December 2016

    Available online 20 April 2017

    Attitude control;

    Base disturbance;

    Dynamics modeling;

    Reaction torque;

    Real time simulator;

    Space robot

    *Corresponding author.

    E-mail address:zhoucheng0818@sina.com(C.ZHOU).

    Peer review under responsibility of Editorial Committee of CJA.

    Production and hosting by Elsevier

    http://dx.doi.org/10.1016/j.cja.2017.02.021

    1000-9361?2017 Chinese Society of Aeronautics and Astronautics.Production and hosting by Elsevier Ltd.

    This is an open access article under the CC BY-NC-ND license(http://creativecommons.org/licenses/by-nc-nd/4.0/).

    ?2017 Chinese Society of Aeronautics and Astronautics.Production and hosting by Elsevier Ltd.This is an open access article under theCCBY-NC-NDlicense(http://creativecommons.org/licenses/by-nc-nd/4.0/).

    久久人人爽人人爽人人片va| 菩萨蛮人人尽说江南好唐韦庄 | 欧美激情在线99| 麻豆成人午夜福利视频| 亚洲国产精品国产精品| 国产精品国产高清国产av| 欧美精品一区二区大全| 高清日韩中文字幕在线| 春色校园在线视频观看| 国产精品麻豆人妻色哟哟久久 | 草草在线视频免费看| 嫩草影院新地址| 少妇人妻精品综合一区二区| 99久久精品一区二区三区| 永久网站在线| 九色成人免费人妻av| av卡一久久| 特大巨黑吊av在线直播| 成人高潮视频无遮挡免费网站| 2022亚洲国产成人精品| 边亲边吃奶的免费视频| 日本色播在线视频| 69av精品久久久久久| 久久精品国产99精品国产亚洲性色| 啦啦啦韩国在线观看视频| 波野结衣二区三区在线| 国产成人91sexporn| 欧美激情在线99| 亚洲欧美精品综合久久99| 亚洲av中文字字幕乱码综合| 亚洲,欧美,日韩| 成人毛片60女人毛片免费| 国产精品久久久久久精品电影| 女人十人毛片免费观看3o分钟| 韩国高清视频一区二区三区| 天堂av国产一区二区熟女人妻| 亚洲电影在线观看av| 久久久成人免费电影| 91精品伊人久久大香线蕉| 极品教师在线视频| 欧美最新免费一区二区三区| 国产成人91sexporn| 亚洲高清免费不卡视频| 日本三级黄在线观看| 99久久九九国产精品国产免费| 日本午夜av视频| 综合色丁香网| 精品久久久久久成人av| 国产av码专区亚洲av| 美女国产视频在线观看| 色播亚洲综合网| 亚洲久久久久久中文字幕| 国内精品一区二区在线观看| 91av网一区二区| 在线观看美女被高潮喷水网站| 青春草国产在线视频| 久久久欧美国产精品| 久久精品人妻少妇| 亚洲av一区综合| 麻豆av噜噜一区二区三区| 日日摸夜夜添夜夜添av毛片| 亚洲自拍偷在线| 偷拍熟女少妇极品色| 国产亚洲一区二区精品| 不卡视频在线观看欧美| 国内揄拍国产精品人妻在线| 99国产精品一区二区蜜桃av| 国产又黄又爽又无遮挡在线| 在线观看av片永久免费下载| 日本猛色少妇xxxxx猛交久久| 99在线视频只有这里精品首页| 欧美一级a爱片免费观看看| 大话2 男鬼变身卡| 噜噜噜噜噜久久久久久91| 国产中年淑女户外野战色| 欧美高清成人免费视频www| 男人舔奶头视频| 舔av片在线| 美女高潮的动态| 国产伦在线观看视频一区| 亚洲乱码一区二区免费版| 我要搜黄色片| 99久久人妻综合| 国产精品1区2区在线观看.| 久久99热这里只有精品18| 久久午夜福利片| 尤物成人国产欧美一区二区三区| 永久免费av网站大全| 国产精品三级大全| 久久久国产成人精品二区| 亚洲在久久综合| www日本黄色视频网| 久久韩国三级中文字幕| 99久久成人亚洲精品观看| 欧美日韩综合久久久久久| 美女黄网站色视频| 国产中年淑女户外野战色| 国产午夜精品论理片| 国产真实伦视频高清在线观看| 精品免费久久久久久久清纯| 天堂影院成人在线观看| 亚洲精品乱码久久久久久按摩| 国产亚洲最大av| 国产精品一二三区在线看| 久久久久免费精品人妻一区二区| 久久这里只有精品中国| 久久久国产成人免费| 小蜜桃在线观看免费完整版高清| 国产真实伦视频高清在线观看| 久久6这里有精品| 亚洲国产精品国产精品| 精品一区二区三区人妻视频| 一级毛片久久久久久久久女| 亚洲国产精品国产精品| 1000部很黄的大片| 乱人视频在线观看| 国产乱人偷精品视频| 国产精品久久久久久久电影| 国产一区二区在线av高清观看| 亚洲人与动物交配视频| 白带黄色成豆腐渣| 九色成人免费人妻av| 日韩一区二区三区影片| 国产精品.久久久| 99热网站在线观看| 韩国av在线不卡| 搞女人的毛片| 日韩三级伦理在线观看| 国产乱人视频| 丝袜美腿在线中文| 免费人成在线观看视频色| 亚洲精品日韩在线中文字幕| 六月丁香七月| 久久久精品大字幕| 亚洲成人精品中文字幕电影| 插阴视频在线观看视频| 一级毛片aaaaaa免费看小| 国产三级在线视频| 国产精华一区二区三区| 国产成人一区二区在线| 亚洲精品乱久久久久久| 麻豆av噜噜一区二区三区| 成年版毛片免费区| 国国产精品蜜臀av免费| 高清毛片免费看| 永久网站在线| 国产精品精品国产色婷婷| 免费黄网站久久成人精品| 亚洲中文字幕日韩| 搡老妇女老女人老熟妇| 午夜精品一区二区三区免费看| 国产三级中文精品| 最后的刺客免费高清国语| 大又大粗又爽又黄少妇毛片口| 午夜日本视频在线| 丰满少妇做爰视频| 看片在线看免费视频| 亚洲国产精品sss在线观看| 97热精品久久久久久| 大又大粗又爽又黄少妇毛片口| 日本五十路高清| 日日摸夜夜添夜夜爱| 男女那种视频在线观看| 久久国产乱子免费精品| 国产精品嫩草影院av在线观看| 日韩欧美在线乱码| 日日撸夜夜添| 国产精品伦人一区二区| 日本免费在线观看一区| 亚洲av一区综合| 天堂影院成人在线观看| 日韩欧美精品v在线| 亚洲婷婷狠狠爱综合网| 久久久欧美国产精品| 日韩精品有码人妻一区| 少妇人妻精品综合一区二区| 全区人妻精品视频| 天美传媒精品一区二区| 欧美潮喷喷水| 99热网站在线观看| 国产一级毛片在线| 中文字幕av成人在线电影| 久久99热这里只频精品6学生 | 成人三级黄色视频| 插阴视频在线观看视频| 伦精品一区二区三区| www.色视频.com| 国产黄色视频一区二区在线观看 | 亚洲欧美日韩高清专用| 人人妻人人澡人人爽人人夜夜 | 日本免费a在线| 国产亚洲av片在线观看秒播厂 | 精品久久久久久电影网 | 中文欧美无线码| 97热精品久久久久久| 小蜜桃在线观看免费完整版高清| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 男女边吃奶边做爰视频| 精品人妻一区二区三区麻豆| 波多野结衣高清无吗| 国产成人a∨麻豆精品| 免费人成在线观看视频色| 婷婷色综合大香蕉| 在线免费观看不下载黄p国产| 亚洲av二区三区四区| 九色成人免费人妻av| 午夜福利成人在线免费观看| 午夜日本视频在线| 老司机影院成人| 国产精品99久久久久久久久| 男女边吃奶边做爰视频| 国产高清三级在线| 国产伦精品一区二区三区视频9| 亚洲人与动物交配视频| 亚洲第一区二区三区不卡| 综合色丁香网| 看黄色毛片网站| 亚洲怡红院男人天堂| 三级毛片av免费| 一个人看的www免费观看视频| 成年女人永久免费观看视频| 久久久久久久久久久丰满| 亚洲激情五月婷婷啪啪| 大香蕉97超碰在线| 搞女人的毛片| 噜噜噜噜噜久久久久久91| 亚洲天堂国产精品一区在线| 午夜a级毛片| 少妇人妻精品综合一区二区| 亚洲中文字幕一区二区三区有码在线看| 一级二级三级毛片免费看| 中文在线观看免费www的网站| 人人妻人人澡欧美一区二区| 久99久视频精品免费| 色综合色国产| 天美传媒精品一区二区| 国产精品人妻久久久久久| 日本wwww免费看| 韩国av在线不卡| 欧美日韩一区二区视频在线观看视频在线 | 又黄又爽又刺激的免费视频.| 村上凉子中文字幕在线| 欧美xxxx性猛交bbbb| 一级黄色大片毛片| 狠狠狠狠99中文字幕| 免费人成在线观看视频色| 搡女人真爽免费视频火全软件| 午夜激情福利司机影院| 纵有疾风起免费观看全集完整版 | 简卡轻食公司| 国产一级毛片七仙女欲春2| 69人妻影院| 亚洲丝袜综合中文字幕| 国产成年人精品一区二区| 欧美成人免费av一区二区三区| 在线观看av片永久免费下载| 好男人在线观看高清免费视频| 人妻制服诱惑在线中文字幕| 欧美一区二区国产精品久久精品| 日韩视频在线欧美| 在线免费观看的www视频| 桃色一区二区三区在线观看| 婷婷色麻豆天堂久久 | 只有这里有精品99| 亚洲人与动物交配视频| 边亲边吃奶的免费视频| 日韩亚洲欧美综合| 欧美日韩一区二区视频在线观看视频在线 | av天堂中文字幕网| 亚洲国产精品成人综合色| 欧美日本视频| 久久久精品94久久精品| 亚洲在线观看片| 国产不卡一卡二| 麻豆精品久久久久久蜜桃| av在线蜜桃| av国产久精品久网站免费入址| 欧美日韩国产亚洲二区| 久久久久久大精品| 午夜福利高清视频| 午夜精品在线福利| 国产av码专区亚洲av| 少妇裸体淫交视频免费看高清| 中文字幕久久专区| 亚洲av免费在线观看| 成人亚洲精品av一区二区| 亚洲精品自拍成人| 欧美xxxx性猛交bbbb| 午夜福利在线观看免费完整高清在| 日韩av在线大香蕉| 伦精品一区二区三区| 国产精品久久久久久精品电影| 热99re8久久精品国产| 黄色一级大片看看| 亚洲av福利一区| 如何舔出高潮| av.在线天堂| 欧美精品国产亚洲| 青春草国产在线视频| 欧美激情久久久久久爽电影| 久久久久久久久大av| 国产真实乱freesex| 91狼人影院| 亚洲欧美成人综合另类久久久 | 永久网站在线| 激情 狠狠 欧美| 午夜精品国产一区二区电影 | 日韩欧美精品免费久久| 国产淫语在线视频| 97人妻精品一区二区三区麻豆| 99在线视频只有这里精品首页| 人体艺术视频欧美日本| 一区二区三区乱码不卡18| 国产毛片a区久久久久| av在线天堂中文字幕| 成人一区二区视频在线观看| 九色成人免费人妻av| 国产一区有黄有色的免费视频 | 淫秽高清视频在线观看| 国产又色又爽无遮挡免| or卡值多少钱| 中国美白少妇内射xxxbb| 久久久久久久午夜电影| 久久久久久久久中文| 色哟哟·www| 成人综合一区亚洲| 国产精品,欧美在线| 国产精品永久免费网站| 在线观看一区二区三区| 午夜福利成人在线免费观看| 免费播放大片免费观看视频在线观看 | 国产成人一区二区在线| 少妇裸体淫交视频免费看高清| 免费av观看视频| 国产成人福利小说| av专区在线播放| 国产av不卡久久| 国产色爽女视频免费观看| 美女国产视频在线观看| 国产 一区 欧美 日韩| 日韩制服骚丝袜av| 亚洲精品一区蜜桃| 精品一区二区免费观看| 国产综合懂色| 亚洲av免费在线观看| 国产大屁股一区二区在线视频| 久久久成人免费电影| 亚洲图色成人| av免费在线看不卡| 国产成人a区在线观看| 在线免费观看的www视频| 中文在线观看免费www的网站| 久久久成人免费电影| 国产精品永久免费网站| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲成人精品中文字幕电影| 久热久热在线精品观看| 色噜噜av男人的天堂激情| 嘟嘟电影网在线观看| 国产一区二区三区av在线| 国产视频内射| 亚洲成人中文字幕在线播放| 亚洲久久久久久中文字幕| 日韩精品有码人妻一区| av黄色大香蕉| 久久久久久久久久黄片| 一级毛片aaaaaa免费看小| АⅤ资源中文在线天堂| 国产伦精品一区二区三区四那| 国产亚洲午夜精品一区二区久久 | 18禁动态无遮挡网站| 欧美日韩精品成人综合77777| 极品教师在线视频| 欧美又色又爽又黄视频| 一本一本综合久久| 国产精品久久久久久久电影| 男人的好看免费观看在线视频| 欧美日本视频| 久久久久久久午夜电影| av在线亚洲专区| 视频中文字幕在线观看| 国产精品.久久久| 欧美成人a在线观看| .国产精品久久| videossex国产| 欧美最新免费一区二区三区| 久久国内精品自在自线图片| 国产女主播在线喷水免费视频网站 | 三级国产精品片| 欧美日本亚洲视频在线播放| 亚洲欧美精品综合久久99| 少妇丰满av| 精品久久久久久成人av| 99久久中文字幕三级久久日本| 国产又黄又爽又无遮挡在线| 少妇的逼好多水| 99热精品在线国产| 欧美性猛交黑人性爽| 91精品一卡2卡3卡4卡| 免费看美女性在线毛片视频| 亚洲在久久综合| 最近中文字幕高清免费大全6| 精品酒店卫生间| 欧美xxxx黑人xx丫x性爽| 人体艺术视频欧美日本| 精品一区二区三区视频在线| 国产极品天堂在线| 亚洲精品456在线播放app| 久久久久精品久久久久真实原创| 亚洲美女搞黄在线观看| 日韩大片免费观看网站 | 久久婷婷人人爽人人干人人爱| 亚洲图色成人| 五月玫瑰六月丁香| 少妇的逼水好多| 亚洲av电影在线观看一区二区三区 | 亚洲人与动物交配视频| 九九在线视频观看精品| 日韩制服骚丝袜av| 午夜亚洲福利在线播放| av专区在线播放| 在线免费观看的www视频| 能在线免费观看的黄片| 老司机影院毛片| 五月玫瑰六月丁香| 国产91av在线免费观看| 国产精品无大码| 久久国内精品自在自线图片| 夫妻性生交免费视频一级片| 亚洲婷婷狠狠爱综合网| 成年版毛片免费区| 听说在线观看完整版免费高清| 99热这里只有是精品50| 少妇被粗大猛烈的视频| 又爽又黄无遮挡网站| 午夜福利在线在线| 99久国产av精品国产电影| 免费无遮挡裸体视频| 一级二级三级毛片免费看| av专区在线播放| 国产成人福利小说| 尤物成人国产欧美一区二区三区| 国产亚洲5aaaaa淫片| 在线免费十八禁| 久久精品91蜜桃| 成人性生交大片免费视频hd| 亚洲欧美精品专区久久| 直男gayav资源| 国产男人的电影天堂91| 97超碰精品成人国产| av在线播放精品| 波野结衣二区三区在线| 国产真实乱freesex| 亚洲精品,欧美精品| 亚洲人成网站在线播| av国产免费在线观看| 国产av一区在线观看免费| 五月玫瑰六月丁香| 国产午夜精品论理片| 一个人观看的视频www高清免费观看| 国产老妇女一区| 成人性生交大片免费视频hd| 亚洲精品日韩av片在线观看| 高清在线视频一区二区三区 | 国产精品,欧美在线| 别揉我奶头 嗯啊视频| 亚洲国产欧洲综合997久久,| 亚洲无线观看免费| 天天躁日日操中文字幕| 三级国产精品欧美在线观看| 三级国产精品片| 国内揄拍国产精品人妻在线| 亚洲三级黄色毛片| 国产黄色小视频在线观看| 女人久久www免费人成看片 | 欧美日韩在线观看h| 国产精品久久电影中文字幕| 久久精品影院6| 欧美日韩国产亚洲二区| 最后的刺客免费高清国语| 22中文网久久字幕| 欧美成人精品欧美一级黄| 丰满人妻一区二区三区视频av| 日韩一区二区视频免费看| 偷拍熟女少妇极品色| 老司机影院毛片| 水蜜桃什么品种好| 在线观看一区二区三区| 国产探花在线观看一区二区| 国产黄色小视频在线观看| 插逼视频在线观看| 麻豆乱淫一区二区| 日韩制服骚丝袜av| 赤兔流量卡办理| 日日干狠狠操夜夜爽| 精品熟女少妇av免费看| 久久精品久久久久久久性| 免费av毛片视频| 午夜视频国产福利| 女人十人毛片免费观看3o分钟| 国产极品精品免费视频能看的| 国产视频内射| 99久久中文字幕三级久久日本| 赤兔流量卡办理| 一二三四中文在线观看免费高清| av在线蜜桃| 国产精品伦人一区二区| av又黄又爽大尺度在线免费看 | 久久久久久国产a免费观看| 久久久久国产网址| 麻豆成人午夜福利视频| 看非洲黑人一级黄片| 一级爰片在线观看| 欧美一区二区亚洲| 夜夜看夜夜爽夜夜摸| 少妇熟女欧美另类| 久久婷婷人人爽人人干人人爱| 国产伦精品一区二区三区四那| 成人午夜高清在线视频| 日韩欧美 国产精品| 国产精品久久久久久久久免| 亚洲熟妇中文字幕五十中出| 少妇熟女aⅴ在线视频| 亚洲最大成人手机在线| 国产一区二区三区av在线| 国产高清三级在线| 午夜亚洲福利在线播放| 大话2 男鬼变身卡| 老司机福利观看| 国产亚洲91精品色在线| 不卡视频在线观看欧美| 国产老妇伦熟女老妇高清| 五月伊人婷婷丁香| 欧美潮喷喷水| 毛片一级片免费看久久久久| 91午夜精品亚洲一区二区三区| 中文字幕制服av| av.在线天堂| 日本爱情动作片www.在线观看| 久久国产乱子免费精品| 国产av不卡久久| 日本一本二区三区精品| 高清av免费在线| 亚洲欧美成人精品一区二区| 久久久精品94久久精品| 国产 一区 欧美 日韩| 国产精品,欧美在线| 欧美变态另类bdsm刘玥| 女人十人毛片免费观看3o分钟| 亚洲成色77777| 午夜免费激情av| 男女啪啪激烈高潮av片| 18禁动态无遮挡网站| 国产精品国产三级国产av玫瑰| 久久精品熟女亚洲av麻豆精品 | 毛片女人毛片| 国产精品国产三级专区第一集| 高清毛片免费看| av又黄又爽大尺度在线免费看 | 变态另类丝袜制服| 看免费成人av毛片| 午夜精品在线福利| 亚洲精品日韩av片在线观看| 成人美女网站在线观看视频| 午夜亚洲福利在线播放| 国产免费福利视频在线观看| 日韩一区二区三区影片| 国产精品一及| 日韩成人伦理影院| 最近手机中文字幕大全| 亚洲自拍偷在线| 亚洲av成人精品一区久久| 69人妻影院| 国产一级毛片七仙女欲春2| 久久久国产成人精品二区| 亚洲av免费在线观看| 特级一级黄色大片| 最近手机中文字幕大全| 亚洲成人中文字幕在线播放| 国产色婷婷99| 69av精品久久久久久| 午夜福利高清视频| 亚洲国产精品久久男人天堂| 综合色av麻豆| 欧美最新免费一区二区三区| 免费在线观看成人毛片| 国产私拍福利视频在线观看| 国产成人一区二区在线| 99久国产av精品国产电影| 97人妻精品一区二区三区麻豆| 男女那种视频在线观看| 搡老妇女老女人老熟妇| 午夜福利网站1000一区二区三区| 亚洲av电影不卡..在线观看| 男女边吃奶边做爰视频| 免费无遮挡裸体视频| 中文字幕精品亚洲无线码一区| 亚洲综合色惰| 麻豆成人av视频| 一本一本综合久久| 亚洲欧美日韩无卡精品| 大又大粗又爽又黄少妇毛片口| 丝袜美腿在线中文| 日韩精品青青久久久久久| 男女那种视频在线观看| 卡戴珊不雅视频在线播放| 亚洲精品自拍成人| 国产乱人偷精品视频| 美女xxoo啪啪120秒动态图| 国产极品精品免费视频能看的| 欧美成人免费av一区二区三区| 成年女人永久免费观看视频| 一级黄色大片毛片| 国产69精品久久久久777片| 国产精品一区二区三区四区久久|