Journal of Control Science and Engineering

Journal of Control Science and Engineering / 2012 / Article
Special Issue

Adaptive Control Theory and Applications

View this Special Issue

Research Article | Open Access

Volume 2012 |Article ID 365067 |

Satoshi Suzuki, Katsuhisa Furuta, "Adaptive Impedance Control to Enhance Human Skill on a Haptic Interface System", Journal of Control Science and Engineering, vol. 2012, Article ID 365067, 10 pages, 2012.

Adaptive Impedance Control to Enhance Human Skill on a Haptic Interface System

Academic Editor: Lili Ma
Received05 Dec 2011
Revised09 Mar 2012
Accepted30 Mar 2012
Published29 May 2012


Adaptive assistive control for a haptic interface system is proposed in the present paper. The assistive control system consists of three subsystems: a servo controller to match the response of the controlled machine to the virtual model, an online identifier of the operator’s control characteristics, and a variable dynamics control using adaptive mechanism. The adaptive mechanism tunes an impedance of the virtual model for the haptic device according to the identified operator’s characteristics so as to enhance the operator’s control performance. The adaptive law is derived by utilizing a Lyapunov candidate function. Using a haptic interface device composed by a xy-stage, an effectiveness of the proposed control method was evaluated experimentally. As a result, it was confirmed that the operator’s characteristics can be estimated sufficiently and that performance of the operation was enhanced by the variable dynamics assistive control.

1. Introduction

An impedance control is a key technology of the force/motion control for any mechanical systems such as an active vehicle suspension, a power steering system, a machining and handling by manipulators, and a tele-operation system. Since dynamics of such controlled mechanism can be adjusted by changing the virtual impedance model, this method is effective to adapt to an ever-changing environment and conditions. Also a biological system has acquired similar strategy of the variable impedance control in the course of an evolution. It is known that an impedance of a musculoskeletal system is changed dynamically during walking, running, and moving the hand [1]. Therefore, variable impedance methods have been studied for artificial legs/orthosis [2] and material-handling machines [3, 4]. Parameters of those impedance control methods are, however, often tuned empirically and intuitively; hence, the system designers have to adjust them according to individuals. Due to this issue, users sometimes have to adapt themselves to the controlled machine when the tuning condition given by the designer is not adequate for the user. To resolve this paradox, the following approach is ideal; the control characteristics of each user are identified during the operation, and then control of the machine is adjusted adaptively according to the identified user’s characteristics. While several similar approaches concerning online variable impedance control are reported, troublesome processes such as a training phase [5] or an empirical tuning for different types of motion [6] are required. Since these approaches are not real adaptive control, a realization of an automatic tuning mechanism without intervention of system designers is expected. Therefore, the present paper presents a design procedure of true adaptive variable impedance control for an assisting system which is considered with the following properties based on the previous method presented in [7]:(a)adaptivity to control characteristics of individual user,(b)derivation of adaptive law for variable impedance tuning based on an adaptive control theory.

Main purpose of the control design proposed here is a development of an adaptive tuning law of the machine dynamics, of which parameters are fixed in an ordinary mechatronics system, in order to enhance manipulation performance of whole of a human-machine system. And, the purposes of the present study are as follows:(c) experimental evaluation,(d) confirmation of the benefit and issue.

Item (a) is realized with an on-line identification based on an assumed model of the human controller. Concerning item (b), an adaptive control law to adjust impedance parameters of the virtual model is derived to ensure the stability and performance of the whole human-machine system. Item (c), evaluation, was performed using a haptic interface device through a point-to-point operation task. Issues and analysis denoted at item (d) are discussed based on the results of the experiment.

This paper is organized as follows. In Section 2, a concept of the adaptive impedance control and its background are mentioned. The haptic interface system and task which were used in the experiment are explained there. Section 3 explains a procedure of the presented assistive control, and its theoretical proof is given there. Section 4 shows results of the experiment and analysis to confirm the effectiveness of the presented method. Last Section 5 presents the conclusion.

2. Human Assistive System

2.1. Concept of Adaptive Impedance Control

