ADAPTIVE FORCE AND MOTION CONTROL OF ROBOT MANIPULATORS IN CONSTRAINED MOTION WITH DISTURBANCES

Similar documents
Adaptive Robust Tracking Control of Robot Manipulators in the Task-space under Uncertainties

THE control of systems with uncertain nonlinear dynamics

Robust Control of Cooperative Underactuated Manipulators

Robot Manipulator Control. Hesheng Wang Dept. of Automation

Adaptive Jacobian Tracking Control of Robots With Uncertainties in Kinematic, Dynamic and Actuator Models

458 IEEE TRANSACTIONS ON CONTROL SYSTEMS TECHNOLOGY, VOL. 16, NO. 3, MAY 2008

Observer Based Output Feedback Tracking Control of Robot Manipulators

Control of industrial robots. Centralized control

Optimal Control of Uncertain Nonlinear Systems

CONTROL OF ROBOT CAMERA SYSTEM WITH ACTUATOR S DYNAMICS TO TRACK MOVING OBJECT

OVER THE past 20 years, the control of mobile robots has

PERIODIC signals are commonly experienced in industrial

q 1 F m d p q 2 Figure 1: An automated crane with the relevant kinematic and dynamic definitions.

AROTORCRAFT-BASED unmanned aerial vehicle

Lecture Schedule Week Date Lecture (M: 2:05p-3:50, 50-N202)

Robot Dynamics II: Trajectories & Motion

Force/Position Regulation for Robot Manipulators with. Unmeasurable Velocities and Uncertain Gravity. Antonio Loria and Romeo Ortega

Adaptive set point control of robotic manipulators with amplitude limited control inputs* E. Zergeroglu, W. Dixon, A. Behal and D.

An Adaptive LQG Combined With the MRAS Based LFFC for Motion Control Systems

A Sliding Mode Controller Using Neural Networks for Robot Manipulator

Video 8.1 Vijay Kumar. Property of University of Pennsylvania, Vijay Kumar

Robust Control of Robot Manipulator by Model Based Disturbance Attenuation

(W: 12:05-1:50, 50-N202)

Force Tracking Impedance Control with Variable Target Stiffness

A Novel Finite Time Sliding Mode Control for Robotic Manipulators

Tracking Control of Robot Manipulators with Bounded Torque Inputs* W.E. Dixon, M.S. de Queiroz, F. Zhang and D.M. Dawson

A Discrete Robust Adaptive Iterative Learning Control for a Class of Nonlinear Systems with Unknown Control Direction

An Adaptive Full-State Feedback Controller for Bilateral Telerobotic Systems

Lecture «Robot Dynamics»: Dynamics 2

Force/Impedance Control for Robotic Manipulators

Mobile Manipulation: Force Control

Real-time Motion Control of a Nonholonomic Mobile Robot with Unknown Dynamics

Neural Network-Based Adaptive Control of Robotic Manipulator: Application to a Three Links Cylindrical Robot

Lecture «Robot Dynamics»: Dynamics and Control

Nonlinear Adaptive Robust Control. Theory and Applications to the Integrated Design of Intelligent and Precision Mechatronic Systems.

A Cascaded-Based Hybrid Position-Force Control for Robot Manipulators with Nonnegligible Dynamics

H State-Feedback Controller Design for Discrete-Time Fuzzy Systems Using Fuzzy Weighting-Dependent Lyapunov Functions

A composite adaptive output feedback tracking controller for robotic manipulators* E. Zergeroglu, W. Dixon, D. Haste, and D.

1030 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 56, NO. 5, MAY 2011

Control of constrained spatial three-link flexible manipulators

GAIN SCHEDULING CONTROL WITH MULTI-LOOP PID FOR 2- DOF ARM ROBOT TRAJECTORY CONTROL

Robotics 2 Robot Interaction with the Environment

Combined NN/RISE-based Asymptotic Tracking Control of a 3 DOF Robot Manipulator

Grinding Experiment by Direct Position / Force Control with On-line Constraint Estimation

Neural Network Control of Robot Manipulators and Nonlinear Systems

EML5311 Lyapunov Stability & Robust Control Design

Inverse differential kinematics Statics and force transformations

Design On-Line Tunable Gain Artificial Nonlinear Controller

Seul Jung, T. C. Hsia and R. G. Bonitz y. Robotics Research Laboratory. University of California, Davis. Davis, CA 95616

AS A POPULAR approach for compensating external

Lyapunov Stability of Linear Predictor Feedback for Distributed Input Delays

Gain Scheduling Control with Multi-loop PID for 2-DOF Arm Robot Trajectory Control

Delay-Independent Stabilization for Teleoperation with Time Varying Delay

Distributed Coordinated Tracking With Reduced Interaction via a Variable Structure Approach Yongcan Cao, Member, IEEE, and Wei Ren, Member, IEEE

Adaptive Robust Control for Servo Mechanisms With Partially Unknown States via Dynamic Surface Control Approach

Simultaneous State and Fault Estimation for Descriptor Systems using an Augmented PD Observer

Exponential Controller for Robot Manipulators

I. MAIN NOTATION LIST

Model Reference Adaptive Control of Underwater Robotic Vehicle in Plane Motion

Stability Analysis and Robust PID Control of Cable-Driven Robots Considering Elasticity in Cables

Design Artificial Nonlinear Controller Based on Computed Torque like Controller with Tunable Gain

On-line Learning of Robot Arm Impedance Using Neural Networks

Case Study: The Pelican Prototype Robot

Trajectory tracking & Path-following control

FINITE HORIZON ROBUST MODEL PREDICTIVE CONTROL USING LINEAR MATRIX INEQUALITIES. Danlei Chu, Tongwen Chen, Horacio J. Marquez

Trajectory-tracking control of a planar 3-RRR parallel manipulator

Control Systems I. Lecture 2: Modeling. Suggested Readings: Åström & Murray Ch. 2-3, Guzzella Ch Emilio Frazzoli

1348 IEEE TRANSACTIONS ON SYSTEMS, MAN, AND CYBERNETICS PART B: CYBERNETICS, VOL. 34, NO. 3, JUNE 2004

Nonlinear PD Controllers with Gravity Compensation for Robot Manipulators

Robust Adaptive Attitude Control of a Spacecraft

NONLINEAR PATH CONTROL FOR A DIFFERENTIAL DRIVE MOBILE ROBOT

TODAY, the developments in technology and the requirements

