[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
Next Article in Journal
An Ink-Jet Printed Eddy Current Position Sensor
Next Article in Special Issue
Flexible PZT Thin Film Tactile Sensor for Biomedical Monitoring
Previous Article in Journal
Performance Evaluation of a Modular Detector Unit for X-Ray Computed Tomography
Previous Article in Special Issue
Research on Bell-Shaped Vibratory Angular Rate Gyro’s Character of Resonator
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Experimental Investigation on Adaptive Robust Controller Designs Applied to Constrained Manipulators

by
Samuel L. Nogueira
1,*,
Tatiana F. P. A. T. Pazelli
2,
Adriano A. G. Siqueira
1 and
Marco H. Terra
3
1
Department of Mechanical Engineering, University of São Paulo, 400 Trabalhador São-Carlense av.,13566-590, São Carlos, Brazil
2
Department of Electrical Engineering, Federal University of São Carlos,Washington Luís, km 235-SP-310, São Carlos, Brazil
3
Department of Electrical Engineering, University of São Paulo, 400 Trabalhador São-Carlense av.,13566-590, São Carlos, Brazil
*
Author to whom correspondence should be addressed.
Sensors 2013, 13(4), 5181-5204; https://doi.org/10.3390/s130405181
Submission received: 5 April 2013 / Revised: 16 April 2013 / Accepted: 16 April 2013 / Published: 18 April 2013
(This article belongs to the Special Issue Piezoelectric Sensors and Actuators)

Abstract

: In this paper, two interlaced studies are presented. The first is directed to the design and construction of a dynamic 3D force/moment sensor. The device is applied to provide a feedback signal of forces and moments exerted by the robotic end-effector. This development has become an alternative solution to the existing multi-axis load cell based on static force and moment sensors. The second one shows an experimental investigation on the performance of four different adaptive nonlinear control methods applied to a constrained manipulator subject to uncertainties in the model and external disturbances. Coordinated position and force control is evaluated. Adaptive procedures are based on neural networks and fuzzy systems applied in two different modeling strategies. The first modeling strategy requires a well-known nominal model for the robot, so that the intelligent systems are applied only to estimate the effects of uncertainties, unmodeled dynamics and external disturbances. The second strategy considers that the robot model is completely unknown and, therefore, intelligent systems are used to estimate these dynamics. A comparative study is conducted based on experimental implementations performed with an actual planar manipulator and with the dynamic force sensor developed for this purpose.

1. Introduction

The first robotic manipulators were developed in order to perform positioning tasks. They were then specifically designed to be robust enough so as to not be affected by external disturbances. This physical robustness of robot manipulators has enabled researchers to obtain accurate positioning systems based on simple control laws. Decades later, the popularization of industrial robotics has heightened researchers' interest in creating a much wider range of applications for robotic manipulators in various environments.

Nowadays, many applications demand robotic manipulators to perform tasks subject to force and motion constraints. For example, the process of milling a piece requires accurate incidence angles, paths, forces and moments exerted by the drill in the milled material. Additionally, in industrial assembly lines, the objects must be assembled along certain paths with predetermined forces and moments. In sheet metal cutting, the cutting angles, paths and forces exerted on the material are also important. Moreover, on surfaces where polishing disks must always be perpendicular to the surface being polished, predetermined force must be applied. Consequently, new concepts of position and force control for lighter and more flexible robots have been created [1,2].

The problem defined in these applications involves three stages: the approach phase, the impact moment and the sustained contact tracking. The approach phase has been addressed in many works and defines the problem of positioning the tool without, or before, touching the environment. The second phase requires controlling the initial impact and damping out the vibrations generated during the event. After the initial impact, sustained contact is desired in many operations. In these cases, not only the motion of the end-effector is required to follow a prescribed path, but also the force exerted by the end-effector is required to follow a predefined reference. In these constrained systems, forces and moments generated between the end-effector and the target must be controlled, rather than being treated as disturbances and rejected. Addressing manipulators subject to model uncertainties and disturbances, the work considered in this paper is concentrated on the sustained contact tracking phase of the problem.

The great majority of solutions presented in the literature for the problem of constrained systems control require knowing the forces and moments of interaction between the robot end-effector and the environment where it is acting. The problem of force control undertaken in this paper was also considered in [1,37]. However, only simulation results were presented.

As a differential and important contribution, this paper addresses an experimental investigation on robust force control as a result of the development of a modular sensor device. The proposed device is designed and built to measure dynamic forces and moments in three orthogonal axes based on unidirectional force sensor units. As a consequence of its independent architecture on the type of sensitive material, static or dynamic force sensors can be applied. Piezoelectric or piezoresistive force sensors are effective solutions for the applications involved in this study due to their inherent dynamic response characteristics.

This paper is organized as follows: next section presents preliminary concepts and relevant results found in the literature; Section 3 introduces the model description of the constrained robot manipulator; Section 4 presents the problem formulation; Section 5 describes the solutions for the nonlinear control problems based on the linear parametrization property of the model, neural networks and fuzzy systems; Section 6 demonstrates the 3D dynamic force and moment sensor; and Section 7 presents the experimental results for a three-link manipulator.

2. Preliminary Concepts

The concept of stiffness control was introduced by Salisbury [8]. It is based on the resistance of the environment in which the robotic end-effector applies the force. The problem is modeled as a mass-spring system and this method made possible the simultaneous position/force control. However, it considers constant desired position and force. In many robotic applications, such as when milling a piece, the end-effector must follow a trajectory along the surface of an object while applying a desired force, which is not necessarily constant. In this case, the stiffness control application does not work properly.

To address this type of limitation, Raibert and Craig [9] partitioned the control problem into two subtasks: one task is for controlling the position trajectory and the other task for controlling the desired force. This approach has been evolutionary for controllers, as proposed by Paul et al. [10], and became the conceptual basis of the hybrid trajectory of position and force control currently found in the literature.

It was shown by McClamroch and Wang [11], that when a manipulator is in contact with a surface, the position degrees of freedom are reduced. In this case, force constraint is added to motion equations through Lagrange multipliers. Thus, the order of the state vector is reduced in the dynamic equations of the manipulator.

Since the phase of controlling robotic systems with restricted position and force has been overcome, researchers began to focus on variables that could degrade the stability of the proposed models. The load on the end-effector that can fluctuate while the manipulator performs various tasks, friction and parameter uncertainties, are some examples that have required much research effort [1215]. However, only a few studies have dealt with adaptive and robust control of robotic systems subject to constraints [1620]. Using the basis developed by [9,11], the work performed by Chang and Chen [21] submitted an adaptive controller with robust performance criterion for robotic systems with position and force constraints.

On the development of adaptive robust controllers, neural networks, combined with a nonlinear control law, were applied by Chang and Chen [22] to control robotic systems subject to parametric uncertainties and external disturbances. A smooth response was achieved with a simple and computationally efficient implementation. Neural networks were applied to estimate the unknown dynamics of the system, thereby not requiring mathematical modeling knowledge. A variable structure controller (VSC) is added to this formulation by Chang and Chen [23]. The inclusion of VSC in the control law weakens the hypothesis used by the authors in [22] that the estimation error should be integrable, and limits it to be only a state-dependent function. Following this approach, the authors in [24] developed an adaptive controller based on fuzzy systems and VSC for robots with position and force constraints.

Karayiannidis and Doulgeri proposed adaptive controllers in [5,6] for force and position trajectory tracking in environments with little or no constraint knowledge. In these papers, the main purpose is to explore the workspace using measurements of position, speed and force in a robotic end-effector, to attenuate impacts caused by unknown environments. In order to improve estimates of these models, the use of a camera was proposed by Cheah et al. [4] to provide a better estimate of the contact constraint through computer vision. However, these works are more interested in exploring the environmental constraint, rather than a repetitive task control with great tracking precision.

To satisfy the need of measuring the forces and moments of interaction between the robot end-effector and the environment, devices known as multi-axis load cell or multi-axis force and moment sensor are used. There are several patents on devices whose purpose is to measure forces and moments in three axes.

Force sensors, such as strain gages, have been generally used as the basic unit of measurement, as shown by the illustrative examples in Figure 1. The multi-axis load cell developed by Meyer and Lowe [25] was built in one piece with internal and external parts that are connected by a pair of axially spaced beams. The beams are fixed in the center of the piece, where the strain gages are attached, and an outside tunnel. The loads are measured by the curvature of the connectors. The load transducer developed by Meyer et al. [26] measures linear forces in three axes and moments of about two axes. The transducer has encapsulations connected by internal and external arms-sensitive loading. The device described by Sommerfeld et al. [27] consists of an external annular body, a hub and four beams that hold the hub to the radially outward portion. Strain gages are fixed to the faces of the beams, with a 90 degree lagging. Forces and moments exerted on the hub are transmitted to the four beams, and consequently, to the strain gages. Differently from [25,26], this device described by Sommerfeld et al. [27] is able to measure forces and moments in three orthogonal axes. However, all the devices mentioned above have an extrusion completely dependent upon the type of sensitive material being used.

3. Model Description

Let a constrained robot manipulator be defined by an n-link serial-chain rigid manipulator whose end-effector is in contact with a rigid environmental constraint, according to Figure 2.

The links of the manipulator are numbered from 1 to n: Ji is the joint connecting the (i − 1)-th and i-th links, li is the vector connecting Ji and Ji+1, qi is the angle of the i-th link around joint Ji, Ci is the center of mass of the i-th link and lci is the vector connecting Ji and Ci. The constraint surface is represented by S and α is the angle of contact between the end-effector and S.

Remark 1 Assume that the end-effector is already in contact with the constraint surface, and the control exerted over the constraint force is such that the force will always maintain the end-effector in contact with the constraint surface.

3.1. Robot Dynamics

The dynamic equations of a constrained robot is given from Lagrange theory as

M ( q ) q ¨ + C ( q , q ˙ ) q ˙ + G ( q ) = τ + f + τ d
where M(q) ∈ ℜn × n is the symmetric positive definite inertia matrix, C(q, q˙) ∈ ℜn × n is the Coriolis and centripetal matrix, G(q) ∈ ℜn is the vector of the gravitational torques, τ ∈ ℜn is the torque vector acting upon the manipulator joint, f ∈ ℜn denotes the vector of joint-space generalized forces on the environmental constraint exerted by the end-effector, and τd defines finite energy unknown disturbances.

Model uncertainties in Equation (1) can be introduced dividing the matrices M(q), C(q, ), G(q), and f into a nominal and a perturbed part:

M ( q ) = M 0 ( q ) + Δ M ( q ) C ( q , q ˙ ) = C 0 ( q , q ˙ ) + Δ C ( q , q ˙ ) G ( q ) = G 0 ( q ) + Δ G ( q ) f = f 0 + Δ f

where M0(q), C0(q, q˙), G0(q), and f0 are nominal matrices and ΔM(q), ΔC(q, q˙), ΔG(q), and Δf represent the uncertainties.

3.2. Constraint Modelling

Considering Remark 1, the m-dimensional surface constraint is described by the holonomic relationship

ϕ ( q ) = 0
where ϕ(q) : ℜnm is a smooth function.

Constraint forces are given by

f = J c T ( q ) λ
where J c ( q ) = δ ϕ ( q ) δ q m × n is the Jacobian matrix that relates the constraint to the controlled variables of the robot and λ ∈ ℜm is a vector of generalized Lagrangian multipliers associated with the constraint.

In this paper, it is considered that parametric uncertainties may also be included into the constraint model since the constraint surface may be not perfectly rigid, frictionless or even that its geometric description may not be exactly known. Thus, consider

ϕ ( q ) = Δ ϕ ( q ) J c ( q ) = J c ( q ) + Δ J c ( q )
and assume that Δϕ and ΔJc are implicit in Δf, described in Equation (2).

3.3. Reduced Order Model

The presence of m constraints causes the manipulator to lose m degrees of freedom, and therefore, nm linearly independent coordinates are sufficient to characterize the constrained movement. Therefore, to formulate a reduced order dynamics for the constrained system, the formulation presented in this paper follows the assumptions made by McClamroch and Wang [11] and posteriorly by Chang and Chen [24].

Define

q = [ q 1 q 2 ]
where q2 = σ(q1) and ϕ(q1, σ(q1)) = 0.

The following reduced model formulation is obtained for the constrained manipulator as the model proposed by Chang and Chen [24]:

A L ( q 1 ) q ¨ 1 + L T ( q 1 ) C L ( q 1 , q ˙ 1 ) q ˙ 1 + L T ( q 1 ) G ( q 1 ) = L T ( q 1 ) ( τ + τ d )
where
L ( q 1 ) = [ I ( n m ) σ ( q 1 ) q 1 ]
and
C L ( q 1 , q ˙ 1 ) q ˙ 1 = M ( q 1 ) L ˙ ( q 1 ) + C ( q 1 , q ˙ 1 ) L ( q 1 )

The necessary structure and properties of the model for controller formulation are maintained since AL(q1) = LT(q1)M(q1)L(q1) is symmetric positive definite and the matrix L(q1) − 2LT(q1)CL(q1, 1) is skew symmetric.

4. Problem Formulation

Let qd(t) ∈ ℜn and d(t) ∈ ℜn be the desired reference trajectory and the corresponding velocity for the joints, respectively Assume that qd(t) and its derivatives d(t) and d(t) are bounded. Define a bounded fd ∈ ℜn as the desired reference contact force. To be consistent with the imposed restrictions, we need to assure that ϕ(qd) = 0 and f d = J c T ( q d ) λ d.

Since q2 = σ(q1), it is only necessary to find a control law that makes q 1 q d 1 when t → ∞. Therefore, define the position tracking error 1(t) and the filtered link tracking error 2(t), as in [24]:

x ¯ ( t ) = ˙ [ x ¯ 1 ( t ) x ¯ 2 ( t ) ] = ˙ [ q 1 ( t ) q d 1 ( t ) q ˙ 1 ( t ) q ˙ d 1 ( t ) + p ( q 1 ( t ) q d 1 ( t ) ) ]
for some constant p > 0.

From Equations (6) and (7), the error dynamic equations can be obtained as

x ¯ ˙ = [ x ¯ ˙ 1 x ¯ ˙ 2 ] = A x ¯ + B u + B w
where
A = [ p I I 0 A L 1 ( q 1 ) C L ( q 1 , q ˙ 1 ) ] , B = [ 0 ( M ( q 1 ) L ( q 1 ) ) 1 ] u = F ( x e ) τ , ω = τ d , and F ( x e ) = ˙ M ( q 1 ) L ( q 1 ) ( q ¨ d 1 p x ¯ ˙ 1 ) + C L ( q 1 , q ˙ 1 ) ( q ˙ d 1 p x ¯ 1 ) + G ( q 1 )

Within this problem formulation, the torques applied to the joints to guarantee the task execution are given by

τ = F ( x e ) u
where the term F(xe) refers to the dynamics of the controlled variables and u is the control law provided by the adaptive controller proposed by Chang and Chen [24].

The term F(xe) can be divided in a nominal and uncertainty part, that is,

F ( x e ) = F 0 ( x e ) + Δ F ( x e )
where term ΔF(xe) contains the uncertainties from parametric Equations (2) and (5).

In this paper, the control problem is solved based on two different approaches. In the first approach, an adaptive intelligent system is applied to estimate only the term ΔF(xe), considering that the nominal model of the robot is well known. In the second approach, it is considered that the system model or the term F0(xe) + ΔF(xe) is completely unknown, then the adaptive intelligent system is applied to estimate it. The nonlinear control and VSC are applied in both approaches to attenuate the effects of estimation errors and external disturbances.

5. Adaptive Nonlinear Controller Based on Intelligent Systems

The adaptive control laws presented as follows are based on two different learning methods to estimate uncertain parameters and also the behavior of unmodeled dynamics. Neural networks and fuzzy systems based on the Takagi–Sugeno model are considered. The structure of these intelligent systems can be seen in details in [28].

Let (xe, Θ) = Ξ(xe)Θ be the output of the adaptive intelligent system, where xe is the input vector and Θ is a vector of adjustable parameters such that ΘTΘ ≤ Mθ|Mθ > 0. Two different cases are considered in the following:

Case 1

Estimation of Model Uncertainties Based on Intelligent Systems

In this case, the nominal model of the robot is considered as well known, the intelligent system then estimates only the model uncertainties, such as

F ^ 1 ( x e , Θ ) Δ F ( x e )

Therefore u and ω in Equation (8) will be rewritten as

u = F 0 ( x e ) + F ^ 1 ( x e , Θ ) τ ω = ( Δ F ( x e ) F ^ 1 ( x e , Θ ) ) + τ d
so that ω includes the estimation error from the intelligent system.

Case 2

Estimation of Complete Model Based on Intelligent Systems

In this case, however, the nominal model of the robot is considered completely unknown, so the intelligent system estimates the complete model, such as

F ^ 2 ( x e , Θ ) F 0 ( x e ) + Δ F ( x e )

Therefore u and ω in Equation (8) is rewritten as

u = F ^ 2 ( x e , Θ ) τ ω = ( F ( x e ) F ^ 2 ( x e , Θ ) ) + τ d
so that ω includes the estimation error from the intelligent system.

Regarding the nonlinear control solution proposed by Chang and Chen [24] for constrained systems, define u = where = uP + uF, such that

u P = k 0 E x ¯ 2 k ( x e ) sgn ( L x ¯ 2 )
u F = J c T λ c
where uP is the +VSC control term for the position enforcement and uF is the control law for the force tracking procedure, with
E : = [ I ( n m ) 0 m × ( n m ) ] and λ c = ˙ λ d k λ 0 T ( λ λ d ) d t
for some constant gain k0, k(xe) > 0 and kλ > 0.

Theorem 1 presented in the following is a variation of the results presented in [21,23,24] with the difference that the nominal model and an intelligent system are considered rather than the linear parametrization of the robot. Thus, this theorem defines an adaptive controller to solve the same problem's position/force tracking control problems presented in these papers but with the advantage of using the known model of the robot.

Theorem 1 Given a desired disturbance attenuation level γ > 0, a weighting matrix Q and the Lyapunov candidate function V(t), the following performance criterion

0 T || x ¯ ( t ) || Q 2 V ( 0 ) + γ 2 0 T | | ω ( t ) | | 2 , T 0
is satisfied, for ω(t) ∈ L2 [0,∞), if there exists a dynamic state feedback controller
Θ ˙ = { ρ T Ξ T L x ¯ 2 if Θ < M θ or ( Θ = M θ and x ¯ 2 T L T Ξ Θ 0 ) ρ T Ξ T L x ¯ 2 + ρ T x ¯ 2 T L T Ξ Θ Θ 2 Θ if Θ = M θ and x ¯ 2 T L T Ξ Θ < 0
τ = F 0 ( x e ) + Ξ Θ k 0 E x ¯ 2 + k ( x e ) sgn ( L x ¯ 2 ) J T λ c
which is the solution of the adaptive nonlinear control problem subject to Equation (6) for the cases 1 and 2.

The stability proof of this dynamic controller follows the line of the proof presented by Chang and Chen [24].

6. Force/Moment Sensor Development

The controllers proposed in Section 5 require the inferred force and moment measures in the robotic end-effector. The most accurate way to obtain such information is by measurements obtained through a multi-axes load cell.

This section presents a detailed description on the modular 3D dynamic force/moment sensor device developed for this application purpose. Its modular architecture is shown in Figure 3.

6.1. Mechanical Description

The 3D dynamic forces/moments sensor essentially comprised of two parts: a moving part and a fixed part. The moving part is shown in yellow in Figure 3(1). It is composed of a base (outside the sensor body) and a body of force transmission (inside the sensor body), which is responsible for transmitting forces and moments to the unidirectional force sensor units mounted on the fixed part, Figure 3(2). The device has an aluminum-made sensor board where twelve force sensor units are disposed in order to measure force and moment in all directions. The sensor units are arranged in pairs, for example, sensors 1A and 1B as shown in Figure 4, which illustrates the body of forces transmission (4) and identifies the sensor units.

Each pair of sensor units measures the force in one direction. This setting is tensioned proportionally to the movement performed by (4) in that axis. Therefore, forces and moments normal to the plane of the base are calculated by the composition of measured forces in the three axes. Figure 5 demonstrates the sensor units operation schematic, and identifies the force measurements.

The formulation of the resultant forces/moments composition is obtained as follows.

Three pairs of sensor units, (1A,1B), (2A,2B) and (3A,3B), define, respectively, the force measurements (F1A,F1B), (F2A,F2B) and (F3A,F3B). The resultant forces of these three points can be expressed as

F 1 = F 1 A F 1 B
F 2 = F 2 A F 2 B
F 3 = F 3 A F 3 B

It can be seen in Figure 5 that F1, F2, and F3 form an equilateral triangle, whose side is given by d1. In [29], Doebelin demonstrated that if three load cells, measuring force in one direction, are arranged triangularly as the vertices of the triangle (the measurement direction normal to the plane formed by the triangle), the following physical quantities can be calculated:

F x = F 1 + F 2 + F 3
M y = d 1 F 1 F 2 ( F 3 F 1 ) 2 3
M z = d 1 ( F 3 F 2 ) 2
where Fx, My, and Mz are the force along x-axis and the moments along axes y and z, respectively.

In the same way, F4, F5, and F6 also form an equilateral triangle, whose side is given by d2. Thus, the resultant forces of these three points can be expressed as

F 4 = F 4 A F 4 B
F 5 = F 5 A F 5 B
F 6 = F 6 A F 6 B
and therefore,
F y = F 5 F 4 ( F 4 F 6 ) 2
F z = 3 ( F 5 F 6 ) 2
M x = d 2 ( F 4 + F 5 + F 6 ) 2 3
where Fy, Fz, and Mx are the forces along axes y and z and the moment along axis x, respectively. The sensor parts and the assembled device is shown in Figure 6.

The applied sensor unit provides a measurement range of 0 to 1.5 kgf, or a maximum load of 14.7 N. By means of Equations (18)(20) and (24)(26), the built sensor presents measurements characteristics as shown in Table 1.

6.2. Electronic Description

The signal from each pair of unidirectional sensors are conditioned by an electronic circuit based on Instrumentation Amplifiers (INA). The electrical schematic of the signal conditioning circuit is summarized in Figure 7.

The sensors are powered by a symmetric supply (+ Vcc and −Vcc) and provide differential output signals (for example, F1Ax and F1Ay).

Six INA comprise the printed circuit board (PCB). Each INA is responsible for conditioning output signals from a pair of force sensors. With this configuration it became possible to measure the differential forces applied on the moving part of the proposed device. The resultant force is then given by the following equation as described in Figure 7.

F I = k F I A k F I B

The signals generated by the INAs are directed to a data acquisition board that interfaces with the computer.

7. Experimental Results

The designed 3D dynamic forces/moments sensor operation is validated in a set of experimental applications, where a three-link planar manipulator is subject to a holonomic constraint defined by a straight rule. The proposed intelligent adaptive robust controllers are applied to this plant for position and force trajectory tracking. The experimental manipulator UARM (UnderActuated Robot Manipulator), whose nominal parameters are given in Table 2, is composed of a DC motor in each joint, a break and an optical encoder with quadrature decoding used to measure joint positions, Figure 8(a). Joint velocities are obtained by numerical differentiation and filtering [30,31]. Modeling matrices for the UARM, M(q), C(q, q˙), and G(q), can be seen in [30]. In Figure 8(b), the designed force sensor is coupled to the UARM end-effector which is constrained in a straight rule. A graphical user interface, whose software was developed using MATLAB® platform, is used to implement control laws and interface with the experimental setup. It was implemented in a modular form, since any new controller can be easily developed and applied.

7.1. Implementation of Control Law

The constraint surface for the robot end-effector is a segment of a straight line on the X–Y plane. The angle α between the end-effector and the constraint line is defined in this application example as π/2. It means that the orientation must remain in a constant value c given through the line inclination β and α, where c = π + atan(β) − α. Hence, the equation of the m = 2 constraints is given by

ϕ ( q ) = [ l 1 s 1 l 2 s 12 l 3 s 123 + β [ l 1 c 1 + l 2 c 12 + l 3 c 123 ] + b q 1 + q 2 + q 3 c ] = [ 0 0 ]
where b is the linear coefficient of the constraint line, sij = sin(qi + qj), cij = cos(qi + qj). Hence, ϕ: ℜ32, and the Jacobian matrix, Jc(q) = ∂ϕ/∂q, is given by
J c = [ J c 11 J c 12 J c 13 J c 21 J c 22 J c 23 ]
with
J c 11 = l 1 c 1 + l 2 c 12 + l 3 c 123 + β [ l 1 s 1 + l 2 s 12 + l 3 s 123 ] J c 12 = l 2 c 12 + l 3 c 123 + β [ l 2 s 12 + l 3 s 123 ] J c 13 = l 3 c 123 + β [ l 3 s 123 ] J c 21 = J c 22 = J c 23 = 1

Defining q1 = [q1] and q2 = [q2 q3], the matrix L(q) of the constraint line is

L ( q ) = [ 1 [ l 1 cos ( q 1 ) + l 2 cos ( q 1 + q 2 ) + β ( l 1 sin ( q 1 ) + l 2 sin ( q 1 + q 2 ) ) ] l 2 [ cos ( q 1 + q 2 ) + β sin ( q 1 + q 2 ) ] [ l 1 cos ( q 1 ) + l 2 cos ( q 1 + q 2 ) + β ( l 1 sin ( q 1 ) + l 2 sin ( q 1 + q 2 ) ) ] l 2 [ cos ( q 1 + q 2 ) + β sin ( q 1 + q 2 ) ] 1 ]

Initial and final coordinates of the movement are (x0,y0) = (0.46,0.38) m and (x(T),y(T)) = (0.53, 0.13) m, respectively. In this case, β = −3.57, b = 2.02, and c = 15.6°. The reference trajectory for the joint variables q d ( t ) = q d 1 ( t ) is a fifth-degree polynomial, with trajectory duration time T = 4 s. It is desired that no force acts on the normal direction of the constrain line and no moment acts on the z-direction, that is, λd = [(Fx)d (Mz)d]T = [0 0]T.

We use as benchmarking for comparison of the proposed controllers, a controller with no intelligent adjust adaptation, that is, only the nominal parameters are used. In other words, we use the controller of Equation (14) without the ΞΘ term.

During the experiment, a limited disturbance was introduced at ti = 1.0s in the following form:

τ d = [ 0 , 01 e ( t t d ) 2 2 μ 2 sen ( 3 , 6 π t ) 0 , 01 e ( t t d ) 2 2 μ 2 sen ( 2 , 7 πt ) 0 , 01 e ( t t d ) 2 2 μ 2 sen ( 1 , 8 π t ) ]

If compared with the nominal torque, the disturbance τd is approximately 64% of its peak value, see Figure 9.

The selected gains are defined in Table 3.

7.1.1. Neural Networks Configuration

For the nonlinear controllers via neural network, Cases 1 and 2, proposed in Section 5, let n = 3 be the number of joints of manipulator. Define

F ^ NNi ( x e , Θ ) = [ F ^ NNi 1 ( x e , Θ 1 ) F ^ NNi 2 ( x e , Θ 2 ) F ^ NNi 3 ( x e , Θ 3 ) ] T = Ξ Θ
with p1 = p2 = p3 = 7 neurons in the hidden layer and the bias vector b1 = b2= b3 = [−3 −2 −1 0 1 2 3]. Using the reduced model, the input vector xe is defined by x e = [ q 1 q ˙ 1 q d 1 q ˙ d 1 q ¨ d 1 ]. In order to provide an input related to the error values, the weighting matrix for the first layer of the neural networks is defined by W i 1 = W i 2 = W i 3 = [ W i j 1 ] = [ W i j 2 ] = [ W i j 3 ] = [ 1 1 1 1 1 ]. The activation function of the neurons in the hidden layer is chosen as the hyperbolic tangent, G(.) = tanh(.). The uncertain vector Θ is defined as Θ = [ Θ 1 T Θ 2 T Θ 3 T ] T and the matrix Ξ = [ ξ 1 T ξ 2 T ξ 3 T ] T as Pazelli et al. and Siqueira et al. showed in [28,30]. Θ values are updated at each control iteration by the adaptive law given in Equation (13). The neural networks outputs are given by
F ^ NNi k ( x e , Θ k ) = i = 1 p k θ k i G ( j = 1 q k w ij k x e j + b i k ) = ξ k T Θ k

Experimental results for the Case 1 (NN1) with neural network are shown in Figures 10(a), 11(a), and 12(a). For the Case 2 with neural network plus nominal model (NN2), experimental results are shown in Figures 10(b), 11(b), and 12(b).

Remark 2 Although NN1 and NN2 use the same structure and the same input data, their outputs are different. The online adaptation law provides the update of Θ values in distinct ways for each case. The contribution of NN1 output to the final torque signal is lower than the contribution of NN2 output, since NN1 estimates the value of ΔF(xe) and NN2 estimates the value of F(xe) + ΔF(xe). Thus, Θ values must achieve that requirement for each case.

7.1.2. Fuzzy Systems Configuration

For the proposed adaptive fuzzy nonlinear controllers, Cases 1 and 2 (Section 5), a set of n = 3 fuzzy systems may be defined by

F ^ F Si ( x e , Θ ) = [ F ^ F Si 1 ( x e , Θ 1 ) F ^ F Si 2 ( x e , Θ 2 ) F ^ F Si 3 ( x e , Θ 3 ) ] T = Ξ Θ
where FSi1(.,.), FSi2(.,.), and FSi3(.,.) correspond to the estimate of the uncertain part of the dynamic behavior of joints 1, 2, and 3, respectively. The input vector xe is defined as x e = [ q ˜ 1 q ˜ ˙ 1 ]. The fuzzy sets A ( x e ) = [ A 1 ( q ˜ 1 ) A 2 ( q ˜ ˙ 1 ) ] are defined according to Figure 13 and the number of linguistic variables in U1 and U2 are defined as r1 = r2 = 3. They are applied to both FSi1(.,.), FSi2(.,.) and Fsi3(.,.) for the universe of discourse of position errors, u 2 1 = u 1 2 = u 1 3 = q ˜ 1 U 1 1 = U 1 2 = U 1 3 = U 1, and for the universe of discourse of velocity errors, u 2 1 = u 2 2 = u 2 3 = q ˜ ˙ 1 U 2 1 = U 2 2 = U 2 3 = U 2.

During the control loop, these graphics are used to determine the grade of membership associated to xe, which defines the value of μ in the output of the Takagi–Sugeno fuzzy model; for more details see [28],

F ^ F Si k ( x e , Θ k ) = i = 1 p k μ i k Y i k i = 1 p k μ i k = i = 1 p k μ i k ( θ i 0 k + θ i 1 k u 1 + θ i 2 k u 2 ) i = 1 p k μ i k = ξ k T Θ k

A fuzzy rule base is defined with pk = r1r2 = 9 rules as

R i : IF ( u 1 is A r 1 i 1 ) and ( u 2 is A r 2 i 2 ) THEN Y i

The vector of adjustable components Θ is defined as Θ = [ Θ 1 T Θ 2 T Θ 3 T ] T, with

Θ 1 T = [ θ 10 1 θ 11 1 θ 12 1 θ 90 1 θ 91 1 θ 92 1 ] Θ 2 T = [ θ 10 2 θ 11 2 θ 12 2 θ 90 2 θ 91 2 θ 92 2 ] Θ 3 T = [ θ 10 3 θ 11 3 θ 12 3 θ 90 3 θ 91 3 θ 92 3 ]
and are updated at each control iteration by the adaptive law in Equation (13).

The matrix Ξ = [ ξ 1 T ξ 2 T ξ 3 T ] T is computed with

ξ 1 T = [ ξ 10 1 ξ 11 1 ξ 12 1 ξ 90 1 ξ 91 1 ξ 92 1 ] ξ 2 T = [ ξ 10 2 ξ 11 2 ξ 12 2 ξ 90 2 ξ 91 2 ξ 92 2 ] ξ 3 T = [ ξ 10 3 ξ 11 3 ξ 12 3 ξ 90 3 ξ 91 3 ξ 92 3 ]

Figures 10(c), 11(c), and 12(c) show the experimental results for Case 3. Figures 10(d), 11(d), and 12(d) show the experimental results of fuzzy system plus nominal model for Case 4.

7.2. Results Discussion

Figure 10 shows the applied torques in the robotic manipulator joints. Note that for the range of the disturbance inserted, the controllers act strongly on the system, reversing the directions of the torques of joints 1 and 2 and increasing the torque of joint 3. In this analysis, it appears that the controllers using neural networks, NN1 and NN2, oscillate less, therefore, they have a slightly slower response time.

The charts of Figure 11 are important because they show the behavior of the end effector during the controller actions, and represent the measurements of force and moment in the robotic end-effector. As expected, it was observed that during the period in which the disturbance appeared in the system, there was a higher intensity of forces and moments which decreased gradually, tending to zero, until the final time (4 s). The proposed configuration for the device sensor proves its effectiveness in providing these data with the required accuracy and without delay to the implementation of the control laws under study.

Figure 12 shows the tracking angles of three joints of the robotic manipulator. Notice that the largest deviation is in the range where the disturbance appeared. Nonetheless, the experiment showed good tracking of the trajectories desired. Again, we can see that for the controller without the intelligent system, NOM, the control actuation is slower, providing greater tracking error of reference signals.

Three performance indexes are used to compare the nonlinear controllers: the 2 norm of the state vector

2 [ x ˜ ] = ( 1 ( t r t 0 ) t 0 t r | | x ˜ ( t ) | | 2 2 d t ) 1 / 2
where ‖ · ‖2 is the Euclidean norm, the sum of the applied torques
E [ τ ] = i = 1 3 ( t 0 t r | τ i ( t ) | d t )
and the sum of areas of contact forces
E [ λ ] = i = 1 3 ( t 0 t r | λ i ( t ) | d t )
where λi(t) the ith component of the contact forces. As the desired values of contact forces are zero, the lower the value of E [λ], the better the controller will be with respect to contact force control.

The results that are shown in Table 4 present the mean value of five experiments.

From Table 4 we can conclude that the controller based only on the nominal model presents a higher state error and greater forces. Although the fuzzy controllers exhibit the best performance, this performance difference with respect to controllers based on neural networks deserves more in-depth study, which will be carried out in future works.

8. Conclusions

In this work, the control problem of trajectory tracking with guaranteed performance was considered for manipulators with force and position constraint. Five controllers were evaluated:

  • The first controller, called nominal, considers that the term F(xe) is completely known, i.e., it does not take into account parametric uncertainties of the model.

  • The second and third controllers have neural network-based estimators, where the second considers the known nominal model of the manipulator and estimates only parameter uncertainties, and the third estimates the full model.

  • The fourth and fifth controllers have fuzzy logic-based estimators which estimate only parametric uncertainties and the full model, respectively.

According to experimental results presented in Section 7, controllers NN1 and FS1 have better tracking trajectory compared, respectively, with the NN2 and FS2. Moreover, it also appears they are more stable because of the lower oscillations after the introduction of disturbance, and this feature occurs for the trajectory tracking in Cartesian space in order to follow the trajectory of joint angles. Comparing the results of intelligent systems it is noted that the controllers via neural networks (NN2 and NN1) have a smoother control action and, consequently, have smaller oscillations than the fuzzy systems (FS1 and FS2), which can be verified by comparing the graphics of torque (Figure 10). On the other hand, the controllers NN1 and NN2 are slower than the FS1 and FS2, hence they tend to have a greater tracking error. In addition, when the system is subject to disturbances, all controllers respond well.

Analyzing the proposed performance indices, it was shown that the state error in the fuzzy system-based controllers (FS1 and FS2) tend to be smaller than those based on neural networks (NN2 and NN1). A possible explanation for this is due to the fact that fuzzy systems acted faster than neural networks. It was observed that choosing values of kλ more directly influence the adjustment of the forces, while the values of ρ influence the tracking errors in the state variables.

The development of the proposed sensor was of crucial importance for the realization of the experimental analysis, a differential of this paper. The measurement results shown proved the effectiveness of the proposed structure. Still, its modular feature enables the application of sensor units of various kinds, such as piezoelectric sensors.

As future work, a comparative study will be performed between the developed sensor and commercially available ones. Furthermore, the use of genetic algorithms to adjust controller gains is a current and ongoing project.

Acknowledgments

This work was supported by Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) and Fundação de Apoio à Pesquisa do Estado de São Paulo (FAPESP), under grant 2012/05552-9.

References

  1. Kilicaslan, S.; Özgören, M.K.; Ider, K. Hybrid force and motion control of robots with flexible links. Mech. Mach. Theory 2010, 45, 91–105. [Google Scholar]
  2. Lewis, F.L.; Abdallah, C.T.; Dawson, D.M. Control of Robot Manipulators; MacMillan Publishing Company: New York, NY, USA, 1993. [Google Scholar]
  3. Bechlioulis, C.P.; Doulgeri, Z.; Rovithakis, G.A. Neuro-adaptive force/position control with prescribed performance and guaranteed contact maintenance. IEEE Trans. Neural Netw. 2010, 21, 1857–1868. [Google Scholar]
  4. Cheah, C.C.; Hou, S.P.; Zhao, Y.; Slotine, J.J.E. Adaptive vision and force tracking control for robots with constraint uncertainty. IEEE/ASME Trans. Mechatron. 2010, 15, 389–399. [Google Scholar]
  5. Karayiannidis, Y.; Doulgeri, Z. Adaptive control of robot contact tasks with on-line learning of planar surfaces. Automatica 2009, 45, 2374–2382. [Google Scholar]
  6. Karayiannidis, Y.; Doulgeri, Z. Robot contact tasks in the presence of control target distortions. Robot. Auton. Syst. 2010, 58, 596–606. [Google Scholar]
  7. Kumar, N.; Panwar, V.; Sukavanam, N.; Borm, S.P.S.J.H. Neural network based hybrid force/position control for robot manipulators. Int. J. Precis. Eng. Manuf. 2011, 12, 419–426. [Google Scholar]
  8. Salisbury, J.K. Active Stiffness Control of a Manipulator in Cartesian Coordinates. Proceedings of the 19th IEEE Conference on Decision and Control including the Symposium on Adaptive Processes, Albuquerque, NM, USA, 10–12 December 1980; Volume 19, pp. 95–100.
  9. Raibert, M.; Craig, J. Hybrid position/force control of manipulators. ASME J. Dyn. Syst. Meas. Control 1981, 102, 126–133. [Google Scholar]
  10. Paul, R.; Shimano, B. Compliance and Control. Proceedings of the Joint Automatic Control Conference, Purdue, IN, USA, 27–30 July 1976; pp. 694–699.
  11. McClamroch, N.H.; Wang, D. Feedback stabilization and tracking of constrained robots. IEEE Trans. Autom. Control 1988, 33, 419–426. [Google Scholar]
  12. Ortega, R.; Spong, M.W. Adaptive motion control of rigid robots: A tutorial. Automatica 1989, 25, 877–888. [Google Scholar]
  13. Slotine, J.J.E.; Li, W. Composite adaptive control of robot manipulators. Automatica 1989, 25, 509–519. [Google Scholar]
  14. Slotine, J.J.E.; Li, W. Applied Nonlinear Control; Prentice-Hall: Englewoods Cliffs, NJ, USA, 1991. [Google Scholar]
  15. Spong, M.W.; Vidyasagar, M. Robot Dynamics and Control; John Wiley and Sons: New York, NY, USA, 1989. [Google Scholar]
  16. Carelli, R.; Kelly, R. An adaptive impedance/force controller for robot manipulators. IEEE Trans. Autom. Control 1991, 36, 967–971. [Google Scholar]
  17. Jean, J.H.; Fu, L.C. Adaptive hybrid control strategies for constrained robots. IEEE Trans. Autom. Control 1993, 38, 598–603. [Google Scholar]
  18. Panteley, E.V.; Stotsky, A.A. Adaptive trajectory/force control scheme for constrained robot manipulators. Int. J. Adapt. Control Signal Process. 1993, 7, 489–496. [Google Scholar]
  19. Su, C.Y.; Stepanenko, Y.; Leung, T.P. Combined adaptive and variable structure control for constrained robots. Automatica 1995, 31, 483–488. [Google Scholar]
  20. Su, C.Y.; Leung, T.P.; Zhou, Q.J. Force/motion control of constrained robots using sliding mode. IEEE Trans. Autom. Control 1992, 37, 668–672. [Google Scholar]
  21. Chang, Y.C.; Chen, B.S. Adaptive tracking control design of constrained robot systems. Int. J. Adapt. Control Signal Process. 1998, 12, 495–526. [Google Scholar]
  22. Chang, Y.C.; Chen, B.S. A nonlinear adaptive tracking control design in robotic systems via neural networks. IEEE Trans. Control Syst. Technol. 1997, 5, 13–28. [Google Scholar]
  23. Chang, Y.C. Neural network-based tracking control for robotic systems. IEEE Proc. Control Theory Appl. 2000, 147, 303–311. [Google Scholar]
  24. Chang, Y.C.; Chen, B.S. Robust tracking designs for both holonomic and nonholonomic constrained mechanical systems: Adaptive fuzzy approach. IEEE Trans. Fuzzy Syst. 2000, 8, 46–66. [Google Scholar]
  25. Meyer, R.A.; Lowe, A.E. Multiple Axis Load Sensitive Transducer. U.S. Patent 4640138, 1987. [Google Scholar]
  26. Meyer, R.A.; Olson, D.J. Load Transducer. U.S. Patent 4821582, 1989. [Google Scholar]
  27. Sommerfeld, J.L.; Meyer, R.A.; Larson, B.A.; Olson, D.J. Multi-Axis Load Cell. U.S. Patent 5969268, 1999. [Google Scholar]
  28. Pazelli, T.F.P.A.T.; Terra, M.H.; Siqueira, A.A.G. Experimental investigation on adaptive robust controller designs applied to a free-floating space manipulator. Control Eng. Pract. 2011, 9, 395–408. [Google Scholar]
  29. Doebelin, E.O. Measurement Systems Application and Design; McGraw-Hill Science: New York, NY, USA, 1989; pp. 410–413. [Google Scholar]
  30. Siqueira, A.A.G.; Terra, M.H.; Bergerman, M. Robust Control of Robots: Fault Tolerant Approaches; Springer: London, UK, 2011. [Google Scholar]
  31. Siqueira, A.A.G.; Terra, M.H. Nonlinear and Markovian controls of underactuated manipulators. IEEE Trans. Control Syst. Technol. 2004, 12, 811–826. [Google Scholar]
Figure 1. (a) Image taken from [25]; (b) Image taken from [26]; (c) Image taken from [27].
Figure 1. (a) Image taken from [25]; (b) Image taken from [26]; (c) Image taken from [27].
Sensors 13 05181f1 1024
Figure 2. Constrained robot manipulator.
Figure 2. Constrained robot manipulator.
Sensors 13 05181f2 1024
Figure 3. Design of the 3D dynamic forces/moments sensor: (1) Moving part; (2) Fixed part.
Figure 3. Design of the 3D dynamic forces/moments sensor: (1) Moving part; (2) Fixed part.
Sensors 13 05181f3 1024
Figure 4. Body of force transmission (4) and arrangement of sensor units: (a) Top view; (b) Bottom view.
Figure 4. Body of force transmission (4) and arrangement of sensor units: (a) Top view; (b) Bottom view.
Sensors 13 05181f4 1024
Figure 5. Schematic of sensor operation.
Figure 5. Schematic of sensor operation.
Sensors 13 05181f5 1024
Figure 6. 3D dynamic forces/moments sensor.
Figure 6. 3D dynamic forces/moments sensor.
Sensors 13 05181f6 1024
Figure 7. Electronic diagram.
Figure 7. Electronic diagram.
Sensors 13 05181f7 1024
Figure 8. (a) UnderActuated Robot Manipulator; (b) Force sensor device coupled to the UARM end-effector.
Figure 8. (a) UnderActuated Robot Manipulator; (b) Force sensor device coupled to the UARM end-effector.
Sensors 13 05181f8 1024
Figure 9. Torque disturbance.
Figure 9. Torque disturbance.
Sensors 13 05181f9 1024
Figure 10. Joint Torque. (a) NN1; (b) NN2; (c) FS1; (d) FS2; (e) NOM.
Figure 10. Joint Torque. (a) NN1; (b) NN2; (c) FS1; (d) FS2; (e) NOM.
Sensors 13 05181f10a 1024Sensors 13 05181f10b 1024
Figure 11. Force and Moment. (a) NN1; (b) NN2; (c) FS1; (d) FS2; (e) NOM.
Figure 11. Force and Moment. (a) NN1; (b) NN2; (c) FS1; (d) FS2; (e) NOM.
Sensors 13 05181f11a 1024Sensors 13 05181f11b 1024
Figure 12. Joint Angles. (a) NN1; (b) NN2; (c) FS1; (d) FS2; (e) NOM.
Figure 12. Joint Angles. (a) NN1; (b) NN2; (c) FS1; (d) FS2; (e) NOM.
Sensors 13 05181f12a 1024Sensors 13 05181f12b 1024
Figure 13. Fuzzy sets A1(1) and A2(1).
Figure 13. Fuzzy sets A1(1) and A2(1).
Sensors 13 05181f13 1024
Table 1. Maximum load.
Table 1. Maximum load.
Axis XAxis YAxis Z
Force (N)44.129.425.5
Moment (N/m)d212.7d117d114.7
Table 2. UARM Parameters.
Table 2. UARM Parameters.
BodymiIililci

(kg)(kgm2)(m)(m)
Link 10.8500.00750.2030.096
Link 20.8500.00750.2030.096
Link 31.7000.09000.2400.177
Table 3. Selected Gains.
Table 3. Selected Gains.
Pk0kλρ
[ 2.50 0 0 0 2.50 0 0 0 2.50 ] [ 0.35 0 0 0 0.35 0 0 0 0.35 ] [ 2.00 0 0 0 2.00 0 0 0 2.00 ] [ 0.75 0 0 0 0.75 0 0 0 0.75 ]
Table 4. Performance Indexes.
Table 4. Performance Indexes.
2[x]E[τ]E[λ]

(m)(Nms)(Ns)
Nominal0.10310.65080.1506
NN10.07800.65640.1435
NN20.10260.43180.1156
FS10.07590.58650.1136
FS20.08990.51820.1061

Share and Cite

MDPI and ACS Style

Nogueira, S.L.; Pazelli, T.F.P.A.T.; Siqueira, A.A.G.; Terra, M.H. Experimental Investigation on Adaptive Robust Controller Designs Applied to Constrained Manipulators. Sensors 2013, 13, 5181-5204. https://doi.org/10.3390/s130405181

AMA Style

Nogueira SL, Pazelli TFPAT, Siqueira AAG, Terra MH. Experimental Investigation on Adaptive Robust Controller Designs Applied to Constrained Manipulators. Sensors. 2013; 13(4):5181-5204. https://doi.org/10.3390/s130405181

Chicago/Turabian Style

Nogueira, Samuel L., Tatiana F. P. A. T. Pazelli, Adriano A. G. Siqueira, and Marco H. Terra. 2013. "Experimental Investigation on Adaptive Robust Controller Designs Applied to Constrained Manipulators" Sensors 13, no. 4: 5181-5204. https://doi.org/10.3390/s130405181

APA Style

Nogueira, S. L., Pazelli, T. F. P. A. T., Siqueira, A. A. G., & Terra, M. H. (2013). Experimental Investigation on Adaptive Robust Controller Designs Applied to Constrained Manipulators. Sensors, 13(4), 5181-5204. https://doi.org/10.3390/s130405181

Article Metrics

Back to TopTop