In order to design an assistive mechanism for user’s manipulation in a human-in-the-loop system, an adequate human modeling and a feasible assistive control are required. Human modeling has been studied in the field of control engineering from its early beginnings such as a linear servo control model [8], a PID-base time-variant model [9], and an optimal control model [10]. Those models can express a human behavior well for each assumed situation; however, it is inadequate to explain the learning process of a user from the beginner phase to the expert phase. To find adequate model which can treat a human adaptability, it is adequate to refer the voluntary motion control of a musculoskeletal system. The reason is that such model is formulated to explain a process of a human development, and most popular model is a feedback-error-learning model [11, 12]. This model can be utilized to explain the learning of an external unknown dynamics since such an external system can be thought as an extension of our body. On the learning process of the external dynamics, a delay has to be considered because it concerns the stability and performance of whole human-machine system. The delay arises certainly at the visual processing and at the neural transmission between a brain and sensory receptor/muscles. Such undesirable effect given by the delay is compensated by an internal feedback compensation using an efferent neuron and by a delay compensation mechanism which is explained by the Smith predictor theory [13]. Additionally, as shown in Kleinman’s research of the dynamics of a pilot [10], a human (i.e., pilot) has a high ability to compensate the delay in the response of a vehicle. Hence, if a time-delay effect inside the machine side system does not change, it is expected that the delay in the human-machine system can be compensated relatively easily. Therefore, the block diagram shown in Figure 1 appears adequate for designing of a human-machine assist system. The model mentions that human learns the unknown dynamics of the machine and uses the identified model as an inverse model for the manipulation of the machine.

The concept shown in Figure 1 indicates that difference of the machine’s dynamics affects indirectly to the learning of the operator. If the operated machine can change its own dynamics characteristics so as to be learned by individual operator without difficulty, performance of the operation would be enhanced. Therefore, in order to enhance the operator’s performance, as shown in Figure 2, an original dynamics of the machine is replaced to a virtual dynamics model from the operator’s side by making a local loop feedback with a virtual internal model control. In short, the impedance of the virtual dynamics model is modified so as to decrease an error which relates with each task performance. To summarize this discussion, the following three functions are required to realize aforementioned adaptive impedance control.

Step 1. Virtual internal model (VIM) control.

Step 2. Online identification of the operator.

Step 3. Adaptive mechanism to tune the VIM.

The VIM control for Step 1 is realized by making a local servo system that tracks the output of a virtual impedance model. The servo control input law is designed using Linear-Quadratic Regulator (LQR). Identification for Step 2 is performed by assuming a parametric model of operator’s control characteristics. Concerning Step 3, the adaptive mechanism is designed by changing the impedance parameters of the VIM obtained at the Step 3 after derivation of the adaptive law of the VIM model based on a Lyapunov-like function. Details of this process are explained in Section 3.

2.2. Experimental System

A haptic interface system, which is shown in Figure 3, was used to evaluate the adaptive assistive control presented in Section 3. The haptic device consists of a two degree-of-freedom planer 𝑥𝑦-stage, produced by NSK corporation, and a real-time CG monitor programmed by visual C++. The 𝑥𝑦-stage is driven by two linear direct drive motors. The operator moves a grip attached to the 𝑥𝑦-stage and a pointer displayed on the CG monitor (see Figure 4) is also moved according to the position of the grip. Operator’s hand force is measured by a 6-axis force sensor embedded between the grip and the stage. The 𝑥- and 𝑦-axes of the stage do not effect each other because of a mechanically independent design. Computations of the control of the stage and the CG displaying are executed by a PC/AT 3 GHz computer under real-time scheduling control. The control interval is 2 ms, and the movable range is about 62mm in both 𝑥 and 𝑦 directions.

2.3. Task of Manipulation

Since a point-to-point (PTP) manipulation is a popular task both in a daily life and in an industrial situation, the PTP-task was adopted for verification of the presented method. The PTP tasks were repeated by changing the target’s position at random so as to keep the distance constant from each last target to the next one. As soon as the target is displayed on the monitor, the operator moves the pointer to the target by manipulating the grip of the 𝑥𝑦-stage. When position of the pointer is kept inside the target circle for 3 seconds, one PTP motion (one trial) is finished, and then a new next target circle is displayed at random. To reduce the fatigue of the participant, ten-second rest was given to the participant after every five trials.

3. Design of Adaptive Impedance Control

3.1. Virtual Internal Model Control (Step 1)

A procedure to apply the virtual internal model (VIM) control [14] to the 𝑥𝑦-stage is mentioned in this section. Any mechanical mechanism includes nonlinearity caused by friction, variances of viscosity, and unknown dynamics; hence, it is difficult to apply a linear system control theory to actual machine without any nonlinear compensation. Since VIM is effective to suppress inherent characteristics of mechanical components such as frictions, an adaptive control for linear systems can produce an effect. Although the haptic system used in the experiment has two degree-of-freedom motions, controllers of the 𝑥- and 𝑦-axes can be designed separately thanks to the mechanical independence; hence, subscripts of 𝑥 and 𝑦 are omitted in later explanation. Variables and parameters of the haptic device model are shown in Figure 5 and Table 1.