, respectively to the inverse and the inverse differential problem. Check the correctness of the obtained results. Exercise 2 y P 2 P 1.

A new large projection operator for the redundancy framework

Asymptotic Tracking for Uncertain Dynamic Systems Via a Multilayer Neural Network Feedforward and RISE Feedback Control Structure

Adaptive Control of a Class of Nonlinear Systems with Nonlinearly Parameterized Fuzzy Approximators

Robotics & Automation. Lecture 25. Dynamics of Constrained Systems, Dynamic Control. John T. Wen. April 26, 2007

RESEARCH SUMMARY ASHKAN JASOUR. February 2016

Model Reference Adaptive Control for Robot Tracking Problem: Design & Performance Analysis

A New Approach to Control of Robot

State and Parameter Estimation Based on Filtered Transformation for a Class of Second-Order Systems

Design of a Nonlinear Observer for a Very Flexible Parallel Robot

Lyapunov Stability Analysis of a Twisting Based Control Algorithm for Systems with Unmatched Perturbations

Balancing of an Inverted Pendulum with a SCARA Robot

RBF Neural Network Adaptive Control for Space Robots without Speed Feedback Signal

Indirect Model Reference Adaptive Control System Based on Dynamic Certainty Equivalence Principle and Recursive Identifier Scheme

LYAPUNOV-BASED ROBUST AND ADAPTIVE CONTROL OF NONLINEAR SYSTEMS USING A NOVEL FEEDBACK STRUCTURE

A Nonlinear Disturbance Observer for Robotic Manipulators

Optimal Design of CMAC Neural-Network Controller for Robot Manipulators

Design and Stability Analysis of Single-Input Fuzzy Logic Controller

Dynamic backstepping control for pure-feedback nonlinear systems

The Design of Sliding Mode Controller with Perturbation Estimator Using Observer-Based Fuzzy Adaptive Network

Integrated Design and PD Control of High-Speed Closed-loop Mechanisms

q HYBRID CONTROL FOR BALANCE 0.5 Position: q (radian) q Time: t (seconds) q1 err (radian)

A Hybrid Control Strategy for Robust Contact Detection and Force Regulation

Passivity-based Control of Euler-Lagrange Systems

Force 9. Force Control

Global robust output feedback tracking control of robot manipulators* W. E. Dixon, E. Zergeroglu and D. M. Dawson

Nonlinear Tracking Control of Underactuated Surface Vessel

Reinforcement Learning and Robust Control for Robot Compliance Tasks

Transcription:

ADAPTIVE FORCE AND MOTION CONTROL OF ROBOT MANIPULATORS IN CONSTRAINED MOTION WITH DISTURBANCES By YUNG-SHENG CHANG A THESIS PRESENTED TO THE GRADUATE SCHOOL OF THE UNIVERSITY OF FLORIDA IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF MASTER OF SCIENCE UNIVERSITY OF FLORIDA 2013 1

2013 Yung-Sheng Chang 2

To my wife, Hsiao-Chia Hu 3

ACKNOWLEDGMENTS I would like to thank my committee chair, Dr. Carl D. Crane III, for his patient guidance and for giving me the great chance to study at the Center for Intelligent Machines and Robotics (CIMAR). I also would like to thank my committee members, Dr. Warren Dixon and Dr. Gloria J. Wiens, for their valuable advice to this study. In addition, I would like to thank my companions at CIMAR. From them I learned lot knowledge about robotics. Thanks also go out to my family and wife for their support and love. 4

TABLE OF CONTENTS page ACKNOWLEDGMENTS... 4 LIST OF FIGURES... 7 ABSTRACT... 9 CHAPTER 1 INTRODUCTION... 11 2 CONSTRAINED ROBOT MANIPULATOR DYNAMICS... 14 2.1 Euler-Lagrange Dynamic Model... 14 2.2 Model of Contact Surface... 16 2.3 Dynamic Model of Constrained Robot Manipulator... 17 3 ADAPTIVE MOTION AND FORCE CONTROL WITH RISE STRUCTURE... 21 3.1 Open-Loop Tracking Error System... 21 3.2 Closed-Loop Tracking Error System... 23 3.3 Stability Analysis... 24 4 SIMULATION RESULTS... 28 4.1 Simulation Environment... 28 4.2 Simulation Results without Disturbances... 31 4.3 Simulation Results with Disturbances... 37 4.4 Discussion... 44 5 CONCLUSIONS AND SUGGESTED FUTURE WORKS... 51 5.1 Conclusions... 51 5.2 Suggested Future Works... 52 APPENDIX A PROOF FOR EQUATION 3-21... 53 B PROOF FOR EQUATION 3-26... 57 C PROOF FOR PROPERTIES 2-4... 59 D PROOF FOR PROPERTIES 2-5... 60 5

REFERENCES... 61 BIOGRAPHICAL SKETCH... 63 6

LIST OF FIGURES Figure page 4-1 A two degree of freedom robot manipulator moving on a semi-circle contact surface.... 28 4-2 Desired position trajectory.... 32 4-3 Desired force trajectory.... 32 4-4 Time response of position trajectory in the presence of unknown system parameters but without disturbance.... 33 4-5 The time response of position tracking errors in the presence of unknown system parameters but without disturbance.... 33 4-6 Enlarged diagram of Figure 4-5.... 34 4-7 The time response of force trajectory in the presence of unknown system parameters but without disturbance.... 34 4-8 The time response of tracking errors in the presence of unknown system parameters but without disturbance.... 35 4-9 Enlarged diagram of Figure 4-8.... 35 4-10 The time response of friction force in the presence of unknown system parameters but without disturbance.... 36 4-11 Estimated system parameters in the presence of unknown system parameters but without disturbance.... 36 4-12 Estimated dry contact surface friction in the presence of unknown system parameters but without disturbance.... 37 4-13 Control torque of the robot manipulator in the presence of unknown system parameters but without disturbance.... 37 4-14 Disturbance used for the simulation.... 38 4-15 Disturbance used for the simulation.... 39 4-16 The time response of position trajectory in the presence of unknown system parameters and disturbance.... 39 4-17 The time response of position tracking errors in the presence of unknown system parameters and disturbance.... 40 7

