Description

A passivity criterion for N-port multilateral haptic systems

Categories

Published

All materials on our website are shared by users. If you have any questions about copyright issues, please report us to resolve them. We are always happy to assist you.

Share

Transcript

Victor Mendez and Mahdi Tavakoli
Abstract
1
.
This paper presents a criterion for passivity of
n
-port networks, which can be used to model multilateral systems involving haptic information sharing between a number of users. Such systems have recently found interesting applications in both cooperative haptic teleoperation and haptics-assisted training. The criterion presented in Theorem 1, which is necessary and sufficient for passivity of the
n
-port network, imposes 2
n
conditions on the immitance parameters of the network and on the residues of the immitance parameters at their imaginary-axis poles. It is shown that when
n
= 2, the proposed conditions reduce to the well-known Raisbeck’s passivity criterion for two-port networks. Another special case for which the proposed criterion has been simplified corresponds to three-port networks. Finally, the passivity of a dual-user haptic system for control of a single teleoperated robot is investigated.
I.
I
NTRODUCTION
Two-port networks are overwhelmingly the method of choice for modeling a bilateral teleoperation system, which consists of a slave robot and a master user interface. The human operator controls the slave and is provided with haptic feedback concerning slave/environment contact forces through the master. Fig.1 shows the equivalent circuit representation of a teleoperation system. Usually, only the linear dynamics of the master and slave are considered:
ssesmmhm
x M f f x M f f
=−=+
, (1) where the hand/master interaction is denoted by
f
h
and the slave/environment interaction is denoted by
f
e
.
M
m
, M
s
, x
m
, x
s
, f
m
,
and
f
s
are the master and slave inertias, positions, and control signals, respectively. Impedances
Z
h
and
Z
e
denote the dynamics characteristics of the human operator’s hand and the remote environment, respectively.
F
h
*
and
F
e
*
are the operator’s and environment’s exogenous input forces, which are independent of the teleoperation system behavior [1].
Fig. 1: Two-port model of a bilateral teleoperation system
Fig. 2 shows a general 2-port network. Depending on which combination of these four quantities (
I
1
, I
2
, V
1
, V
2
) are chosen as independent and dependent variable pairs, six different ways for modeling the 2-port network exist. For instance, using the impedance parameters, the 2-port network can be modeled as
1
This research was supported by the Natural Sciences and Engineering Research Council (NSERC) of Canada under grants RGPIN-372042 and EQPEQ-375712. The authors are with the Department of Electrical and Computer Engineering, University of Alberta, Edmonton, Alberta, T6G 2V4 Canada. (e-mail: vmendez@ualberta.ca, tavakoli@ece.ualberta.ca). Fig. 2: A general 2-port network
=
212221121121
I I z z z zV V
(2) Accordingly, the impedance model of the bilateral teleoperation system in Fig. 1 is given by
−
=
smeh
X X z z z zF F
22211211
(3) The main goals of teleoperation control are transparency and stability. Transparency is the ability of a teleoperation system to present the undistorted dynamics of the remote environment to the human operator [2], and requires the master and the slave positions and interactions to match regardless of the operator and environment dynamics
smeh
x x f f
==
, (4) Taking precedence to transparency is closed-loop stability, which is crucial for safe teleoperation. For the analysis of closed-loop stability of a teleoperation system, according to Fig. 1, the knowledge of the human operator and the environment dynamics are needed in addition to the teleoperation system
immitance
parameters
(z, y, h,
or
g
). In practice, however, the model for the human operator and environment are usually unknown, uncertain, and/or time-varying. In fact, the dynamic parameters of the human operator change in response to the specific requirements of the task at hand [3], [4], and the dynamic parameters of an environment may also change over time. This makes it impossible to use conventional techniques to investigate the closed-loop stability of a teleoperation system. However, assuming that
Z
h
(
s
) and
Z
e
(
s
) in Fig. 1 are passive, we might be able to draw stability conditions that are independent of the human operator and the environment (see Section II). On the other hand, some tasks can be performed more effectively using two hands rather than one or through collaboration rather than individual operation. Also, by using multiple master interfaces each with a corresponding slave robot, multilateral tele-cooperation systems enable haptic information sharing and collaboration in performing a task in a remote environment by multiple users [5], [6], [7], [8], [9]. The key difference between a multilateral tele-cooperation system and a bilateral teleoperation (i.e., single-master/single-slave) system is that the former cannot be modeled as a 2-port network. Therefore, conventional theories for absolute stability or passivity analysis of 2-port networks will not be adequate for multilateral haptic systems. Therefore, there is a need for tools that can analyze the passivity of multi-user haptic systems modeled as
n
-port networks. The only tool available to date for checking the passivity of an
n
-port is
A Passivity Criterion for N-Port Multilateral Haptic Systems
49th IEEE Conference on Decision and ControlDecember 15-17, 2010Hilton Atlanta Hotel, Atlanta, GA, USA
978-1-4244-7744-9/10/$26.00 ©2010 IEEE274
based on the singular value of the scattering matrix of the network [10]:
))()((
ω ω σ
jS jS S
T
=
∞
(5) The above condition is difficult to verify in the general case and without knowledge of the model parameter values for the robots and the controllers (making it not suitable for control synthesis). The criterion we propose in Section III is necessary and sufficient for passivity of the
n
-port network and is easy to check as it imposes 2
n
conditions directly on the immitance parameters of the network and on the residues of the immitance parameters at their imaginary-axis poles.
II.
A
BSOLUTE
S
TABILITY AND
P
ASSIVITY FOR
2-
PORT
N
ETWORKS
Two well-known methods have been developed to investigate the stability of a 2-port network connected to unknown terminations. These methods are known as Llewellyn’s absolute stability criterion and Raisbeck’s passivity criterion. Both criteria work under the assumption that the operator and the environment are passive. By definition, a 2-port network is absolutely stable if it remains stable under all possible uncoupled passive terminations. Also by definition, a 2-port network is passive if the total energy delivered to the network at its input and output ports is non-negative [11].
A
.
Llewellyn’s absolute stability criterion:
If
p
represents any of the four immitance parameters (z
, y, h, g
) of a 2-port network, the criterion establishes that the network is absolutely stable if and only if [11]
•
p
11
and
p
22
have no poles in the right-half plane (RHP)
•
Any poles of
p
11
and
p
22
on the imaginary axis are simple with real and positive residues
•
For all real values of frequencies
,
we have
0)Re(ReRe2
0Re0Re
211221122211
2211
≥−−≥≥
p p p p p p
p p
(6) where
Re
denotes the real part.
B
.
Raisbeck’s passivity criterion
: The necessary and sufficient conditions for passivity of a 2-port network with the immitance parameter
p
are [11]
•
The
p
-parameters have no RHP poles.
•
Any poles of the
p
-parameters on the imaginary axis are simple, and the residues of the
p
-parameters at these poles satisfy the following conditions:
o
If
ij
k
denotes the residue of
ij
p
and
∗
ij
k
is the complex conjugate of
ji
k
, then
∗
=≥−≥≥
122121122211
2211
,000
k k with k k k k
k k
(7)
•
The real and imaginary part of the
p
-parameters satisfy the following conditions for all real frequencies
( ) ( )
0ImImReReReRe4
0Re0Re
22112221122211
2211
≥−−+−
≥≥
p p p p p p
p p
(8) where
Im
denotes the imaginary part. So far, these criteria have not been extended to
n
-port networks where
>2
. Our objective is to give necessary and sufficient conditions for passivity of a general
n
-port network, which can represent a system where multiple haptic teleoperators are used to perform a physical or virtual task collaboratively. Nonetheless, the proposed criterion can be applied for passivity analysis of any
n
-port network including microwave circuits [12], [13].
Fig. 3: A general
n
-port network
III.
M
AIN
R
ESULT
:
P
ASSIVITY
C
ONDITIONS FOR N
-P
ORT
N
ETWORKS
An
n
-port network can be defined as a network containing
n
pairs of terminals for external connections. Each pair of terminals represents a port to which an external network can be connected (Fig. 3). The external behavior of the
n
-port network can be determined if all the
I
i
currents and
V
i
voltages are known. In this section, we present the necessary and sufficient conditions for passivity of an
n
-port network. By analogy with 2-port networks, an
n
-port network is passive if, for all excitations, the total energy delivered to the network at its input and output ports is non-negative. Mathematically, this passivity definition is expressed as
( )
0)()()()()()()(
2211
≥+++=
∞−
τ τ τ τ τ τ τ
d ivivivt E
t nn
(9) where )(
t E
is the total energy delivered to the
n
-port network. We know that the passivity of a 1-port network is a necessary and sufficient condition for positive-realness of its impedance, which is equivalent to
( )
0)()(Re
≥
∗
s I sV
for
0Re
≥
s
(10) By analogy, (9) is equivalent to the following condition
( )
0)()()()(Re
11
≥++
∗∗
s I sV s I sV
nn
for
0Re
≥
s
(11) where
)(
s I
i
∗
is the complex conjugate of
)(
s I
i
.
Using impedance parameters of the
n
-port network, the relations between voltages and currents are
=
)()()()()()(
)()()(
)()()(
)()()(
2121222211121121
s I s I s I s zs zs z
s zs zs z
s zs zs z
sV sV sV
nnnnn
nnn
(12) which can be compactly described as
ZIV
=
. The following Theorem holds for any of the four immittance parameters, yet for brevity we write it only in terms of impedance parameters.
Theorem 1:
The necessary and sufficient conditions for passivity (defined by (11)) of an
n
-port network are:
•
The
p
-parameters have no RHP poles.
275
•
Any poles of the
p
-parameters on the imaginary axis are simple, and the residues of the
p
-parameters at these poles satisfy the following conditions:
0.0)()()(.30.2,,1 ,0.1
'2111112313211
21122211
13212311
1131133311
1121122211
≥−≥−−−−−≥−=≥
−=
iini innnii
k uk n
k k k k k
k k k k
k k k k
k k k k k
k k k k k
nik
(13) where
ij
k
denotes the residue of
ij
p
and
∗
ij
k
is the complex conjugate of
ji
k
. The terms
ij
u
are the elements of an upper triangular matrix
U
used to diagonalize the matrix
K
according to
KUKU
*
=′
. The coefficients
'
ii
k
are the elements of the diagonal matrix
K
′
.
•
The real and imaginary parts of the
z
-parameters satisfy the following conditions for all real frequencies
0.1
'11
≥
z
0.2
'11'21'12'22'11
≥−
z z z z z
(14) 0)()()(.3
'11'12'31'32'11'21'12'22'11'13'21'23'11'11'31'13'33'11
≥−−−−−
z z z z z
z z z z
z z z z
z z z z z
0.
'211'
≥−
−=
iini innn
zw zn
where
)(
*'
jiijij
z z z
+=
21. The terms
ij
w
are the elements of an upper triangular matrix
W
used to diagonalize the matrix
Z
′
according to
WZWZ
*
′′=′
. The elements
''
ii
z
are the entries of the diagonal matrix
Z
′′
.
Proof:
Eliminating the voltages in (11) by using (12), we find that the
n
-port network passivity is equivalent to
[ ]
0)(Re
≥
sF
for
0Re
≥
s
(15) where
)]()()()()()(
)()()()()()(
)()()()()()(Re[)(Re
11222121
111111
s I s I s zs I s I s z
s I s I s zs I s I s z
s I s I s zs I s I s zsF
nnnnnn
nnnn
∗∗∗∗∗∗
++++++++++=
(16) Also, the rational function )(
sF
is positive real (i.e., (15) holds) if and only if, in addition to being real for real
s
, )(
sF
meets the following conditions:
A.
)(
sF
has no RHP poles
B.
Any poles of )(
sF
on the imaginary axis are simple with real and non-negative residues
C.
[ ]
ω ω
∀≥
,0)(Re
jF
For condition A, we require that none of the
z
-parameters of the
n
-port network have any poles in the RHP. For condition B, assume that
F(s)
has a simple pole at
0
ω
js
=
with a residue
0
k
. Let
nn
k k k k
,,
211211
denote the residues of
nn
z z z z
211211
,, respectively, at this pole. Expanding
F(s)
in a Laurent series about
0
ω
js
=
and keeping only the dominant terms in the neighborhood of the pole, we get
0000000101
001001
00101011
00
)()()()()()(
)()()()()()(
ω ω ω ω ω ω ω ω ω ω ω ω ω ω ω ω ω
js j I j I jk
js j I j I jk
js j I j I jk
js j I j I jk
jsk
nnnnnn
nn
−++−++−++−=−
∗∗∗∗
(17) which is equivalent to
)()()()()()(
)()()()()()(
00000101
010010101011
0
ω ω ω ω ω ω
ω ω ω ω ω ω
j I j I jk j I j I jk
j I j I jk j I j I jk
k
nnnnnn
nn
∗∗∗∗
++++++=
(18) In (18),
0
k
must be a real and non-negative number to satisfy condition B. All the
ii
k
for are real and positive since the impedances )(
s z
ii
are positive real functions. Also,
)()(
00
ω ω
j I j I
ii
∗
is real and positive. Note that in the pairs
)()()(
000
ω ω ω
j I j I jk
i jij
∗
)()()(
000
ω ω ω
j I j I jk
ji ji
∗
+
, since
)()(
00
ω ω
j I j I
i j
∗
and
)()(
00
ω ω
j I j I
ji
∗
are complex conjugates,
ij
k
and
ji
k
are also complex conjugates. Since the right side of (18) is a Hermitian form, it can be diagonalized. To do so, (18) can be written in matrix form as
[ ]
=
nnnnn
nnn
I I I k k k k k k k k k I I I k
21212222111211**2*10
=
IKI
*
(19) The
-matrix is diagonalizable and we want to find a linear transformation such as
KUKU
*
=′
where
K
′
is a diagonal matrix,
U
is an upper triangular matrix, and
*
U
(the transpose complex conjugate of
U
) is a lower triangular matrix. By representing the
U
matrix in the reduced row-echelon form, we arrive to the following system
=
⋅
⋅
∗∗∗∗∗∗
nnnnn
nnnnnnnnnnn
k k k k
k k k k
k k k k
k k k k
uuuuuuk k k k uuuuuu
3213333231
2232221
1131211
322311312''33'22'11321231312
1000
1001010000000000001010010001
(20) Solving for
K
′
and
U
will lead us to expressions for each
'
ii
k
as a function of
ij
k
elements. The left side of (20) is
ni
,,1
=
276
+++
+++++++
=∗
ni iiinnnnnnnnnnnnn
k uuk uk uuk u
uk uk u
uk k uuk uuk uk uuk uk u
k uk uuk uk uuk uk u
uk uk uk uk u
uk uk k
1'23'332'22*231'11*132'221'11*121'11'33'323'22*2'13'11*1'22*212'11*1'11*1'3323'22*2313'11*13*222312'11*13'11*1313'11*12'2212'11*12'11*1213'1112'11'11
(21) The solution to (20) is
'11
k
=
11
k
1121122211
'22
k k k k k
k
−=
1112313211
21122211
13212311
1131133311
'33
)()()(
k k k k k
k k k k
k k k k
k k k k k
k
−−−−−=
'211'
iini innnnn
k uk k
−=
−=
(22) Now, (19) can be rewritten as )()(
*0
UIKUIUIKUIIKI
***
′=′==
k
(23) implying that
0
k
will be non-negative and equivalently condition B holds iff
'
ii
k
in (22) are all non-negative (this also implies 0000
443322
≥≥≥≥
nn
k k k k
,,,
). Therefore, it is established that condition B holds iff (13) holds. Regarding condition C, the real part of )(
ω
jF
can be obtained from
[ ]
[ ]
)()(
21)(Re
*
ω ω ω
jF jF jF
+=
. By using)(21
*'
jiijij
z z z
+=
, the real part of )(
ω
jF
can be written as
[ ]
)]()()()()()(
)()()()()()(
)(Re
'1'11'111'11
ω ω ω ω ω ω
ω ω ω ω ω ω
ω
j I j I j z j I j I j z
j I j I j z j I j I j z
jF
nnnnnn
nn
∗∗∗∗
++++++
=
(24) or equivalently as
[ ]
IZI
*
′=
)(Re
ω
jF
(25) where
=
=′
)()()(
)()()(
)()()(
''2'1'2'22'21'1'12'11
ω ω ω ω ω ω ω ω ω
j z j z j z
j z j z j z
j z j z j z
nnnnnn
Z
++++++
)()()()(
)()()()(
)()()()(
21
**11*22*1221*11*1111
ω ω ω ω
ω ω ω ω
ω ω ω ω
j z j z j z j z
j z j z j z j z
j z j z j z j z
nnnnnn
nnnn
(26) The
z
-parameters have complex values, i.e.,
=
+
where
is the real part and
is the imaginary part of
.It is easy to see that (24) is a Hermitian form. Using a procedure similar to (19)-(23), which was for the residue matrix,
Z
′
can be expressed as
WZWZ
*
′′=′
where
Z
′′
is a diagonal matrix and
W
is an upper triangular matrix.
=
∗∗∗∗∗∗
''3'2'1'3'33'32'31'2'23'22'21'1'13'12'11322311312''''33''22''11321231312
1000
1001010000000000001010010001
nnnnn
nnnnnnnnnnn
z z z z
z z z z
z z z z
z z z z
wwwwww z z z zwwwwww
(27) The solution to (27) is
''11
z
=
'11
z
'11'21'12'22'11''22
z z z z z
z
−=
(28)
'11'12'31'32'11'21'12'22'11'13'21'23'11'11'31'13'33'11''33
)()()(
z z z z z
z z z z
z z z z
z z z z z
z
−−−−−=
'211'''
iini innnnn
zw z z
−=
−=
Now, (25) can be rewritten as
[ ]
I)(WZI)(WIWZWIIZI
****
′′=′′=′=
)(Re
ω
jF
(29) Therefore,
[ ]
ω ω
∀≥
,0)(Re
jF
(condition C holds) iff the
z
''
-parameters in (28) are non-negative (this also implies
0,0,0,0
''44'33'22
≥≥≥≥
nn
z z z z
). Therefore, condition C holds iff (14) holds. In summary, conditions A, B and C are necessary and sufficient for (15) or equivalently (11), which defines the
n
-port network passivity. This concludes the proof.
IV.
C
ASE
S
TUDY
:
2-
PORT AND
3-
PORT NETWORKS
In this section, we consider the special cases resulting from substituting
n
= 2 and
n
= 3 in Theorem 1.
A. Two-port networks
It is easy to see that solving (20) for
n
= 2 results in
=
/
,
∗
=
/
, and
=
,
=
−
(!"
It is straightforward that condition (13) in Theorem 1 is same as condition (7) in the Raisbeck’s criterion. Writing
as
+
where
is the real part and
is the imaginary part of
,
we have
#
("
("
("
("$=%
(
+
"+
(
−
"
(
+
"−
(
−
"
&
which can be diagonalized as
277
#!
∗
$#
!!
$)
!*
where
''11
z
=
'11
z
'11'21'12'22'11''22
z z z z z
z
−=
(31) with
=
and
∗
=
. Using
)(21
*'
jiijij
z z z
+=
, the second condition in (31) becomes
−(
+
"
−(
−
"
-!
(32)
which is same as condition (8) in the Raisbeck’s criterion. Therefore, Theorem 1 is in agreement with the Raisbeck’s criterion for the case where
n
= 2.
B. Three-port networks
For
n
= 3, (20) is
.!!
∗
!
∗
∗
0%
!!!
!!!
&.
!
!!0=.
0
(33) which can be solved to get
=
(34)
=
−
1112313211
21122211
13212311
1131133311
'33
)()()(
k k k k k
k k k k
k k k k
k k k k k
k
−−−−−=
Also, for
n
= 3, (27) is
.!!
∗
!
∗
∗
0%
!!!
!!!
&.
!
!!0=%
&
(35) which renders the following solution
''11
z
=
'11
z
'11'21'12'22'11''22
z z z z z
z
−=
(36)
'11'12'31'32'11'21'12'22'11'13'21'23'11'11'31'13'33'11''33
)()()(
z z z z z
z z z z
z z z z
z z z z z
z
−−−−−=
As a result, a 3-port network is passive if and only if
•
The
-parameters have no poles in the RHP
•
The following residues conditions must be satisfied by the residues of the
z
-parameters at their imaginary-axis poles
0)()()(03,2,10
1112313211
21122211
13212311
1131133311
1121122211
≥−−−−−≥−=≥
k k k k k
k k k k
k k k k
k k k k k
k k k k k
ik
ii
(37)
•
The real and imaginary parts of the
-parameters satisfy the following inequalities
-!, 1=,2,
−
(
+
"
+(
−
"
3
−
(
+
"
+(
−
"
3
−
(
+
"
+(
−
"
3
+(
+
"(
+
"(
+
"
(38)
+(
+
"(
−
"(
−
"
−(
+
"(
−
"(
−
"
+(
+
"(
−
"(
−
"-!
V.
A
PPLICATION OF THE
P
ROPOSED
C
RITERION TO A
D
UAL
-U
SER
H
APTIC
T
ELEOPERATION
S
YSTEM
The goal is to train the trainee to do a task under haptic guidance from the mentor. In such a shared-control haptic environment, a parameter
4
can be adjusted such that the trainee and the mentor collaborate and each contributes to the position command while receiving some force feedback. This provides “hand-over-hand” training using haptic assistance. Consider the “four-channel multilateral shared control architecture” given in [14] and depicted in Figure 4. In this 3-robot cooperative manipulation system, the desired velocity and force of each robot is a function of the velocities and forces of the other two robots.
The two human operators are in contact with the two master devices and the slave is in contact with an environment. In frequency domain, this robot models are represented as
5
6
7
8
=9
8
+9
:6
5
6
7
8
=9
8
+9
:6
(39)
5
;
7
<
=−9
<
+9
:;
where
5
6
=
6
, 5
6
=
6
?@A 5
;
=
;
are the
mass models of the two masters and the single slave, respectively. Also,
9
8
,
9
8
and
9
<
are the contact forces between each master and its human operator, and between the slave and its environment. The controller outputs in the 4-channel architecture are
9
:6
=−B
6
7
8
−B
C6
7
8D
+B
E6
9
8
−B
6
9
8D
9
:6
=−B
6
7
8
−B
C6
7
8D
+B
E6
9
8
−B
6
9
8D
(40)
9
:;
=−B
;
7
<
+B
7
<D
−B
F
9
<
+B
9
<D
where
B
6
=G
6
+
H
IJ
;
and
B
;
=G
;
+
H
K
;
denote local PD position controllers, and
B
E6
,B
F
are local force feedback terms for the two masters (
1=,2"
and the slave, respectively. Controllers
B
, B
C
are position compensators similar to
B
;
and
B
6
, respectively. Lastly,
B
6
,B
are feedforward force terms for the two masters and the slave. In (40),
7
8D
,7
8D
,7
<D
are the desired positions and
9
8D
,9
8D
,9
<D
are the desired forces for the two masters and the slave, respectively. The desired velocity and force of one robot depends on the actual velocities and forces of the other two robots as the following expressions and Fig. 4 show it
7
8D
=47
<
+(−4" 7
8
7
8D
=(−4"7
<
+4 7
8
7
<D
=47
8
+(−4" 7
8
9
8D
=49
<
+(−4"9
8
(41)
9
8D
=(−4"9
<
+4 9
8
9
<D
=49
8
+(−4" 9
8
where
4 L !,3
is the weight or authority parameter specifying the relative authority that each operator has over the slave. For perfect transparency (assuming no time delay), the choices
B
=5
;
+B
;
,B
6
=+B
E6
,B
=+B
F
,B
C6
=−(5
6
+B
6
"
are normally made in the 4-channel architecture. For simplicity, however, let us choose
B
=
278

We Need Your Support

Thank you for visiting our website and your interest in our free products and services. We are nonprofit website to share and download documents. To the running of this website, we need your help to support us.

Thanks to everyone for your continued support.

No, Thanks

SAVE OUR EARTH

We need your sign to support Project to invent "SMART AND CONTROLLABLE REFLECTIVE BALLOONS" to cover the Sun and Save Our Earth.

More details...Sign Now!

We are very appreciated for your Prompt Action!

x