Variables/parameters Unit Meanings

𝑥 𝑝 [m] Position of a grip
𝑓 [N] Force from an operator
𝑓 𝑎 [N] Force from an actuator
𝑚 𝑝 [kg] Mass of a stage
𝑑 𝑝 [Ns/m] Viscosity of a stage

*: 𝑥 or 𝑦 for 𝑥 - and 𝑦 -axis.

The block diagram of the virtual internal model control is shown in Figure 6 and the related variables and parameters are summarized in Table 2. Dynamic equations of the stage and the virtual model are expressed as follows:𝑚𝑝̈𝑥𝑝+𝑑𝑝̇𝑥𝑝=𝑓+𝑓𝑎,𝑚(1)𝑟̈𝑥𝑟+𝑑𝑟̇𝑥𝑟=𝑓.(2) Defining an error as 𝑒𝑟=𝑥𝑟𝑥𝑝, (1) and (2) are transformed into𝑑𝑒𝑑𝑡𝑟̇𝑒𝑟=𝑒0100𝑟̇𝑒𝑟+01𝑢𝑢=𝑑𝑟̇𝑥𝑟+𝑓𝑚𝑟+𝑑𝑝̇𝑥𝑝𝑓𝑓𝑎𝑚𝑝.(3) Minimizing the error defined by (3) makes the stage conform to a response of the virtual model described by (2). To compensate steady-state error, the integral variable 𝑒𝑟 is taken into consideration in the state vector as follows:𝑑𝑑𝑡𝑧=𝐴𝑧+𝐵𝑢,(4) where001𝑒𝐴=010001000,𝐵=,𝑧=𝑟𝑒𝑟̇𝑒𝑟.(5) The control law is calculated using an LQR method with the quadratic criterion:𝐽=0𝑧𝑇𝑄𝑧+𝑢𝑇𝑅𝑢𝑑𝑡,(6) where a positive semi-definite 𝑄𝑅3×3 and a positive definite 𝑅𝑅1×1 are weighting matrices. The input is given as𝑢=𝐹𝑧,𝐹=𝑅1𝐵𝑇𝑃,(7) where 𝑃 is a symmetry positive-definite matrix of a Riccati Algebraic Equation given by𝑃𝐴+𝐴𝑇𝑃𝑃𝐵𝑅1𝐵𝑇𝑃+𝑄=0.(8) Since (7) is expanded as𝐹𝑧=𝑑𝑟̇𝑥𝑟+𝑓𝑚𝑟+𝑑𝑝̇𝑥𝑝𝑓𝑓𝑎𝑚𝑝,(9) a final form of the control law is obtained as follows:𝑓𝑎=𝑚𝑝𝑚𝑟𝑓𝑑𝑟̇𝑥𝑟+𝑑𝑝̇𝑥𝑝𝑓𝑚𝑝𝐹𝑧.(10) For the actual apparatus used in the experiment, parameters of VIM were specified as 𝑚𝑟𝑥=𝑚𝑟𝑦=50[kg] and 𝑑𝑟𝑥=𝑑𝑟𝑦=50[Ns/m] to intentionally obtain a slightly difficult manipulation feeling as a training test with consideration of input range of the actuators. The weighting matrices in (6) for the LQR servo design were decided as 𝑄=diag(150,1.12×107,5100) and 𝑅=1. As a result, the feedback gain matrix was obtained as 𝐹=[10.9,3007.9,100.7].

Variables/parameters Unit Meanings

𝑥 𝑟 [m] Position of a virtual model
𝑒 𝑟 [m] Error ( = 𝑥 𝑟 𝑥 𝑝 )
𝑚 𝑟 [kg] Mass of a virtual model
𝑑 𝑟 [Ns/m] Viscosity of a virtual model

3.2. Online Identification of Human Control Characteristics (Step 2)