4-18 Enlarged diagram of Figure 4-17.... 40 4-19 The time response of force trajectory in the presence of unknown system parameters and disturbance.... 41 4-20 The time response of force tracking errors in the presence of unknown system parameters and disturbance.... 41 4-21 Enlarged diagram of Figure 4-20.... 42 4-22 The time response of friction force in the presence of unknown system parameters and disturbance.... 42 4-23 Estimated system parameters in the presence of unknown system parameters and disturbance.... 43 4-24 Estimated dry contact surface friction in the presence of unknown system parameters and disturbance.... 43 4-25 Control torque of the robot manipulator in the presence of unknown system parameters and disturbance.... 44 4-26 Position tracking error with disturbance and without disturbance.... 45 4-27 Enlarged diagram of Figure 4-26.... 46 4-28 Force tracking error with disturbance and without disturbance.... 46 4-29 Enlarged diagram of Figure 4-28.... 47 4-30 The time response of position tracking errors with different gains.... 47 4-31 The time response of force tracking errors with different gains.... 48 4-32 The time response of position tracking errors with different gains.... 48 4-33 The time response of force tracking errors with different gains.... 49 4-34 The time response of control torque with different gains.... 49 4-35 The time response of position tracking errors with.... 50 4-36 The time response of force tracking errors with.... 50 8

Abstract of Thesis Presented to the Graduate School of the University of Florida in Partial Fulfillment of the Requirements for the Degree of Master of Science ADAPTIVE FORCE AND MOTION CONTROL OF ROBOT MANIPULATORS IN CONSTRAINED MOTION WITH DISTURBANCES Chair: Carl D. Crane, III Major: Mechanical Engineering By Yung-Sheng Chang May 2013 The goal of this study is to design a controller to reject the disturbances caused by sensor noise and unmodeled effects during hybrid adaptive force and motion control of robot manipulators in constrained motion. A continuous robust integral of the sign of the error (RISE) feedback term is incorporated with adaptive force and motion control to yield asymptotic tracking results in the presence of disturbances and unknown system parameters and dry contact surface friction coefficient. The main reason to use the RISE method is that it can enhance disturbance rejection capabilities. It is assumed that the system parameters of the robot manipulator and the dry friction coefficient of contact surface are unknown. These unknown parameters can be updated by the adaptive update law. The suggested controller can guarantee semiglobal asymptotic motion and force tracking results which are supported through Lyapunov-based stability analysis under the condition that the position and velocity of end-effector and the normal contact force between end-effector and contact surface are measurable. The contact surface of the environment is modeled by the set of m rigid and mutually independent hypersurfaces. The dynamic model of constrained robot 9

manipulators is developed through an Euler-Lagrange formulation. Two degree of freedom (DOF) robot manipulator simulation results are given to illustrate the efficacy of the suggested controller. 10

CHAPTER 1 INTRODUCTION Today, more and more robot manipulators are applied in many different areas. Only position control of manipulators is not enough to satisfy various tasks. Thus, in order to broaden the tasks of manipulators, it will be necessary to control force exerted between the manipulator end-effectors and the contact environment. The major tasks requiring force control are deburring, grinding, contour following, and assembly tasks. In order to perform these tasks successfully, it is important to simultaneously control the motion of the end-effectors and the force of the contact surfaces. A way to achieve simultaneously motion and force control is hybrid position/force control which controls the position of end-effectors in certain directions and controls the force in other directions which is orthogonal to the direction of position [1-3]. During the contact with a rigid surface, the contact object may impose some constraints on robot manipulators. There are some directions that the end-effectors cannot move to. This situation is denoted as constrained motion [4-6]. It must be noted, however, that the orthogonality conditions used in some of these hybrid control algorithms have mixed units. This causes the algorithm to not be invariant to changes in the choice of units of length [21]. There are many other challenging problems to perform complex tasks well, such as unstructured environments, unmodeled effects, and external and force sensor disturbances. Many researches have focused on how to compensate for these effects. Different methods have been proposed to deal with these conditions, such as sliding mode control [7] and robust control [8]. To address unstructured environments and unknown parameters of robot manipulators, an adaptive motion and force control method [9-10] is adopted in this study. 11

An adaptive law is updated by position and force tracking errors to yield semiglobal asymptotic tracking results. To address external and force sensor disturbances, a recently developed continuous Robust Integral of the Sign of the Error (RISE) technique is used [11-12]. The reason to choose RISE technique is that this method can deal with the system with sufficiently smooth bounded disturbances and guarantee asymptotic stability results [13-14]. In this study, the dynamic equation of a general rigid link robot manipulator having n degrees of freedom is modeled using the Euler-Lagrange formulation. The contact surface of the environment is modeled by the set of m rigid and mutually independent hypersurfaces [4, 5 17]. Adaptive motion and force control with RISE feedback terms is used to account for disturbances in the sensor and unmodeled effects, and the presence of uncertain parameters in the robot manipulator and environment. Due to the friction force appearing in the contact surface, the dynamic model proposed by [4] cannot be used. A new transformed dynamic model proposed by [6] is adopted. This model is particularly suitable for the presence of friction force in the contact surface. The uncertain parameters of the robot manipulator and the dry friction coefficient of the contact surface are updated by an adaptive law. The proposed controller can yield semi-global asymptotic tracking results. According to the papers [15-16], the experimental results propose that the PI type force feedback control yields the best force tracking result. The controller proposed in this study has the same structure suggested in the papers. Simulation results are given to illustrate the suggested controller. 12

This study is organized as follows. The dynamic model of the robot manipulator is developed using the Euler-Lagrange formulation in Chapter 2. The filtered tracking errors and the suggested adaptive motion and force controller with the RISE feedback term is given in Chapter 3. Simulation results with different gains and external disturbances are presented in Chapter 4 and conclusions and future works are given in Chapter 5. 13

CHAPTER 2 CONSTRAINED ROBOT MANIPULATOR DYNAMICS 2.1 Euler-Lagrange Dynamic Model The nonlinear dynamic system to be controlled is a general rigid link manipulator with n degrees of freedom. This system can be described using the Euler-Lagrange formulation: ( ) ( ) ( ) ( ) ( ) (2-1) where ( ) is the generalized inertial matrix, ( ) is the generalized Centrifugal and Coriolis force, ( ) is the generalized gravitational force, ( ) is the generalized nonlinear disturbances (e.g., unmodeled effects or force sensor noise), ( ) is the applied torque input, and ( ) ( ) ( ) are the generalized joint position, velocity and acceleration, respectively. When the robot manipulator makes contact with its environment, the contact forces occur between the end-effector and contact surface. The dynamic Equation 2-1 can be modified to [4] ( ) ( ) ( ) ( ) ( ) ( ) (2-2) where ( ) ( ) is the Jacobian matrix, is the position and orientation of the end-effector in Cartesian space, and is the forces/moments exerted by the end-effector on the contact surface. The following development is based on the assumption that ( ), ( ), and the normal contact force are measurable, that ( ) is nonsingular in a finite work space that the end-effector is in contact with the rigid surface at first and the contact forces will keep the end-effector on the contact surface, that the parameters of ( ) ( ) ( ) are unknown and that is unknown and that the robot manipulator is nonredundant. In 14

