Jose Luis Huayanay - Robotica

1
G. Legnani
Robotica Industriale
Casa Ed. Ambrosiana - isbn 88-408-1262-8
http://bsing.ing.unibs.it/~glegnani
http://robotics.ing.unibs.it
inviare commenti e suggerimenti a giovanni.legnani@ing.unibs.it
Errata Corrige ed Ampliamenti

2 maggio 2016
pagina riferimento errata corrige

a cavallo delle pagine 19-20 aggiungere:
19-20
... consentono un movimento relativo di traslazione rettilineo mentre i secondi consentono un movimento...
Ellissi di manipolabilità: i diametri di sono proporzionali all’inverso dei valori singolari σi di J.
I valori singolari di J sono la radice quadrata degli autovalori di JJ T .
46
L’area dell’ellisse è inversamente proporzionale al determinante di J.
Le ellissi dS = kJJ T Fs sono di cedevolezza (inverso di quelle di rigidezza).
52 ultima eq. . . . atan2(x, y) . . . . . . atan2(xp , yp ) . . .
52 matrice J δ ∂
T
57 eq.(2.17) S = [x, y, z] S = [x, y, ψ]T
61 eq.(2.32) β = atan2(v, u) . . . β = atan2(u, v) . . .
= atan2(x cos . . . , z − h) − . . . = atan2(z − h, x cos . . .) − . . .
62 eq.(2.34) (J T M Q̇) ˙
(J T M J)
65 3◦ riga §2.17 parallelo all’asse x parallelo all’asse y
√ √
65 ◦
6 riga dal fondo v2 + v2 u2 + v 2
65 eq. (2.44) (2.45) z = l3 cos(γ) sin(β) z = (l2 + l3 cos(γ)) sin(β)
66 eq. (2.47) de- −2l1 z − (x2 + . . . −2l1 z + (x2 + . . .

nominatore eq.
tc 2=
l’equazione in tb è corretta e serve per determinare β,
mentre (per sin(β) 6= 0) γ ! si determina come
66 eq. (2.47) z l2
γ = ± acos −
l3 sin(β) l3
per sin(β) = 0 il robot si comporta come un robot SCARA in cui il primo link è lungo
l1 + l2 ed il secondo è lungo l3 (vedi §2.9.2)
√ √
66 eq. (2.49) γ = 2 atan tc2 + kπ γ = 2 atan ± tc2
66 eq. (2.50) e poco u2 + y 2 u2 + v 2

dopo (3 volte)
66 eq. (2.51) α = (−xy, . . . α = atan2(−xv + yu, xu + yv) singolare
per −xv + yu = xu + yv = 0
2 G.Legnani, Robotica Industriale Ed. CEI. Correzioni e ampliamenti. 2 maggio 2016
67 eq. (2.53) l2 + . . . 0 = l2 + . . .
l1 + . . . 0 = l1 + . . .
T
71 eq.(2.60) . . . − T Fq . . . . . . T T Fq . . .
infatti se se sul robot agisce Fq sulla trasmissione agisce −Fq . (Gli altri segni sono corretti)
79 sotto matrice R XrT Zr = . . . = 1 XrT Xr = . . . = 1

   
91 nel riquadro = T
R10  = T
R01 
96 Rtras(z . . .) − sin(ϑy ) − sin(ϑz )

...−Z
99 ultima formula u = ...−Zz
x ... ux = ... y
⋄
125 2 formula . . . [W01(1) + . . . . . . [W01(0) + . . .
135 aggiungere la didascalia alla tabella di inizio pagina: Descrizione dei termini della eq.(6.2)
136 M34 elemento 1,3 S4 −S4
aggiungere: −ϕ1 = ϕ2 = −ϕ4 = ϕ5 = 90◦ , ϕ3 = ϕ6 = 0
136 §6.2.2 elenco (vedi figura 6.2.2) (vedi figura 6.5)
   
148 M03 =  0  M03 =  −S2 

1 1
152 5◦ eq. s6 = −c5 (...) + s5 (...) c6 = +c5 (...) − s5 (...)
6◦ eq. c6 = ...Xx s6 = ...Xz
ultima eq. c2 = ...∆x − s3 ... c2 = ...∆x + s3 ...
153 fig.6.14 ϑ2 è indicato in maniera erronea (è def. in fig.6.13 ove tutti gli angoli valgono 0).
155 righe 5 e 7 M45 M34
161 Mbp (2 volte) (1 − c) (1 − c3 )
◦
175 3 eq. W23(0) = L23(0) q̇ W23(0) = −L23(0) q̇
   
175 fine pagina M0a = . . . =   M0a = . . . =  

0 1
176 7◦ riga Ṗ = W03(a) P Ṗ = W03(0) P
177 1◦ eq. . . . + L23(0) q̇3 . . . − L23(0) q̇3
178 prima formula H = Lα̈ + L2 α̇2 H = Lq̈i + L2 q̇i2
   
. . . + l1 β̈ cos α . . . + l1 β̈ sin α
178 ultima formula H03(0) =  0  H03(0) =  0 
0 0
181 eq. (6.13) Ṡ = J Q̈ Ṡ = J Q̇
183 fig. 6.33 α2 −α2 (dal link 1 al 2 ma nel verso opposto)
185 nota piè pagina −l1 s1 − l2 s12 = −a −l1 c1 − l2 c12 = −a
ax = −l1 c1 α̇1 α̇2 − ((l1 c1 + l2 c12 )α̇1 + l2 c12 α̇2 ) α̇3
186
ay = −l1 s1 α̇1 α̇2 − ((l1 s1 + l2 s12 )α̇1 + l2 s12 α̇2 ) α̇3
1
207 eq. (7.8) Epi = traccia(...) Epi = traccia(...)
2
215 riga 9, 11 ±25 ±90 ±25◦ ±90◦
G.Legnani, Robotica Industriale Ed. CEI. Correzioni e ampliamenti. 2 maggio 2016 3

223 dopo eq. (8.19) La cinematica inversa . . . La cinematica diretta . . .
251 punto b) . . . ∆qn ] . . . ∆qn ]T
... ...
255 prima eq.
τn . . . τn−1 . . .
262 4◦ riga dal fondo d = P0 − P1 d = P1 − P0
264 eq. calcolo α, β { . . . = −V1T V1 . . . = −V2T V2 { . . . = V1T V1 . . . = V2T V2
(P2 − P1 ) × (P0 − P1 )
264 fondo pagina z = ... z=
||(P2 − P1 ) × (P0 − P1 )||
270 eq.(12.25) amg = bck (amg − bck) cos(α) = 0
soddisfatta se α = ±π/2 per qualsiasi
valore di a, m, b, c, k
oppure per qualsiasi α se amg = bck
271 figura 9.46 v̇x , v̇y vx , vy
Ω s Ω sTi /kp
282 ultima formula G(s) = ≃ G(s) = ≃
Cr ... Cr ...
( (
ϑ 1 ϑ kp′
285 = F ′ (s) = . . . ≃ . . . = F ′ (s) = . . . ≃ . . .
ϑ0 s ϑ0 s
ω2 ωt ωt ωt
286 ω2 = = ω2 = =
2ξ 2 2ξ 2
′
288 sopra eq.(10.5) R (s) = kp (1 + Td s) R (s) = kp′ (1 + Td s)
′
288 eq.(10.5)(10.6) τω = 1/ωω = 1/ωT è l’inverso della banda passante dell’anello di velocità
288 dopo eq.(10.6) kd = kp Td kd = kp′ Td
288 eq.(10.7) . . . 1 + kp Td ) . . . 1 + kp′ Td )
′
T′ kp Td′ ′
kp
ϑ 1 + s k′d J + ′ s
kp ϑ 1 + sTd′ (1 + Td′ s)
p J
288 eq.10.8 = = = = = =
ϑ0 ... ... ϑ0 ... ...
numeratore
1 1
296 inizio prima eq. (s non è un pedice)
Js Js
296 prima eq. J2 Js (s poi si semplifica con quello
che compare nel termine J1 s)
Ω ... Ω ...
296 fine prima eq. = ··· = = ··· = " #
C0 s2 2ξ2 C0 s2 2ξ2
s(J1 + J) 2 + s+1 s(J1 + J) + s+1
ω2 ω2 ω22 ω2
309 figura 10.55 y z
317 fig.11.4 scambiare i nomi “pezzo 1” e “pezzo 2”
353 terzultima riga . . . del rotore. . . . dello statore.
366 J2 = . . . = 4.4210−4 kgm2 J2 = . . . 4.4910−4 kgm2
383 fig.13.18 scambiare r e k
. . . En−1 . . . . . . En . . .
417 1◦ equazione Q Qn (p = n + 1)
n−1
i=0 i=0
425 eq.(14.15) m2 g2 β̇12 c −m2 g2 β̇12 c

425 dopo eq.(14.15) sin β cos β sin β2 cos β2
448 arresto, terzo ca- • categoria 1 • categoria 2 (non ammessa per arresto
so emergenza)
484 tab. A.5 −λSα Cα β̇ + λSα Cβ γ̇ Cα β̇ + Sα Sβ γ̇
485 tab. A.7 λSα Cβ λSα Sβ
486 par. A.5 φ = . . . angolo di precisione φ = . . . angolo di precessione
487 elem.1,2 matrice . . . + sφ cψ . . . − sφ cψ

488 elem.3,1 matrice cφ sθ cψ + . . . −cφ sθ cψ + . . .
495 ultima riga . . . soluzione unica t = −b/a. . . .soluzione t = −b/a; l’altra è α =
π + 2kπ. Se a2 + b2 = c2 le 2 soluzioni
coincidono.
499 prima formula atan2(x, y) = atan . . . atan2(y, x) = atan . . .
509 prima formula f (q1 ) + f ′ (q1 )(q − q1 ) = 0 f (qi ) + f ′ (qi )(q − qi ) = 0
509 ultima riga f ′ (qi ) = 0 f ′ (qi ) 6= 0
510 fig. D.2 nell’esempio centrale scambiare q1 con q2
G.Legnani, Robotica Industriale Ed. CEI. Correzioni e ampliamenti. 2 maggio 2016 999
Bibliografia
Riguardo al metodo matriciale basato su matrici 4 × 4 ci si può anche riferire a
[1] A homogeneous matrix approach to 3D kinematics and dynamics - I. Theory

Mechanism and Machine Theory, Volume 31, Issue 5, July 1996, Pages 573-587
Giovanni Legnani, Federico Casolo, Paolo Righettini, Bruno Zappa
[2] A homogeneous matrix approach to 3D kinematics and dynamics - II. Applications to chains of
rigid bodies and serial manipulators
Mechanism and Machine Theory, Volume 31, Issue 5, July 1996, Pages 589-605
Giovanni Legnani, Federico Casolo, Paolo Righettini, Bruno Zappa
Appendice H
Ampliamenti e chiarimenti
Questa sezione è dedicata a piccole aggiunte che possono rendere più comprensibile il testo
o completarlo.
I riferimenti a pagine, capitoli, figure, equazioni sono relative al libro.
H.1 Ellissi di Manipolabilità, Isotropia (Pagg.45-48)

Le ellissi di manipolabilità sono degli strumenti sintetici per rappresentare le prestazioni cinetostatiche
dei manipolatori nelle diverse zone dell’area di lavoro. Il loro significato verrà spiegato con l’ausilio di
esempi.
H.1.1 Ellissi di manipolabilità in velocità e in forza

Consideriamo un manipolatore cartesiano xy; ciascuno dei due assi avrà una velocità massima ammissibi-
le ±q̇x,max e ±q̇y,max . Le componenti della velocità del dispositivo d’estremità lungo gli assi x, y avranno
quindi modulo minore o uguale a quella del motore corrispondente.qTuttavia in direzioni diagonali esso
potrà raggiungere velocità più elevate ma mai superiori a Vmax = q̇x,max2 2
+ q̇y,max . Il luogo descritto
dai vettori velocità massima nelle varie direzioni è perciò un rettangolo. Sei i due motori hanno pari
prestazioni il rettangolo diventa un quadrato. Il risultato esposto non cambia all’interno dell’area di
lavoro.
Un discorso analogo può essere fatto qualora interessi studiare le forze massime esercitabili nelle diverse
direzioni.
Figura H.1: Manipolatore cartesiano: velocità

(o forze) massime ammissibili in funzione delle
direzioni.
Se si considera un manipolatore polare la velocità “radiale” massima è limitata da quella del motore li-
neare Vr,max = ρ̇max , mentre quella “tangenziale” è limitata dal prodotto della velocità angolare massima
del motore rotativo per la distanza del dispositivo d’estremità dalla torretta Vt,max = ρωmax . Inoltre in
direzioni diverse le due velocità
p si possono sommare e quindi, per particolari direzioni, la velocità può
raggiungere il valore Vmax = (ρ̇max )2 + (ρωmax )2 . Il luogo dei vettori velocità massime è ancora un
rettangolo con lati paralleli alle direzioni radiale e tangenziale. Dato che il valore della distanza ρ varia
1001
1002 APPENDICE H. AMPLIAMENTI E CHIARIMENTI
Vr , max
Fr , max
wmax
Cmax
Figura H.2: Manipolatore polare: velo- Figura H.3: Forze massime esercitabili da
cità massima ammissibile in funzione delle un manipolatore polare in funzione della
direzione e della configurazione. direzione e della configurazione.
nell’area di lavoro, si deduce che la velocità massima tangenziale aumenta con la distanza e il rettangolo
delle velocità si deforma di conseguenza.
Analogamente per le forze possiamo dire che in direzione radiale il manipolatore potrà esercitare una
forza massima costante limitata dalle prestazioni del motore lineare, mentre in direzione tangenziale il
modulo della forza massima esercitabile sarà inversamente proporzionale alla distanza Ft,max = Cmax /ρ
essendo limitata dalla coppia massima del motore rotativo. Il luogo delle forze massime esercitabili è un
rettangolo in cui il lato “tangenziale” decresce con la distanza.
y Fmax y
y
P P
P
b b2
O2 O2
O1 a O2 x O1 x O1 x
b1
Figura H.4: Composizione del-

le velocità in un manipolatore Figura H.5: Determinazione delle forze massime esercitabili
SCARA. da un manipolatore SCARA.
Nel caso di un manipolatore SCARA la componente di velocità fornita da ciascun motore è proporzionale
alla distanza P − Oi tra il dispositivo di estremità P e l’asse del motore Oi ed è limitata dalla velocità
massima di quest’ultimo. La componente di velocità fornita da ogni motore è ortogonale al vettore P −Oi .
In generale, le due componenti di velocità non sono tra loro ortogonali. Il luogo delle velocità massime è
quindi un parallelogramma i cui lati sono proporzionali alle quantità Vi,max = |P − Oi |ωi,max .
La forza massima esercitabile in una direzione dipende dai bracci bi della forza rispetto agli assi dei
motori oltre che dalla loro coppia massima Ci,max . Nessun motore deve essere sovraccaricato per cui:
Fmax = min(C1,max /b1 , C2,max /b2 ). Ovviamente il braccio dipende sia dalla posizione della pinza, sia
dalla direzione della forza. Il parallelogramma delle forze massime generabili in funzione della posizione si
ottiene dall’intersezione di due “corridoi” di semi-larghezza fi = Ci,max /||P − Oi || paralleli alla direzione
P − Oi .
Nel caso di manipolatori a diversi gradi di libertà, la determinazione delle velocità massime raggiungibili
e delle forze massime generabili può essere molto complesso, ma segue i criteri generali qui descritti.
Nel caso di moto xyz (3D), i parallelogrammi sono sostituiti da parallelepipedi.
H.1. ELLISSI DI MANIPOLABILITÀ, ISOTROPIA (PAGG.45-48) 1003
Le ellissi di manipolabilità in velocità ed in forza danno un metodo standard, applicabile in modo automa-
tico anche a manipolatori complicati con diversi gradi di libertà, per indicare in quali direzioni è possibile
avere velocità elevate o esercitare forze più o meno considerevoli. Le ellissi di manipolabilità sono delle
ellissi che approssimano il luogo dei punti di velocità o forza massime. Il disegno di una famiglia di ellissi
relative a diverse posizioni del manipolatore evidenzia il suo diverso comportamento nell’area di lavoro.
Nel caso di manipolatori piani si considera l’ellisse inscritta nei parallelogrammi descritti e tangenti ai
lati nelle loro mezzerie.
Nel caso di manipolatori che si muovono nelle tre direzioni xyz, le ellissi sono sostituite da ellissoidi.
Queste ellissi (o gli ellissoidi) sono determinabili anche nel caso di manipolatori a più gradi di libertà.
La determinazione delle ellissi può essere effettuata sulla base di procedure standard che fanno uso degli
autovettori ed ai valori singolari della matrice jacobiana o di matrici derivate da essa.
Quando i diversi motori di un manipolatore sono tutti dello stesso tipo (rotativo o lineare) spesso per
semplicità si assume che abbiano tutti le stesse prestazioni massime sia in termini di velocità che di coppia
(o forza) massima. Il calcolo delle ellissi si semplifica dato che spesso è sufficiente calcolarne la forma e
non la dimensione. In questo caso la scelta delle unità di misura non influenza la forma dell’ellisse.
Nel caso in cui il manipolatore abbia attuatori dei due tipi, allora è indispensabile tenere conto delle reali
prestazioni di ciascuno altrimenti si ottengono ellissi la cui eccentricità dipende dalle unità di misura
impiegate e non hanno quindi significato fisico.
Le ellissi di velocità e di forza non pesate sono associate alle seguenti equazioni
X
1= q̇i2 = Q̇T Q̇ ⇒ Ṡ T J −t J −1 Ṡ = 1
X
2
1= fqi = FqT Fq ⇒ FsT JJ t Fs = 1
mentre quelle pesate sono associate alle seguenti equazioni
!2
X q̇i
1= = Q̇T DvT Dv Q̇ ⇒ Ṡ T J −t DvT Dv J −1 Ṡ = 1
q̇i,max
!2
X fqi
1= = FqT DfT Df Fq ⇒ FsT JDfT Df J T Fs = 1
fqi,max
con    
1 1
 q̇1,max 0 ...   fq1,max 0 ... 
   
 1   1 
Dv =  0 ...  Df =  0 ... 
 q̇2,max   fq2,max 
   
.. .. .. .. .. ..
. . . . . .
Attenzione: l’ellisse di forza e quella di velocità sono tra loro or-

togonali solo nel caso di ellissi non pesate oppure nel caso di ellissi
pesate quando le matrici dei pesi sono scelte in modo opportuno
per cui Df = kDv−1 ove k 6= 0 è uno scalare di valore arbitrario.
H.1.2 Ellissi d’inerzia (di massa)

Le ellissi di inerzia (o di massa) danno informazioni sulla capacità di accelerazione nelle varie direzioni
dell’area di lavoro quando è sottoposto a forze esterne applicate sul centro utensile. Queste ellissi sono
particolarmente significative nei casi in cui il manipolatore scambia forza con l’ambiente esterno; ad
esempio mentre movimenta un utensile che esegue lavorazioni meccaniche. In congiunzione con le ellissi
di elasticità (o di rigidezza), le ellissi di massa danno informazioni sulla possibilità di insorgenza di
vibrazioni.
Nel §2.8 viene mostrato come le forze motrici Fq necessarie a produrre un’accelerazione dei giunti Q̈
valgono
Fq ≃ MQ̈
ove è stato usato il simbolo “circa uguale” perché vengono qui trascurati i termini quadratici nelle velocità,
cioè il prodotto J˙Q̇, che è trascurabile a basse velocità. Inoltre si considera il caso con forze esterne nulle.
Consideriamo ora il caso che i motori siano spenti in folle (coppia nulla) e ci chiediamo quali forze
applicate all’estremità dal manipolatore avrebbero lo stesso effetto delle coppie motrici appena calcolate.
Dalla cinetostatica sappiamo trasferire forze ad accelerazioni dai giunti al dispositivo d’estremità con le
y
Ellissi d’inerzia: rappresentazione grafica
S̈
Per visualizzare la forza Fs necessaria a produrre una
m2 α certa accelerazione S̈ si può applicare la costruzione
seguente.
Fs x Si considera inizialmente un’accelerazione unitaria
diretta come S̈ = [a1 a2 ]T e la si scompone nelle
m1 direzioni dei due diametri. La componente parallela
a u1 va moltiplicata per m1 , quella parallela a u2 va
moltiplicata per m2 . Le due componenti vanno quindi
moltiplicate per il modulo dell’accelerazione

f1 m1 cos(α)
Fs = = |S̈| =
f m2 sin(α)
2
m1 0 a1
0 m2 a2
Figura H.6: Costruzione grafica ellisse di
massa (inerzia).
relazioni S̈ ≃ J Q̈ e Fq = J T Fs (il segno meno non compare perchè le forze sono quelle equivalenti e non
quelle equilibranti). Si ottiene quindi
Fs ≃ J −t MJ −1 S̈
Questa relazione permette di calcolare quale forza deve essere applicata al dispositivo d’estremità del
manipolatore per generare un’accelerazione desiderata. Possiamo considerare questa relazione come una
→
−
generalizzazione della legge di Newton
f = m− →a.
Alla matrice Ms = J −t MJ −1 può essere associata una ellisse rappresentante l’inerzia (massa) equiva-
lente del manipolatore in una determinata direzione. Infatti consideriamo un’accelerazione S̈ di ampiezza
unitaria e di direzione generica e calcoliamo la forza necessaria a produrla; se facciamo ruotare in tutte
le direzioni questo vettore, il vettore della forza descriverà un ellisse. Questo può essere dimostrato os-
servando che la matrice Ms è simmetrica definita positiva.
Si osserva che forza ed accelerazione sono tra loro parallele solo se la forza è diretta lungo uno degli assi
dell’ellisse.
Le direzioni degli assi del’ellisse sono definite dagli autovettori di Ms e i raggi (semidiametri) dai
corrispondenti autovalori.
Ellissi di rigidezza ed elasticità
L’ellisse di cedevolezza o quella di rigidezza permettono di indicare lo spostamento della pinza del mani-
polatore dovuto alla applicazione di una forza quando i giunti sono cedevoli.
L’ellisse associata all’equazione dS = kJJ T Fs definita al capitolo 2 è di cedevolezza (inversa di quelle di
1 −1
rigidezza). L’ellisse di rigidezza sarebbe quindi associata a Fs = JJ T dS.
k
Si osserva che forza e spostamento sono tra loro paralleli solo se la forza è diretta lungo uno degli assi
dell’ellisse. Per una visualizzazione grafica di questo fenomeno si consideri per analogia il caso di ellisse
di inerzia (massa).
Consideriamo l’ellisse associata alla matrice JJ T . Il diametro in una direzione rappresenta la cedevolezza
equivalente in quella direzione ed è quindi proporzionale allo spostamento dovuto all’applicazione di una
forza in quella direzione.
Le direzioni degli assi del’ellisse di cedevolezza (elasticità) sono definite dagli autovettori ui di JKq J T (o
kJJ T ) e i raggi (semidiametri) dai corrispondenti autovalori λi . Per la rigidezza occorre fare riferimento
alla matrice inversa i cui autovettori sono i medesimi e i cui autovalori valgono 1/λi .
H.1.3 Forza, Velocità massima in una direzione.

La forza massima esercitabile in una direzione si può facilmente determinare esprimendo la forza applicata
all’end-effector come prodotto del modulo per un versore Fs = f us , l’ellissi di forza è cosı̀ definita come
f uT Auf = 1 A = JDfT Df J T
da cui
1
f=√
T
u Au
H.1. ELLISSI DI MANIPOLABILITÀ, ISOTROPIA (PAGG.45-48) 1005
in maniera analoga è calcolabile la velocità massima.
H.1.4 Nota
Si osservi come nei diversi casi considerati le ellissi siano state definite in maniera differente.
In alcuni casi (velocità, forze) si è fatto riferimento a equazioni del tipo v T AT Av = 1 e le direzioni degli
assi dell’ellisse corrispondono agli autovettori della matrice AT A mentre le lunghezze
p ri dei semidiametri
sono legate all’inverso della radice quadrata dei corrispondenti autovalori ri = 1/λi (inverso dei valori
singolari di A).
In altri casi (massa, forza) la relazione di riferimento è del tipo v = Bu e le direzioni degli assi dell’ellisse
corrispondono agli autovettori della matrice B mentre le lunghezze ri dei semidiametri coincidono con il
valore dei corrispondenti autovalori ri = λi .
H.1.5 Esempi applicativi
y
a& (P - O )
l2 b
P
b& (P - A) Nel caso di motori identici o comunque nel caso non
pesato, il robot SCARA (§2.9, pag. 51) può esse-
A re in configurazione isotropa solo se le componenti
delle velocità generate dai due motori sono di egual
l1 modulo |P − O| = |P − A| e tra loro ortogonali
a x −−−−→ −−−−→
P − O⊥P − A. Ciò implica che devono valere le √due
seguenti condizioni: β = ±135◦(3/4π) e l1 = 2l2 ;
O P’ il valore di α è indifferente. In questo caso le ellis-
si di manipolabilità in forza e velocità diventano dei
cerchi.
Figura H.7: Robot SCARA in configu-

razione isotropa (forza e velocità non
pesate).
Figura H.9: Le ellissi di velocità del robot pola-

re (§2.11, pag.58) si allungano in direzione tan-
Figura H.8: Nel caso di motori identici genziale via via che il centro utensile si allonta-
o comunque nel caso non pesato, il ro- na dall’asse di rotazione. Corrispondentemente
bot cartesiano è isotropo in ogni punto le ellissi di forza (non mostrate in figura) si ri-
dello spazio di lavoro. durranno. L’isotropia si potrà avere solo per una
determinata distanza dall’asse di rotazione.
y Ellissi d’inerzia (massa): robot cartesiano

u2 Le ellissi di massa di un robot cartesiano sono estre-
mamente semplici da calcolare. Con riferimento
all’esempio di figura, come è intuitivo, l’ellisse è pa-
u1 rallela agli assi di scorrimento, la massa equivalente
(autovalori) in direzione x vale m1 + m2 mentre quella
in direzione y vale m2 .
m2
Ellissi di elasticità: robot cartesiano
Le ellissi di elasticità hanno diametri paralleli agli
assi e lunghezza proporzionale alla elasticità dei giunti
x corrispondenti. Nel caso di assi tra loro identici il
manipolatore è isotropo. Ovviamente, come per tutti
m1 i manipolatori, questo vale solo se si considera solo
l’elasticità dei giunti e non quella flessionale degli
elementi che compongono i link.
Figura H.10: Robot cartesiano: ellissi di
massa o di elasticità.
y u1 y
u2 u1
u2
m2
Jg 1 m2
b m1 r
m1 g1
x
O1 a
O2
O1 a x
Figura H.12: Ellisse di massa di un robot polare
Figura H.11: Ellisse di massa di un robot SCARA nel caso semplificato in cui la massa del secondo
nel caso semplificato in cui Jg1 = Jg1 = 0, m = 0, link sia concentrata all’estremo h = 0 (confronta
g1 = l1 , g2 = l2 (confronta §2.9.5 pag.55). §2.11.3 pag.59).
Robot SCARA
Immaginiamo per semplicità che le uniche masse significative possano essere considerate puntiformi e
localizzate al termine dei link: m1 al termine del link 1 e m2 al termine del link2.
Le coordinate necessarie allo studio dinamico sarebbero
   
x l1 cos(α) + l2 cos(α + β)
x l1 cos(α) + l2 cos(α + β)  y   l1 sin(α) + l2 sin(α + β) 
S= = Se =   
 xg1  = 


y l1 sin(α) + l2 sin(α + β) l1 cos(α)
yg1 l1 sin(α)
M = diag(m2 , m2 , m1 , m1 )
In questo caso le direzioni degli assi dell’ellisse (autovettori della matrice di massa) ed i relativi semidia-
metri (autovalori) risultano

cos(α + β) − sin(α + β) m1
u1 = u2 = λ1 = m2 + λ2 = m2
sin(α + β) cos(α + β) sin(β)2
L’asse u1 è parallelo al secondo link. Per spostamenti in direzione u2 conta solo la massa m2 , mentre in
direzione u1 occorre sommare una quota derivante dalla massa m1 . Qualora i baricentri non siano nelle
posizioni indicate o i momenti d’inerzia non siano nulli le direzioni (autovettori) e le masse equivalenti
(autovalori) sarebbero differenti.
H.2. BIELLA EQUIVALENTE 1007
0.6 0.6
0.4 0.4
0.2 0.2
y[m]
y[m]
0 0
−0.2
−0.2
−0.4
−0.4
−0.6 −0.4 −0.2 0 0.2 0.4 0.6 −0.6 −0.4 −0.2 0 0.2 0.4 0.6
x[m] x[m]
(a) β > 0 (b) β < 0
Figura H.13: Le ellissi di inerzia del robot ICOMATIC SCARA03.
Robot polare
Immaginiamo per semplicità che le uniche masse significative possano essere rappresentate da due masse
puntiformi: m1 sul link 1, m2 al termine del link2, m la massa trasportata. Siano inoltre Jg1 , Jg2 , eJ i
momenti d’inerzia baricentrali dei due link e della massa trasportata. Le coordinate necessarie allo studio
dinamico sarebbero
   
x ρ cos(α)
 y   ρ sin(α) 
x ρ cos(α)    
S= = Se =   
 xg1  =  g1 cos(α) 

y ρ sin(α)  yg1   g1 sin(α) 
α α
M = diag(m + m2 , m + m2 , m1 , m1 , Jg1 + Jg2 + J)

In questo caso le direzioni degli assi dell’ellisse (autovettori della matrice di massa) ed i relativi semidia-
metri (autovalori) risultano

cos(α) − sin(α) g12 Jg1 + Jg2 + J
u1 = u2 = λ1 = m+m2 λ2 = m+m2 +m1 +
sin(α) cos(α) ρ2 ρ2
L’asse u1 è parallelo ai link, u2 ortogonale ad essi. Per spostamenti in direzione u2 conta solo la massa
m + m2 , mentre in direzione u1 occorre sommare una quota derivante dalla massa m1 e dai momenti
d’inerzia.
A titolo d’esempio si riportano le ellissi d’inerzia del robot Icamatic SCARA03 calcolate per le due con-
figurazioni (a sinistra q2 > 0, a destra q2 < 0). Esse sono tracciate solo per le configurazioni realmente
raggiungibili tenendo conto dei fine corsa reali. Nel caso specifico di questo particolare modello di mani-
polatore (Icomatic SCARA03), le proporzioni e l’orientamento delle ellissi di rigidezza sono molto simili a
quelle delle ellissi di inerzia; tuttavia questo non avviene in generale per altri modelli o altri manipolatori.
H.2 Biella equivalente

Il metodo per l’analisi dinamica esposto nel capitolo 2 è molto adatto per lo studio di sistemi piani e per
qualche sistema 3D particolarmente semplice. Tuttavia esso non si applica facilmente a sistemi 3D in cui
si possono avere rotazioni attorno ad assi di orientamento variabile nel tempo e che determinano quindi
la nascita di effetti dinamici giroscopi.
Il problema può essere risolto facilmente sostituendo masse distribuite (proprietà dinamiche descritte da
tensore (o pseudo-tensore) d’inerzia) con un gruppo di masse puntiformi che nel loro insieme abbiano la
stessa proprietà di massa del corpo distribuito.
Un caso particolarmente semplice è quello di un corpo (ad esempio un’asta, figura H.14(a)) allungato
che ha un momento d’inerzia longitudinale Jz trascurabile e due momenti d’inerzia trasversali tra loro
mb mb
z b
b
x
g a mg
a
y ma ma
(a) corpo allungato (b) modello 2 masse (c) modello 3 masse
Figura H.14: Approssimazione della dinamica di un corpo allungato con un’asta e 2 o 3 masse puntiformi.
identici e molto più grossi del momento longitudinale
Jx = Jy >> Jz
Consideriamo quindi un’asta di massa m, momento d’inerzia baricentrico trasversale Jg , lunghezza ℓ =

a + b il cui baricentro disti delle quantità a e b dagli estremi (figura H.14(b)). In questo caso si può
immaginare di sostituire la massa distribuita con due masse ma e mb poste alle estremità in modo tale
che si conservi la massa e la posizione del baricentro. Deve perciò essere

(  b

 ma = m Jg′ = a2 ma + b2 mb 6= Jg
m = ma + mb a+b
perciò e
ma a = mb b 
 mb = a m
 J ′z = 0
a+b
Tuttavia cosı̀ facendo (a meno di casi fortunati) non si conserva il valore del momento d’inerzia trasversale.
Infatti, ad esempio, considerando un’asta uniforme il suo momento d’inerzia reale Jg e quello del modello
approssimato Jg′ valgono
Z ℓ/2
a = b = ℓ/2 1 1 m 2
Jg′ = mℓ2 Jg = mℓ2 = z dz
ma = mb = m/2 4 12 −ℓ/2 ℓ
Un modello più preciso può essere sviluppato con tre masse, di cui due poste agli estremi ed una mg
posta nel baricentro (figura H.14(c)); in questo caso anche il valore del momento d’inerzia può essere
mantenuto esatto. Si ottiene

 Jg

 ma =
 
 a(a + b)

 m = ma + mb + mg
 
 Jg
ma a = mb b perciò mg = m −

 
 ab
Jg = m a a 2 + m b b 2 


 Jg

 mb =
b(a + b)
Perciò nel caso di barra omogenea si ottiene
m 2 m
ma = mg = m mb =
6 3 6
H.3 Robot 3R (capitolo 2.17)

L’equazione della cinematica inversa necessaria per calcolare α ammette una singolarità quando i due
argomenti della funzione atan2 sono contemporaneamente nulli (−xv + yu = xu + yv = 0 ⇒ u2 + y 2 = 0)
α = atan2(−xv + yu, xu + yv)
tuttavia questa condizione può capitare solo per cos(β) = ±1 se contemporaneamente cos(γ) = ±1 e solo
per particolari valori delle lunghezze per cui capita l1 ± l2 ± l3 = 0. Quindi questo caso rientra nelle già
H.4. COORDINATE ANGOLARI (CAPITOLO 3.3) 1009
citate configurazioni singolari.
Le equazione per la cinematica inversa che permettono il calcolo di β e γ si ottengono con il seguente
procedimento.
L’eq.(2.45) rappresenta il seguente sistema di equazioni

x2 + y 2 = (l1 + (l2 + l3 cos(γ)) cos(β))2 + (l3 sin(γ))2
(H.1)
z = (l2 + l3 cos(γ)) sin(β)
Aggiungendo z 2 alla prima si ottiene

 2
 x + y2 + z 2 = l12 + l22 + l32 + 2l1 (l2 + l3 cos(γ)) cos(β) + 2l2 l3 cos(γ)
z (H.2)
 − l2 = l3 cos(γ)
sin(β)
e per sin(β) 6= 0 si ottiene
 2

 (x + y 2 + z 2 − l12 + l22 − l32 ) sin(β) − 2l1 z cos(β) − 2l2 z = 0
 !
z l2 (H.3)


 γ = ± acos l sin(β) − l
3 3
e quindi si giunge al nuovo sistema che sostituisce eq.(2.47)

 2

 tb (l1 − l2 )z + tb (x2 + y 2 + z 2 − l12 + l22 − l32 ) − (l1 + l2 )z = 0
 !
z l (H.4)

 γ = ± acos −
2
 l3 sin(β) l3
dalla prima equazione di secondo grado si determina tb . In generale si avranno due soluzioni (eventual-
mente coincidenti). Detta equazione può degenerare in una di primo grado se l1 = l2 o se z = 0. Per
ognuna delle soluzioni trovate, l’angolo β si determina come β = 2 atan(tb ).
Noto β, se sin(β) 6= 0, è immediato ricavare γ.
Se sin(β) = 0 allora anche z = 0, ed il manipolatore si comporta come un robot SCARA in cui il primo
link è lungo l1 + l2 ed il secondo è lungo l3 . Il robot SCARA è trattato al §2.9.2.
H.4 coordinate angolari (capitolo 3.3)

La difficoltà nello scegliere un gruppo di tre variabili angolari che descriva “bene” la posizione angolare di
un corpo deriva dal principio di non integrabilità della velocità angolare. Cioè si dimostra che l’integrale
della velocità angolare non ho particolare significato fisico. Come conseguenza di ciò non è possibile
definire tre coordinate angolari la cui derivata rispetto al tempo coincida con il vettore velocità angolare.
Questo concetto è illustrato dalla figura seguente ove si confrontano due situazioni. Nella prima un dado
è inizialmente ruotato di 90◦ intorno a Z, poi di 90◦ intorno a Y . Nella seconda situazione le rotazioni
avvengono in ordine opposto. Si osservi come l’orientamento finale del dado sia diverso nei due casi
mentre l’integrale della velocità angolare è identico.
H.5 Stanford Arm (Pagg.135-137)

Lo Stanford Arm è descritto dai seguenti parametri di Denavit e Hartenberg:
n.link giunto θ d a ϕ
1 R q1 0 0 −π/2 rotazione della torretta attorno all’asse z0
2 R q2 d2 0 π/2 elevazione del braccio attorno all’asse z1
3 P 0 q3 0 0 allungamento del braccio lungo z2
4 R q4 0 0 −π/2 prima rotazione giunto sferico attorno a z3 = z2
5 R q5 0 0 π/2 flessione polso sferico attorno a z4
6 R q6 0 0 0 terza rotazione giunto sferico attorno a z5
wy
t
wz
z
Rot(y)
Rot(z)
t
(b) Caso 1: prima Rot(z), poi Rot(y)
x y wy
Rot(y) Rot(z)
(a) Due diverse sequenze di rotazioni di un corpo rigido attorno agli assi y e z. t
In entrambi
R i casi si ha: R wz
Rot(x) = 0; Rot(y) = ωy dt = 90◦ Rot(z) = ωz dt = 90◦
t
(c) Caso 2: prima Rot(y), poi Rot(z)
Figura H.15: Effetto di due differenti sequenze di rotazione di 90◦ .
H.6 §6.2.3 Sistemi di riferimento (Pag. 136-138)

La figura H.16 illustra il posizionamento delle terne tecniche previste dalla norma UNI EN ISO 9787.
Sono definite:
(0) a terna del laboratorio (non completamente definita) numerata con 0 ed asse z verticale;
(1) la terna di base numerata con 1, asse z diretta verso il manipolatore, asse x verso il centro dello
spazio di lavoro (CW = Center of Workspace, stà nel piano x0 z0 );
(m) dell’interfaccia meccanica (con m = n + 1, n=numero dei gradi di libertà), quando i motori si
trovano circa a metà corsa i piani xz della terna (0) e (m) coincidono o sono paralleli, asse z uscente
dall’interfaccia;
(t) dell’utensile (tool), non completamente definito, generalmente con asse y in direzione del movimento
delle dita (se presenti), ed asse z uscente dalla pinza (utensile).
Il centro dello spazio di lavoro è il punto raggiunto dal punto di riferimento del polso quando gli assi
principali sono (circa) a metà corsa.
H.7 Elbow Arm (Pagg.149-153)

L’Elbow Arm è descritto dai seguenti parametri di Denavit e Hartenberg:
n.link giunto θ d a ϕ
1 R q1 0 0 π/2 rotazione della torretta attorno all’asse z0
2 R q2 0 a2 0 elevazione del braccio attorno all’asse z1
3 R q3 0 a3 0 elevazione dell’ avambraccio attorno z2
4 R q4 0 a4 −π/2 prima rotazione giunto sferico attorno a z3
5 R q5 0 0 π/2 flessione polso sferico attorno a z4
6 R q6 0 0 0 terza rotazione giunto sferico attorno a z5
Enumerazione delle 8 soluzioni dell’Elbow Arm
braccio TCP gomito k1 k2 k3 ϑ1 ϑ2 ϑ3 ϑ234 ϑ4 ϑ5 ϑ6
1 Avanti Fuori giù 0 0 1 ϑa ϑb1 ϑc1 ϑg ϑd1 ϑe ϑf
2 Avanti Fuori su 0 0 -1 ϑa ϑb2 −ϑc1 ϑg ϑd2 ϑe ϑf
3 Avanti Dentro giù 0 1 1 ϑa ϑb3 ϑc2 ϑg ± π ϑd3 -ϑe ϑf − π
4 Avanti Dentro su 0 1 -1 ϑa ϑb4 −ϑc2 ϑg ± π ϑd4 -ϑe ϑf − π
5 Dietro Fuori su 1 0 1 π − ϑa π − ϑb2 ϑc1 ϑg −ϑd2 π − ϑe ϑf − π
6 Dietro Fuori giù 1 0 -1 π − ϑa π − ϑb1 −ϑc1 ϑg −ϑd1 π − ϑe ϑf − π
7 Dietro Dentro su 1 1 1 π − ϑa π − ϑb4 ϑc2 ϑg ± π −ϑd4 −π + ϑe ϑf
8 Dietro Dentro giù 1 1 -1 π − ϑa π − ϑb3 −ϑc2 ϑg ± π −ϑd3 −π + ϑe ϑf
I valori di k1 , k2 , k3 sono quelli utilizzati nella selezione delle soluzioni:
ϑ1 = atan2(Ty , Tx ) + k1 π, k1 = 0, 1 ϑ234 = atan2(Zz , c1 Zx + s1 Zy ) + k2 π, k2 = 0, 1ϑ3 = k3 acos(..), k3 = −1, 1
H.8. METODO DEL GRAFO DI PAUL; METODO SCISSIONE CATENA CINEMATICA (§6.7.1 E §6.7.2 PAGG.154-155)1011
z
y
(m)
z
x
z
y x
Sistema coordinate
generale
(1)
(0)
y
(1) x
(t) z cw
Figura H.16: Esempi di posizionamento delle terne tecniche previste dalla norma UNI EN ISO 9787: del
laboratorio (0), della base (1), dell’interfaccia meccanica (m) e dell’utensile (t).
Considerate le dimensioni geometriche di uno specifico manipolatore ed una posa da raggiungere non
è detto che tutte le otto soluzioni elencate siano matematicamente ammissibili. Da un punto di vista
matematico la non raggiungibiltà di un punto viene verificata se l’argomento dell’arcocoseno ha modulo
maggiore di uno. Per questo tipo di manipolatore si possono avere nessuna, 4 o 8 soluzioni.
Una volta identificata una soluzione, all’interno del gruppo contenete le numero 1, 2, 5 e 6 (oppure del
gruppo 3, 4, 7 e 8), si può affermare che anche le altre del gruppo sono ammissibili. Ciò si verifica imme-
diatamente in maniera geometrica osservando che per ogni soluzione con gomito alto esiste anche quella
corrispondente con gomito basso e viceversa (ad es. soluzioni 1 e 2); inoltre per ciascuna di esse una
volta considerata una soluzione, un’altra si ottiene ruotando la torretta (angolo ϑ1 ) di 180◦ e ribaltando
il braccio da davanti a dietro o viceversa (ad es. da soluzione 2 a 5). Non è detto invece che se esiste
una soluzione con TCP in avanti esista quella in dietro o viceversa; infatti il polso non può avvicinarsi
ne allontanarsi troppo dalla torretta. A questo proposito si osservi che per determinare quali punti siano
raggiungibili dal polso possiamo osservare che i link 2 e 3 movimentati da ϑ2 e ϑ3 definiscono una zona
di piano simile allo spazio di lavoro di un robot SCARA (vedi figura 1.21).
Ovviamente non è detto che tutte le soluzioni matematicamente ammissibili siano poi raggiungibili
fisicamente infatti ogni grado di libertà avrà in generale limitazioni di fine corsa.
H.8 Metodo del grafo di Paul; metodo scissione catena cinema-

tica (§6.7.1 e §6.7.2 Pagg.154-155)
Come esempio si consideri la soluzione analitica del manipolatore Elbow (§6.6 pagg.149-153).
Tutte le volte che una coppia di equazioni è stata sommata dopo aver moltiplicata una di esse per il
seno di un’angolo e l’altra per il coseno dello stesso angolo si è ottenuto lo stesso effetto che si sarebbe
ottenuto trasportando una delle matrici di posizione dal primo al secondo membro. Ad esempio, una volta
determinato l’angolo ϑ1 , la matrice M01 può essere portata al secondo membro ottenendo la relazione
2, 5 2, 5
replacemen
a a
4, 7 TCP 4, 7 TCP
ϑ2 ϑ2
b b
3, 8 3, 8
ϑ1 1, 6 ϑ1 1, 6
Fig. 6.14 Le otto soluzioni del robot 6R (Elbow Arm).

Per ϑ1 = ϑa (braccio in avanti) il manipolatore può raggiungere le configurazioni 1, 2, 3 e 4, mentre per
ϑ1 = π − ϑa (indietro) può raggiungere le altre (5, 6, 7 e 8). Per ciascuno dei due casi precedenti il TCP
in confronto del centro del polso può stare in fuori (a) o in dentro (b). Infine per ciascuna dei 4 casi
citati il gomito può stare in alto o in basso.
matriciale
−1
M12 ...M56 = M01 M06 = M ∗
dalla quale è più facile determinare alcune relazioni tra le quali quelle dedotte dalle eq. 7bis e 8bis o quella
che porta alla determinazione di ϑ234 ove compare c1 Zx + s1 Zy che deriva dal termine di posizione (1,3)
di M ∗ . Inoltre nella soluzione del manipolatore compare il termine c1 Xx + s1 Xy che deriva dall’elemento
di posizione (1,1) di M ∗ mentre c1 Tx + s1 Ty dalla posizione (1,4). Infine il termine s1 Zx − c1 Zy deriva
dalla posizione (3,3).
H.9 Metodo dei vincoli equivalenti (§6.7.1 e §6.7.3 Pagg.155)

Come esempio si consideri la soluzione analitica del manipolatore Elbow (§6.6 pagg.149-153) dove per
determinare l’angolo ϑ3 si è fatto uso di alcune sostituzioni temporanee di variabili denominate ∆x e ∆y.
Si può verificare che ciò si otterrebbe sostituendo la parte di struttura che comprende i giunti 2, 3 e 4, e
quindi i link 2 e 3, con una struttura composta da due accoppiamenti prismatici con movimento pari a
∆x e ∆y seguiti da un accoppiamento rotoidale che effettui la rotazione ϑ234 . In altre parole si immagina
di sostituire la struttura di sinistra di figura 6.18 con quella di destra ponendo ∆x = a2 c2 + a3 c23 e
∆y = a2 s2 + a3 s23 . Quindi i valori di ϑ2 e ϑ3 si ottengono osservando che il movimento dei link 2 e 3
nel piano verticale contenente l’asse z0 ed i link in oggetto ricordano la cinematica di un manipolatore
SCARA. I termini c234 a4 e s234 a4 sono le proiezioni nel piano orizzontale e verticale dell’ultimo link.
H.10 Scelta coordinate, cinematica inversa numerica (Pag.170

§6.9.7)
In pratica alla i-esima iterazione occorre risolvere il seguente sistema lineare
 
  dq1  
bx1 bx2 bxi bxn   b∗x
   dq2
  
 by1 by2 byi byn    
.. b∗y 
 bz1 bz2 bzi bzn    
. b∗z 
 ··· ···     ∆Qi = J −1 dS
 ax 1 ax 2 ax i ax n   dqi  =  a∗x  ⇔ J∆Qi = dS ⇒
     
 ay 1 ay 2 ay i ay n   ..   a∗y 
 . 
az 1 az 2 u zi az n a∗z
dqn
 
bxi  
   byi  0 −azi ayi bxi
   
 bzi   azi 0 −axi byi 
J =  J1 J2 . . . Ji . . . Jn  Ji = 

 = F (Li )
 Li = 



 axi   −ayi axi 0 bzi 
 ayi 
0 0 0 0
azi
H.11. MOMENTI D’INERZIA - ESEMPI (PAG.201) 1013
ove la i-esima colonna della matrice Jacobiana è ottenuta scegliendo ordinatamente i sei elementi non
nulli di Li , mentre la colonna dei termini noti dS è ottenuta in maniera analoga da L∗ .
La matrice inversa può essere proficuamente sostituita con la pseudoinversa per il caso di n 6= 6 o di
posizioni vicine a configurazioni singolari.
Step risolutivi della cinematica inversa
1. leggi M on % posa da raggiungere
2. poni Q = Q0 % valore di primo tentativo
3. poni j = 0 % contatore iterazioni
4. iterazione:
5. poni M = I
6. per i = 1, 2, . . . , n % cinematica diretta (per ogni giunto)
7. calcola Mr = M (qi , Λi ) % matrice di posa relativa del link i (Mi−1,i )
8. calcola L = L(Λi ) % matrice dell’asse del giunto i in terna locale (Li(i−1) )
9. calcola L0 = M LM −1 % matrice dell’asse del giunto in terna 0 (Li(0) )
10. calcola M = M Mr % matrice di posa assoluta del link i (M0,i )
11. calcola J i = F (L0 ) % colonna dello jacobiano
12. fine ciclo i
13. calcola ∆M = M on − M % errore di posa
14. se ||∆M || < ε % convergenza raggiunta (fine algoritmo): Q calcolato correttamente
15. poni j = j + 1 % conta iterazioni
16. se j > jmax % il metodo non converge (fine algoritmo): Q non calcolato correttamente
17. calcola L∗ = ∆M M −1
18. calcola dS = F (L∗ ) % termine noto
19. calcola ∆Q = J −1 dS % incremento variabili di giunto
20. calcola Q = Q + ∆Q % valore di nuovo tentativo
21. ripeti da iterazione
La funzione F (L) trasforma una matrice L in una matrice colonna scegliendo ordinatamente i 6 elementi
significativi; con il simbolo Λi si è indicato l’insieme dei parametri geometrici del giunto i-esimo.
Si osserva che vale:
L∗ = M 0n − M M −1 = ∆M M −1 = M 0n M −1 − I
H.11 Momenti d’inerzia - esempi (pag.201)

Con riferimento alla figura H.17, ed indicando con m = vρ la massa, con v il volume, con ρ la densità i
momenti d’inerzia del prisma con le notazioni J ed I sono (m = vρ, v = ABL)
1 1 1
Jx = m(B 2 + L2 ) Jy = m(A2 + L2 ) Jz = m(A2 + B 2 )
12 12 12
1 1 1
Ixx = mA2 Iyy = mB 2 Izz = mL2
12 12 12
(H.5)
1 1 1 1 1
Jx ′ = mB 2 + mL2 Jy ′ = mA2 + mL2 Jz ′ = m(A2 + B 2 )
12 3 12 3 12
1 1 1
Ix′ x′ = mA2 Iy′ y′ = mB 2 Iz′ z′ = mL2
12 12 3
x x
′ ′
x x
L
L
A R
z = z′ z = z′
B
y′ y y′ y
Figura H.17: Esempi di forme geometriche semplici e notazione per il calcolo dei momenti d’inerzia
rispetto a vari assi (xyz sono assi baricentrali).
2
1
Figura H.18: Esempi di forme geometriche complesse composte da forme semplici.
L L
a B R
r
A
Figura H.19: Esempi di forme geometriche semplici cave.
mentre per il cilindro si ottiene (m = vρ, v = πR2 L)
1 1 1 1 1
Jx = mR2 + mL2 Jy = mR2 + mL2 Jz = mR2
4 12 4 12 2
1 1 1
Ixx = mR2 Iyy = mR2 Izz = mL2
4 4 12
(H.6)
1 1 1 1 1
Jx′ = mR2 + mL2 Jy′ = mR2 + mL2 Jz′ = mR2
4 3 4 3 2
1 1 1
Ix′ x′ = mR2 Iy′ y′ = mR2 Iz′ z′ = mL2
4 4 3
Nel caso di corpi di forma complessa, ma assimilabili alla somma di corpi semplici, è sufficiente sommare
le inerzie dopo averle espresse rispetto alla medesima terna d’assi (ad es. fig. H.18).
Mentre nel caso di corpi cavi fig. H.19), i momenti d’inerzia si calcolano immaginado il corpo pieno e
sottraendo l’inerzia corrispondente al corpo vuoto.
Ad esempio per il prisma si ottiene
1 1
Jz′ = ρABL(A2 + B 2 ) − ρabL(a2 + b2 ) (H.7)
12 12
ed anche
ρ ρ 1 1
Ix′ ′ x′ = L BA3 − ba3 Iy′ ′ y′ = L AB 3 − ab3 Iz′ ′ z′ = mL2 = ρABL3 (H.8)
12 12 3 3
ed analogamente per gli altri momenti d’inerzia.
H.12. LEGGE DI MOTO “AD ACCELERAZIONE COSTANTE” (PAG.238 FIG.9.2, PAG.249 FIG.9.18) 1015
Nel caso del cilindro pieno il momento d’inerzia Jz vale
1 1 π
Jz = mR2 = πR2 LρR2 = ρLR4 (H.9)
2 2 2
per cui il momento d’inerzia del cilindro cavo vale
π
Jz′ = ρL R4 − r4 (H.10)
2
e quindi anche
π 1 π
′
Ixx ′
= Iyy = ρL R4 − r4 ′
Izz = mL2 = ρR2 L3 (H.11)
4 12 12
H.12 Legge di Moto “ad accelerazione costante” (Pag.238 fig.9.2,

pag.249 fig.9.18)
Consideriamo una legge di moto “ad accelerazione costante” che preveda uno spostamento ∆s da effet-
tuarsi nel tempo T .
Con riferimento alle figure ed alle formule delle pagine citate, possiamo scrivere le seguenti equivalenze
∆t1 = t1 − t0 λ1 = ∆t1 /T primo tratto;

∆t2 = t2 − t1 secondo tratto;
∆t3 = t3 − t2 λ2 = ∆t3 /T terzo tratto;
T = t3 − t0 durata totale;
A, D valore assoluto dell’accelerazione e della decelerazione;
vm velocità massima, quella cioè nel secondo tratto.
Valgono i vincoli 0 < λ1 < 1, 0 < λ2 < 1, 0 < λ1 + λ2 ≤ 1.
Si ottiene
∆s vm ∆s 1 vm ∆s 1
Vm = ! A= = 2 ! D= = 2 !
λ1 + λ2 λ1 T T λ1 + λ2 λ2 T T λ1 + λ2
T 1− λ1 1− λ2 1−
2 2 2
quindi i coefficienti di velocità nonché quelli di accelerazione positiva e negativa valgono
1 1 1
cv = ! Ca+ = ! Ca− = !
λ1 + λ2 λ1 + λ2 λ1 + λ2
1− λ1 1− λ2 1−
2 2 2
Nel caso di leggi di moto simmetriche si ha ovviamente λ = λ1 = λ2 , Ca = Ca+ = Ca− .

Le definizioni date di A accelerazione, D decelerazione, Ca+ e Ca− sono corrette nel caso di movimento
positivo (∆s > 0); nel caso di movimento negativo (∆s < 0) le formule di A e D si scambiano tra loro
come anche quelle di Ca+ e Ca− .
H.13 Traiettoria circolare per 3 punti (Pag. 264)

La rotazione è sempre positiva e l’asse di rotazione è calcolato con la formula
(P2 − P1 ) × (P0 − P1 )
z=
||(P2 − P1 ) × (P0 − P1 )||
L’ampiezza della rotazione da P0 a P2 passando per P1 si calcola in due passi. Innanzitutto, considerato
il triangolo P0 P1 P2 , si calcola l’angolo in P1 indicato in figura come φ′ dal prodotto scalare tra i vettori
(P2 − P1 ) e (P0 − P1 ). Infine l’angolo di rotazione φ si ottiene ricordando che, considerato un arco di
circonferenza, l’angolo alla circonferenza è la metà dell’angolo al centro. In questo caso interessa l’angolo
explementare !
(P 2 − P 1 ) · (P0 − P1 )
φ = 2π − 2φ′ = 2π − 2 acos
||(P2 − P1 )|| ||(P0 − P1 )||
P0 P0
P1 P0 P0
P1
f’
P2 P2
2f’ f=2p-2f’ P2 Da2 a0 P2 Da2 a0
a1 Da1
Pc a2 a1
2f’ x a2 x
f=2p-2f’ f’ Pc Pc
P1 Da1
P1
Figura H.20: Segmento circolare per 3 punti: Figura H.21: Segmento circolare per 3 punti:
caso generale 3D. L’asse z è entrante nel foglio caso piano xy. A sinistra rotazione antioraria,
nell’esempio di sinistra, uscente in quello di a destra oraria.
destra.
P2
P1
Figura H.22: Movimento in traiettoria assegnata da P1 a P2 .
Casi singolari
L’arco da percorrere non è definito e la procedura non è applicabile se i tre punti risultano allineati o
se due di essi coincidono. In questo caso si può generare un segmento rettilineo. Queste condizioni si
identificano osservando che in questi casi si verifica almeno una delle due identità ||V1 ×V2 || = 0, V1 = ±V2 .
Cerchio nel piano xy

Nel caso si debba generare un arco di cerchio nel piano xy non è necessario utilizzare la procedura indicata,
ma si può più semplicemente implementare la procedura seguente. Innanzitutto si calcolano gli angoli αi
che rappresentano l’inclinazione dei vettori Pi − Pc rispetto all’asse delle x nonché la loro differenza con
α0
αi = atan2 (Pyi − Pyc , Pxi − Pxc ) i = 0, 1, 2
∆αi = mod (αi − α0 , 2π) i = 1, 2
dove l’operatore “mod” (modulo) indica che l’angolo va espresso nell’intervallo 0÷2π; quindi, ad esempio,
è necessario aggiungergli 2π se è negativo.
L’angolo effettivo di rotazione φ si calcola quindi con due diverse modalità assecondo che l’arco debba
essere percorso con verso orario (angolo positivo) o antiorario (negativo)

 ∆α2 ∆α2 > ∆α1 orario
φ= se

∆α2 − 2π ∆α2 < ∆α1 antiorario
risulta (
x = ρ cos(α + α0 ) 0≤α≤φ
y = ρ sin(α + α0 ) α = λ/ρ
H.14. MOVIMENTO IN TRAIETTORIA §9.8 1017
P2
y V P
y s
α P1
r
U yc β
α
Pc
x
xc x
(a) Rotazione (b) Traiettoria circolare
Figura H.23: Movimento in traiettoria circolare.
H.14 Movimento in traiettoria §9.8

Questa sezione riordina alcuni discorsi già discussi al paragrafo §9.8 del libro.
H.14.1 Movimento in traiettoria

Nel movimento in traiettoria il manipolatore deve muoversi da un punto iniziale P1 a uno finale P2
garantendo di rimanere su un percorso predefinito (fig H.22); ad es. segmento di retta, arco di cerchio,...
Secondo i casi la posizione di un punto generico P sarà indicato da due o tre coordinate cartesiane
 
x
x
P = oppure P = y  (H.12)
y
...
Il movimento è quindi caratterizzato dalla forma della traiettoria e dalla legge di moto con la quale la si
percorre.
Come esempio di paragone consideriamo una strada con rettilinei e curve (la traiettoria) su cui una
automobile si muove con una certa legge di moto indicata dal variare del contachilometri e dalle letture
del tachimetro.
Nel caso del robot il problema consiste nel calcolare la legge di moto da dare ai motori in modo che
la pinza del robot segua la traiettoria prevista con la legge di moto desiderata. Il problema può essere
scomposto nei sotto problemi:
• scelta della forma della traiettoria P (λ);
• scelta della legge di moto di percorrenza della traiettoria λ(t), λ̇(t), λ̈(t);
• calcolo del movimento della pinza del robot S(t), Ṡ(t), S̈(t);
• calcolo del corrispondente movimento dei motori Q(t), Q̇(t), Q̈(t) tramite la cinematica inversa.
Consideriamo il caso semplice in cui la traiettoria sia esprimibile in forma parametrica, ovvero la posizione
di ogni punto della traiettoria sia esprimibile in funzione di un parametro λ. Particolarmente semplice è
il caso in cui λ rappresenta la distanza percorsa; in questo caso talvolta indicheremo questo parametro
con s.
Una volta che sia nota la variazione di λ con il tempo, la posizione, la velocità e l’accelerazione del punto
saranno:
P = P (λ)
∂P
Ṗ = λ̇ λ = f (t) (H.13)
∂λ
∂P ∂2P
P̈ = λ̈ + λ̇
∂λ ∂λ2
ove f (t) indica la legge di moto di percorrenza che definisce il valore di λ in funzione del tempo t.
Nel nostro caso siccome ci interessa il movimento della pinza del robot porremo P = S. In generale la
legge di moto da assegnare ai motori si ottiene
• assegnando la traiettoria P (λ);
• scegliendo una legge di moto f (t) che nel tempo scelto faccia variare λ da λ1 (cioè punto P1 ) a λ2
(cioè punto P2 );
• scegliendo un opportuno intervallo temporale ∆T , campionando nel tempo la legge λ = f (t) e per
ogni punto calcolando posizione, velocità ed accelerazione della pinza;
λ = f (t), S = P (λ), Ṡ = . . ., S̈ = . . .
• infine calcolando la cinematica inversa per ciascuno

dei punti.

Q = G(S) Q̇ = J −1 Ṡ Q̈ = J −1 S̈ − J˙Q̇
Verranno analizzati i casi in cui esiste una semplice rappresentazione parametrica delle traiettorie in
funzione dell’ascissa curvilinea: retta, arco di cerchio, elicoide.
H.14.2 Traiettoria rettilinea

Scegliamo come parametro l’ascissa curvilinea λ = s, e cioè la distanza da P1 . La lunghezza della
traiettoria L è ovviamente L = ||P1 − P2 ||. Occorre quindi scegliere una legge di moto s = f (t, T, L) che
nel tempo assegnato T faccia variare s da 0 a L. Si può utilizzare ad esempio una legge cicloidale, tre
tratti, o altro.
Applicando le equazioni (H.13) alla retta, si verifica facilmente che risulta
P2 − P1
U=
S = P (s) = P1 + U s ||P2 − P1 ||
Ṡ = Ṗ = U ṡ (H.14)
0 ≤ s = f (t, T, L) ≤ L
S̈ = P̈ = U s̈
0≤t≤T
H.14.3 Traiettoria circolare in un piano coordinato

Rotazione
Consideriamo in via propedeutica il semplice caso in cui si voglia ruotare di un angolo α un vettore U
giacente nel piano xy ottenendo V (fig. H.23(a)). Si dimostra che questa operazione può essere fatta
moltiplicandolo per una matrice di rotazione R

Vx cos(α) − sin(α) Ux
V = RU = (H.15)
Vy sin(α) cos(α) Uy
Questa relazione può essere spiegata indicando con β la direzione del vetture U , ed osservando che la
direzione di V è ovviamente α + β. A questo punto le componenti di V sono proporzionali al coseno
ed al seno di α + β e quindi ricordando le funzioni trigonometriche della somma di angoli si giunge alla
dimostrazione.
Generazione del movimento circolare
Con riferimento alla figura H.23(b), scegliamo come parametro l’ascissa curvilinea λ = s, e cioè la distanza
da P1 ed assumiamo che il moto avvenga nel piano xy. Indichiamo con Pc il centro del cerchio, con r il
raggio. Vale ovviamente s = αr, α = s/r.
La posizione di un punto generico si ottiene ruotando il vettore P1 − Pc dell’opportuno angolo α
P = P1 + R(P1 − Pc ) (H.16)
H.14. MOVIMENTO IN TRAIETTORIA §9.8 1019
h = 2πp
Figura H.24: Movimento in traiettoria elicoidale.
Quindi scegliendo una legge di moto opporuna che faccia variare s da 0 a L (distanza tra i due punti P1
e P2 misurata sulla circonferenza) nel tempo prefissato T , si ottiene
P = Pc + R(P1 − Pc )
R = R(α)
∂
Ṗ = R (P1 − Pc )α̇ α = s/r (H.17)
∂α
∂ ∂2 s = f (t, T, L)
Ṗ = R (P1 − Pc )α̈ + R (P1 − Pc )α̇2
∂α ∂α2
con
cos(α) − sin(α)
R= (H.18)
sin(α) cos(α)

∂ sin(α) cos(α) ∂2 − cos(α) sin(α)
R= R= (H.19)
∂α − cos(α) − sin(α) ∂α2 − sin(α) − cos(α)
H.14.4 Traiettoria elicoidale

Una traiettoria elicoidale (figura H.24) è derivabile da una traiettoria circolare aggiungendo un avanza-
mento nel piano ortogonale a quello della rotazione. Se il movimento circolare avviene nel piano xy e
quello di avanzamento è in direzione z, allora le generiche equazioni si possono derivare da quelle del
moto circolare aggiungendo la terza coordinata con spostamento proporzionale alla rotazione
     
x xc + r cos(α + α0 ) 0
P =  y  =  yc + r sin(α + α0 )  = Pc + R(P1 − Pc ) +  0  (H.20)
z zc + pα αp
con  
cos(α) − sin(α) 0
R =  sin(α) cos(α) 0  (H.21)
0 0 1
ove p è il passo dell’elica che si ottiene dividendo lo spazio percorso in un giro lungo z per l’angolo
corrispondente (2π ragionando in radianti) cioè p = h/(2π), e Pc è un punto qualsiasi dell’asse dell’elicoide.
La lunghezza di un tratto di elicoide si ottiene dal teorema di Pitagora sommando lo spostamento nel
piano αr, con quello verticale αp ottenendo
p
s = α r 2 + p2 (H.22)
questa relazione si dimostra immaginando che l’elicoide sia avvolto su un cilindro di raggio r, e immagi-
nando di srotolare il cilindro stesso su di un piano (figura H.25)
H.14.5 Traiettorie circolari ed elicoidali con giacitura generale

Qualora le traiettorie circolari o elicoidali avvengano secondo piani o direzioni non parallele a piani coordi-
nati, è possibile sviluppare il movimento in un sistema di riferimento locale con assi coordinati opportuni e
poi effettuare una rototraslazione dei vari punti calcolati per riportarli nel sistema di riferimento generale.
αp
s
αr
Figura H.25: Calcolo della lunghezza di un’elica.
H.15 Effetto elasticità nelle trasmissioni (Pag.294 §10.5.5)

H.15.1 Introduzione, misura diretta e indiretta di posizione
Si considera il caso in cui si voglia controllare il movimento di un organo meccanico che si trova al ter-
mine di una catena cinematica che include elementi elastici. In questo caso il movimento dell’organo
in questione può essere significativamente differente da quello che si potrebbe dedurre immaginando la
trasmissione perfettamente rigida.
Consideriamo il caso in cui l’elasticità si possa ritenere concentrata in un punto.
Facciamo riferimento alle figure H.26, H.27 e H.28
La retroazione di velocità è sempre ottenuta misurando il movimento dell’albero motore, mentre quella
di posizione può basarsi sulla misura del movimento del motore o su quella del carico. Nel primo caso
si parla di misura indiretta di posizione, nel secondo di misura diretta. In caso di trasmissione ideale
e rigida, con misura diretta o indiretta di posizione, si otterrebbero risultati identici tuttavia nei casi
pratici si possono ottenere risultati molto diversi.
La retroazione con misura diretta permette di avere maggiore precisione di movimento perché la retroazio-
ne tiene conto del movimento effettivo del carico compensando eventuali errori nella trasmissione dovuti
a imperfezione geometrica, giochi o deformazione dovute dai carichi. Tuttavia rischia di rendere instabile
il sistema nel caso in cui la trasmissione abbia cedevolezza significativa o giochi e sia presente un elevato
attrito statico (o a basse velocità) che potrebbero determinare il fenomeno dello stick-slip.
Inoltre, nel caso di movimenti rotatori con misura tramite encoder, il posizionamento del trasduttore
sull’albero motore da una migliore risoluzione di misura a causa della riduzione e consente l’impiego di
trasduttori con meno tacche al giro. Infatti un encoder con un numero Ns di passi al giro garantisce una
risoluzione angolare pari a ∆ϑ = 2π/(4Ns ). Il fattore 4 è valido per i normali encoder incrementali con
canali in quadratura. Se tra encoder e carico è presente un riduttore di velocità con riduzione pari a η
si ottiene una risuluzione più elevata e pari a ∆ϑ = 2πτ /(4Ns ) (generalmente τ < 1). Quindi per avere
una prefissata risoluzione occorre un numero di passi Ns ≥ 2π̇/(4τ ∆ϑ).
Infine nel caso di moderni motori brushless con resolver, l’adozione di una misura indiretta consente il
risparmio di un trasduttore in quando il resolver consente sia la misura della posizione che della velocità.
In ogni caso, per motivi costruttivi, nei manipolatori industriali non è quasi mai possibile realizzare una
misura diretta di posizione, mentre lo è quasi sempre negli assi a geometria cartesiana delle macchine
utensili ove una riga ottica (o altro dispositivo di misura) può essere agevolmente accostato al carro
mobile (figura H.28).
La decisione sull’utilizzo della misura diretta o indiretta di posizione va fatta di volta in volta.
H.15.2 Modello di controllo con trasmissione elastica

Lo schema di riferimento è quello di figura H.26 che come vedremo può essere trasformato in quello
equivalente di figura H.27. I simboli hanno il seguente significato
H.15. EFFETTO ELASTICITÀ NELLE TRASMISSIONI (PAG.294 §10.5.5) 1021
J J
C Cr
ϑ0 + + r, k
R′ (s) R(s) motore D E τ E2
− −
Ω1
Ω Ω2
Ω
ϑ ϑ1 ϑ2
1/τ
Figura H.26: Schema di un azionamento ad un grado di libertà con riduttore di velocità ed elasticità
nella trasmissione e possibili schemi di retroazione.
J J1
C Cr
ϑ0 + Ω 0 +
r, k
R′ (s) R(s) motore D E E2
− −
Ω Ω1
Ω
ϑ ϑ1
ϑ ϑ1
Figura H.27: Schema equivalente di un azionamento ad un grado di libertà senza riduttore di velocità
(τ = 1) ed elasticità nella trasmissione e possibili schemi di retroazione.
ϑ posizione angolare del motore;

ϑ1 = ϑ2 /τ posizione angolare equivalente del carico per τ = 1;
ϑ2 posizione angolare del carico;
Ω = ϑ̇ velocità del motore;
Ω1 = ϑ̇1 velocità del carico dopo la connessione elastica equivalente;
Ω2 = ϑ̇2 velocità del carico;
Ω2 ϑ2
τ= = rapporto di trasmissione del riduttore di velocità (supposto idealmente rigido);
Ω1 ϑ1
J momento di inerzia di motore e trasmissione;
J momento di inerzia del carico (braccio del robot) a valle del riduttore;
J1 = τ2 J inerzia equivalente del carico per trasmissione unitaria;
k rigidezza equivalente della trasmissione (albero + riduttore);
r costante di smorzamento;
C coppia motrice;
C2r coppia resistente;
Cr = τ C2r coppia resistente equivalente;
In assenza di elasticità (trasmissione rigida) si avrebbe ϑ = ϑ1 , Ω = Ω1 , Ω2 = τ Ω.
Per semplificare i calcoli è possibile eliminare il rapporto di trasmissione τ pur di sostituire i valori di
coppia resistente, inerzia, posizione angolare, velocità e accelerazione con quelli equivalenti; si giunge cosı̀
allo schema equivalente di figura H.27. Nell’esempio seguente viene trattato il caso di carico con moto
rotatorio, ma è evidente l’estensione al caso di moto traslatorio ove compaiono masse al posto di momenti
d’inerzia, forze al posto di coppie, traslazioni al posto di rotazioni (figura H.29).
Con riferimento alle figure H.27 e H.29, lo sviluppo delle equazioni necessarie allo studio può essere
ottenuto a partire dall’equilibrio alla rotazione delle masse solidali al motore e quelle del carico i cui
movimenti sono tra loro accoppiati dalla trasmissione tramite la quale si scambiano coppie elastiche e
viscose

 C = k (ϑ1 − ϑ) + r ϑ̇1 − ϑ̇ + J ϑ̈ masse su albero motore

 Cr = k (ϑ − ϑ1 ) + r ϑ̇ − ϑ̇1 + J1 ϑ̈1 masse su asse del carico
Passando al dominio delle trasformate di Laplace, con semplici calcoli, si può giungere allo schema a
blocchi ed alle funzioni di trasferimento di figura H.30. Si tratta di uno schema con due ingressi (la
coppia motrice C e quella resistente Cr ) e di due uscite (la posizione del motore ϑ e quella equivalente del
carico ϑ1 ). Questo schema può essere facilmente convertito in quello ad esso equivalente di figura H.31
nel quale sono state aggiunte linee tratteggiate che indicano un classico schema di controllo con anelli di
x=pa
coppia coppia
w D M E a k
motrice resistente
convertitore riga (ottica) J J1
i, C C Cr
(azionamento)
J ,J&,J&& r J1 ,J&1 ,J&&1
wo
diretta
CNC o plc o
controllore Figura H.29: Schema equivalente lineare di una
indiretta trasmissione con elasticà.
Figura H.28: Retroazione tramite misura diretta o indiretta di

posizione (schema di principio).
retroazione di velocità e posizione annidati. Viene considerato sia il caso di misura diretta che quello di
misura indiretta di posizione.
Le funzioni di trasferimento sulle quali è necessario indagare sono quelle che legano il movimento del
carico ϑ1 ai comandi ϑo (posizione desiderata) ed alla coppia resistente Cr .
H.15.3 Retroazione di velocità

Analizziamo innanzitutto la regolazione di velocità. Con riferimento allo schema di figura H.31 calcoliamo
la funzione di trasferimento d’anello L nel caso di regolatore proporzionale R = kp nonchè la funzione
che lega tra loro le rotazioni di motore e carico (figure H.32 e H.33)
s2 2ξ1 r
+ s+1 s+1
ϑ̇ ϑ̇ 1 ω12 ω1 ϑ1 ϑ̇1 ϑ̈1 rs + k k
L= = kp = kp H(s) = = = = =
Eϑ̇ C s (J + J1 ) s2 2ξ2 ϑ ϑ̇ ϑ̈ J1 s2 + rs + k s2 2ξ1
+ s+1 + s+1
ω22 ω2 ω12 ω1
ove le pulsazioni ω e i relativi smorzamenti ξ valgono
s
k r
ω1 = ξ1 = √
J
r 1 2kJ 1
k r J J1
ω2 = ξ2 = √ J∗ =
J∗ 2 kJ ∗ J + J1
si riconosce che ω1 è la pulsazione propria di vibrazione del carico immaginando che l’albero motore sia
bloccato, mentre ω2 è la pulsazione propria di vibrazione del sistema immaginando albero motore e carico
liberi (gli alberi vibrerebbero in controfase tra di loro).
Osserviamo le funzioni di trasferimento di figura H.32; notiamo una antirisonanza (modulo quasi nullo)
alla pulsazione ω1 ed un picco di risonanza (modulo molto elevato) alla pulsazione ω2 . Al crescere
della rigidezza k i due picchi si spostano verso destra. Per frequenze basse e inferiori alla pulsazione di
antirisonanza (ω << ω1 ) il sistema si comporta come se fosse rigido e quindi con inerzia totale pari alla
somma di quella dell’albero motore più quella equivalente del carico ϑ̇/C ≃ ϑ̇1 /C ≃ 1/(s(J + J1 )). Per
pulsazioni elevate e superiori a quella di risonanza (ω >> ω1 ) il sistema si comporta come se motore e
carico fossero sconnessi. Il motore si muove secondo la funzione di trasferimento ϑ̇/C ≃ 1/(sJ), ma il
carico rimane quasi fermo (vedi anche figura H.33).
La figura H.34 descrive lo schema retroazionato di regolazione della velocità. Utilizzando un regolatore
puramente proporzionale kp si ottiene la funzione di trasferimento di figura H.37. A sinistra è quello che
accade per bassi guadagni, a destra la situazione per alti guadagni.
Nel grafico in alto è riportata la funzione di trasferimento d’anello e quella che si ottiene per il sistema
retroazionato con ingresso il comando di velocità e uscita la velocità del motore.
Il grafico centrale rappresenta la funzione di trasferimento da albero motore a carico, l’ultimo grafico è la
funzione di trasferimento da comando a carico. Quest’ultima funzione di trasferimento è il prodotto delle
due precedentemente citate e graficamente si ottiene come somma dei due diagrammi che lo precedono
nella figura (i diagrammi sono in scala doppio logaritmica).
Osserviamo che per bassi guadagni abbiamo banda passante ridotta, ma il sistema si comporta come se
fosse rigido. Per alti guadagni, la banda passante non può comunque superare ω1 al cui valore si verifica
risonanza.
H.15. EFFETTO ELASTICITÀ NELLE TRASMISSIONI (PAG.294 §10.5.5) 1023
C
J1s 2 + rs + k + J&& = W
& 1 J& = W 1 J
JJ1s 2 + (J + J1)rs + k (J + J1 ) s s
+
rs + k
JJ1s 2 + (J + J1)rs + k (J + J1 )
rs + k
JJ1s 2 + (J + J1)rs + k (J + J1 )
+
Js 2 + rs + k J&&1 = W
& &
1 J1 = W1 J1
+ 1 1
Cr JJ1s 2 + (J + J1)rs + k (J + J1 ) s s
Figura H.30: Schema a blocchi e funzioni di trasferimento da coppia motrice C e coppia resistente Cr a
posizioni angolari dell’albero motore ϑ e del carico ϑ1 .
selezione J, J1
J0 EJ J&0 Ew C J1s 2 + rs + k + 1 1
R’ R JJ1s 2 + (J + J1)rs + k (J + J1 ) J&& = W
& s
J& = W
s J
+
rs + k
J1s 2 + rs + k
rs + k
J1s 2 + rs + k
+ J&&1 = W
& &
Js 2 + rs + k 1 1 J1 = W1 1 J1
JJ1s 2 + (J + J1)rs + k (J + J1 ) s s
Cr +
Figura H.31: Schema a blocchi equivalente a quello di figura H.30; tratteggiato reatroazione (diretta) di
velocita e (in)diretta di posizione.
1 J&
sJ C
w2 aumentando k
J1
w J
w
k w1
w1
1
s (J + J1 )
Figura H.32: La funzione di trasferimento da coppia Figura H.33: La funzione di trasferimento da

motrice a velocità del motore ϑ̇/C. movimento motore a movimento del carico.
J&0 EJ& C J&& 1 J& J&1

R A(s) H(s)
+ - s
Figura H.34: Anello di retroazione di velocità.
J0 EJ
&
J 0 J& 1 J J1 J0 EJ J& 0 J& J& 1 1 J1
R’ H R’ N H
s s
Figura H.35: Anello di retroazione di posizione con Figura H.36: Anello di retroazione di posizione con
misura indiretta. misura diretta.
H.15.4 Retroazione di posizione

L’anello di posizione può essere chiuso tramite misura diretta o indiretta di posizione. In ogni caso si
vuole evitare che la funzione di trasferimento includa delle risonanze.
Consideriamo innanzitutto il caso in cui l’anello di posizione venga chiuso sull’albero motore (misura
indiretta) con un regolatore proporzionale (figura H.35). Per evitare di avere poli complessi e coniugati
(risonanza) è necessario che l’anello di velocità non ne abbia ed è necessario tenere guadagni limitati
anche nell’anello di posizione (vedi caso trattato al §10.4.2 pag. 285-286).
In figura H.36 viene illustrato uno schema che può essere ottenuto nel caso di misura diretta di posizione
e che consente di avere una banda passante più elevata. Il blocco indicato con N è un notch filter (filtro
a gola) che ha la seguente funzione di trasferimento
s2 2ξn
2 + s+1
ωf ωf
N= ξn < 1 ξd > 1
s2 2ξd
+ s+1
ωf2 ωf
Numeratore e denominatore sono del secondo ordine con la stessa pulsazione ωf ; il numeratore ha smor-
zamento piccolo (< 1) mentre il denominatore ha smorzamento elevato (> 1). Si ottiene una funzione di
trasferimento quasi unitaria con una gola di antirisonanza che se opportunamente tarata in pulsazione e
smorzamento può compensare la risonanza dell’anello di velocità (figura H.38). Ciò consente di chiudere
l’anello di posizione sull’albero del carico, includendo al suo interno la trasmissione, e tenere una banda
passante più elevata.
Le caratteristiche più salienti dei due schemi illustrati sono evidenziate in tabella H.38.
H.15.5 Effetto delle coppie resistenti

???
???
misura di posizione indiretta diretta

retroazione su prima massa seconda massa
realizzabilità sempre non sempre
trasmissione, elasticità, giochi fuori dall’anello di retroazione all’interno dell’anello di retroa-
zione
effetto dell’attrito riduce la precisione pericolo di instabilità
banda passante limitata da pericolo di risonanza maggiore
limitazione rischio risonanza occorre ridurre banda passante eliminabile con notch filter o
regolatori evoluti
precisione di regolazione minore maggiore
sensibilità a coppie resistenti maggiore minore
Tabella H.1: Caratteristiche della regolazione di posizione con misura diretta ed indiretta di posizione.
H.16. CONTROLLO CENTRALIZZATO E DECENTRALIZZATO (PAG.298 §10.6.1) 1025
J&
kp alto J& 1
J& kp basso C
C J& 0
J&
J&0 w
w w1
J&
J& 0 wt su prima massa wt
banda elevata N w f = w1
J&1 J&1 wf
J& J&
w1 w1
J&1 J&
J&0 N &1
J&1 J0
J&0
w
w1 w1
wt il sistema si 2° massa può
comporta andare in
-1 -2
come rigido risonanza
Figura H.38: Retroazione di po-
sizione con misura diretta di po-
Figura H.37: Retroazione di velocità con R = kp . Dall’alto: fun- sizione. Il notch filter può com-
zione di trasferimento d’anello L = kp ϑ̇/C e dei comandi ϑ̇/ϑ̇o pensare l’effetto della risonanza
(velocità motore); da velocità motore a velocità carico; da coman- presente Nell’anello di velocità
do a velocità carico. A sinistra guadagno basso, a destra guadagno a causa del guadagno e banda
elevato. elevati.
x1o - e1 u1 x1
x10 - e1 u1 x1
R1 motore1
R1 motore1 +
+
x20 x2
+ e2 u2
R2 motore2 x2o x2
- + u2
e2 R2 motore2
-
Figura H.39: Regolatore decentralizzato (ad

Figura H.40: regolatore Centralizzato (ad assi
assi indipendenti).
interdipendenti).
H.16 Controllo centralizzato e decentralizzato (Pag.298 §10.6.1)

Vengono illustrati meglio i concetti che stanno alla base di controlli centralizzati e brevemente presentati
degli esempi.
Nei controlli decentralizzati viene generato un movimento di riferimento per ogni giunto e ciascuno di essi
è poi controllato in maniera indipendente dagli altri. In altre parole l’azione di controllo su ogni motore
dipende solo dal suo errore di inseguimento e non da quello degli altri attuatori.
Nei controlli centralizzati invece gli anelli di retroazione di due o più giunti interferiscono tra di loro.
Ad esempio in un sistema a 2 gradi di libertà indichiamo con x1 e x2 le grandezze da controllare, con x1o
e x2o i loro valori desiderati e con u1 e u2 le due azioni di controllo. Nel caso di controllore decentralizzato
u1 dipende solo da x1 e x1o (fig. H.39), mentre in un controllo centralizzato esso può dipendere anche da
x2 e x2o (ed analogamente per u2 ) (fig. H.40); cioè
controllo decentralizzato controllo centralizzato
u1 = f1 (x1 , x1o ) u1 = f1 (x1 , x1o , x2 , x2o )

u2 = f2 (x2 , x2o ) u2 = f2 (x1 , x1o , x2 , x2o )
Esempi di controlli centralizzati sono
• controlli a coppia precalcolata;
• controlli a dinamica inversa;

• controlli ad assi virtuali;

• controlli nello spazio di lavoro.
Il controllo ad assi virtuali e nello spazio di lavoro posso essere sostanzialmente molto simili tra loro; la
loro diversa denominazione deriva da diversi ambiti applicativi in cui sono utilizzati.
Analizziamo ora brevemente i controlli dinamica inversa e quelli a coppia precalcolata. Con riferimento
al §2.8 pag. 48, si assume che il comportamento dinamico del robot sia descritto dall’equazione

Fq = M Q̈ + V Q, Q̇
ove la matrice delle masse M dipende da Q ed il termine V (.) che dipende da Q e da Q̇ congloba anche
gli effetti dovuti al peso ed agli attriti.
Si immagina di avere a disposizione delle stime del modello dinamico del manipolatore e cioè delle matrici
c e Vb tali che
M
c≃M
M Vb (.) ≃ V (.)
è pertanto possibile stimare la coppia necessaria a fare avvenire un certo movimento

Fbq = McQ̈ + Vb Q, Q̇
Coppia precalcolata nello spazio dei giunti

Si potrebbe cosı̀ ipotizzare di controllare il manipolatore in anello aperto chiedendo ai motori di generare
la coppia Fbq . Se la stima dello stato iniziale del manipolatore (Q e Q̇) fosse perfetta e se il modello
dinamico fosse perfetto il manipolatore si muoverebbe esattamente come previsto. Tuttavia è inevitabile
la presenza di qualche imperfezione e quindi il manipolatore devierà dalla traiettoria prevista. Sarà quindi
necessario introdurre una retroazione di posizione e di velocità (vedi schema 1 figura H.41). Indichiamo
con Eq = Q0 − Q l’errore di posizione e con sEq la sua derivata che coincide con l’errore di velocità. Il
termine (Kp + sKv ) Eq costituisce quindi un regolatore di tipo PD che garantisce stabilità e limitato errore
di inseguimento. Se il modello dinamico è buono, questo schema garantisce un errore di inseguimento
ridotto rispetto a quello che si potrebbe ottenere con un PD senza il termine di coppia precalcolata.
Lo schema può anche essere ridisegnato in forma perfettamente equivalente come in schema 1bis (figura
H.43) che evidenzia come il termine non lineare di predizione di coppia possa essere visto come un termine
di compensazione in avanti (Feed Forward).
In queste condizioni l’equazione differenziale che descrive la dinamica del sistema si può scrivere sommando
i termini di coppia calcolata ed eguagliandoli a quelli entranti nel manipolatore

Mc(Q0 )Q̈0 + Vb Q0 , Q̇0 + Kp (Q0 − Q) + Kv Q̇0 − Q̇ = M (Q)Q̈ + V Q, Q̇

c(Q0 ) ≃ M (Q) e Vb Q0 , Q̇0 ≃ V Q, Q̇ giungiamo all’equazione
nelle ipotesi che M

M Q̈0 − Q̈ + Kv Q̇0 − Q̇ + Kp (Q0 − Q) = 0
Questa è un’equazione differenziale di ordine 2n ove n il numero dei gradi di libertà del manipolatore.
L’equazione non è lineare perché M dipende da Q in maniera non lineare. Inoltre i movimenti dei vari
gradi di libertà sono tra loro accoppiati perché la matrice M non è diagonale.
Dinamica inversa nello spazio dei giunti

Le ipotesi semplificative indicate nel regolatore a coppia precalcolata possono essere in parte rilassate
utilizzando lo schema a dinamica inversa (schemi 2 o equivalentemente 2bis figure H.42 e H.44) ove il
termine Vb è calcolato utilizzando posizioni e velocità misurate e non quelle desiderate.
In questo ca-
b b
so la condizione su V è meno stringente essendo sufficiente avere V Q, Q̇ ≃ V Q, Q̇ e non anche

b
∂ V /∂Q (Q − Q0 ) < ε1 e ∂ Vb /∂ Q̇ (Q̇ − Q̇0 ) < ε2 con ε1 e ε2 costanti piccole.
Disaccoppiamento dei giunti

c in posizione differente in modo che
Gli schemi precedenti possono essere migliorati spostando il blocco M
anche il termine del regolatore PD lo attraversi (vedi schemi 1b e 2b figure H.45 e H.46). In questo caso,
sotto le medesime ipotesi semplificative appena viste l’equazione differenziale diventa

M Q̈0 − Q̈ + M Kv Q̇0 − Q̇ + M Kp (Q0 − Q) = 0
Fq S S
F̂q + Fq
Mˆ Q
&& + Vˆ Q , Q& ( ) robot Mˆ Q
&& robot
Qo 0 0 0
+ Q0 0
+ Q + + Q
K p + sK v K p + sK v ( )
Vˆ Q, Q&
EQ EQ
+ - + -
Figura H.41: Coppia precalcolata nello spazio Figura H.42: Dinamica inversa nello spazio dei
dei giunti (schema 1). giunti (schema 2).
&&
Q 0
(
Vˆ Q0 , Q& 0 ) M̂
&&
Q + + EQ + Fq S
0 +
M̂ K p + sK v robot
EQ F̂q S Q0 + +
+ + -
K p + sK v robot Q
Q0
- + Fq Q
( )
Vˆ Q, Q&
Figura H.43: Controllo a coppia precalco-

lata nello spazio dei giunti (schema 1bis Figura H.44: Controllo a dinamica inver-
perfettamente equivalente a schema 1). sa nello spazio dei giunti (schema 2bis
perfettamente equivalente a schema 2).
&&
Q 0
(
Vˆ Q0 , Q& 0 )
&&
Q 0 EQ Fq S
+ + +
K p + sK v M̂ robot
EQ Fq S + +
+ + + - Q
K p + sK v M̂ robot
Q0 +
- +
Q
( )
Vˆ Q, Q&
Figura H.45: Controllo a coppia precalcolata Figura H.46: Controllo a dinamica inversa
nello spazio dei giunti a giunti disaccoppiati nello spazio dei giunti a giunti disaccoppiati
(schema 1b). (schema 2b).
S&&0 + S
Es + S&&' Fs' Fs Fq
R M̂ s JT robot
S0 + - + +
( )
Vˆs S , S&
Figura H.47: Controllo a dinamica inversa nello spazio operativo - teorico (schema 3a).
S&&0 + S
+ Es + +
&&'
Q Fq
R J-1 M̂ robot
- S&&' - + +
S0 Q
S Q& ( )
Vˆ Q, Q&
J&
cinematica diretta
Figura H.48: Controllo a dinamica inversa nello spazio operativo - pratico (schema 3b).
ricordando che la matrice delle masse è sempre definita positiva, e quindi invertibile, la relazione appena
scritta può essere semplificata moltiplicando a sinistra per M −1 ottenendo

Q̈0 − Q̈ + Kv Q̇0 − Q̇ + Kp (Q0 − Q) = 0
Questa relazione è di tipo lineare, è scomparsa la dipendenza dalla dinamica del manipolatore, e se
le matrici dei guadagni Kp e Kv sono diagonali, i vari giunti risultano tra loro disaccoppiati e quindi
facilmente regolabili agendo sulle matrici dei guadagni.
Indicando con kp e kv i guadagni di un asse, e ricordando che il denominatore di una funzione di trasfe-
rimento del secondo ordine è esprimibile come s2 + 2ξωn s + ωn2 , la pulsazione naturale e lo smorzamento
dell’asse valgono
p kv
ωn = kp ξ= p
2 kp
e sembrerebbe quindi possibile fissare arbitrariamente questi valori agendo sui guadagni. Il sistema
avrebbe funzione di trasferimento unitaria con errore di inseguimento tendente asintoticamente a zero
con transitorio del secondo ordine.
E’ questo un assunto teorico al quale ci si può solo avvicinare qualora si disponga di un modello veramente
buono infatti saranno necessariamente presenti inevitabili errori di misura e approssimazioni di modelliz-
zazione. Inoltre il modello considera gli attuatori ideali con tempo di risposta nullo e senza limitazioni
di potenza.
Dinamica inversa nello spazio operativo
Gli schemi di controllo indicati possono anche essere realizzati nello spazio operativo per un controllo
più diretto dell’errore d’inseguimento (vedi ad es. schemi 3a, 3b figure H.47 e H.48) ove il regolatore
può essere di tipo PI ponendo R = Kp + sKv . Per fare ciò è necessario esprimere tutte le grandezze
dinamiche in funzione delle coordinate utensile utilizzando le relazioni cinetostatiche S = F (Q), Ṡ = J Q̇,
S̈ = J Q̈ + J˙Q̇, Fq = J T Fs le quali permettono anche di trasformare la matrice delle masse
Ms = J −T M J −1
si noti a questo proposito l’analogia tra le relazione seguenti
Fq = M Q̈ ⇔ Fs = Ms S̈
Convertendo lo schema 2b (figura H.46 dinamica inversa nello spazio dei giunti) si perviene cosı̀ allo
schema 3a (H.47 dinamica inversa nello spazio di lavoro - schema teorico). Tuttavia la realizzazione
pratica di questo schema presenta diversi problemi:
1. normalmente non è disponibile una misura della posa dell’utensile ed è necessario stimarla a partire
dalle coordinate ai giunti;
2. normalmente le relazioni dinamiche (jacobiano, matrice delle masse, termini Vb , . . .) sono più fa-
cilmente esprimibili in funzione delle coordinate ai giunti Q piuttosto che in quelle nello spazio
operativo;
3. la sola conoscenza delle coordinate nello spazio operativo non dà un’informazione completa sullo
stato del sistema (ad esempio la stessa posizione S può essere raggiunta con diversi valori di Q).
Da un punto di vista pratico è quindi inevitabile dover modificare lo schema ottenendo lo schema “pratico”
3b (figura H.48); il blocco tratteggiato è quello che esegue l’operazione Q̈′ = J −1 (S̈ ′ − J˙Q̇).
Anche in questo caso, con regolatore PI (R = Kp +sKv ), si perviene ad un sistema di equazioni differenziali
lineari del secondo ordine in cui, se le matrici dei guadagni sono diagonali, le coordinate sono tra loro
disaccoppiate
S̈0 − S̈ + Kv Ṡ0 − Ṡ + Kp (S0 − S) = 0
x1 j
- x1
e1 u1
PID1 motore1
x0 +
x l x2
+ e2 u2
PID2 motore2
-
Figura H.50: Controllo decentralizzato a un

x2 grado di libertà.
Figura H.49: Sistema a 2 g.d.l. per la

movimentazione di un’asta.
x +
1
2
+
x0 - ex + u1
x1
x0 + x1 - e1 u1 x1
PIDx motore1
PID1 motore1 +
-
+ -
x2 j0 + u x2
j0 l + x2 + e2 u2 ej 2
PID2 motore2 PIDj motore2
2 + - + - +
-
1
l
j +
Figura H.51: Controllo decentralizzato a due
gradi di libertà.
Figura H.52: Controllo centralizzato a due

gradi di libertà (assi virtuali).
Assi virtuali Si consideri un sistema di movimentazione che deve movimentare un’asta. Possono veri-
ficarsi due casi. Nel primo caso l’asta deve semplicemente traslare e non ruotare (ϕ = 0), ma essendo
molto sottile è necessario utilizzare due motori che si devono muovere in maniera sincrona per evitare
che essa si infletta.
Alternativamente può essere necessario spostare e ruotare l’asta ottenendo quindi un sistema a due gradi
di libertà.
Nel primo caso è possibile utilizzare due motori retroazionati in posizione i quali ricevono lo stesso segnale
di posizione desiderata x0 . Si ottiene un sistema decentralizzato ad assi indipendenti.
Qualora si intenda anche regolare l’orientamento dell’asta un primo schema risolutivo prevede di calcolare
il set point di posizione di ogni motore xi a partire dal valore desiderato di posizione x0 ed orientamento
ϕ0

 l

 x1 = x0 − ϕ0
2

 x2 = x0 + l ϕ

2 0
I due motori sono poi controllati in maniera indipendente. Si ottiene cosı̀ uno schema di controllo
decentralizzato a 2 gradi di libertà.
Alternativamente si può decidere di controllare direttamente la posizione e l’orientamento dell’asta. Si
procede quindi al calcolo della sua posizione ed orientamento

 x1 + x2

 x=
2

 ϕ = 2 x1
 x −
l
questi valori vengono confrontati con quelli desiderati, gli errori di inseguimento vendono processati da
dei blocchi PID e le azioni dei due regolatori vengono ricombinate per essere inviate ai due motori. Gli
anelli di retroazione dei due gradi di libertà interagiscono tra di loro: siamo in presenza di un controllore
centralizzato a due gradi di libertà.
In pratica regoliamo due grandezze x e ϕ alle quali non corrisponde direttamente un attuatore; x e ϕ
sono talvolta detti assi virtuali.
Questo schema può anche essere considerato un controllo nello spazio operativo in cui S = [x ϕ]T
e Q = [x1 x2 ]T ; questo sarebbe particolarmente interessante qualora S fosse misurabile direttamente
permettendo cosı̀ di ottenere precisioni di controllo elevate.
Si osservi che nel caso decentralizzato i PID dei due assi saranno presumibilmente tarati in maniera
identica, mentre nel caso centralizzato i due PID devono controllare grandezze diverse con diverse inerzie
equivalenti e quindi andranno tarati in maniera differente.
Inoltre il caso centralizzato, a fronte di una maggiore complessità, permette di gestire meglio alcune si-
tuazioni. Immaginiamo per esempio che l’inclinazione debba essere garantita in maniera precisa mentre
la regolazione di posizione possa essere regolata in maniera più grossolana. In questo caso i parametri
del PID di posizione potranno avere guadagni inferiori rispetto a quelli di orientamento. Cosı̀, se ad
esempio, uno dei motori dovesse incontrare una coppia resistente molto alta ed imprevista, si otterrebbe
un rallentamento anche dell’altro motore mentre rimarrebbe sempre garantita la corretta inclinazione.
Viceversa se si fosse impiegato uno schema decentralizzato uno dei due motori avrebbe rallentato molto
più dell’altro producendo un’elevato errore di inclinazione.
La necessità di garantire una corretta inclinazione anche in presenza di disturbi asimmetrici potrebbe
giustificare l’uso del sistema centralizzato anche nel caso ad un grado di libertà ove si porrebbe sempre
ϕ0 = 0, cioè inclinazione desiderata nulla.
Figura H.53: Schema semplificato di robot SCARA con azionamento cilindrico dell’end-effector.
H.17 Giunto cilindrico attuato

In alcuni manipolatori è presente un giunto cilindrico motorizzato in entrambi i gradi di libertà. Ciò può
avvenire, per esempio, in un robot SCARA (fig. H.16) ove i due gradi di libertà servono per conferire il
moto verticale della pinza nonchè la rotazione attorno al proprio asse. L’attuazione può avvenire grazie a
due motori rotativi montati all’interno dell’ultimo link del robot. In alternativa i motori possono essere
posizionati nel piantone “a terra” ed il loro moto trasmesso ove necessario tramite due cinghie.
Il moto cilindrico viene generato tramite un dispositivo meccanico simile alle viti a ricircolo di sfere, ma
con due gradi di libertà (fig. H.54). Il dispositivo comprende un nocciolo con due scanalature di cui
una longitudinale (detta spline) ed una elicoidale (detta screw). Sono inoltre presenti due madreviti a
ricircolo di sfere ciascuna delle quali si accoppia con una delle due scanalature. Chiamiamo x e α le
coordinate che indicano il movimento traslatorio e rotatorio del nocciolo, e con x1 , α1 , x2 , α2 quello
delle due madreviti; p è il passo della vite. Per semplicità di rappresentazione, nelle figure H.54(a) e
H.54(b) le due scanalature e viti sono indicate separate, ma nella realtà esse sono contemporaneamente
presenti (H.54(c)).
Ciascuno dei due accoppiamenti produce un vincolo che lega il movimento del nocciolo a quello delle due
H.18. TRASMISSIONI A CAVI §12.3.2 PAG.332 1031
x1 x2 x2
α2 α2
α1
x1
p α1
x p
x
x
α α
α
α − α2 x − x2
(b) Ball-Screw = (c) Ball-Screw-Spline
(a) Ball-Spline α = α1 2π p
Figura H.54: Dispositivo con vite (screw) a spline a ricircolo di sfere per la movimentazione di
un’accoppiamento cilindrico (2 gradi di libertà).
madreviti: 

 α = α1
α − α2 x − x2 (H.23)

 =
2π p
che opportunamente riscritte portano a risolvere tanto la cinematica diretta che inversa:
 
 α = α1 
 α1 = α
x2 − x (H.24)
 x = p (α1 − α2 ) + x 
 α2 = α + 2π
2
2π p
Considerato che normalmente il movimento di traslazione delle due madreviti è bloccato x1 = x2 = 0 e
che esse possono solo ruotare, si ottiene
 
 α = α1 
 α1 = α
x (H.25)
 x = p (α1 − α2 ) 
 α2 = α − 2π
2π p
Quindi una pura rotazione si ottiene per α = α1 = α2 , e una pura traslazione si ottiene per α1 = 0 e
quindi x = −pα2 /(2π).
H.18 Trasmissioni a cavi §12.3.2 pag.332

La figura H.55 rappresenta l’uso di cavi per la trasmissione di moto tra pulegge ad assi paralleli o
ortogonali. Questo tipo di trasmissioni simula l’accoppiamento tra ingranaggi (fig. H.55(b)), ma permette
di eliminare giochi ed attriti. Per evitare slittamenti le estremità dei cavetti sono resi solidali alle pulegge.
Dato che i cavi possono trasmettere forze di trazione, ma non di compressione l’utilizzo di un solo cavo (fig.
H.55(a)) permette la trasmissione della coppia in un solo verso. Dato che le trasmissioni sono usualmente
bidirezionali occorre montare i cavetti in coppia (fig. H.55(c)); nella figura uno dei cavetti della coppia è
rappresentato tratteggiato. E’ anche possibile realizzare rapporti di trasmissione non unitari impiegando
ruote di diametro differente (fig. H.55(d)).
E’ possibile realizzare anche trasmissioni tra albero ad assi concorrenti anche tra più alberi coassiali (fig.
H.55(e)).
H.19 Polso e trasmissioni ad alberi cavi (pag. 339, figura 12.14)

Una foto della trasmissione del moto nel polso del robot DOGHI è riportata in figura H.56.
motrice
(a) trasmissione monodirezionale a cavi tra due pulegge (b) ingranaggi equivalenti
(c) trasmissione bidirezionale a cavi tra due pulegge

(d) trasmissione bidirezionale con τ 6= 1
(e) trasmissioni ad assi ortogonali e alberi coassiali
Figura H.55: Trasmissioni a cavi tra pulegge ad assi paralleli o concorrenti.

H.19. POLSO E TRASMISSIONI AD ALBERI CAVI (PAG. 339, FIGURA 12.14) 1033
Figura H.56: Il polso del robot DOGHI con trasmissione ad alberi cavi, ingranaggio conici a denti
elicoidali.
Accuratezza
Buona Scadente
Buona
Ripetibilità
Scadente
Figura H.57: Significato di ripetibilità (errori casuali) ed accuratezza (errori ripetitivi) di un manipolatore.
Gli errori ripetitivi possono essere compensati tramite calibrazione.
q1
q3
q2 φ3
φ1
φ2
Figura H.58: Calibrazione di una PKM: per migliorare le prestazioni, in aggiunta alle misure di errore di
posa, si possono utilizzare anche le misure dei movimenti dei giunti non motorizzati φi .
potenziometro o encoder
rocchetto
∆α anello calibrato ∆l
Figura H.59: Sistema di misura economico di
posa di un manipolatore tramite 6 rocchetti
Figura H.60: Rocchetto strumentato che
strumentati. La posa del manipolatore può
misura la quantità di filo estratto ∆l = r∆α.
essere stimata conoscendo la lunghezza dei sei
fili e la posizione dei rocchetti.
3 3
4
2
1
4
5
2
5
1
(b) struttura cinematica
(a) Braccio di misura
Figura H.61: Braccio seriale di misura a 5 gradi di libertà.
α, β, d ⇒ x, y, z
α
Figura H.62: Dima per forzare il manipolatore

ad una posa nota; qui è utilizzato per forzare
la posa del braccio di figura H.61(a). Figura H.63: Il laser tracker consente di
misurare la posizione assoluta del centro
di uno speciale specchio (corner cube).
H.20. CALIBRAZIONE(§14.1 PAG. 403 E SEGUENTI) 1035
Pa
Pr
gradiente termico
Figura H.64: A causa di un gradiente termico il raggio laser viene deviato. Il laser tracker crede che lo
specchio si trovi in Pa (posizione apparente) mentre in realtà esso si trova in Pr (posizione reale).
H.20 Calibrazione(§14.1 Pag. 403 e seguenti)

H.20.1 Concetti generali
La figura H.57 richiama i concetti di ripetibilità ed accuratezza di un manipolatore (vedi anche capitoli
relativi alle misura delle prestazioni §17).
La calibrazione di un manipolatore richiede una strumentazione in grado di misurare la posa reale rag-
giunta da un manipolatore per calcolare poi l’errore per differenza tra posa reale e teorica. La misura
assoluta di posa è molto complessa e richiede strumentazione complessa e costosa. A puro esempio in-
dicativo in figura H.59 viene schematizzato un sistema di misura economico di posa di un manipolatore
tramite 6 rocchetti strumentati (figura H.60). La posa del manipolatore può essere stimata conoscendo
la posizione dei rocchetti, la lunghezza dei sei fili e utilizzando equazioni simili a quelle per la cinematica
diretta di manipolatori paralleli (vedi §8.15).
Quando non è conveniente misurare la posa effettivamente raggiunta da un manipolatore si può usare un
metodo “duale” che prevede di forzare la posa a valori predefiniti (esempio figura H.62). Poi si stima la
posa tramite la cinematica diretta utilizzando come ”movimento motori“ Q la lettura dei trasduttori dei
giunti. La differenza tra posa imposta e posa stimata serve ad innescare il procedimento di calibrazione.
La figura H.63 rappresenta un laser tracker ovvero una strumentazione laser per la misura della posi-
zione di uno speciale specchio di tipo corner cube da porsi sul dispositivo d’estremità del robot. La
sorgente laser può venire diretta in ogni direzione ruotando due specchi (interni allo strumento). Se il
raggio laser incontra lo speciale specchio esso viene riflesso verso la sorgente. Opportuni sensori sono in
grado di misurare la distanza d tra strumento e specchio. Per la misura di distanze esistono metodolo-
gie interferometriche di tipo incrementali molto precise che consentono di arrivare a precisioni di alcuni
micron, mentre metologie di tipo assoluto arrivano a frazioni di decimo di millimetro. La conoscenza
delle rotazioni α e beta degli specchi e della distanza d permette di determinare le coordinate x, y, z dello
specchio. Appositi sensori misurano la temperatura ambiente per compensare gli effetti sulla strumen-
tazione, tuttavia se la temperatura ambientale varia lungo il percorso del raggio laser, esso non viaggia
più in linea retta e la strumentazione viene ingannata (figura H.64); occorre quindi effettuare le misure
in un ambiente adeguato. Per misurare l’orientamento del dispositivo di estremità è necessario misurare
la posizione di almeno tre diversi punti (specchi) posti sul dispositivo d’estremità.
La figura H.58 si riferisce a una possibile tecnica calibrazione di un manipolatore a struttura cinematica
parallela (PKM Parallel Kinematic Machine).
H.20.2 Calibrazione non parametrica (§14.1.5 pag.407)

La calibrazione parametrica prevede che l’errore di posa di un manipolatore venga misurato in punti
appartenenti ad una griglia regolare (es. figura H.66). Normalmente l’errore varierà gradualmente da
una posa all’altra e questo consentirà di prevedere l’errore per ogni posa intermedia tramite interpolazione
dai valori adiacenti nella griglia.
Conoscendo l’errore di posa sarà possibile compensarlo producendo uno spostamento uguale e contrario.
Tuttavia questo metodo non consente di capire quale sia la fonte dell’errore e quindi la calibrazione può
essere applicata solo nei punti interni alla griglia mentre non si è in grado di calibrare il manipolatore
quando si trova fuori da essa. Inoltre se il root deve essere calibrato per orientamento qualsiasi della pinza,
per ogni punto della griglia è necessario misurare l’errore di posa per parecchi orientamenti del dispositivo
d’estremità. Questa operazione può essere molto complessa e richiedere parecchi tempo tuttavia si è
sgravati dall’onere di sviluppare un modello sofisticato del manipolatore.
2
R
6
4
Figura H.65: Interpretazione geometrica della scelta dei parametri di calibrazione: un modello completo,
minimo, proporzionale.
H.21 Calibrazione: scelta dei parametri: un modello completo,

minimo, proporzionale (§14.1.14 pag.416)
La natura dei parametri strutturali considerata in questo approccio è descritta dalla figura H.65.
H.22 Modello dinamico con elasticità concentrate (pag. 425

§14.2.4, §2.9 pag. 56 e §2.8.1 pag. 49)
Le equazioni del modello dinamico con elasticità concentrata si possono ottenere sostituendo nell’equa-
zione della dinamica Fq = MQ̈ + V(Q, Q̇) + G1 (Q, F ) il valore delle forze trasmesse dalla elasticità delle
trasmissioni. Risulta

Kq (Qm − Q) + Kr Q̇m − Q̇ = MQ̈ + V(Q, Q̇) + G1 (Q, F )
ove Qm rappresenta le rotazioni dei motori “equivalenti” (uscita riduttori supposti rigidi) e Q le rotazioni
dei link; le matrici diagonali Kq e Kr indicano la rigidezza e le perdite di deformazione viscosa delle
trasmissioni
Kq = diag (k1 , k2 , . . .) Kq = diag (r1 , r2 , . . .)
nell’esempio del manipolatore SCARA si ottiene

α1 β1 k1 0 r1 0
Qm = Q= Kq = Kr =
α2 β2 0 k2 0 r2
Definiamo con ∆Q la differenza tra movimento dei motori e dei link e con ∆S la differenza tra posizione
teorica dell’organo terminale nel caso teorico rigido e quella del caso reale con elasticità
∆Q = (Qm − Q) ∆S = (Sm − S)
Si può operare un cambio di coordinate dallo spazio dei giunti in quello di lavoro ponendo ∆S = J∆Q
ottenendo
Ms ∆S̈ + Krs ∆Ṡ + Ks ∆S = J −T MQ̈ + V(Q, Q̇) + G1 (Q, F )
con
Ms = J −T MJ −1 Krs = J −T Kr J −1 Ks = J −T Kq J −1
Le equazioni scritte sono quelle tipiche di sistemi vibranti forzati. In questo caso la forzante è il movi-
mento dei motori e l’integrazione della equazione consente di determinare il movimento del manipolatore
reale. Nel caso che i termini quadratici in Q̇ siano trascurabili e nell’ipotesi di piccole deformazioni si
ottiene un’equazione differenziale lineare; i modi principali di vibrare si possono dedurre dagli autovettori
−1
di M−1s Ks o di M Kq mentre le pulsazioni proprie si deducono dai corrispondenti autovalori.
H.23. STIMA PARAMETRI DEL MODELLO DINAMICO (PAG. 427 §14.2.7, §2.9 PAG. 56, §2.21 PAG. 70 E PAG. 50 §2.8.3)10
Figura H.66: Griglia di calibrazione (non parametrica). L’errore di posa ∆S viene misurato su una griglia
regolare di punti. Quando è richiesto di raggiungere una posa S, in realtà il manipolatore raggiungerà la
posa S +∆S. L’errore di posa può venire annullato chiedendo al robot di raggiungere la posa S −∆S. Per
i punti per i quali non è disponibile un valore di errore di posa questo valore si stima per interpolazione
dai punti adiacenti.
H.23 Stima Parametri del modello dinamico (pag. 427 §14.2.7,

§2.9 pag. 56, §2.21 pag. 70 e pag. 50 §2.8.3)
Nota La funzione segno dovrebbe essere indicata come sgn() ma altre volte compare erroneamente come
sign(); essa è definita al §B.3.1 come sgn(x) = 1 se x > 0, 0 se x = 0, −1 se x < 0.
Consideriamo il modello a parametri dinamici minimi

T
C = B Pe Pe = PdT PaT
La simbologia usata nei due capitoli non è omogenea. Inoltre nel paragrafo sulla calibrazione dinamica
non sono stati esplicitati la massa m ed il momento d’inerzia J del corpo trasportato mentre sono stati
introdotti i momenti d’inerzia dei motori e i rapprti di riduzione τ1 e τ2 . Le varie equazioni non sono cosı̀
direttamente confrontabili. Ciò non influenza la validità dei risultati sperimentali descritti perché, come
dichiarato, non ci interessa (ne è completamente possibile) trovare il valore dei singoli parametri, ma solo
quello di parametri globali equivalenti in un modello rigido a parametri dinamici minimi.
Inoltre per omogeneità conviene eliminare il coefficiente l1 dalla matrice B ed inglobarlo nel vettore dei
termini dinamici.
E’ quindi possibile rifinire il vettore dei parametri P e la matrice B nella seguente maniera:
 
    a01
p1 I1z + (m + m2 )l12 + Jm1 /τ12  a11 
 p2   l1 (ml2 + m2 sx )   
     a21 
 
Pd =  p3  =   l 1 m2 s y  Pa =    (H.26)
 a02 
 p4   I2z   
 a12 
p5 Jm2 /τ22
a22
 
τ1 q̈1 0
 τ1 (c2 (2q̈1 + q̈2 ) − s2 q̇2 (2q̇1 + q̇2 )) τ2 c2 q̈1 + s2 q̇12 
 
 τ1 (−s2 (2q̈1 + q̈2 ) − c2 q̇2 (2q̇1 + q̇2 )) τ2 −s2 q̈1 + c2 q̇12 
 
 τ1 (q̈1 + q̈2 ) τ2 (q̈1 + q̈2 ) 
 
 0 τ2 q̈2 
 
T
B =  sgn (q̇1 ) 0 

 q̇1 0 
 
 sgn (q̇1 )q̇12
0 
 
 0 sgn(q̇2 ) 
 
 0 q̇2 
2
0 sgn(q̇2 )q2
I tre blocchi evidenziati nella matrice B T evidenziano la parte relativa alle masse (azioni d’inerzia) e
quelle relative alle perdite nei due riduttori.
Per compare le equazioni del §14 con quelle del §2 occorre osservare che quest’ultime non tengono conto
delle trasmissioni per cui è necessario considerare le seguenti eguaglianze
Jm1 = 0 Jm2 = 0 τ1 = 1 τ2 = 1
infine la corrispondenza tra simboligie si completa osservando che vale
I1z = Jg1 + m1 l12 I2z = Jg2 + J + m2 g22 + ml22 sx = g 2 sy = 0
sy è nullo perché si assume che il baricentro stia sulla mezzeria del link.
Avendo spostato il coefficiente l1 = 0.33 dalla matrice B al vettore P occorre moltiplicare i parametri
p2 e p3 per questo stesso valore ed aggiornare la sua unità di misura. Nell’esempio sperimentale citato il
parametro p3 vale esattamente zero perché la stima è stata fatta con riferimento ad un modello ridotto
che assume sy = 0.
Risulta
cerchio OJD
P1 16.3286514074 16.2606169659 kgm2
P2 -0.0550159053 0.837423654 kgm2
P3 0.0 0.0 kgm2
P4 -0.3813491029 0.6642731627 kgm2
P5 4.8753219665 4.2045177593 kgm2
a01 0.1846270328 0.2100604969 Nm
a11 0.2945835590 0.0940617391 N m/(rad/s)
a21 -0.1448462244 -0.0168580957 N m/(rad/s)2
a02 0.0841700456 0.0876724328 Nm
a12 0.0251936773 0.0508717965 N m/(rad/s)
a22 0.0009255059 -0.0073197797 N m/(rad/s)2
Si osservi che il parametro p5 non è nullo evidenziando l’importanza del riduttore sull’inerzia vista dal
motore. A questo riguardo si richiamano i discorsi fatti al §12.6. Per comparare i valori numerici del
§12.6 con quelli del §14.2 si osservi che in un caso l’inerzia è valutata prima del riduttore, nell’altra dopo.
Inoltre si ricorda che nel caso in esame τ = 1/100 da cui τ 2 = 10−4 .
Appendice I
Ellissi: proprietà e rappresentazioni
In questo paragrafo trattiamo alcuni argomenti relativi ad ellissi in un piano xy, ma tutti i risultati
generali basati sui concetti di autovalori ed autovettori sono validi anche per ellissoidi nello spazio xyz
(o a più dimensioni), anche se alcune formule particolari vanno ovviamente adattate.
Un’ellisse con diametri paralleli agli assi coordinati x, y e con centro nell’origine può essere espressa in
forma canonica come
x2 y 2
+ =1
r12 r22
ove r1 e r2 sono i raggi (semidiametri) dell’ellisse stessa. La stessa ellisse può essere espressa in forma
matriciale come
   
1 1
0  0 
 r2 x x  r2
v T Av = x y 
1  =1 v= A= 1 
1  y y  1 
0 0
r22 r22
ed in forma parametrica come

x r1 cos(α)
=
y r2 sin(α)
L’area dell’ellisse è
a = πr1 r2
Ruotando l’ellisse di un angolo α0 , l’area e la lunghezza dei semidiametri non cambiano. Le equazioni che
rappresentano l’ellisse ruotata possono essere ricavate dalla relazione che esprime le coordinate ruotate
x′ , y ′ in funzione di quelle originali x, y
′
′ x cos(α0 ) − sin(α0 ) x
v = Rv =
y′ sin(α0 ) cos(α0 ) y
y
y’
u2
r2 a x u1
r1
r1 a0
r2
x’
Figura I.1: Costruzione di un’ellisse ispi-

Figura I.2: Ellisse ruotata di un angolo α0 .
rata alla sua definizione tramite equazione
parametrica.
1039
1040 APPENDICE I. ELLISSI: PROPRIETÀ E RAPPRESENTAZIONI
v t Bv = 1 v t B-1v = 1
Figura I.3: Ellissi associate ad una matrice B ed alla sua inversa B −1 .

x′ cos(α0 ) − sin(α0 ) r1 cos(α)
=
y′ sin(α0 ) cos(α0 ) r2 sin(α)
ove le colonne della matrice ortogonale “di rotazione” R possono essere interpretate come i versori degli
assi originali x, y (su cui giacciono i diametri) visti da x′ , y ′
 

cos(α0 ) − sin(α0 )
R= =  u1 u2  RT = R−1
sin(α0 ) cos(α0 )
ricordando la relazione inversa v = RT v ′ , l’ellisse può essere riscritta come

 
  1
0
′ ′  r2  uT1 x′
T
1 = v Av = v ′T T ′T ′ ′
RAR v = v A v = x y ′  u1 u2    1 
1  uT2 y′
0
r22
con  
 
1
 r2 0 
uT1
′
A = RAR T
=  u1 u2  

1 
1  uT2
0
r22
ove ui sono gli autovettori di A′ ed i corrispondenti autovalori sono λi = 1/ri2 .
Si deduce quindi che data un’ellisse (o un’ellissoide) nella forma v T A′′ v = 1, gli autovettori ui della matrice
A′′ rappresentano le direzioni dei suoi assi, mentre l’inverso della radice quadrata dei corrispondenti
1
autovalori λi rappresentano i corrispondenti semidiametri ri = √ .
λi
Inoltre dato che il determinate di una matrice è pari al prodotto degli autovalori, si deduce che l’area
dell’ellisse è proporzionale all’inverso della radice quadrata del suo determinante
π π
a = πr1 r2 = √ = p
λ1 λ2 det(A′′ )
Talvolta la matrice A′′ (che risulta simmetrica e definita positiva) è ottenuta dal prodotto di una matrice
B per la sua trasposta: A′′ = B T B. In questo caso si definiscono valori singolari σi di B la radice
quadrata degli autovalori di A′′
p
σi = λi det(A′′ ) = det(B)2
per cui si può affermare che i semidiametri sono l’inverso dei valori singolari di B, l’area dell’ellisse è
inversamente proporzionale al determinante di B. La matrice A′′ è chiaramente definita positiva. Si
ottiene
1 π π
σi = a = πr1 r2 = =
ri σ1 σ2 det(B)
Comparando una matrice B con la sua inversa B −1 si osserva che i loro determinanti (come pure i loro
autovalori) sono uno l’inverso dell’altro mentre gli autovettori coincidono
1 1
det (B) = λi (B) = ui (B) = ui B −1
det (B −1 ) λi (B −1 )
1041
quindi l’asse maggiore e quello minore delle ellissi si scambiano di posto ed ad un’ellisse di area “piccola”
(minore di 1) corrisponde una di area “grande” (maggiore di 1) o viceversa.
Nel caso di ellissoidi nello spazio xyz in luogo di area occorre parlare di volume calcolabile dal prodotto
dei tre semidiametri Vol = 4/3πr1 r2 r3 .
In generale si può verificare che un’ellisse con generico orientamento e centro nell’origine degli assi può
essere scritta come
a b/2 x
ax2 + bxy + cy 2 = [x y] = v T Av = 1
b/2 c y
con a, c > 0 e det A > 0.
Dimetro in una direzione.

Consideriamo un ellisse associata alla matrice A. Il diametro in una direzione u si può facilmente deter-
minare rappresentandolo con un vettore v di modulo d moltiplicato per il versore associato alla direzione
(v = ud), l’ellisse è cosı̀ definita come
dT uT Aud = 1
da cui
1
d= √
T
u Au
1042 APPENDICE I. ELLISSI: PROPRIETÀ E RAPPRESENTAZIONI
Appendice J
Appendice richiami matematici -

proprietà matrici
J.1 rango
Alcune proprietà del rango di una matrice generica A
rank(A) = rank(AT ) = rank(AAT ) = rank(AT A)
J.2 proiezione
proiezione x′ di un vettore x sullo spazio definito dalle colonne di A
−1 T
x′ = P x con P = A AT A A

valida se AT A risulta invertibile.
Il residuo xr (cioè la parte ortogonale allo spazio) risulta
−1 T
xr = x − P x = Bx con B = I − A AT A A
ove I è la matrice identità; risulta ovviamente x = x′ + xr . B e P sono matrici proiezione; tutte le matrici
proiezione sono singolari tranne la matrice identità.
Per dimostrare la relazione per il calcolo della matrice P , si possono considerare le proprietà della matrice
−1 T
pseudo inversa (vedi appendice C) A+ = AT A A
J.3 derivate di funzioni lineari e quadratiche

Siano c e x delle matrici colonna ed A una matrice quadrata. Valgono le seguenti proprietà
∂cT x ∂xT c ∂Ax ∂xT Ax

= =c = AT = A + AT
∂x ∂x ∂x ∂x
1043
1044 APPENDICE J. APPENDICE RICHIAMI MATEMATICI - PROPRIETÀ MATRICI
Appendice K
Pentalatero (§8, pag.213) -

Manipolatore parallelo a 2 gradi di
libertà
Esistono più versioni del meccanismo denominato pentalatero. In questa sede si esamina la versione
classica azionata da motori rotativi (fig. K.1, K.4) e quella azionata da motori lineari (fig. K.5). Si
tratta di manipolatori piani a due gradi di libertà destinati a posizionare un utensile nel piano xy (pol-
so del robot) tramite due braccia motorizzate. A scopo propedeutico è preliminarmente studiato l’arco
a tre cerniere ed un manipolatore che si può ottenere facendo variare la lunghezza delle sue aste (fig. K.2).
Una realizzazione commerciale del pentalatero classico è raffigurata in fig. K.1.
K.1 L’arco a 3 cerniere

In figura K.2 è rappresentato un arco a tre cerniere. Indichiamo con xa , ya e xb , yb le coordinate delle
cerniere a terra e con x, y quelle della cerniera mobile; le lunghezze delle due aste sono ρa e ρb .
Questa struttura può essere utilizzata per generare meccanismi a uno o due gradi di libertà. Ciò può
avvenire facendo variare la lunghezza di una, o entrambe le aste, oppure movimentando una o entrambe
cerniere a terra. Al §K.2 viene descritto come realizzare un pentalatero classico facendo muovere le
cerniere su traiettorie circolari, mentre al §K.3 viene descritto il caso di cerniere mobili su guide lineari.
In questa sezione studiamo invece il caso in cui la cerniera mobile venga azionata facendo variare le
lunghezze delle due aste tramite azionamenti lineari (o pistoni). In questo caso le coordinate nello spazio
dei giunti Q e nello spazio di lavoro S saranno

ρa x
Q= S= (K.1)
ρb y
Figura K.1: Un pentalatero commerciale RP Mitsubishi a 4 gradi di libertà; apposite trasmissioni nascoste
all’interno dei bracci permettono la rotazione del polso attorno al suo asse ed una movimentazione in
direzione verticale.
1045
1046 APPENDICE K. PENTALATERO (§8, PAG.213) - MANIPOLATORE PARALLELO A 2 GRADI DI LIBERTÀ
xa ya
ρa xb yb
ρb
xy
Figura K.2: L’arco a tre cerniere: la base

cinematica del pentalatero. Figura K.3: Configurazioni singolari dell’arco
a 3 cerniere.
xa ya
la xb yb xa ya
la xb yb
l1 l2
lb
α1 xy ρ1 lb
α2 x1 y1 xy ρ2
x1 y1
α1 α2
x2 y2
x2 y2
Figura K.4: Pentalatero tradizionale ad Figura K.5: Pentalatero ad azionamenti

azionamenti angolari. lineari.
La posizione dell’utensile si trova nell’intersezione di due circonferenze di centro xa ya e xb yb e raggio,

rispettivamente, ρa e ρb (fig. K.2)
 2
 ρa = (x − xa )2 + (y − ya )2
(K.2)
 2
ρb = (x − xb )2 + (x − xb )2
É evidente che potranno esistere due possibili configurazioni corrispondenti alle due intersezioni. Questo
procedimento risolve il problema cinematico diretto. Le due soluzioni potranno essere coincidenti qualora
le due cerniere a terra distino esattamente quanto il valore d pari alla somma delle lunghezze due aste
d = ρa + ρb o alla loro differenza d = |ρa − ρb | (fig. K.3); in questo caso avremo una configurazione
singolare per la cinematica diretta (cerchi tangenti). In questo caso i due motori non saranno in grado
di controllare il movimento dell’utensile in direzione normale a quella delle due aste allineate; una forza
applicata all’utensile in direzione ortogonale alle aste provoca in esse forze infinite.
La cinematica inversa è invece risolta da una semplice trasformazione delle eq.(K.2)
 p
 ρa = ± (x − xa )2 + (y − ya )2
p (K.3)

ρb = ± (x − xb )2 + (x − xb )2
Esistono 4 soluzioni (combinazioni di lunghezze positive o negative); la configurazioni singolari sono quel-
le con ρa = 0, oppure con ρb = 0; infatti in questi casi 2 soluzioni coincidono. In queste configurazioni
singolari l’asta a lunghezza nulla può liberamente ruotare su se stessa senza vincoli.
Questo manipolatore è uno di quelli in cui è più semplice risolvere la cinematica diretta, che quella inversa.
K.2 Pentalatero classico

In questo pentalatero i due motori rotativi azionano le due aste incernierate nei punti x1 y1 e x2 y2 (spalle)
la cui posizione angolare è indicata dagli angoli α1 e α2 . Le coordinate ai giunti Q e nello spazio di lavoro
K.2. PENTALATERO CLASSICO 1047
S sono
α1 x
Q= S= (K.4)
α2 y
La posizione delle due cerniere intermedie (gomiti) è indicata dalle coordinate xa ya e xb yb che si deter-
minano facilmente come
xa = x1 + l1 cos(α1 ) xb = x2 + l2 cos(α2 )
(K.5)
ya = y1 + l1 sin(α1 ) yb = y2 + l2 sin(α2 )
La posizione dell’utensile si trova nell’intersezione di due circonferenze di centro xa ya e xb yb e raggio,

rispettivamente, la e lb (fig. K.6)
 2
 la = (x − xa )2 + (y − ya )2
(K.6)

lb2 = (x − xb )2 + (x − xb )2
E’ evidente che si troveranno due soluzioni corrispondenti alle due intersezioni. Questo procedimento
risolve il problema cinematico diretto. Le due soluzioni potranno essere coincidenti qualora i due gomiti
distino esattamente quanto il valore d pari alla somma delle due aste d = la + lb o alla loro differenza
d = |la −lb | (fig. K.8); in questo caso avremo una configurazione singolare per la cinematica diretta (cerchi
tangenti). In questo caso i due motori non saranno in grado di controllare il movimento dell’utensile in
direzione normale a quella delle due aste allineate; una forza applicata all’utensile in direzione ortogonale
alle aste provoca in esse forze infinite. Ovviamente non è possibile azionare il pentalatero in modo da
allontanare maggiormente i due gomiti (cerchi non intersecanti).
Per passare da una delle due configurazioni all’altra è necessario che il polso attraversi la configurazione
singolare, operazione difficile da compiere per le difficoltà di controllo.
Per risolvere il problema cinematico inverso si osserva che una volta nota la posizione xy dell’utensile,
quella del gomito xa ya è determinata dall’intersezione di due cerchi. Il primo cerchio ha centro xy e
raggio la , mentre il secondo ha centro x1 y1 e raggio l1 (fig. K.7); l’angolo α1 si determina quindi con un
arcontangente
 2
 la = (x − xa )2 + (y − ya )2
α1 = atan2(ya − y1 , xa − x1 ) (K.7)
 2
l1 = (x1 − xa )2 + (x1 − xa )2
Analogamente si ragiona per l’altro braccio determinando α2 . In generale, per ogni braccio, si avranno
due soluzioni corrispondenti alle due intersezioni dei cerchi. Le due soluzioni potranno coincidere qualora
le due aste di un braccio si allineino (cerchi tangenti); questa situazione corrisponde ad una configurazione
singolare per la cinematica inversa (fig. K.9). In questa configurazione un movimento infinitesimo di un
motore non determina cambiamenti della posizione dell’utensile che quindi localmente perde un grado
di libertà. La mancata intersezione dei cerchi indica che il punto xy è fuori dallo spazio di lavoro del
manipolatore.
Si hanno due possibili soluzioni per ciascun braccio e quindi, in totale, 4 soluzioni per il manipolatore.
Lo studio delle velocità ed accelerazioni può essere fatta risolvendo la cinematica diretta e/o inversa
S = F (Q), Q = G(S) ed effettuando la derivazione rispetto al tempo previo calcolo dello jacobiano
J = ∂F/∂Q, J −1 = ∂G/∂S. In alternativa si può utilizzare la metodologia descritta al §8.8 pag.218
basata su una diversa forma dell’equazione cinematica descritta dall’eq. (8.2) che nel nostro caso si
ottiene combinando le equazioni K.5 e K.6
   
(x − (x1 + l1 cos(α1 )))2 + (y − (y1 + l1 sin(α1 )))2 − la2 0
F (S, Q) =  =  (K.8)
2 2 2
(x − (x2 + l2 cos(α2 ))) + (x − (y2 + l2 sin(α2 ))) − lb 0
per cui
   
x − (x1 + l1 cos(α1 )) y − (y1 + l1 sin(α1 )) x − xa y − ya
∂F
Js = = 2  = 2  (K.9)
∂S
x − (x2 + l2 cos(α2 )) x − (y2 + l2 sin(α2 )) x − xb x − xb
 
A 0
∂F
Jq = − = −2   (K.10)
∂q
0 B
con
A = ((x − x1 − l1 cos(α1 )) sin(α1 ) − (y − y1 − l1 sin(α1 )) cos(α1 )) la
= ((x − xa ) sin(α1 ) − (y − ya ) cos(α1 )) la

(K.11)
B = ((x − x2 − l2 cos(α2 )) sin(α2 ) − (x − y2 − l2 sin(α2 )) sin(α2 )) lb
= ((x − xb ) sin(α2 ) − (x − yb ) cos(α2 )) lb

È abbastanza facile dimostrare che il determinante dello jacobiano Js si annulla quando sono allineati i
tre punti xa ya , xy, xb yb (configurazione singolare per la cinematica diretta); invece il determinante di Jq
si annulla quando si allineano i punti x1 y1 , xa ya , xy oppure i punti x2 y2 , xb yb , xy (configurazioni singolari
per la cinematica inversa).
K.3 Pentalatero con azionamenti lineari

In questa versione del pentalatero, la posizione dei gomiti è determinata dalla movimentazione degli
azionamenti lineari. Le coordinate ai giunti Q e nello spazio di lavoro S sono

ρ1 x
Q= S= (K.12)
ρ2 y
La posizione delle due cerniere intermedie (gomiti) è indicata dalle coordinate xa ya e xb yb che si deter-
minano facilmente come
xa = x1 + ρ1 cos(α1 ) xb = x2 + ρ2 cos(α2 )
(K.13)
ya = y1 + ρ1 sin(α1 ) yb = y2 + ρ2 sin(α2 )
ove x1 , y1 sono le coordinate di un punto arbitrario, assunto come origine del movimento, che si trova
sulla retta di inclinazione α1 che contiene il movimento della cerniera di coordinate xa ya ; analogamente
x2 , y2 è l’origine della retta di inclinazione α2 su cui scorre xb yb .
La posizione dell’utensile, si trova dall’intersezione di due circonferenze di centro xa ya e xb yb e raggio,
rispettivamente, la e lb come già visto per il pentalatero tradizionale (eq. K.6), risolvendo cosı̀ il problema
cinematico diretto. Anche in questo caso si possono avere due soluzioni, oppure una (configurazione
singolare), oppure nessuna.
La soluzione della cinematica inversa invece si ottiene per ciascun braccio come intersezione tra le rette
degli azionamenti lineari (eq. K.13) e i due cerchi di centro xy, raggio pari a la (oppure lb per l’altro
braccio) descritti dalle eq. K.6. Si ottengono due equazioni di secondo grado in ρ1 e ρ2 . Si hanno due
possibili soluzioni per ciascun braccio e quindi, in totale, 4 soluzioni per il manipolatore.
Le configurazioni singolari per la cinematica inversa accadono quando uno (o entrambi) i bracci sono
ortogonali alle rispettive guide lineari e le equazioni di 2◦ grado hanno soluzioni reali coincidenti (retta
tangente al cerchio). Ovviamente non si hanno soluzioni reali se il polso è troppo lontano dalle guide.
Lo studio delle velocità, accelerazioni e configurazioni singolari può essere svolto anche con l’impiego degli
jacobiani Js e Jq analogamente a come visto per il pentalatero tradizionale.
K.4 Spazio di lavoro

Lo spazio di lavoro del pentalatero tradizionale si determina osservando che ciascuno dei due bracci del
pentalatero tradizionale ha la stessa struttura cinematica del robot scara. Pertanto, considerando il pri-
mo braccio, si deduce che il polso è costretto a stare nella corona circolare di centro x1 y1 , raggio minore
r1 = l1 + la , e raggio maggiore R1 = |l1 − la |. Analogamente, considerando il secondo braccio, il polso
è costretto a stare nella corona circolare di centro x2 y2 , raggio minore r2 = l2 + lb , e raggio maggiore
R2 = |l2 − lb |. Secondo le dimensioni delle aste e le posizioni delle cerniere a terra si possono ottenere
diverse forme dello spazio di lavoro. Lo spazio di lavoro può essere unico o diviso in due zone. In questo
caso se il polso si trova in una delle due zone non può essere spostato nell’altra a meno di smontare e
rimontare il manipolatore (fig. K.10).
I bordi dell’area di lavoro corrispondono alla posizione del polso quando si trova nelle configurazioni
singolari per la cinematica inversa. Il luogo delle posizioni singolari per la cinematica diretta è di forma
complessa ed è una curva di 4◦ grado (verificare).
Lo spazio di lavoro si riduce se l’ampiezza del movimento dei motori è limitato da fine corsa.
Lo spazio di lavoro del pentalatero ad azionamenti lineari è un parallelogramma che si determina osser-
vando che il polso deve stare in una fascia di semiampiezza la rispetto alla retta 1 su cui si muove il gomito
K.4. SPAZIO DI LAVORO 1049
Figura K.7: Costruzione geometrica per la ci-

Figura K.6: Costruzione geometrica per la cinematica inversa del pentalatero tradizionale.
nematica diretta del pentalatero tradizionale.
Figura K.8: Configurazioni singolari per la ci-

nematica diretta del pentalatero tradizionale;
nei 2 casi di sinistra e centro sono possibili
Figura K.9: Configurazioni singolari per
movimenti infinitesimi, nel caso di destra sono
la cinematica inversa del pentalatero ad
possibili movimenti finiti.
azionamenti lineari.
a e in una fascia di semiampiezza lb rispetto alla retta 2 su cui scorre il gomito b (fig. K.11 a sinistra).
I bordi del parallelogramma rappresentano il luogo dei punti (posizione del polso) nelle configurazioni
di singolarità per la cinematica inversa; invece si dimostra che il luogo dei punti di singolarità per la
cinematica diretta sone due ellissi (fig. K.11 a destra).
Figura K.10: Esempi di area di lavoro del

pentalatero tradizionale.
la
2
lb
Figura K.11: Area di lavoro del pentalatero ad azionamenti lineari.

Jose Luis Huayanay - Robotica

Caricato da

Copyright:

Formati disponibili

Jose Luis Huayanay - Robotica

Caricato da

Informazioni sul documento

Descrizione originale:

Titolo originale

Copyright

Formati disponibili

Condividi questo documento

Condividi o incorpora il documento

Opzioni di condivisione

Hai trovato utile questo documento?

Questo contenuto è inappropriato?

Copyright:

Formati disponibili

Jose Luis Huayanay - Robotica

Caricato da

Copyright:

Formati disponibili

1

Errata Corrige ed Ampliamenti

pagina riferimento errata corrige

66 eq. (2.47) de- −2l1 z − (x2 + . . . −2l1 z + (x2 + . . .

66 eq. (2.50) e poco u2 + y 2 u2 + v 2

pagina riferimento errata corrige

79 sotto matrice R XrT Zr = . . . = 1 XrT Xr = . . . = 1

96 Rtras(z . . .) − sin(ϑy ) − sin(ϑz )

148 M03 =  0  M03 =  −S2 

175 fine pagina M0a = . . . =   M0a = . . . =  

pagina riferimento errata corrige

425 eq.(14.15) m2 g2 β̇12 c −m2 g2 β̇12 c

pagina riferimento errata corrige

[1] A homogeneous matrix approach to 3D kinematics and dynamics - I. Theory

H.1 Ellissi di Manipolabilità, Isotropia (Pagg.45-48)

H.1.1 Ellissi di manipolabilità in velocità e in forza

Figura H.1: Manipolatore cartesiano: velocità

Figura H.4: Composizione del-

Attenzione: l’ellisse di forza e quella di velocità sono tra loro or-

H.1.2 Ellissi d’inerzia (di massa)

H.1.3 Forza, Velocità massima in una direzione.

in maniera analoga è calcolabile la velocità massima.

H.1.5 Esempi applicativi

Figura H.7: Robot SCARA in configu-

Figura H.9: Le ellissi di velocità del robot pola-

y Ellissi d’inerzia (massa): robot cartesiano

(a) β > 0 (b) β < 0

Figura H.13: Le ellissi di inerzia del robot ICOMATIC SCARA03.

M = diag(m + m2 , m + m2 , m1 , m1 , Jg1 + Jg2 + J)

H.2 Biella equivalente

identici e molto più grossi del momento longitudinale

Consideriamo quindi un’asta di massa m, momento d’inerzia baricentrico trasversale Jg , lunghezza ℓ =

H.3 Robot 3R (capitolo 2.17)

α = atan2(−xv + yu, xu + yv)

citate configurazioni singolari.

L’eq.(2.45) rappresenta il seguente sistema di equazioni

Aggiungendo z 2 alla prima si ottiene

e quindi si giunge al nuovo sistema che sostituisce eq.(2.47)

Noto β, se sin(β) 6= 0, è immediato ricavare γ.

H.4 coordinate angolari (capitolo 3.3)

H.5 Stanford Arm (Pagg.135-137)

Figura H.15: Effetto di due differenti sequenze di rotazione di 90◦ .

H.6 §6.2.3 Sistemi di riferimento (Pag. 136-138)

H.7 Elbow Arm (Pagg.149-153)

H.8 Metodo del grafo di Paul; metodo scissione catena cinema-

Fig. 6.14 Le otto soluzioni del robot 6R (Elbow Arm).

H.9 Metodo dei vincoli equivalenti (§6.7.1 e §6.7.3 Pagg.155)

H.10 Scelta coordinate, cinematica inversa numerica (Pag.170

1. leggi M on % posa da raggiungere

2. poni Q = Q0 % valore di primo tentativo

3. poni j = 0 % contatore iterazioni

6. per i = 1, 2, . . . , n % cinematica diretta (per ogni giunto)

7. calcola Mr = M (qi , Λi ) % matrice di posa relativa del link i (Mi−1,i )

8. calcola L = L(Λi ) % matrice dell’asse del giunto i in terna locale (Li(i−1) )

9. calcola L0 = M LM −1 % matrice dell’asse del giunto in terna 0 (Li(0) )

10. calcola M = M Mr % matrice di posa assoluta del link i (M0,i )

11. calcola J i = F (L0 ) % colonna dello jacobiano

12. fine ciclo i