Human control characteristics are complex because various kinds of compensators, such as an oculomotor control, a proprioceptive control, and a neuromuscular control, are related to each other [13]. There is, however, a fairly large body of data that can be explained by a linear model plus time delay [15] when an operation condition is limited. One of most famous models supporting such linear model assumption is a crossover model. This model insists that a frequency transfer function of a skilled operator in a man-machine system adapts to make the total system kept unchanged under a variation of the controlled system dynamics. In other words, human changes own control characteristics so that a closed-loop transfer function of whole human-machine system becomes a first-order system at a wide frequency band and the human plays a role of simple linear model. Also in previous study of the present authors, an identification analysis of the skilled operator’s frequency characteristics showed an existence of the cross-over model through a juggling task using a haptic test device [16]. Hence, whole system relating a voluntary motion is simplified into the three components in this study: a linear controller inside a brain, a neuromuscular dynamics, and reaction time delay. After a learning of the machine dynamics is sufficiently finished, the human can be considered as a simple feedback controller which moves the grip to the target position by watching the monitor in case of the PTP task. Finally, a block diagram of a visual voluntary motion control is assumed as a feedback model as shown in Figure 8. In the figure, 𝑟 is a reference position for a pointer, 𝑒 is an error between the target and the present pointer, and 𝑢 is an input computed by a brain controller. Here, the plant block is a virtual 𝑥𝑦-stage of which impedance property is adjusted by the VIM control. The neuromuscular dynamics can often be approximated by a first-order lag [17] and a simplest human controller is a PD controller [9]; hence, the human transfer function, 𝐺(𝑠), is assumed in this study as𝐺𝐾(𝑠)=𝑑𝑠+𝐾𝑝𝑒𝑇𝑠+1𝐿𝑠,(11) where 𝐾𝑝,𝐾𝑑,𝑇, and 𝐿 are a proportional gain of the human brain controller, the differential gain of it, a time constant of the neuromuscular system, and reaction time delay, respectively. As discussed in Section 2.1, compensation of the delay factor is necessary not only for voluntary motion control in a human but also for an adequate human-machine system, and a human has an excellent ability to compensate the delay effect. And an influence of the delay to the control characteristics of the whole human-machine system depends on a response speed of the machine and the task condition. Therefore, in the present study, the response delay of participants was investigated as a preliminary experiment using the VIM control which was designed at Section 3.1. Participant aged 22 years was requested to execute the PTP manipulation hundred times. The time that the pointer begins moving just after the new target circle was displayed on the monitor was counted as the response delay. Figure 7 shows the change of the measured time delay. The dots represent measured values, and the solid lines express an approximated third-order polynomial fitting curve from the measured data. This graph shows no conclusive relationship between time delay and the number of trials, and the value is almost constant at about 0.4 second. Additional nine participants showed same tendency, and significant difference between individuals was not confirmed. For this reason, it was expected that the simple data shift would be sufficient to compensate the delay effect in the identification for the present study. Therefore, the time delay factor described in (11) was omitted for the identification by shifting the measured data for the 0.4 second as a rest time, and the following model was considered for later process:𝐺𝐾(𝑠)=𝑑𝑠+𝐾𝑝.𝑇𝑠+1(12)

Applying a bilinear transformation2𝑠Δ1𝑧11+𝑧1(13) to (12) yields the following discrete impulse transfer function 𝐺[𝑧]:𝐺[𝑧]=𝑏1𝑧1+𝑏0𝑎1𝑧1,𝑎+1(14)1=2𝑇+Δ𝑏2𝑇+Δ,(15)0=2𝐾𝑑+𝐾𝑝Δ𝑏2𝑇+Δ,(16)1𝐾=𝑝Δ2𝐾𝑑,2𝑇+Δ(17) where Δ is a sampling interval. From (15)–(17), following equations are derived:Δ𝑇=21𝑎11+𝑎1,𝐾𝑝=2𝑇+Δ𝑏2Δ0+𝑏1,𝐾𝑑=2𝑇+Δ4𝑏0𝑏1.(18) If 𝑎0,𝑏0, and 𝑏1 are identified from the input/output response data, characteristic parameters of the human controller can be derived using (18). These parameters are used in a design of the next variable dynamics assistive controller.

3.3. Variable Dynamics Assist Control (Step 3)

An assistive control proposed in this paper changes dynamics of the internal model on-line depending on operator’s characteristics. A block diagram of the assistive control is shown in Figure 9(a). In the figure, 𝑟,𝑦,𝑒,𝑣, and 𝑓 are a positional reference, a position of the stage, the error, an output of a brain controller, and a force generated by the hand, respectively. It is assumed that (a) 𝐾𝑝,𝐾𝑑, and 𝑇 are time-slowing changing parameters and that (b) parameters of a virtual machine 𝑚 and ̃𝑏 can be tuned, because the assumption (b) is realized by the VIM control, that is, 𝑚 and ̃𝑏 are adjustable parameters in this scheme. Figure 9(a) expresses a general human-machine system that includes a human controller (𝐾𝑝+𝐾𝑑)/(𝑇𝑠+1)=𝐶 and a plant ̃1/(𝑚𝑠+𝑏)𝑠=𝑃 for virtual haptic interface device. It can be considered conversely that the system consists of a plant 𝐶 changing slowly the parameters (𝐾𝑝,𝐾𝑑, and 𝑇) and the controller 𝑃 having directly variable coefficients (𝑚 and ̃𝑏), as shown in Figure 9(b).