addition, the following assumptions will be made to facilitate the subsequent Lyapunovbased stability analysis. Assumption 2-1: The inertial matrix ( ) is symmetric, positive definite, and satisfies the following inequality ( ) : ( ) (2-3) where are known positive constant and denotes the standard Euclidean norm. Assumption 2-2: The matrix ( ) ( ) ( ) is a skew-symmetric matrix satisfied the following relationship. ( ) ( )/ (2-4) Assumption 2-3: The desired trajectory is designed as ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ), which exist, and are bounded. Assumption 2-4: The nonlinear disturbance is satisfied such that ( ) ( ) ( ) are bounded by known constants. Assumption 2-5: ( ) ( ) ( ) can be linearly parameterized in terms of suitably selected parameters for robot manipulator. ( ) ( ) ( ) ( ) where is regression matrix and represents the parameters of the robot manipulator. Assumption 2-6: If ( ) ( ), then ( ) and ( ) are bounded. In addition, if ( ) ( ), then the first and second partial derivatives of the elements of ( ) ( ) and ( ) with respect to ( ) exist and are bounded, and the first and 15

second partial derivatives of the elements of ( ) with respect to ( ) exist and are bounded. 2.2 Model of Contact Surface Suppose that the contact surface of the environment can be described by the following set of m rigid and mutually independent hypersurfaces [4, 5, 17] ( ) ( ), ( ) ( )- (2-5) where ( ) is assumed to be twice differentiable with respect to x. At the point of contact on the surface, the contact force including the normal contact force and friction force can be described by the following equation ( ) ( ), ( ) ( )- ( ) ( ) ( ) ( ) (2-6) where is a vector of normal contact force components, is normal contact force in Cartesian space, is friction force in Cartesian space and ( ) is the direction of, is the direction of. The magnitude of is depended on and the friction coefficient, and the direction of is the opposite direction of the end-effector velocity. is differentiable with respect to x except at the point when changes direction. When the robot manipulator is in contact with environment, the end-effector of the robot manipulator is constrained to be on that contact surface and only ( ) degree of freedom can be control independently. Thus, motion control can be described by the following set of ( ) mutually independent curvilinear coordinates ( ), ( ) ( )- (2-7) 16

where ( ) is assumed to be twice differentiable. Because the position control cannot be done in the m rigid and mutually independent hypersurfaces ( ), ( ) is independent of ( ). ( ) and ( ) can uniquely determine the configuration of the robot manipulator. 2.3 Dynamic Model of Constrained Robot Manipulator According to the above definition, the configuration of the robot manipulator can be defined by [4], - (2-8) where, ( ) ( )-, ( ) ( )- Differentiating Equation 2-8 with respect to time, gives (2-9) where ( ), ( ) - ( ) ( ) ( ( )) ( ( )) ( ) (2-10) Substituting Equation 2-8 and 2-9 into the joint space dynamic model Equation 2-2 and multiplying both sides by, gives the operational space dynamic model [1] with the constraints Equation 2-5 and the contact force Equation 2-6 as ( ) ( ) ( ) ( ) ( ) [ ] 0 1 ( ) (2-11) Equation 2-11 can be expressed as ( ) ( ) ( ) ( ) ( ) 17

( ) ( ) ( ) ( ) (2-12) where ( ) ( ) ( ) ( ) [ ( ) ( ) ( ) ( ) ] ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) [ ( ) ( ) ( ) ( ) ] ( ) ( ) [ ( ) ( ) ] ( ) [ ] ( ) ( ) [ ] ( ) [ ] In Equation 2-11, the constraint set by the environment can be simply presented by. The position of the robot manipulator is described by. Note that in operational space dynamic model, normal contact force has a simple structure, i.e.,, -. Also, the contact surface friction force only appears in the second equation of Equation 2-12. Due to this condition, motion and force control cannot be dealt with separately. In other words, they are coupled. The strategy developed by [6] is adopted to deal with this condition. According to [6], ( ) is added and subtracted to the left hand side of the second equation of Equation 2-12, where { }, and is added to the both sides of the first equation of Equation 2-12, where { }. In this way, Equation 2-12 can be reformed as the following 18

( ) ( ) ( ) ( ) ( ) (2-13) where [ ] ( ) [ ( ) ( ) ( ) ] ( ) [ ] ( ) ( ) ( ) ( ) [ ( ) ] [ ] According to the assumption 2-1, 2-2, and 2-5, the following properties can be obtained for Equation 2-13 in the Appendix. Properties 2-1: The matrix ( ) is symmetric, positive definite, and satisfies the following inequality ( ) : ( ) (2-14) where are known positive constant and denotes the standard Euclidean norm. Properties 2-2: The matrix ( ) ( ) ( ) is a skew-symmetric matrix which satisfies the following relationship. ( ) ( )/. (2-15) Properties 2-3: ( ) ( ) ( ), and ( ) can be linear parameterized in terms of a suitably selected set of parameters for robot manipulators. 19

Properties 2-4 1 : The matrix ( ) is a symmetric positive definite matrix on the assumption that the maximum eigenvalue of is small enough. Properties 2-5 2 : The matrix ( ) ( ) ( ) is a skew-symmetric matrix which satisfies the following relationship ( ( ) ( )). (2-16) The controller proposed in this study is based on Equation 2-13 which designs the control torque ( ) to let the position and force follow the desired trajectory. The controller is also incorporated with a RISE feedback structure to reject the disturbance term ( ) in Equation 2-13. 1 See Appendix for the proof of matrix ( ) is a symmetric positive definite matrix 2 See Appendix for the proof of ( ) is skew-symmetric matrix. 20