Note that the output 𝑦 of new controller 𝑃 cannot be changed arbitrarily and that only tuning of the controller’s coefficients is possible. Moreover, transformation of the block diagram shown in Figure 9(b) yields a general feedback form, as shown in Figure 9(c). In the following, in order to avoid misunderstanding owing to habits, characters for variables 𝑥 and 𝑢 are used instead of 𝑓 and 𝑒, respectively. Then, the following equations are obtained:𝑥𝐾(𝑠)=𝑝+𝐾𝑑𝑠𝑢1𝑇𝑠+1(𝑠),(19)𝑢(𝑠)=̃𝑏𝑠𝑒𝑚𝑠+(𝑒𝑠),(20)(𝑠)=𝑟𝑟(𝑠)𝑥(𝑠),(21)̃𝑏(𝑠)=𝑚𝑠+𝑠𝑟(𝑠).(22)

The purpose of the PTP task is a tracking such that 𝑦𝑟 in the original block diagram shown in Figure 9(a). This means that 𝑒0 (in Figure 9(b)), that is, 𝑢0 (in Figure 9(c)); then (20) indicates that 𝑒0 as 𝑡. First, choosing a Lyapnov candidate 𝑉 as 𝑉=(1/2)𝑒(𝑡)2, the condition of convergence is investigated. It can be considered that a closed-loop system shown in Figure 9(c) is almost stable under the assumptions of (a) and (b); hence, it is not always necessary that 𝑑𝑉/𝑑𝑡<0 holds for keeping the stability. Second, an update law for 𝑚 and ̃𝑏 is derived using a Lyapunov-like analysis. If a step input is chosen for 𝑟(𝑡) for the PTP motion, the response of 𝑟(𝑡) defined by (22) becomes almost impulse shape. The moment of 𝑡=0 is, however, not important practically because the purpose of the control is an enhancement of the performance of the motion by making the tracking error small which occurs mainly by the positioning near the target position. In addition, the impulsive response converges into zero rapidly, hence, an approximation as 𝑑𝑟(𝑡)/𝑑𝑡0 holds if 𝑡0. Then, the time-derivative of 𝑉 can be approximated and can be transformed as follows:𝑑𝑑𝑡𝑉(𝑡)=𝑒𝑑(𝑡)𝑒𝑑𝑡(𝑡)𝑒𝑑(𝑡)𝑑𝑡𝑥(𝑡)(𝑡>0)=𝑒𝑑(𝑡)𝑑𝑡1𝐾𝑝+𝐾𝑑𝑠1𝑇𝑠+1̃𝑏𝑠𝑒𝑚𝑠+(𝑠)=𝑒(𝑡)1𝐾𝑝+𝐾𝑑𝑠1𝑇𝑠+1̃𝑏𝑚𝑠+𝑒(𝑠)=𝑒(𝑡)1𝐾𝑝𝐾𝑑/𝑇̃1𝑏𝑇𝑚𝑒𝑠+1/𝑇+𝐾(𝑠)𝑝𝐾𝑑̃𝑏/𝑚̃𝑏1𝑚𝑇̃𝑒𝑠+𝑏/𝑚(𝑠)=𝑒𝐾(𝑡)𝑝𝐾𝑑/𝑇̃1𝑏𝑇𝑚𝜙𝑡,𝑇+𝐾𝑝𝐾𝑑̃𝑏/𝑚̃𝑏̃𝑏𝑚𝑇𝜙𝑡,,𝑚(23) where the function 𝜙(𝑡,𝛼) is defined as𝜙(𝑡,𝛼)=𝑡0𝑒𝛼(𝑡𝜏)𝑒(𝜏)𝑑𝜏.(24) Since it is necessary for each term in (23) to be negative in order to satisfy 𝑑𝑉/𝑑𝑡<0 as long as possible, the following conditions are considered:𝐾𝑝𝐾𝑑/𝑇̃𝑏𝑇𝑚>(<)0if𝑒1(𝑡)𝜙𝑡,𝑇𝐾>(<)0,(25)𝑝𝐾𝑑̃𝑏/𝑚̃𝑏𝑚𝑇>(<)0if𝑒̃𝑏(𝑡)𝜙𝑡,𝑚>(<)0.(26)