CHAPTER 3 ADAPTIVE MOTION AND FORCE CONTROL WITH RISE STRUCTURE 3.1 Open-Loop Tracking Error System The goal of this study is to design an adaptive motion and force controller with a RISE feedback term to track the desired time-varying motion and force trajectory with disturbance, based on the system described by Equation 2-14. Let the position and force tracking error be defined as ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) (3-1) where ( ) is the desired robot manipulator trajectory and ( ) is the desired constrained force trajectory. To facilitate the following Lyapunov-based stability analysis, the filtered tracking errors are defined as ( ) ( ) ( ) ( ) (3-2) ( ) ( ) ( ) ( ) (3-3) (3-4) (3-5) where denote positive constants. The proposed controller is based on the assumption that position, velocity, and normal contact force are measurable and acceleration is not measurable, so the filtered tracking error ( ) and are not measurable. Multiplying both sides of Equation 3-5 by ( ) and utilizing Equation 2-13, 3-1, 3-2, 3-3 and 3-4 gives the open-loop tracking error system as ( ) ( ) ( ) ( ) ( ) ( ) 21

( ) ( )( ). (3-6) By adding and subtracting ( ) ( ) ( ) ( ) to Equation 3-6, gives ( ) ( ) ( ) (3-7) where is defined as ( ) ( ) ( ) ( ). (3-8) In Equation 3-8, represents the unknown parameters needed to be updated online and is the desired regression matrix containing the desired robot manipulator position, velocity, acceleration, and force, respectively. In Equation 3-7, the auxiliary function is defined as ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ). (3-9) Based on the Equation 3-7, the control torque is designed to be ( ) (3-10) where is the estimated system parameter vector and is the RISE feedback term. is defined as ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ), ( ) ( ). ( )/ ( ) ( ). ( )/- (3-11) where are positive constant control gains. The estimated system parameter vector is generated by the following update law (3-12) 22

where is a known, constant, diagonal, positive-definite gain matrix. Because Equation 3-12 contains the unknown signal, acceleration and derivative of force error, can be integrated by parts as follows to get rid of the immeasurable signal ( ) ( ) ( ) [ ( )], (3-13) ( ( ) ( )). (3-14) According to the above equations, the terms in the control torque ( ) are all measurable and thus the control torque is implementable. 3.2 Closed-Loop Tracking Error System The closed-loop tracking error system can be obtained by putting Equation 3-10 into Equation 3-7 as ( ) ( ) ( ) (3-15) where ( ) denotes the estimated system parameter error vector defined as. (3-16) According to the subsequent Lyapunov-stability analysis, the time derivative of Equation 3-15 is obtained as ( ) ( ) ( ) ( ) ( ) (3-17) where the auxiliary term ( ) is defined as ( ) ( ) (3-18) and ( ) is defined as ( ) ( ). (3-19) The time derivative of ( ) is obtained as ( ) ( ) ( ) ( ) ( ). (3-20) 23

According to [18], the Mean Value Theorem can be applied to Equation 3-18 to obtain the upper bound 1 of ( ) as ( ) ( ) (3-21) where ( ) is defined as ( ) [ ]. (3-22) The upper bound of disturbance term ( ) and its time derivative ( ) can be obtained based on assumption 2-4 as ( ) ( ) (3-23) where are known positive constants. 3.3 Stability Analysis Theorem 3-1: The proposed controller given in Equation 3-10, 3-11, and 3-13 ensures that the signals of the constrained robot manipulator described by Equation 2-13 are bounded under closed-loop operation and the position and force tracking error is regulated as ( ) ( ) (3-24) on the condition that gain and are chosen large enough based on the initial conditions of the constrained robot manipulator system, and are chosen according to the subsequent proof as (3-25) and and are chosen according to the Lemma 2 in the Appendix as. (3-26) 1 See Lemma 1 of the Appendix for the proof of the upper bound of ( ). 24

Proof: A positive definite function is selected as (3-27) where the auxiliary function ( ) ( ) are defined as ( ) ( ) ( ) ( ) ( ) (3-28) ( ) ( ) ( ) ( ) ( ) (3-29) and the subscript denotes the th element of the vector. In Equation 3-28 and 3-29, the auxiliary function ( ) ( ) are defined as ( ) ( ( ) ( )) (3-30) ( ) ( ( ) ( )). (3-31) The derivative of ( ) ( ) with respect to time can be obtained as ( ) ( ) ( ( ) ( )) (3-32) ( ) ( ) ( ( ) ( )). (3-33) If the sufficient condition in Equation 3-26 is satisfied, the following inequality can hold 2 ( ) ( ) ( ) ( ) (3-34) ( ) ( ) ( ) ( ). (3-35) Thus, ( ) ( ) are obtained. Taking the time derivative of, the following equation can be obtained (3-36) ( ) ( ) ( ( ) ( ) ( ) ( ) ( ) 2 See Lemma 2 of the Appendix for the proof of the inequality in Equation 3-26 25

( ) ) ( ). ( ) ( )/. ( ) ( )/ (3-37) ( ) ( ) ( ) (3-38) Using Equation 3-22 and the following inequality (3-39) can be upper bounded in terms of ( ) as ( ) ( ) ( ( ) ). ( ) / (3-40) where * + { } and ( ) is a positive globally invertible nondecreasing function. According to Equation 3-27 and 3-40, and. Thus, ( ) ( ) ( ) ( ) ( ). Since ( ) ( ) ( ) ( ) and ( ), then according to Equation 3-1, 3-2, 3-3 and 3-4,. Since and is unknown constant vector, then according to Equation 3-16,. Since ( ) ( ) ( ) ( ) ( ) and assumption 2-3 that ( ) ( ) ( ) ( ) ( ) exist and are bounded, then ( ) ( ) ( ) ( ) ( ). Since ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) exist and are bounded, then 26

( ) ( ) ( ). Since ( ) ( ) ( ) ( ) ( ), then input torque ( ) and is implementable. Since ( ) ( ) ( ) ( ), then ( ) and ( ) are uniformly continuous. Moreover, since ( ) ( ). ( ) / ( ) ( ) 2 (3-41) applying Barbalat s lemma, if ( ), then ( ) and ( ) semiglobally asymptotically regulates to zero, ( ) ( ) as. 27

controller. CHAPTER 4 SIMULATION RESULTS 4.1 Simulation Environment The simulation environment depicted in Figure 4-1 was used to test the proposed Figure 4-1. A two degree of freedom robot manipulator moving on a semi-circle contact surface. A two degree of freedom robot manipulator is chosen for the simulation. The dynamic system matrices, ( ) ( ), forward kinematics and Jacobian matrix are given by ( ) [ ] ( ) [ ( ) ] 0 1 [ ] ( ) ( ) [ ] 28

( ) [ ( ) ( ) ] where ( ) ( ) ( ) ( ) ( ) ( )., - is system parameter calculated by [ ] [ ] The actual values of system parameters are, -, and, It is assumed that the value of is unknown in simulation and its initial estimate value is, -. As seen in Figure 4-1, the contact surface between the robot manipulator and the environment is a semi-circle surface S with a dry friction coefficient. The surface is described as ( ) The task space of the robot manipulator is described as, - ( ) where is orthogonal to the curvilinear coordinate of. According to Equation 2-6, the interaction force between the contact surface and the end-effector are given by ( ) ( ) [ ( ) ( ) ] ( ) [ ( ) ] ( ) ( ) ( ) 29

[ ( ) ( ) ( ) ( ) ] ( ( ) ( ) ( ) [ ( ) ( ) ] ( ) [ ( ) ( ) ] ( ) where, normal contact force component, can be measured by the force sensor. Based on the above equations, Operational space dynamics can be derived and the dynamics Equation 2-13 controlled by the proposed controller can be obtained. The following simulation is based on the assumption that the position and velocity of the end-effector and the normal contact force are measurable, that ( ) is nonsingular in a finite work space that the end-effector is in contact with the rigid surface with dry friction force at first, and the end-effector will not leave the contact surface. The control torque is given by Equation 3-10, 3-11, and 3-13. The parameter estimated vector is defined as [ ]. (4-1) The parameter estimated vector Equation 4-1 is generated by the updated law Equation 3-13. The regression matrix can be obtained by Equation 3-8. Parameter values of the controller are, - and * +. The desired position and force trajectory are given by ( ( )) Sampling time is 0.001 sec ( ) 30

4.2 Simulation Results without Disturbances The desired position trajectory is shown in Figure 4-2 and desired force trajectory is shown in Figure 4-3. The time response of position trajectory and its position tracking errors in the presence of unknown system parameters but without disturbance are shown in Figure 4-4 and Figure 4-5, respectively. Figure 4-6 is the enlarged diagram of Figure 4-5. The time response of force trajectory and its force tracking errors in the presence of unknown system parameters but without disturbance are shown in Figure 4-7 and Figure 4-8, respectively. Figure 4-9 is the enlarged diagram of Figure 4-8. The time response of friction force is shown in Figure 4-10. The estimated system parameters and dry contact surface friction are shown in Figure 4-11 and Figure 4-12, respectively. Figure 4-13 shows the control torque of the robot manipulator. As shown in the Figure 4-5 and Figure 4-6, the position tracking error falls between -0.001m and 0.001m. The proposed controller has a good position tracking ability. According to Figure 4-8 and Figure 4-9, the proposed controller also has good force tracking ability. The force tracking error falls between -1N and 1N. In Figure 4-11 and Figure 4-12, the estimated parameters do not converge to their actual value. This phenomenon may be addressed by using composite adaptive update law which may be studied in the future work. 31

Figure 4-2. Desired position trajectory. Figure 4-3. Desired force trajectory. 32

Figure 4-4. Time response of position trajectory parameters but without disturbance. in the presence of unknown system Figure 4-5. The time response of position tracking errors in the presence of unknown system parameters but without disturbance. 33

Figure 4-6. Enlarged diagram of Figure 4-5. Figure 4-7. The time response of force trajectory parameters but without disturbance. in the presence of unknown system 34

Figure 4-8. The time response of tracking errors parameters but without disturbance. in the presence of unknown system Figure 4-9. Enlarged diagram of Figure 4-8. 35

Figure 4-10. The time response of friction force in the presence of unknown system parameters but without disturbance. Figure 4-11. Estimated system parameters parameters but without disturbance. in the presence of unknown system 36

Figure 4-12. Estimated dry contact surface friction parameters but without disturbance. in the presence of unknown system Figure 4-13. Control torque of the robot manipulator in the presence of unknown system parameters but without disturbance. 4.3 Simulation Results with Disturbances The desired position and force trajectory is the same as shown in Figure 4-2 and Figure 4-3, respectively. The disturbance is shown in Figure 4-14 and Figure 4-15. As per assumption 2-4, the nonlinear disturbance is bounded by known constants. The time response of the position trajectory and its position tracking errors in the 37

presence of unknown system parameters and disturbance are shown in Figure 4-16 and Figure 4-17, respectively. Figure 4-18 is the enlarged diagram of Figure 4-17. The time response of force trajectory and its force tracking errors in the presence of unknown system parameters and disturbance are shown in Figure 4-19 and Figure 4-20, respectively. Figure 4-21 is the enlarged diagram of Figure 4-20. The time response of friction force is shown in Figure 4-22. The estimated system parameters and dry contact surface friction are shown in Figure 4-23 and Figure 4-24, respectively. Figure 4-25 shows the control torque of the robot manipulator. As shown in the Figure 4-17 and Figure 4-18, the position tracking error falls between -0.001m and 0.001m like Figure 4-8 and Figure 4-9. The proposed controller has good position disturbance rejection ability. According to Figure 4-20 and Figure 4-21, the proposed controller also has good force disturbance rejection ability. The force tracking error falls between -1N and 1N. Like the estimated system parameter vector in Figure 4-11 and Figure 4-12, the estimated parameters of system with disturbances do not converge to their actual value. Figure 4-14. Disturbance used for the simulation. 38

Figure 4-15. Disturbance used for the simulation. Figure 4-16. The time response of position trajectory system parameters and disturbance. in the presence of unknown 39

Figure 4-17. The time response of position tracking errors unknown system parameters and disturbance. in the presence of Figure 4-18. Enlarged diagram of Figure 4-17. 40

Figure 4-19. The time response of force trajectory parameters and disturbance. in the presence of unknown system Figure 4-20. The time response of force tracking errors system parameters and disturbance. in the presence of unknown 41

4 Figure 4-21. Enlarged diagram of Figure 4-20. Figure 4-22. The time response of friction force in the presence of unknown system parameters and disturbance. 42

Figure 4-23. Estimated system parameters parameters and disturbance. in the presence of unknown system Figure 4-24. Estimated dry contact surface friction parameters and disturbance. in the presence of unknown system 43

Figure 4-25. Control torque of the robot manipulator in the presence of unknown system parameters and disturbance. 4.4 Discussion Figures 4-26 and 4-27 show the position tracking error without disturbance and with disturbance in the same figure. From these figures, the position tracking error with disturbance is almost the same as the position tracking error without disturbance. It shows that, due to the RISE feedback term, the disturbance rejection ability is good. Figures 4-28 and 4-29 show the force tracking error without disturbance and with disturbance in the same figure. From these figures, the force tracking error with disturbance is a little bigger than the force tracking error without disturbance. Despite this situation, the proposed controller also has a good force disturbance rejection ability. Noted that in these figures, there are abrupt changes at. Because the direction of the friction force depends on the direction of end-effector velocity, the abrupt direction change of the end-effector velocity can be seen in Figures 4-4 and 4-16. To make comparison, the simulation is run with the same parameters as before except to show how the different gain will affect the result in Figure 4-30 and 4-31. 44

The position tracking errors look alike, but the force tracking errors have much difference. It agrees with the Properties 2-4 that has to be small enough. The simulation is also run with the same parameters as before except in Figures 4-32, 4-33, and 4-34. The transient position tracking error improves a lot, but the force tracking error doesn t ameliorate in this condition and the joint torque become bigger. In the last simulation, the system is run to test the validity of Theorem 3-1. Figures 4-35 and 4-36 are given to show that the position and force tracking error don t converge to zero under the condition that. These figures show that if do not satisfy Theorem 3-1 conditions, the system would not be stable. Figure 4-26. Position tracking error with disturbance and without disturbance. 45

Figure 4-27. Enlarged diagram of Figure 4-26. Figure 4-28. Force tracking error with disturbance and without disturbance. 46

Figure 4-29. Enlarged diagram of Figure 4-28. Figure 4-30. The time response of position tracking errors with different gains. 47

Figure 4-31. The time response of force tracking errors with different gains. Figure 4-32. The time response of position tracking errors with different gains. 48

Figure 4-33. The time response of force tracking errors with different gains. Figure 4-34. The time response of control torque with different gains. 49

Figure 4-35. The time response of position tracking errors. with Figure 4-36. The time response of force tracking errors. with 50

CHAPTER 5 CONCLUSIONS AND SUGGESTED FUTURE WORKS 5.1 Conclusions This work focused on designing an adaptive controller for constrained robot manipulator with unknown parameters both in the robot manipulator and the contact surface with disturbance caused by force sensor noise and unmodeled effects. First, the dynamic equation of a general rigid link robot manipulator having n degrees of freedom modeled by Euler-Lagrange formulation is used in this study. Because of the friction force occurring between contact surface and end-effector, the position and force control cannot be controlled separately. The Euler-Lagrange dynamic model is modified to suit this condition. Second, the contact surface of the environment is modeled by the set of m rigid and mutually independent hypersurfaces. When the robot manipulator is in contact with environment, the end-effector of robot manipulator is constrained to be on that contact surface. Thus, only (n-m) position coordinates and m force coordinates need to be controlled. Third, the proposed controller is an adaptive hybrid position/force controller with RISE feedback structure. By using the Lyapunov-stability analysis, the suggested controller can guarantee semi-global asymptotic position and force tracking result even if in the presence of the disturbance. The only assumption about the disturbance is that it has to be bounded by known constants. In conclusion, the suggested force feedback term in the controller is in accordance with PI type. Intensive simulation results were given to show the validity of the proposed controller. 51

5.2 Suggested Future Works The contact surface model in this study is assumed to be a rigid contact surface, but in most application areas, the contact surface may not be rigid. The stiffness of the contact surface needs to be considered in the contact model [19]. In addition, the robot manipulator is assumed to be non-redundant in this study. The future work can focus on how to extend the controller proposed in this study to the redundant robot manipulator which has some advantages. For example, a problem occurs, when the Jacobian matrix becomes linearly dependent. In this case, a nonredundant robot manipulator is at a singular configuration and loses the degree of freedom in some direction [1]. Redundant robot manipulator can deal with this situation. Another suggested future work is to improve the estimated system parameters and dry contact surface friction coefficient. In this study, the estimated vector doesn t converge to their actual values. According to [20], the composite adaptive control skill may be incorporated with the controller proposed in this study. The main idea of composite adaptive control skill is that it includes the prediction error of the system parameters estimation in the adaptive updated law. The prediction error is defined as the difference between the estimated system parameters and actual system parameters. Thus, including prediction error in the Lyapunov-stability analysis can prove that estimated system parameters converge to their actual values as time goes to infinity. 52

APPENDIX A PROOF FOR EQUATION 3-21 Lemma 1. The Mean Value Theorem can be applied to prove the upper bound of ( ) in Equation 3-21 ( ) ( ) where ( ) is defined as ( ) [ ] Proof: ( ) can be written in the following form ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) 53