Conversely, if parameters do not fulfill the previous inequality conditions, variable parameters 𝑚 and ̃𝑏 are tuned so as the unsatisfied condition will be recovered. Now, the following intermediate variables are introduced:𝛿1=𝜂1̃,𝛿𝑏𝑇𝑚2=𝜂2̃𝑏,𝜂𝑚𝑇1𝐾=sgn𝑝𝐾𝑑𝑇𝑒sgn1(𝑡)𝜙𝑡,𝑇,𝜂2𝐾=sgn𝑝𝐾𝑑̃𝑏𝑒𝑚sgñ𝑏(𝑡)𝜙𝑡,.𝑚(27) By checking signs of a numerator and a denominator of (25) and signs of 𝑚 and ̃𝑏, the following update law can be considered:̃𝑏[]̃𝑏[𝑡]𝑡+Δ+𝑘1𝜎𝛿1𝜂1[][𝑡]|𝑒|,𝑚𝑡+Δ𝑚𝑘2𝜎𝛿1𝜂1|𝑒|,(28) where 𝑘1 and 𝑘2 are positive constant parameters, Δ is a control interval, brackets in previous equations mean a discrete-time point, and a function 𝜎 is defined as||𝛿||𝜎(𝛿)=0,𝛿>0,,𝛿<0.(29) The other update law is derived from (26) in same manner as follows:̃𝑏[]̃𝑏[𝑡]𝑡+Δ𝑘3𝜎𝛿2𝜂2[][𝑡]|𝑒|,𝑚𝑡+Δ𝑚+𝑘4𝜎𝛿2𝜂2|𝑒|,(30) where 𝑘3 and 𝑘4 are positive constants. Equations (28)–(30) are summarized into the following parameter update law.̃𝑏𝑚[𝑡+Δ]=̃𝑏𝑚[𝑡]+𝑘1𝑘3𝑘2𝑘4𝜎𝛿1𝜂1𝜎𝛿2𝜂2.(31) On the implementation, these parameters are updated under the following practical limit to avoid an input saturation of actual actuators: 𝑏<̃𝑏<𝑏,𝑚<𝑚<𝑚,(32) where 𝑏,𝑏,𝑚, and 𝑚 are constant. Here parameters 𝑘𝑖 are chosen as they satisfy 𝑘1𝑘4𝑘2𝑘30. Integral computation described in (24) is executed by using the following alternative online recursive computation:𝜙[]𝑡,𝛼=𝑒𝛼Δ𝜙[]𝑡Δ,𝛼+𝑒[𝑡]Δ.(33) Since (22) cannot be computed directly, an approximation as ̃̃(𝑚𝑠+𝑏)𝑠(𝑚𝑠+𝑏)𝑠/(0.01𝑠+1)2 is used, and the response is computed by the Eular integration with the state-space model which is derived with a controllable canonical form. 𝐾𝑝,𝐾𝑑, and 𝑇 are identified on every PTP motion and are updated according to an appropriateness of the identification result.

4. Experimental Result and Analysis

4.1. Online Identification of Human Controller

For a design of the VIM control of the 𝑥𝑦-stage, the initial parameters were chosen as 𝑚[0]=50[kg], and ̃𝑏[0]=50[Ns/m]. Input information for the identification was chosen as an error between the current position and the target one. Measurement value of force filtered through a 36 Hz LPF was used as output information for the identification. The time-delay effect was compensated by sifting the measured input signal at every PTP motion. Measured data was decimated by a factor of 10 for an identification; in short, an identification sampling time is 20ms to suppress oscillation in the identified parameters. One result of the identification is shown in Figure 10. The solid curve is a simulated step response that was computed using an identified human controller model and virtual dynamics model of the stage. Those identified parameters were 𝐾𝑝𝐾=779.0,𝑑=288.0, and 𝑇=0.18, and the time delay was treated as 𝐿=0.406 in the simulation. Since the response of the identified model resembles to the actual response, it can confirmed that the identification process was reasonable.

4.2. Verification of Assistive Effect

Since a key point of the proposed assistive method is to increase performance of the operator’s manipulation by adjusting the machine dynamics, we investigated whether the performance of an operator who was used to the PTP operation without the adaptive control could be increased with the proposed adaptive control. From this aim, before the presented assistive control was applied to a participant, sufficient training was given to become a skilled operator using the haptic device tuned with fixed parameters which were same initial values on the assistive control. As a result of this preliminary training, it was confirmed that the performance of the participant became good and did not indicate no further improvement by checking the settling time on the PTP operation.

For the assistive control, parameters of the update law in (31) were chosen as 𝑘1=𝑘4=1×104, and 𝑘2=𝑘3=2×104. Since even the expert showed perturbation in the performance at the beginning of several trials, the normal VIM control was executed from the first trial and the adaptive impedance control was activated after 50 seconds (about 10 trials). Figure 11 shows an evolution of the settling time from when the new target was displayed on the monitor till when the pointer was reached into the target circle. When three seconds passed after the pointer was kept staying inside the target circle on the monitor, it was judged that the pointer had been moved to the target by the operation. Values of the 𝑦-axis in the figure show the settling time that does not include three seconds. The solid line shows the result of the adaptive impedance control, and the dotted line shows the other result obtained by nonadaptive impedance control before the participant did not yet try the adaptive impedance control. Each line shows the evolution of trend computed by the moving average computation against five PTP tasks. While the nonadaptive impedance control case shows roughly steady state of 1.8 seconds after 50 trials, the other adaptive impedance control case shows a decrease to about 1.1 seconds. In short, speed of the PTP motion was improved by the adaptive impedance control.

Figure 12 shows an evolution of the accumulation errors |𝑒(𝑡)|𝑑𝑡 till each settling time at each trial. Similarly, the solid and dotted lines show the results of the adaptive impedance control and nonadaptive control, respectively. It can be confirmed that the accumulation error was also decreasing in case of the adaptive control. Both Figures 12 and 11 demonstrated effects of the presented method.

Figure 13 shows change of tuned parameters 𝑚 and ̃𝑏. At the beginning of trials, these values were constant because the assistive control was activated after 50 seconds. After 150 seconds, 𝑚 was saturated at the lower limit that was specified for the safety. The reason of this nonconvergence is that the update law (31) cannot guarantee to stop the update of the parameters since the law was designed so as to make the tracking error be zero. This practical issue can be avoided by introducing a dead-zone against small error against the update law.

Finally, transitions of the identified parameters of the operator’s control model, 𝑇,𝐾𝑝, and 𝐾𝑑, are shown in Figures 14, 15, and 16, respectively. Identified parameter varies during till 30 trials. Transitions of their moving averages are comparatively flat at period of 30–80 trials except rapid change due to large outlier in 53rd trial. Although it is difficult to find tendency of change of the identified characteristics, it was confirmed that their moving averages of parameters 𝑇,𝐾𝑝, and 𝐾𝑑 are almost constant after 60 trials of when the tracking error keeps small in Figures 11 and 12. Their constant values do not differ much from their initial values which are ones before the activation of the adaptive control law. In short, it can be considered that total performance was increased by changing the machine side mainly without imposition of large change in human side. This supposition is not authentic since it is not demonstrated by statistical analysis with sufficient number of participants. These are future work.

Results of the experiment, however, showed that the proposed assistive control approach works well, and it can be said that the total performance of whole human-machine system can be enhanced by changing the dynamics of the machine itself.

5. Conclusion

For a force-feedback haptic interface system, an adaptive impedance assistive control to enhance the manipulation performance was proposed. The strategy consists of an identification of an operator’s control characteristics and an adaptive online tuning of the dynamic property of the machine. The tuning is executed by changing impedance parameters of a virtual internal model for the machine. The adaptive law of the tuning was derived by utilizing a Lyapunov stability concept. Using a haptic 2-DOF interface device, it was demonstrated that proposed adaptive impedance assistive control worked effectively. The tuning law, however, cannot guarantee a convergence to a steady state without reaching to the safety limit yet since the presented algorithm was designed by only focusing on an enhancement of the performance of human manipulation without consideration of the machine limit such as an input saturation and frequency bandwidth. This practical issue will be resolved by introducing a tradeoff computation between a performance and requirement of the machine side such as an energy consumption. In the present study, however, a basic strategy for the design of human assistive system could be shown; hence, we would like to treat such practical issues in future.