The Mean Value Theorem can be applied to ( ) and rewritten as ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) where ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) Thus, ( ) can be upper bounded as ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) Noted that the partial derivatives can be upper bounded as 54

( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) Using the following inequality and ( ), ( ) can be further upper bounded as ( ). ( ) ( ) ( )/ ( ( ) ( )) ( ( ) ( )) ( ) ( ) ( ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( )) ( ) Thus, 55

( ) ( ) where ( ) is some positive globally invertible nondecreasing function. 56

APPENDIX B PROOF FOR EQUATION 3-26 Lemma 2 The function ( ) ( ) is defined as ( ) ( ( ) ( )) ( ) ( ( ) ( )) If the following sufficient conditions is satisfied Then ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) Proof: Integrating both sides of Equation 3-30 ( ) ( ) ( ( ). ( )/) (A-1) ( ( )) ( ( )). ( ) ( )/ Using integration by part, the first term of Equation A-1 can be expressed as ( ( )) ( ) ( ) ( ) ( ) ( ) Using the following property, ( ) ( ) ( ) Equation A-1 can be rewritten as ( ) ( ) ( ) ( ) ( ) ( ( ) ) ( ) ( ( ) ( ) ) 57

Thus, according to Equation 3-23, if, then ( ) is positive definite function. In the similar way as above proof, integrating both sides of Equation 3-31 ( ) ( ) ( ( ). ( )/) (A-2) ( ( )) ( ( )). ( ) ( )/ Using integration by part, the first term of Equation A-2 can be expressed as ( ( )) ( ) ( ) ( ) ( ) ( ) Using the following property, ( ) ( ) ( ) Equation A-2 can be rewritten as ( ) ( ) ( ) ( ) ( ) ( ( ) ) ( ) ( ( ) ( ) ) Thus, according to Equation 3-23, if, then ( ) is positive definite function. 58