This work is supported by the Grant-in-Aid for 21st Century (Center of Excellence) COE Program in Ministry of Education, Culture, Sports, Science and Technology. Preparation of the experimental system and the experiment were supported by Keiichi Kurihara and other participants who embraced the authors requests kindly. The authors are grateful to all of them for supporting this work.


  1. T. A. McMahon, “Mechanics of locomotion,” International Journal of Robotics Research, vol. 3, no. 2, pp. 4–28, 1984. View at: Google Scholar
  2. A. B. Joaquin and H. Herr, “Adaptive control of a variable-impedance ankle-foot orthosis to assist drop-foot gait,” IEEE Transactions on Neural Systems and Rehabilitation Engineering, vol. 12, no. 1, pp. 24–31, 2004. View at: Publisher Site | Google Scholar
  3. Y. Yamada, H. Konosu, T. Morizono, and Y. Umetani, “Proposal of skill-assist: a system of assisting human workers by reflecting their skills in positioning tasks,” in Proceedings of the IEEE International Conference on Systems, Man, and Cybernetics, vol. 4, pp. 11–16, October 1999. View at: Google Scholar
  4. H. Konosu, I. Araki, and Y. Yamada, “Practical Development of Skill- Assist,” Journal of The Robotics Society of Japan, vol. 22, no. 4, pp. 508–514, 2004 (Japanese). View at: Google Scholar
  5. M. Uemura, K. Kanaoka, and S. Kawamura, “Power assist system for sinusoidal motion by passive element and impedance control,” in Proceedings of the IEEE International Conference on Robotics and Automation (ICRA '06), pp. 3935–3940, Orlando, Fla, USA, May 2006. View at: Publisher Site | Google Scholar
  6. V. Duchaine and C. M. Gosselin, “General model of human-robot cooperation using a novel velocity based variable impedance control,” in Proceedings of the 2nd Joint EuroHaptics Conference and Symposium on Haptic Interfaces for Virtual Environment and Teleoperator Systemsv (WHC '07), pp. 445–451, March 2007. View at: Publisher Site | Google Scholar
  7. S. Suzuki, K. Kurihara, K. Furuta, F. Harashima, and Y. Pan, “Variable dynamic assist control on haptic system for human adaptive mechatronics,” in Proceedings of the 44th IEEE Conference on Decision and Control, and the European Control Conference (CDC-ECC '05), pp. 4596–4601, Seville, Spain, December 2005. View at: Publisher Site | Google Scholar
  8. A. Tustin, “The Nature of the Operator’s Response in Manual Control and its Implications for Controller Design,” Journal of the Institution of Electrical Engineers, vol. 94, no. 2A, pp. 190–202, 1947. View at: Google Scholar
  9. J. R. Ragazzini, “Engineering aspects of the human being as a servo-mechanism,” in Proceedings of the Meeting of the American Psychological Association, 1948. View at: Google Scholar
  10. S. Baron, D. L. Kleinman, and W. H. Levison, “An optimal control model of human response part II: prediction of human performance in a complex task,” Automatica, vol. 6, no. 3, pp. 371–383, 1970. View at: Google Scholar
  11. D. M. Wolpert and M. Kawato, “Multiple paired forward and inverse models for motor control,” Neural Networks, vol. 11, no. 7-8, pp. 1317–1329, 1998. View at: Publisher Site | Google Scholar
  12. M. Kawato, “Internal models for motor control and trajectory planning,” Motor Systems, vol. 9, no. 6, pp. 718–727, 1999. View at: Publisher Site | Google Scholar
  13. R. C. Miall, D. J. Weir, D. M. Wolpert, and J. F. Stein, “Is the cerebellum a smith predictor?” Journal of Motor Behavior, vol. 25, no. 3, pp. 203–216, 1993. View at: Google Scholar
  14. K. Kosuge, K. Furuta, and T. Yokoyama, “Virtual model following control of robot arms,” IEEE Robotics and Automation, pp. 1549–1554, 1987. View at: Google Scholar
  15. D. L. Kleinman, S. Baron, and W. H. Levison, “An optimal control model of human response part I: theory and validation,” Automatica, vol. 6, no. 3, pp. 357–369, 1970. View at: Google Scholar
  16. K. Furuta, Y. Kado, S. Shiratori, and S. Suzuki, “Assisting control for pendulum-like juggling in human adaptive mechatronics,” Journal of Systems and Control Engineering IMechE, vol. 225, no. 6, pp. 709–720, 2011. View at: Google Scholar
  17. A. Phatak, H. Weinert, I. Segall, and C. N. Day, “Identification of a modified optimal control model for the human operator,” Automatica, vol. 12, no. 1, pp. 31–41, 1976. View at: Google Scholar

Copyright © 2012 Satoshi Suzuki and Katsuhisa Furuta. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

Article of the Year Award: Outstanding research contributions of 2020, as selected by our Chief Editors. Read the winning articles.