APPENDIX C PROOF FOR PROPERTIES 2-4 Properties 2-4 The matrix ( ) is a symmetric positive definite matrix on the assumption that the maximum eigenvalue of is small enough. Proof: ( ) can be rewritten as [ ( ) ( ) ] ( ) (A-3) If the maximum eigenvalue of, then the minimum eigenvalue of ( ). According to assumption 2-1, ( ) ( ) is a symmetric positive semidefinite matrix. Thus, ( ) is a symmetric positive definite matrix on the assumption that the maximum eigenvalue of is small enough. 59

APPENDIX D PROOF FOR PROPERTIES 2-5 Properties 2-5 The matrix ( ) ( ) ( ) is a skew-symmetric matrix satisfied the following relationship ( ( ) ( )) Proof: According to Equation A-3, ( ) can be rewritten as the following ( ) ( ) [ ( ) ( ) ] From assumption 2-2 and the above equation lead to the properties 2-5. 60

REFERENCES [1] O. Khatib, "A Unified Approach for Motion and Force Control of Robot Manipulators - the Operational Space Formulation," IEEE Journal of Robotics and Automation, vol. 3, pp. 43-53, Feb 1987. [2] M. H. Raibert and J. J. Craig, "Hybrid Position-Force Control of Manipulators," Journal of Dynamic Systems Measurement and Control-Transactions of the ASME, vol. 103, pp. 126-133, 1981. [3] D. E. Whitney, "Historical-Perspective and State-of-the-Art in Robot Force Control," International Journal of Robotics Research, vol. 6, pp. 3-14, Spr 1987. [4] N. H. Mcclamroch and D. W. Wang, "Feedback Stabilization and Tracking of Constrained Robots," IEEE Transactions on Automatic Control, vol. 33, pp. 419-426, May 1988. [5] T. Yoshikawa, T. Sugie, and M. Tanaka, "Dynamic Hybrid Position Force Control of Robot Manipulators - Controller-Design and Experiment," IEEE Journal of Robotics and Automation, vol. 4, pp. 699-705, Dec 1988. [6] B. Yao and M. Tomizuka, "Adaptive control of robot manipulators in constrained motion - Controller design," Journal of Dynamic Systems Measurement and Control-Transactions of the ASME, vol. 117, pp. 320-328, Sep 1995. [7] J. J. E. Slotine and W. P. Li, "On the Adaptive-Control of Robot Manipulators," International Journal of Robotics Research, vol. 6, pp. 49-59, Fal 1987. [8] Z. Qu, Robust Control of Nonlinear Uncertain Systems. New-York: Wiley, 1998. [9] J. J. E. Slotine and W. P. Li, "On the Adaptive-Control of Robot Manipulators," International Journal of Robotics Research, vol. 6, pp. 49-59, Fal 1987. [10] J. J. Craig, P. Hsu, and S. S. Sastry, "Adaptive-Control of Mechanical Manipulators," International Journal of Robotics Research, vol. 6, pp. 16-28, Sum 1987. [11] P. M. Patre, W. MacKunis, C. Makkar, and W. E. Dixon, "Asymptotic tracking for systems with structured and unstructured uncertainties," IEEE Transactions on Control Systems Technology, vol. 16, pp. 373-379, Mar 2008. [12] B. Xian, D. M. Dawson, M. S. de Queiroz, and J. Chen, "A continuous asymptotic tracking control strategy for uncertain nonlinear systems," IEEE Transactions on Automatic Control, vol. 49, pp. 1206-1211, Jul 2004. [13] Z. Cai, M. S. de Queiroz, and D. M. Dawson, "Robust adaptive asymptotic tracking of nonlinear systems with additive disturbance," IEEE Transactions on Automatic Control, vol. 51, pp. 524-529, Mar 2006. 61

[14] B. Xian, M. S. de Queiroz, D. M. Dawson, and M. L. McIntyre, "A discontinuous output feedback controller and velocity observer for nonlinear mechanical systems," Automatica, vol. 40, pp. 695-700, Apr 2004. [15] E. Paljug, T. Sugar, V. Kumar, and X. Yun, "Some important considerations in force control implementation," in Robotics and Automation, 1992. Proceedings., 1992 IEEE International Conference on, 1992, pp. 1270-1275 vol.2. [16] R. Volpe and P. Khosla, "A Theoretical and Experimental Investigation of Explicit Force Control Strategies for Manipulators," IEEE Transactions on Automatic Control, vol. 38, pp. 1634-1650, Nov 1993. [17] R. M. Murray, Z. Li, and S. S. Sastry, A Mathematical Introduction to Robotic Manipulation, 1st ed. CRC Press, Mar. 1994. [Online]. Available: http://www.worldcat.org/isbn/0849379814 [18] B. Xian, D. M. Dawson, M. S. de Queiroz, and J. Chen, "A continuous asymptotic tracking control strategy for uncertain nonlinear systems," IEEE Transactions on Automatic Control, vol. 49, pp. 1206-1211, Jul 2004. [19] R. Volpe and P. Khosla, "An experimental evaluation and comparison of explicit force control strategies for robotic manipulators," in Robotics and Automation, 1992. Proceedings., 1992 IEEE International Conference on, 1992, pp. 1387-1393 vol.2. [20] P. M. Patre, W. MacKunis, M. Johnson, and W. E. Dixon, "Composite adaptive control for Euler-Lagrange systems with additive disturbances," Automatica, vol. 46, pp. 140-147, Jan 2010. [21] J. Duffy, "The fallacy of modern hybrid control theory that is based on orthogonal complements of twist and wrench spaces," Journal of Robotic Systems, vol. 7, pp. 139-144, 1990. 62

BIOGRAPHICAL SKETCH Yung-Sheng Chang was born in Taipei, the capital city of Taiwan. He received his bachelor s degree in the Department of Electrical and Control Engineering from the National Chiao Tung University in June 2009. He then began working in Ministry of Nation Defense in August 2009. After that, he continued his education for pursuing master s degree in the Department of Mechanical Engineering from the University of Florida and joined the Center for Intelligent Machines and Robotics under the advisement of Dr. Carl Crane in August 2011. 63