Dispense Avanzini Meccanica
Dispense Avanzini Meccanica
Dispense Avanzini Meccanica
Giulio Avanzini
Universit`a del Salento
Facolt`a di Ingegneria Industriale
April 16, 2013
Indice
I
Meccanica Newtoniana
3
3
4
7
8
9
9
10
11
12
12
13
15
16
18
2 Leggi di Newton
2.1 Cinematica . . . . . . . . . . . . . . . . . . . . . . . . . . .
2.1.1 Legge oraria del moto . . . . . . . . . . . . . . . . .
2.1.2 Terna intrinseca . . . . . . . . . . . . . . . . . . . .
2.1.3 Accelerazione nella terna intrinseca . . . . . . . .
2.2 Le tre leggi della Dinamica . . . . . . . . . . . . . . . . . .
2.2.1 Definizioni preliminari . . . . . . . . . . . . . . . .
2.2.2 Enunciati . . . . . . . . . . . . . . . . . . . . . . . .
2.2.3 Impulso lineare . . . . . . . . . . . . . . . . . . . . .
2.2.4 Momento di una forza rispetto a un polo fisso O
19
19
20
21
21
23
23
23
24
24
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
ii
INDICE
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
51
51
51
53
55
56
57
INDICE
iii
II
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
Meccanica Lagrangiana
85
85
85
86
87
87
88
89
90
91
92
94
97
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
99
100
100
101
101
102
104
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
107
107
107
109
110
110
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
12 Principi fondamentali
113
12.1 Il Principio dei Lavori Virtuali . . . . . . . . . . . . . . . . . . . 113
12.2 Il Principio di dAlembert . . . . . . . . . . . . . . . . . . . . . . 114
iv
INDICE
12.3 Il Principio di Hamilton . . . . . . . . . . . . . . . . . . . . . . . 115
12.4 Il Principio di Hamilton per sistemi continui . . . . . . . . . . 117
13 Equazioni di Lagrange
119
13.1 Formulazione delle equazioni di Lagrange . . . . . . . . . . . . 119
13.2 Rappresentazione dei fenomeni dissipativi . . . . . . . . . . . . 122
III
127
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
129
130
131
131
131
131
131
132
133
134
135
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
137
138
138
139
140
141
142
142
142
143
.
.
.
.
.
.
.
.
.
149
149
149
150
151
151
152
152
152
155
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
INDICE
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
157
157
157
159
159
160
161
161
161
161
162
163
IV
Stabilit`
a
19 Stabilit`
a: prime definizioni
19.1 Moti ed equilibri L-stabili . . . .
19.1.1 Definizioni preliminari .
19.1.2 Requisito per la stabilit`a
19.1.3 Altre definizioni . . . . .
19.2 Disturbi e incertezze . . . . . . .
19.2.1 Stabilit`
a totale . . . . . .
19.2.2 Stabilit`
a strutturale . . .
19.3 Stabilit`
a di moti periodici . . . .
171
. . . . . . . . . . . .
. . . . . . . . . . . .
secondo Lyapunov
. . . . . . . . . . . .
. . . . . . . . . . . .
. . . . . . . . . . . .
. . . . . . . . . . . .
. . . . . . . . . . . .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
20 Stabilit`
a di un equilibrio
20.1 Linearizzazione nellintorno di un equilibrio . . . . . . . . . . .
20.2 Soluzione di unequazione differenziale lineare a coefficienti
costanti . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
20.3 Condizioni per la stabilit`
a . . . . . . . . . . . . . . . . . . . . .
20.4 Stabilit`
a per sistemi nonlineari . . . . . . . . . . . . . . . . . . .
20.5 Secondo metodo di Lyapunov . . . . . . . . . . . . . . . . . . .
20.5.1 Funzioni definite in segno . . . . . . . . . . . . . . . . .
20.5.2 Teoremi principali per il metodo diretto . . . . . . . . .
20.6 Sistemi conservativi . . . . . . . . . . . . . . . . . . . . . . . . .
173
174
174
174
175
175
176
176
177
179
179
180
180
181
183
183
183
184
21 Stabilit`
a di sistemi meccanici lineari
187
21.1 Sistema massivo . . . . . . . . . . . . . . . . . . . . . . . . . . . 187
21.2 Sistema giroscopico . . . . . . . . . . . . . . . . . . . . . . . . . 188
vi
INDICE
21.3 Richiamo elastico . . . . . . . . . . . . . . . . . . . . . . . . . . . 188
21.4 Smorzamento . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 188
21.5 Frequenze proprie . . . . . . . . . . . . . . . . . . . . . . . . . . 189
22 Esercizi
22.1 Dinamica Newtoniana . . . . .
22.2 Moti relativi . . . . . . . . . . .
22.3 Equazioni di Lagrange . . . . .
22.4 Geometria delle masse . . . . .
22.5 Stabilit`
a e piccole oscillazioni
22.6 Moti rigidi . . . . . . . . . . . .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
191
191
193
194
195
197
199
Parte I
Meccanica Newtoniana
Capitolo 1
Notazione e concetti
introduttivi
1.1
Quando questultima ipotesi non `e valida, i fenomeni appaiono contraddire alcuni dei
pi`
u elementari e intuitivi postulati della Meccanica (e delle Geometria!). Occorre in questo
caso rinunciare a nozioni apparentemente scontate, quali quella di evento contemporaneo,
e fare posto alla Meccanica Relativistica. Tutto questo, per ora, `e lavoro per i Fisici
sperimentali e teorici.
2
Un caso esemplare, in questo ambito, sar`
a quello della dinamica del corpo rigido, idealizzazione di un corpo massivo di dimensioni non trascurabili, che verr`
a assunto come in-
1.2
Per quanto i postulati alla base della Meccanica Classica appaiano intuitivi,
per non dire ovvi, essi sono stati una conquista ottenuta alla fine di un
lungo processo speculativo e conoscitivo che, a partire dalla Fisica di Aristotele, ha condotto fino alla grande rivoluzione galileiana e alla nascita della
scienza moderna, basata sulle evidenze sperimentali e la modellizzazione
matematica. Le intuizioni del grande pisano, suffragate dalla verifica del
comportamento di oggetti reali, portarono alla formulazione, ancorche embrionale, di concetti chiave della Meccanica, quale quello di inerzia (ovvero
di massa) di un corpo o di accelerazione uniforme dovuto alla forza peso.
Lisocronismo delle piccole oscillazioni di un pendolo, dimostrato empiricamente dal Galilei, o le fasi delle lune di Giove, osservate da Arcetri e portate
a supporto della teoria copernicana del moto dei pianeti intorno al Sole,
non avevano ancora alcuna giustificazione teorica. Basteranno per`o pochi
decenni (e il genio di Isaac Newton) per formalizzare in modo matematicamente rigoroso il comportamento di pendolo e pianeti e ottenere una teoria
del moto basato sui tre principi della dinamica e la legge di gravitazione
universale: la fisica segnava il suo primo grande successo nella descrizione
dei meccanismi alla base dei fenomeni osservati. Era nata la Meccanica.
Ci furono anche insuccessi, come quello dellinterpretazione delle forze
aerodinamiche e dei fenomeni luminosi quali la riflessione e la rifrazione. Nel
primo caso linsuccesso fu parziale, in quanto la teoria newtoniana sulla generazione di una resistenza allavanzamento nei mezzi fluidi, non applicabile al
caso di mezzi fluidi continui (gas a pressione ambiente e liquidi) risult`o invece
adeguata alla modellizzazione del moto di corpi in un gas rarefatto, sebbene
il problema venne posto quasi 200 anni dopo la formulazione della teoria.
Per altro la fluidodinamica risulta comunque fondata sui principi base della
deformabile, sotto lazione delle forze e dei momenti esterni che ne determinano lequilibrio
o lo stato di moto.
sempre pi`
u sullanalisi dellinfinitamente piccolo: nascono la teoria atomica e i modelli delle interazioni nucleari. Nasce, soprattutto, la Meccanica
Quantistica che dimostra come i passaggi di energia su scala atomica non
possono avvenire con continuit`a, ma in base a pacchetti di energia discreti
chiamati quanti. Ma questa volta bastarono pochissimi anni dai due fondamentali lavori di Einstein sulla Teoria della Relativit`a generale e ristretta per
una ulteriore, profondissima, rivoluzione. Il principio di indeterminazione
di Heisemberg, sancendo che non si possono misurare simultaneamente posizione e quantit`
a di moto di una particella, minava, in un certo senso, le
fondamenta stesse della Fisica sperimentale. In questo caso non `e pi`
u lo
stato di moto dellosservatore, ma il suo stesso osservare a rendere impossibile unosservazione neutra rispetto al fenomeno studiato.
I fenomeni su scala atomica dimostravano di sfuggire al meccanicismo
causalistico che aveva generato la fiducia positivistica di riuscire, prima o
poi, a interpretare qualunque fenomeno fisico sulla base di un opportuno
modello matematico verificato attraverso prove sperimentali. Da questo
punto di vista la Teoria della Relativit`a era inserita con continuit`a sulla
Meccanica Classica, riuscendo semplicemente a fornire un modello valido
per condizioni in cui le teorie di Newton perdevano di validit`a, senza rinunciare per`
o al fondamento teorico di una conoscenza certa e assoluta dei
fenomeni fisici a partire da condizioni misurabili. Ed Einstein, cos` come a
suo tempo Newton, dimostr`o i limiti del suo genio nel non accettare la sfida
dellindeterminazione.3 Per altro, allinizio del XX secolo, seri problemi di
prevedibilit`
a si manifestarono anche nellambito di fenomeni deterministici. Il meccanismo della transizione al caos fa s` che piccole variazioni nelle
condizioni iniziali portino a unevoluzione di alcuni sistemi completamente
diversa e, di fatto, non prevedibile, data la precisione (inevitabilmente finita)
con cui `e possibile determinare lo stato iniziale del sistema stesso.
Il Principio di Indeterminazione di Heisemberg e, sebbene in misura minore, la scoperta di fenomeni caotici, hanno consentito il salto epistemologico che ha portato a rinunciare alla pretesa (che ora sappiamo infondata)
di poter descrivere analiticamente in modo esatto con mezzi matematici e
misure i fenomeni fisici su scala atomica e subatomica. Il modello matematico diviene a questo punto non pi`
u la verit`a in senso conoscitivo, ma pi`
u
modestamente, uno strumento descrittivo che consente di prevedere in modo
anche accurato landamento di fenomeni fisici che avvengano in condizioni
note e misurabili. Il cambiamento di approccio al problema `e radicale: la
matematica rimane s` il linguaggio con cui leggere il mondo fisico, come
gi`
a indicato dal Galilei ne Il Saggiatore,4 ma il modello matematico diviene
3
La famosa frase attribuita ad Einstein, Dio non gioca a dadi con il mondo, testimonia
la difficolt`
a di accettare aspetti fondamentali (e a quel punto ampiamente comprovati
sperimentalmente) della moderna fisica delle particelle.
4
La filosofia `e scritta in questo grandissimo libro che continuamente ci sta aperto innanzi
a gli occhi (io dico luniverso), ma non si pu intendere se prima non simpara a intender
` DI MISURA
1.3. UNITA
1.3
Unit`
a di misura
1.4
Per molto tempo `e stato utilizzato il sistema pratico degli ingegneri, che aveva come
grandezze fondamentali la lunghezza, il tempo e la forza, ragion per cui era la massa
ad essere una grandezza derivata, avente come dimensioni quelle di una forza diviso
unaccelerazione, ovvero [F L1 T 2 ], dove il chilogrammo era lunit`
a di misura della forza,
valutata sul peso di un campione, e la massa si misurava quindi in kg m1 s2 .
6
Si osservi come in questultimo caso lunit`
a risultante risulti molto piccola, e pertanto
si usano spesso dei multipli (e.g. il kPa). Per altro la pressione `e talvolta espressa in unit`
a
diverse, frutto dei diversi procedimenti di misura (come i millimetri di mercurio, mmHg) o
di diversi campioni, (come latmosfera, atm, il cui valore equivale allincirca alla pressione
atmosferica standard a livello del mare).
7
Non ci si soffermer`
a sul procedimento di misura delle grandezze meccaniche, argomento
che esulta dalle finalit`
a di un corso di Meccanica Razionale, presupponendo che, per
tutte le grandezze considerate in seguito, sia disponibile uno strumento che permetta di
misurarle in modo preciso. La teoria della misura, ovvero la descrizione degli strumenti
e dei procedimenti tecnici che permettono di effettuare la misura delle grandezze fisiche
e le tecniche statistiche che ne stimano lincertezza, sono oggetto dei corsi di Misure
Meccaniche.
1.5
1.5.1
(ovvero u
v + w)
v+u
La differenza di due vettori (v. Fig. 1.3.a) `e data dalla somma del
= u
v
= u
+ (
primo pi`
u lopposto del secondo vettore, w
v ), dove il
, ma verso opposto, tale che
vettore
v ha la stessa direzione e modulo di v
+ (
v
v ) = 0.
= ku
, fornisce
La moltiplicazione di un vettore per uno scalare, v
10
b)
v
u +v
v
u
Fig. 1.1: Somma di due vettori: a) regola del parallelogramma; b) concatenazione dei vettori.
v
u
u + v + w +r
1.5.2
Prodotto scalare
a)
-v
b)
ku
u -v
u
b)
u xv
Fig. 1.4: a) angolo convesso compreso fra due vettori e b) prodotto vettoriale.
(1.1)
dove `e langolo fra le direzioni dei due vettori. Si osservi che il prodotto
v
=v
u
e u
(
=
scalare `e commutativo e distributivo, ovvero u
v + w)
v
+u
w.
u
Il prodotto scalare di due vettori perpendicolari `e nullo, mentre se due
vettori sono paralleli, il prodotto scalare `e dato dal prodotto dei loro moduli.
Se un vettore viene moltiplicato scalarmente per se stesso, si ottiene quindi
il quadrato del suo modulo,
u
= u2
u
1.5.3
Prodotto vettoriale
v
`e dato da un vettore w
diretto
Il prodotto vettoriale di due vettori u
lungo la normale al piano individuato dai due vettori, avente come modulo il
prodotto dei moduli dei vettori per il seno dellangolo acuto (ovvero inferiore
a rad) compreso fra essi (Fig. 1.4.b). Posto che
=u
v
(1.2)
risulta
= uv sin , w
u
=w
v
=0
w
lungo la perpendicolare al piano identificato da
La scelta del verso di w
ev
viene effettuata in modo tale che i tre vettori, u
, v
ew
siano diretti
u
come il pollice, lindice e il medio della mano destra, disposti allincirca
ev
sono paralleli,
perpendicolarmente fra loro.9 Si pu`o osservare che, se u
9
In alternativa, si pu`
o dire che il prodotto vettoriale `e diretto lungo la direzione del
pollice della mano destra che si chiude a pungo nella direzione in cui il primo vettore ruota
ev
dalla direzione di w,
12
(
u+v
1.6
(1.4)
1 if i = j
0 if i j
1.7
Componenti di un vettore
(1.5)
v3
e3
13
v2
e2
e1
v1
data una terna ordinata di numeri reali, ad essa corrisponde uno e un solo
vettore.10
Per le grandezze vettoriali la freccia comparir`a solo quando le grandezze
corrispondenti sono evocate indipendentemente da un particolare riferimento.
Una volta scelto un opportuno sistema di riferimento `e possibile per`o rappresentare una grandezza vettoriale per componenti tramite un vettore colonna
indicato dallo stesso simbolo del grassetto usato per il vettore, senza freccia,
ma con un pedice che indica il particolare riferimento in cui sono rappresentate le componenti.
1.8
Trasformazione di coordinate
i sono rappreSi osservi che nello spazio delle terne ordinate (x, y, z)T R3 i versori e
3 = (0, 0, 1)T .
1 = (1, 0, 0)T , e
2 = (0, 1, 0)T e e
sentati dalle terne e
11
Per semplicit`
a si ipotizza che il riferimento G abbia la stessa origine di F.
14
= [
e1G
e2G
e3G ]
Propriet`
a della matrice di trasformazione
Si osservi che la matrice T GF `e formata ordinando per colonne le compoi espressi nella terna G, versori mutuamente ortogonali.
nenti dei versori e
Pertanto le colonne della matrice risultano essere tutte di norma unitaria e
rappresentano versori di R3 perpendicolari fra loro. Una matrice cos` formata `e detta matrice ortogonale e gode di alcune interessanti propriet`a:
1. linversa di una matrice ortogonale T coincide con la sua trasposta:
T 1 = T T ;
2. il determinante di una matrice ortogonale `e det(L) = 1; se la matrice
`e una matrice di trasformazione di coordinate il determinante `e 1;
3. se T 1 e T 2 sono matrici ortogonali, il loro prodotto righe per colonne
T 1 T 2 `e ancora una matrice ortogonale.
...
T
T GF = e
2
. . I.
T
e
3I
Il lettore dimostri come esercizio che valgono le relazioni
T GF = [
g 1F
g 2F
g 3F ] e T F G
1.9
T1F
g
...
T2F
g
...
T3F
g
Il prodotto scalare di due vettori pu`o essere facilmente espresso per componenti, tenendo presente lassociativit`a delloperazione. Consideriamo un
= u1 e
1 + u2 e
2 + u3 e
3 e
riferimento cartesiano ortogonale e due vettori u
= v1 e
1 + v2 e
2 + v3 e
3 . Il loro prodotto scalare `e dato da
v
v
= (u1 e
1 + u2 e
2 + u3 e
3 ) (v1 e
1 + v2 e
2 + v3 e
3 )
u
1 (v1 e
1 + v2 e
2 + v3 e
3 ) +
= u1 e
2 (v1 e
1 + v2 e
2 + v3 e
3 ) +
+ u2 e
3 (v1 e
1 + v2 e
2 + v3 e
3 )
+ u3 e
i e
j = ij , si ha che
Tenendo conto delle relazioni e
v
= u1 v1 + u2 v2 + u3 v3 = uTF v F
u
Procedendo in modo analogo per il prodotto vettoriale si ha che
v
= (u1 e
1 + v2 e
2 + v3 e
3 )
1 + u2 e
2 + u3 e
3 ) (v1 e
u
2 + v3 e
3 ) +
1 (v1 e
1 + v2 e
= u1 e
1 + v2 e
2 + v3 e
3 ) +
2 (v1 e
+ u2 e
3 (v1 e
1 + v2 e
2 + v3 e
3 )
+ u3 e
(1.7)
16
=u
v
pu`o essere scritto
Il vettore w
minante della matrice
RRR e
RRR 1
=u
v
= RRR u1
w
RRR
RR v1
2 e
3 RRRR
e
R
u2 u3 RRRR
R
v2 v3 RRRR
sono date dai minori principali della
Di conseguenza, le componenti di w
matrice
u u2 u3
1
v1 v2 v3
con
wF = (u2 v3 u3 v2 , u3 v1 u1 v3 , u1 v2 u2 v1 )T
` possibile anche una rappresentazione matriciale del prodotto vettoriale
E
v F , dove la matrice antisimmetrica
per componenti, dove wF = uF v F = U
U `e data da
0 u3 u2
0 u1
U = u3
u2 u1
0
1.10
(1.8)
Non verranno considerati operatori lineari fra spazi di dimensione diversa. Per
una esposizione completa della teoria degli operatori lineari il lettore pu`
o consultare un
qualunque manuale di Algebra Lineare.
(1.9)
dove ciascuno degli n vettori Abi ammette una sua rappresentazione per
componenti. Scriviamo tale rappresentazione nella forma
n
AB =
18
1.11
Capitolo 2
Leggi di Newton
2.1
Cinematica
Lanalisi `e per limitata ai corpi per i quali sia definita la sola posizione P
nello spazio S, ovvero corpi la cui dimensione sia trascurabile rispetto al (o
ininfluente sul) loro spostamento. Supporremo il corpo dotato di massa m e
lo definiremo per questo motivo elemento materiale1 . Lelemento materiale occupa dunque la posizione P (t) identificata dal vettore posizione r (t)
rispetto allorigine O di un riferimento cartesiano opportunamente definito
(Fig. 2.1). La curva C S definita come
t0
r (t + t) r (t)
r d
r
= lim
=
t0 t
t
dt t
(t) = lim
a
1
In accordo con il testo di Piero Giorgio Bordoni, Lezioni di Meccanica Razionale, Ed.
Veschi, Roma, ristampa 1982, si `e preferita la dizione elemento materiale in luogo di punto
materiale, pure tipica della tradizione italiana, o quella di particella, tipica della tradizione
anglosassone (particle). Con il termine elemento si indica un corpo dotato di massa che
occupa una posizione P nello spazio. La dizione punto materiale `e considerata infatti
fuorviante dal Bordoni, in quanto il punto come ente geometrico `e per definizione privo di
dimensione, mentre i corpi considerati dalla Meccanica possono avere dimensioni tuttaltro
che trascurabili, anche nel caso in cui si consideri come rilevante la loro sola posizione. Un
esempio `e costituito dai moti orbitali, in cui la dimensione dei corpi celesti rende difficile
assimilarli a punti nello spazio. Considerazione analoga vale per il termine particella,
che suggerisce una dimensione limitata se non proprio puntiforme delloggetto studiato.
19
20
v (t)
P(t+t)
r (t)
z
e3
r (t+t)
e2
C
O
e1
dx dy dz T
, , ) = (x,
y,
z)
T
dt dt dt
aF (t) = (
d2 x d2 y d2 z
,
,
) = (
x, y, z)T
dt2 dt2 dt2
2.1.1
s(t) =
v(t)dt
t0
2.1. CINEMATICA
21
v(t) = x 2 + y 2 + z 2
2.1.2
Terna intrinseca
lungo la direzione n.
= t n
`e detta binormale, ed `e perpendicolare
La direzione del versore b
formano la
eb
al piano su cui giace il cerchio osculatore. I tre versori t, n
rispetto ad s
terna intrinseca (Fig. 2.2). Il modulo della derivata di b
rappresenta la rapidit`
a con cui il piano osculatore ruota intorno a t ed `e
detto torsione,
= db/ds
2.1.3
d
v
dt ds
= st + s
.
dt
ds dt
Questo risultato vale per qualunque vettore di modulo costante, la cui derivata risulta
(t) e assumiamo che
perpendicolare al vettore dato. Consideriamo il generico vettore u
u
, il che vuol dire
22
P(t0)
s=0
s
Rc
C
t^
v
^
n
^
b
at
v (t)
an
a (t)
a
da cui risulta evidente come laccelerazione abbia due sole componenti in tale
riferimento, lungo la tangente e la normale alla traiettoria, dette, appunto,
componenti tangenziale e normale dellaccelerazione (Fig. 2.3).
Dallultima espressione `e inoltre possibile determinare il valore della curvatura,
c=
a st
s 2
2.2
2.2.1
23
, si
Dato un elemento materiale di massa m che si muove con velocit`a v
il prodotto della sua massa per la velocit`a,
definisce quantit`
a di moto p
= m
p
v . Considerando un generico polo Q (ad esempio, ma non necessariamente, lorigine O del riferimento cartesiano ortogonale rispetto al quale
viene descritto il moto), si definisce momento della quantit`
a di moto (o
il prodotto vettoriale h
= r p
= m
.
momento angolare) h
rv
2.2.2
Enunciati
Postulando lesistenza di un riferimento assoluto, detto riferimento inerziale (o galileiano), `e possibile formulare i tre principi fondamentali della
meccanica newotniana:
In assenza di forze esterne agenti su di esso, un elemento materiale si
muove di moto rettilineo a velocit`a costante.
La forza agente su un elemento materiale `e pari alla derivata temporale
della sua quantit`
a di moto.
Quando due corpi si scambiano una mutua azione meccanica, le forze
agiscono lungo la congiungente i due corpi, sono uguali in modulo e
opposte in verso.
In termini formali il secondo principio della dinamica indica che
d
p d(m
v)
f =
=
= m
v + m
a
dt
dt
Per sistemi a massa costante, m
= 0, e si ottiene la classica formula forza
uguale massa per accelerazione:
f = m
a
valida solo ed esclusivamente in un riferimento inerziale. Si osservi come il
primo principio, gi`
a enunciato dal Galilei, ricada come caso particolare del
secondo, in quanto in assenza di forze la derivata della quantit`a di moto
risulta nulla e, di conseguenza, la quantit`a di moto risulta vettorialmente
costante (principio di conservazione della quantit`a di moto).
Infine, il terzo principio `e enunciato sopra nella sua formulazione forte, in
cui la retta dazione della forza scambiata fra due corpi A e B mutuamente
interagenti `e diretta lungo la loro congiungente AB. Nella formulazione
debole si postula invece solo che, detta f AB la forza che il corpo A esercita
24
(2.1)
m
z = fz
2.2.3
Impulso lineare
t2
f dt =
t1
t2
t1
d
p
(t2 ) p
(t1 )
dt = p
dt
2.2.4
Si consideri lespressione del Secondo Principio della Dinamica e se ne moltiplichino vettorialmente ambo i membri per r da sinistra:
r f = r
d
p
dt
t1
(
r f )dt =
t2
(
r
t1
d
p
) dt
dt
(2.2)
25
= d
Ricordando che v
r /dt, si ha che
d
p d
r
d(
r p)
= r
+
p
dt
dt dt
`e per definizione parallelo a v
=
dove lultimo termine `e nullo, in quanto p
d
r /dt, quando il polo O `e fisso. Pertanto lequazione (2.2) pu`o essere scritta
nella forma
t2
t1
(
r f )dt =
t2
t1
t2
t1
)
d(
rp
dt
dt
dh
2 ) h(t
1)
dt = h(t
dt
t1
2 ) h(t
1)
= h(t
mdt
26
Capitolo 3
Lavoro
3.2
Energia cinetica
f d
r =
= m
=
d2 r d
r
dt
2
dt
tA dt
tB d
d 1 d
r d
r
1
(v 2 ) dt
(
) dt = m
dt 2 dt dt
2
tA dt
m
a d
r = m
tB
tA
1
2
2
)
m (vB
vA
2
27
tB
28
$"!!$
#
!"!!$
%"!!$
#&
!"!!"#!$
#'
!
"
#%
WAB =
f d
r = TB TA
(3.1)
3.3
Forze conservative
WAB =
f d
r
(3.2)
29
&%
!
!!
!#
!!
!#
!"
"
"
` evidente che forze che non cambiano segno lungo una traiettoria chiusa
E
non possono essere conservative. Un esempio di tali forze `e dato da tutte le
forze di attrito, sempre diretto in verso opposto a quello del moto, che sono
forze dissipative, tali cio`e da ridurre (dissipare) lenergia del sistema.
La verifica della condizione espressa da Eq. (3.2) risulta difficile nella
maggior parte dei casi pratici, ma, grazie al teorema si Stokes si pu`o sostituire la condizione espressa sullintegrale di una generica traiettoria chiusa
con una condizione locale sul rotore del campo di forza. Loperatore rotore,
indicato con il simbolo () o, talvolta, come rot(), applicato a un generico
in tre dimensioni `e dato da
campo vettoriale v
=(
v
vy vx
vz vy
vx vz
) e1 + (
) e2 + (
) e3
y
z
z
x
x
y
il che `e possibile, per qualunque curva chiusa e campo di forza sufficientemente regolare, se e solo se
f (
r) = 0
Daltra parte, condizione necessaria e sufficiente affinch`e il rotore di un
campo vettoriale sia nullo ovunque `e che il campo vettoriale sia il gradiente
30
V = (
WAB =
f d
r=
V d
r = VA VB
(3.3)
WAB =
f d
r = TB TA = VA VB .
V V V T
,
,
)
x y
z
)=
(
)+
(
)=0
y
z
y z
z y
per il terorema di Schwarz sullinversione nellordine delle derivate miste. Analogamente
si annullano tutte le altre componenti di f e dunque, se il campo vettoriale ammette
un potenziale, allora il suo rotore `e nullo. La necessit`
a della condizione `e pi`
u ostica da
dimostrare e a questo proposito il lettore `e rimandato ai testi di Analisi Matematica.
31
V(0) V(x) =
f (x)dx =
kxdx = (1/2)kx2
a = 2E/k
vmax = a k/m
Tutte queste considerazioni permettono di analizzare molte propriet`a del moto
delloscillatore, prescindendo dalla soluzione dellequazione differenziale che ne descrive
il moto. In questo semplice esempio, per`o, `e disponibile una soluzione analitica per il
32
E,V
b)
k
m
0
x
-xmax
xmax
Capitolo 4
Sistemi di N elementi
materiali.
4.1
Definizioni
mi r i = f i + (1 ij )f ji
(4.1)
j=1
4.2
N N
i=1
i=1
i=1 j=1
mi r i = f i + (1 ij )f ji
33
34
#%
#%$!
%%%
!"
"& #$
!#
#!
&
%%%
##
#$'"
""
!
%%%
#"$!
#"
#"
!!
dove il primo termine del secondo membro pu`o essere posto pari alla risul
tante delle forze esterne agenti sul sistema, f tot = N
i=1 f i , mentre il secondo termine del secondo membro risulta nullo per il terzo principio della
dinamica, dal momento che f ji = f ij . Moltiplicando e dividendo il
primo membro per la massa totale M e ricordando la definizione di centro
di massa si ha che
M r C = f tot
tot pari alla derivata della quantit`a di moto totale del sistema,
con M r C = p
N
i=1
i=1
tot = p
i = mi v
i = M r C .
p
(4.2)
Dunque, la posizione del centro di massa del sistema di N elementi materiali si comporta come un elemento fittizio di massa pari alla massa totale
M del sistema, mosso da una forza pari alla risultante delle forze esterne.
Le forze interne non hanno effetto sul moto del centro di massa e pertanto,
in assenza di forze esterne o per forze esterne a risultante nulla, il centro di
massa si muove di moto rettilineo uniforme e la quantit`a di moto totale del
sistema risulta costante, il che comporta la dimostrazione del seguente
4.3
35
i=1
i=1
i=1
A = r Ai mi r i + r Ai mi r i
h
tot
` possibile scrivere i vettori r A e r i come somma del vettore posizione del
E
i
centro di massa (contato rispetto ad A ed O nei due casi, rispettivamente)
i dellelemento iesimo contato a partire dal centro
pi`
u il vettore posizione
di massa stesso, ovvero (Fig. 4.2)
i
i , r i = r C +
r Ai = r AC +
Si noti che, dalla definizione di centro di massa, risulta che1
N
i = 0.
mi
i=1
Sostituendo queste espressioni nella definizione della derivata del momento angolare totale e ricordando la definizione di centro di massa si ha
che
N
A = (r AC +
i ) mi (r C +
i ) + r Ai f i + (1 ij )f ji
h
tot
i=1
j=1
i=1
A
= M r AC r C + m
tot
dove
A
Ai f i
m
tot = r
i=1
`e il momento totale delle forze esterne rispetto al polo A, mentre il momento totale delle forze interne, che agiscono a coppie lungo la stessa direzione e hanno pertanto lo stesso braccio, risulta nullo. La cosa pu`o
1
36
'
"'
!"
!#
#%
##
"'& &
%%%
#%$!
#
%%%
"
#"$!
#
"&
$ "'"
#"
""
$"
!!
Ai (N
e composta da coppie di termini del tipo
N
i=1 r
j=1 (1 ij )f ji ) `
r Ai f ji + r Aj f ij = (
r Ai r Aj ) f ji = 0
in quanto, assumendo la validit`a del terzo principio della dinamica nella
sua formulazione forte, i vettori (
r Ai r Aj ) e f ji sono paralleli, entrambi
diretti lungo la congiunte gli elementi i e j, implicando con ci`o la nullit`a del
prodotto vettoriale.
Vi sono due scelte per il polo A che portano a equazioni di bilancio del
momento angolare di forma particolarmente semplice. Se scegliamo come
polo A lorigine del riferimento cartesiano ortogonale cui riferiamo il moto,
risulta r AC = r C e dunque il primo elemento nel membro di destra si annulla
perche `e nullo il prodotto vettoriale di un vettore per se stesso. Abbiamo
cos` lequazione di bilancio espressa nella forma
O = m
O
h
tot
tot
Se invece scegliamo come polo A il centro di massa C, il vettore r AC
risulta nullo e lequazione di bilancio del momento angolare `e data da una
forma analoga alla precedente:
C = m
C
h
tot
tot
Se dunque scegliamo come polo di riduzione dei momenti lorigine del
riferimento cartesiano ortogonale O o il centro di massa del sistema di elementi materiali C, la derivata del momento angolare calcolata rispetto al
polo scelto risulta pari al momento totale delle forze esterne valutato rispetto
allo stesso polo, il che permette di dimostrare banalmente il seguente
37
4.4
A = r Ai mi r i
h
tot
i=1
A = (
i ) mi (r C +
i )
h
r AC +
tot
i=1
N
i mi
i
= r AC M r C +
i=1
rel
tot + h
= r AC p
tot
dove
N
rel =
i mi
i
h
tot
i=1
i=1
i=1
i=1
rel
C
i mi
i mi (
i mi
i ) =
i h
h
r i =
r C +
tot =
tot .
38
4.5
1 N
i )]
i ) (r C +
[mi (r C +
2 i=1
N
1
1 N
i
i
r C r C + r C mi
i + mi
2
2 i=1
i=1
4.6
Si assuma ora che tutte le forze in gioco, tanto quelle esterne, f i , quanto
quelle interne, f ji , siano conservative, ovvero esprimibili come il gradiente
di un potenziale. Ci`
o significa che le forze esterne sono date da
f i = V ext i
dove il pedice i indica che il gradiente `e calcolato in r i . Le forze interne
sono date da
int
j
f ji = Vijint i = f ij = Vji
40
Capitolo 5
Fa eccezione il pianeta Mercurio che, per la sua vicinanza al Sole, risente di effetti
relativistici legati alla curvatura dello spaziotempo indotta dal campo gravitazionale del
Sole stesso, effetto trascurabile per tutti gli altri corpi che si muovono nel sistema solare.
Lanomalia orbitale di Mercurio e la sua interpretazione nellambito della teoria della
relativit`
a ristretta costitu` una delle prime prove empiriche per la teoria di Einstein.
m1
e3
f21= -f
r1
e2
f12= f
m2
rC
r2
e1
41
42
r1
e2
rC
m2
r2
O
e1
5.1
1
(m1 r 1 + m2 r 2 )
m1 + m2
si trova lungo la congiungente i due punti P1 e P2 che individuano la posizione dei due elementi materiali (Fig. 5.2). Questo lo si pu`o dimostrare
facilmente scrivendo
1 ; r 2 = r C +
2
r 1 = r C +
2 le espressioni
1 e
Sostituendo la definizione di r C , si ottengono per
1
(m1 r 1 + m2 r 2 )
m1 + m2
m2
m2
(
r 1 r 2 ) =
r
m1 + m2
m1 + m2
1 = r 1
2 , si ottiene
Procedendo analogamente per
2 =
m1
r
m1 + m2
43
(m1 + m2 )r C = 0
ovvero, in assenza di forze esterne il centro di massa C ha accelerazione
nulla e la quantit`
a di moto totale del sistema `e costante. Sottraendo la
prima equazione divisa per m1 dalla seconda divisa per m2 si ottiene invece
1
1
r 2 r 1 = r =
f 12
f
m2
m1 21
Ponendo f 12 = f e ricordando il terzo principio della dinamica si ottiene
lequazione del moto per il vettore r , che fornisce la posizione relativa di m2
rispetto ad m1 :
r = m1 + m2 f
m1 m2
da cui si deduce che il vettore che esprime la posizione relativa del secondo
elemento rispetto al primo si comporta come il vettore posizione di un elemento materiale isolato di massa m = (m1 m2 )/(m1 + m2 ) sottoposto alla
sola forza f .
Dal momento che il centro di massa C non subisce accelerazione, esso
pu`o essere assunto come origine del riferimento inerziale, senza perdita di
generalit`
a, e levoluzione nel tempo di r `e sufficiente a descrivere lintero
sistema di due elementi.
Osservazione: Per m1 m2 il centro di massa del sistema coincide, di
fatto, con la posizione della massa pi`
u grande, che risulta quindi praticamente fissa nello spazio. Il vettore posizione dellelemento di massa pi`
u
piccola diviene dunque r 2 r e ai fini pratici si torna al caso di una singola
massa m2 che si muove rispetto a un punto assunto come origine immobile
del riferimento inerziale. Inoltre, per m2 /m1 0 la massa dellelemento
equivalente m = m1 m2 /(m1 + m2 ) m2 .
5.2
Dal momento che la forza f agisce lungo la congiungente i due elementi, essa
ha braccio nullo rispetto al centro di massa e dunque il momento totale delle
forze agenti sul sistema (in assenza di altre forze esterne) rimane identica `e costante.
mente nullo. Questo implica che il momento angolare totale h
Scegliendo come origine il centro di massa C (e assumendo dunque che
si pu`
r C = r C = 0)
o dimostrare che
C = r 1 m1 r 1 + r 2 m2 r 2
h
tot
m2
m1 m2
m1
m1 m2
r ) +
r
= (
r ) (
r
m1 + m2
m1 + m2
m1 + m2
m1 + m2
m1 m2
=
r r
m1 + m2
44
.
r
e^2
m2
r
e^3
m1
v
u^r
e^1
Ci`
o implica che i vettori r e r giacciono sempre sullo stesso piano e dunque
il moto relativo di un elemento rispetto allaltro (o quello di entrambi gli
elementi rispetto al comune centro di massa) `e vincolato ad appartenere a
un piano invariante, detto piano del moto P.
Si osservi che momento angolare totale e relativo, in questo caso, co e con h il suo
C = h
incidono. Nel seguito si indicher`a per semplicit`a h
tot
modulo.
5.3
Coordinate polari
Date le caratteristiche ora evidenziate del moto, confinato su un piano invariante, risulta conveniente esprimere il vettore r in coordinate polari, (r, ),
dove r =
r `e, come al solito, il modulo del vettore r (ovvero la distanza
fra i due elementi), mentre rappresenta una distanza angolare contata
rispetto a una direzione di riferimento arbitraria, fissata sul piano P, che,
stante quanto sopra discusso, risulta fisso nel riferimento inerziale (Fig. 5.3).
Senza perdita di generalit`a, `e possibile associare a tale direzione il primo ver1 , in modo tale che e
2 giaccia su P ed e
3 sia ad
sore della terna inerziale, e
3 = h/h.
esso perpendicolare, in modo tale che e
r e u
i versori radiale (parallelo a r ) e trasversale (perpenSiano u
r e tale che u
r > 0). Stante la scelta fatta per i versori del
dicolare a u
riferimento cartesiano ortogonale, risulta
(
ur )F = (cos , sin , 0)T ; (
u )F = ( sin , cos , 0)T
Il terzo versore della terna polare, normale al piano del moto, risulta dato
n = u
r u
e
3 .
da u
r = r u
r = r u
e laccelerazione nella forma
u
r = ru
r + (2r + r)
r + r u
+ ru
+ ru
r + r u
= (
r r2 ) u
5.4
46
v
e^2
m2
r
r d
m1
e^1
5.5
5.5.1
d
r = dx
e1 + dy
e2 = dru
il gradiente di F `e dato dal vettore
F =
F
F
F
1 F
1 +
2 =
r +
e
e
u
u
x
y
r
r
5.5.2
Condizione di conservativit`
a
47
r
La condizione `e necessaria in quanto, per avere una forza f = V = f u
diretta lungo la direzione radiale, deve essere nulla la componente trasversale. Osservando che la componente trasversale di forza `e data da
f = (V) =
1 V
=0
r
V(r) = V0 +
5.5.3
f (r)dr
(5.1)
r0
d 2
(r ) = 0 m r2 = h = cost
dt
Questo significa che nella prima equazione in (5.2) il secondo termine nella
parentesi tonda pu`
o essere riscritto nella forma r2 = h2 /(m 2 r3 ). Inoltre
r =
d2 r d dr
d dr dr dr
1 d 2
=
( )=
( )
=
r =
(r )
2
dt
dt dt
dr dt dt dr
2 dr
1 d 2
h2
(r ) 2 3 ] f (r) = 0
2 dr
m r
48
v = 2 [E V(r)] /m
quanto la sua componente radiale
2
h2
vr = r =
[E
V(r)
]
m
2m r2
dipendono cos` solo dalla distanza r fra i due elementi materiali del sistema.
5.5.4
Traiettoria
2
m r 2
dr dr dt r
2 (E V(r) h )
=
= =
d dt d
h
m
2m r2
la cui inversa `e integrabile formalmente come3
R
h
dr
(R) = 0
2
2m r0 r E V(r) h2 /(2m r2 )
U
h
du
2m u0
E V(u) u2 h2 /(2m )
(5.3)
che pu`
o essere integrata in alcuni casi di interesse pratico e, in particolare,
quando la forza `e funzione di una potenza di r, ovvero
f (r) = Krn = Kun
2
49
In questo caso
V(r) =
K (n+1)
K n+1
r
=
u
n+1
n+1
du
u0
a bu2
+ cu(n+1)
che risulta integrabile in modo analitico in diversi casi di interesse, fra cui
quello delle forze gravitazionali, per cui n = 2.
50
Capitolo 6
Moti Kepleriani
6.1
Attrazione gravitazionale
K
r2
dove il segno + vale per le forze repulsive, come quelle elettrostatiche fra
cariche dello stesso segno, mentre il vale per forze attrattive, come le forze
elettrostatiche fra cariche di segno opposto o la forza gravitazionale. In
questultimo caso fu Newton a proporre la legge che porta il suo nome, per
cui la forza di attrazione fra due corpi di massa m1 e m2 `e proporzionale
al prodotto delle masse e inversamente proporzionale al quadrato della distanza, cos` che risulta
m1 m2
f (r) = G 2
r
11
2
2
dove G = 6.67 10
N m kg `e la costante di gravitazione universale.
6.2
(6.1)
52
(6.2)
r r2 = 2
r
Si noti che a questo punto levoluzione del sistema non dipende pi`
u ne
dalle singola masse m1 e m2 degli elementi, ne dalla massa dellelemento
equivalente, m = m1 m2 /(m1 + m2 ), ma solo dal valore del parametro gravi` quindi possibile definire la forza di attazionale, = G(m1 + m2 ). E
trazione gravitazionale per unit`a di massa equivalente, f /m = /r2 e, di
conseguenza, una energia potenziale per unit`
a di massa e unenergia
cinetica per unit`
a di massa date rispettivamente da
r
V(r) =
dr =
2
r
r
(r ) 3 = 2
2 dr
r
r
che pu`
o essere integrata rispetto ad r. Tale operazione fornisce
1 2 k2
(r + 2 ) = cost
2
r
r
1`
6.3
dr d k dr
d 1
du
= 2
= k ( ) = k
d dt r d
d r
d
r =
du
d2 u
dr d k d
= 2
(k ) = k 2 u2 2
d dt r d
d
d
Di conseguenza
d2 u
k 2 u3 = u2
d2
k 2 u2
+u= 2
2
d
k
Lequazione differenziale ordinaria lineare del secondo ordine non `e altro
che lequazione del moto delloscillatore armonico di pulsazione unitaria e
forzante costante. La sua soluzione generale pu`o essere espressa nella forma
u() =
+ C cos( 0 )
k2
54
si ottiene
C2 =
2k 2 E
2
(1
+
)
k4
2
In realt`
a si `e soliti considerare come terzo integrale del moto il valore della
costante
1/2
2k 2 E
e = (1 + 2 )
= 2 [1 + e cos( 0 )]
r() k
ovvero
r() =
k 2 /
1 + e cos( 0 )
Lequazione sopra scritta rappresenta inoltre lespressione di una generica conica rappresentata in coordinate polari2 di parametro pari a p = k 2 /
ed eccentricit`
a pari ad e. Al variare di e le orbite possono essere circolari (e = 0), ellittiche (0 < e < 1), paraboliche (e = 1) o iperboliche (e > 1).
Quanto sopra conferma e dimostra su basi teoriche la validit`a della prima
legge di Keplero, che affermava come le orbite dei pianeti fossero delle ellissi
in cui il Sole occupa uno dei due fuochi.3
Figura 6.1 riporta la geometria di unorbita per il caso ellittico. La
distanza minima fra i due corpi viene ottenuta al periasse, per = 0 ,
dove si ha rP = rmin = (k 2 /)/(1 + e). Pertanto langolo 0 rappresenta
la direzione in cui si annulla la derivata r,
cambiando segno da negativa
a positiva. Per orbite chiuse ellittiche si ha anche una distanza massima
ottenuta allapoasse4 per 0 = , pari a rA = rmax = (k 2 /)/(1 e).
Periasse e apoasse sono entrambi paralleli allasse maggiore (detto anche
linea degli apsidi), ragion per cui il semiasse maggiore risulta pari ad
a = (rP + rA )/2 = (k 2 /)/(1 e2 ). Data la relazione fra parametro p = k 2 /,
semiasse maggiore, apoasse e periasse, questi ultimi possono essere ottenuti
come
rA = a(1 + e) e rP = a(1 e).
2
55
/2
rP
rA
2a
6.4
Periodo orbitale
A = ak a/.
5
Si noti che per e = 0 si ha r = a = b =cost e si ottiene la ben nota formula per larea
del cerchio, A = r2 .
56
2.5
2.0
2.0
1.5
1.5
1.0
1.0
0.5
0.5
0.0
0.0
- 0.5
- 0.5
- 1.0
- 1.0
- 1.5
- 1.5
- 2.0
- 2.0
a)
- 2.5
b)
- 3.0
- 2.5
2.5
- 2.0 - 1.5 - 1.0
- 0.5
0.0
0.5
1.0
1.5
- 2.5
2.0 - 3.0
- 2.5
- 2.0
- 1.5
- 1.0
- 0.5
0.0
0.5
1.0
1.5
2.0
2.0
1.5
1.0
0.5
0.0
- 0.5
- 1.0
- 1.5
- 2.0
c)
- 2.5
-5
-4
-3
-2
-1
A
T=
= 2
dA/dt
a3
6.5
Energia totale
2k 2 E
2
(1 e2 )2
2k 2
Ricordando che p = k 2 / = a(1 e2 ), risulta
E =
E =
2a
6.6
2 2r2 r
Grazie alla conservazione del momento angolare, in tale espressione non
Si consideri fissato
compaiono ne la variabile angolare ne la sua derivata .
il valore del momento angolare k (e dunque quello del parametro orbitale
p = k 2 /). I termini dipendenti dalla sola r costituiscono una energia
potenziale equivalente,
V (r) =
k2
2r2 r
58
E>0
E=0
E<0
Emin
rA
rc
rP
r / rc
Orbite circolari
Lenergia totale associata al minimo di V `e pari a E = 2 /(2k 2 ), da cui
si evince che a = k 2 / = p, ragion per cui leccentricit`a `e nulla e lorbita `e
circolare. Dunque lorbita circolare rappresenta lorbita di minima energia
totale E per un dato valore del momento angolare k.
Essendo V = E deve risultare r = 0, il che conferma che lorbita `e circolare, in quanto la derivata del raggio `e identicamente nulla e dunque r
`e costante. La velocit`
a a cui viene percorsa tale orbita `e detta velocit`
a
circolare ed `e pari a
vc = /r0 = /k.
Orbite ellittiche
Per 2 /(2k 2 ) < E < 0 le rette E = cost intersecano la curva V (r) in due
punti, corrispondenti ai valori del periasse rP e dellapoasse rA .
=0
2
rP
k2
=0
2rP2 rP
Orbite iperboliche
Per E > 0 si ha una sola intersezione per r = rp , per cui la traiettoria risulta
di nuovo aperta, come nel caso precedente e il comportamento `e analogo a
quello osservato per traiettorie paraboliche. Questa volta per`o, come gi`a osservato, si ha un eccesso iperbolico di velocit`a con cui m2 percorre lasintoto
uscente delliperbole.7
6
60
e`
1 + e cos
P
d
F
l
p
1 + e cos
circonferenza;
ellisse;
parabola;
iperbole.
dove le prime due sono curve chiuse,8 mentre le ultime due sono aperte.9
La distanza minima dal fuoco la si ottiene sempre per = 0 ed `e pari a
rmin = p/(1 + e).
Nel caso dellellissi, la distanza massima la si ottiene per = ed `e pari
a rmax = p/(1 e). Il semiasse maggiore dellellisse `e dunque pari a
a = (rmin + rmax )/2
valida solo per le ellissi. Nel caso pi`
u generale si ha che
a = p/(1 e2 )
da cui si vede chiaramente che per e 1 si ha a , mentre per e > 1
il semiasse maggiore a diviene negativo (e perde pertanto il suo significato
geometrico).
Geometria dellellissi
Preso il centro dellellisse, C, posto allintersezione fra gli assi maggiore e
minore, si consideri il secondo fuoco G come il punto simmetrico rispetto a
F posto sullasse maggiore dalla parte opposta di C. La distanza fra i due
fuochi F e G `e pari a
2ep
2c = rmax rmin =
= 2ea.
1 e2
(ovvero la distanza fra uno dei due fuochi e il centro C `e pari a c = ea).
Lellisse pu`
o essere definita anche come il luogo dei punti per cui `e
` facile evidenziare per i punti
costante la somma delle distanze dai fuochi. E
posti alle estremit`
a dellasse maggiore come tale somma sia pari a 2a, risultato che deve valere dunque anche per i punti posti allestremit`a dellasse
`
minore. Il semiasse minore dunque
dato da
b = a2 c2 = a 1 e2 .
Tenuto contro del fatto che larea dellellisse A = ab si ottiene per A
lespressione
A = a2 1 e2 = a ap.
8
62
Capitolo 7
Premessa
7.2
1 , e
2 , e
3 } e FB =
Si considerino due terne cartesiane ortogonali FA = {O; e
1, g
2, g
3 }, dove si `e assunto, per semplicit`a, che abbiano la stessa origine
{O; g
O. Nei richiami iniziali si `e gi`a mostrato come sia possibile associare a
uno stesso vettore diverse scomposizioni per componenti in terne diverse e
come tali scomposizioni siano legate da una matrice di trasformazione di
coordinate, T BA . La matrice T BA , che `e ortogonale e gode pertanto di
alcune propriet`
a notevoli, pu`o essere costruita utilizzando le componenti dei
versori di una terna espressi per componenti nellaltro riferimento.
Si osservi che le componenti dei versori di una terna FB espressi in nel
riferimento cartesiano assunto come fisso o assoluto FA rappresentano
tutta linformazione necessaria per conoscere la posizione di FB rispetto a
FA , definita assetto di FB . Questo richiede per`o la conoscenza dei 9 elementi della matrice T BA , coefficienti fra loro non indipendenti, date le 6
7.2.1
Una sequenza di tre rotazioni elementari porta alla definizione di due terne
j della
3 }, tali che lasse e
2 , e
1 , e
3 } e F2 = {O; e
2 , e
1 , e
intermedie, F1 = {O; e
terna di partenza intorno al quale si effettua la prima rotazione coincide con
j della prima terna intermedia. Dopo le prime due
il corrispondente asse e
i della seconda terna intermedia deve coincidere
rotazioni, uno dei tre assi e
i della terna di arrivo, in modo che gli assi e
j e e
k (con j, k
con lasse g
i . In questo modo la terza
i g
i) giacciano sul piano perpendicolare a e
i e pu`o completare la sovrapposizione
rotazione avverr`
a proprio intorno a g
delle due terne. Daltro canto `e noto che bastano due rotazioni per allineare
un asse con una direzione arbitraria e dunque 3 rotazioni sono sufficienti per
portare due riferimenti cartesiani a coincidere.
Nel seguito gli indici i, j e k potranno assumere ciascuno uno dei valori
1, 2 o 3. Come verr`
a discusso nel seguito in dettaglio, non `e necessario che
le rotazioni avvengano intorno ad assi coordinati tutti diversi, per cui i tre
` evidente,
indici i, j e k non devono essere necessariamente tutti diversi. E
per`o, che due rotazioni elementari successive non possono avvenire intorno
allo stesso asse, altrimenti le due rotazioni potrebbero essere sostituite da
ununica rotazione di ampiezza pari alla somma delle due.
i lungo g
i vi sono due possibilit`a. La prima consiste
Per puntare e
i
nelleffettuare due rotazioni corrispondenti allazimut e allelevazione di g
i . La proiezione di g
i
rispetto a uno dei piani coordinati ij che contiene e
j porta e
i su e
i g
i , dove j `e lelevazione di g
i rispetto al piano
intorno a e
7.2.2
Stante il vincolo che due rotazioni successive non possono aver luogo intorno
ad un medesimo asse, `e possibile definire le 12 sequenze che consentono la
sovrapposizione della terna FA con una qualunque terna FB . Le sequenze
sono identificate dagli indici i, j e k che indicano la successione degli assi
attorno ai quali si effettuano le tre rotazioni elementari:
1-2-3
1-3-2
2-1-3
2-3-1
3-1-2
3 - 2 - 1 (angoli di Cardano o di Bryan)
1-2-1
1-3-1
2-1-2
2-3-2
3 - 1 - 3 (angoli di Eulero classici)
3-2-3
7.2.3
Sequenze famose
e3
g3
g2
e2
`
e1
e1' = n
e2'
g1
e2
,
3 `e parallelo al momento angolare h = r v
vettore posizione r , il terzo asse g
2 completa una terna destra.
eg
In questo caso la direzione dellintersezione fra i piani perpendicolari a
3 e g
3 , data da n
= (
3 )/
3 , `e detta linea dei nodi. Langolo
e
e3 g
e3 g
ed e
1 `e detto ascensione retta del nodo ascendente. Langolo i fra
fra n
3 e g
3 `e detto inclinazione (in quanto rappresenta linclinazione del piano
e
3 , rispetto al piano coordinato, perpendicolare
orbitale, perpendicolare a g
3 ), e infine lanomalia = + `e data dalla somma fra la distanza
a e
e la direzione del periasse, detta argomento del periasse,
angolare fra n
e lanomalia vera , che determina la posizione del veicolo spaziale lungo
lorbita, contata a partire dal periasse.
Angoli di Cardano (o di Bryan)
Gli angoli di Cardano (noti nella letteratura anglosassone come angoli di
Bryan) sono molto utilizzati in meccanica del volo e in ingegneria navale per
indicare lassetto di una terna associata a un veicolo (sia esso un natante
o un aeromobile) rispetto a una terna di assi verticali locali, in cui lasse
3 `e parallelo alla verticale locale e glia assi e
1 ed e
2
della terna fissa e
e3=e3'
e3
g3
g2
`
e1
e1'
e2'=e2
e2
g1
Fig. 7.2: Angoli di Cardano.
giacciono quindi nel piano orizzontale. Le rotazioni associate alla sequenza
3-2-1 (Fig. 7.2) sono le seguenti:
, detta imbardata, corrispondente alla prima rotazione intorno alla
3 , che porta lasse e
1 lungo la proiezione e
1 dellasse
verticale locale e
1 sul piano orizzontale;
longitudinale del mezzo g
2 e
, detta angolo di beccheggio, che viene effettuata intorno allasse e
7.2.4
Commenti
Una volta scelta la sequenza da utilizzare, i tre angoli associati alle tre
rotazioni elementari permettono di identificare la posizione della terna mobile rispetto a quella assunta come fissa. Levoluzione temporale di tali
angoli descrive quindi lassetto della terna mobile in funzione del tempo. Si
osservi per`
o che la sequenza di rotazioni utilizzata per descrivere lassetto
non costituisce di per s`e la rappresentazione di alcun moto. Le tre rotazioni
elementari vengono descritte in modo sequenziale solo per fini didattici, al
fine di evidenziare il significato geometrico di ciascuna rotazione elementare,
7.2.5
Singolarit`
a
7.3
7.3.1
y1
os
y 1c
os
x 1c
y2
x2
in
y 1s
x1
sin
-x 1
Fig. 7.3: Trasformazione di coordinate per una rotazione elementare intorno allasse z.
R3 () = sin cos 0
0
0
1
1
0
R1 () = 0 cos sin ; R2 () = 0
0 sin cos
sin 0 cos
7.3.2
Composizione
La matrice di trasformazione di coordinate pu`o ora essere espressa combinando le diverse matrici associate alle rotazioni elementari. Consideriamo i
due casi discussi prima, relativi agli angoli di Eulero e agli angoli di Cardano.
Matrice di trasformazione di coordinate per gli angoli di Eulero
Siano r A = (xA , yA , zA )T le componenti del vettore r in FA . Le componenti
nella prima terna intermedia sono date da
r = R3 ()r A
Procedendo in modo analogo si ha che le componenti nella seconda terna
intermedia sono date da
r = R1 ()r = R1 () [R3 ()r A ]
e quelle nella terna FB sono quindi pari a
r B = R3 ()r = R3 () {R1 () [R3 ()r A ]}
= [R3 ()R1 ()R3 ()] r A
dove lultimo passaggio `e reso possibile dalla propriet`a associativa della
moltiplicazione righe per colonne fra matrici. Dallultima relazione `e possibile ricavare che
T BA = R3 ()R1 ()R3 ()
cos cos
sin cos cos sin sin
sin cos
sin sin
sin sin
sin cos
cos
cos cos
cos sin
sin
+ cos cos
T BA = cos sin
+ sin sin
sin cos
7.3.3
Trasformazioni inverse
Le matrici di trasformazione delle coordinate risultano univocamente definite dallassetto relativo delle due terne FA e FB , ma sono (ovviamente)
indipendenti dalla particolare sequenza di rotazioni elementari scelta. Una
volta scelta la sequenza `e possibile valutare le ampiezze delle tre rotazioni elementari a partire da funzioni goniometriche inverse di alcuni elementi della
matrice o di loro rapporti.
In particolare, indicati come ti,j gli elementi della matrice T BA , per gli
angoli di Eulero risultano valide le seguenti relazioni:
= acos(t3,3 )
= atan(t3,1 , t3,2 )
= atan(t1,3 , t2,3 )
7.3.4
Propriet`
a e problemi
T BA
0
0
1
cos( + ) sin( + ) 0
(sin + ) cos( + ) 0 = R3 ( + )
0
0
1
cos cos sin sin
T BA
0
0
1
sin( ) cos( )
0
cos( ) sin( ) 0
0
T BA = 1
T BA = I3
4
La noncommutativit`
a delle rotazioni pu`
o essere dimostrata anche in modo pratico,
ma efficace, prendendo un libro chiuso poggiato sul tavolo e considerando due rotazioni
di 90 intorno a un asse perpendicolare alla copertina e un asse parallelo al dorso. Se
si effettua prima la rotazione intorno allasse perpendicolare alla copertina e poi quella
intorno al dorso, il libro finir`
a con il dorso poggiato sul tavolo (o il lato dal quale si aprono
le pagine). Se loperazione viene ripetuta ruotando prima intorno al dorso e poi intorno
allasse perpendicolare alla copertina il libro alla fine `e poggiato sul piano del tavolo come
fosse su uno scaffale. Quindi le rotazioni non commutano!
5
Se ne lascia la verifica al lettore.
0 1
= 3
2 1
TBA
= TBA
= I3 +
1
+ 0
1
T BA = ( + )
Capitolo 8
La velocit`
a angolare
T B0 B = T 1
BB0 = T BB0 = I3 +
t0 )r B
(t
d
r
B
) = lim
= r
dt B0 t0
t
dove
0 3 2
= 3
0
2 1
0
8.1.1
`
Ma
e veramente un vettore?
(t t )r
d
r
0 B
r
) = lim
=
C
dt C0 t0
t
B r B = T B C (
C r C ) = (T B C
T C B ) rB
vB =
C
0 0
0 0
0 0
B = TB C
T C B , ovvero B = T B C
C . Ma allora le
Dunque
C
0 0
0 0
0 0
componenti B e C sono collegate da una relazione di trasformazione di
coordinate, dunque esprimono le componenti di un vettore
= 1 g
1 + 2 g
2 + 3 g
3
` RELATIVA
8.2. VELOCITA
8.1.2
79
8.1.3
Formule di Poisson
Se consideriamo come vettore r fisso nel riferimento mobile uno dei suoi
versori, si ha
d
g
( i ) = B (
g i )B
dt B
ma allora
d
g1
g
1 = 3 g
2 2 g
3
=
dt
d
g2
g
2 = 1 g
3 3 g
1
=
dt
d
g3
g
1 1 g
2
3 = 2 g
=
dt
Tali formule, note come formule di Poisson, sono alla base della
definizione delle propriet`
a della cinematica relativa.
8.2
Velocit`
a relativa
r = r O +
e
dt
Proiettando le componenti di tutti i vettori in FB si ha che
=v
O +
v
v B = T BA v OA + B B + B
8.3
Accelerazione relativa
d2
2
dt
ottenendo
B
aB = T BA aOA + B B + B ( B B ) + 2 B B +
(8.1)
dove T BA aOA + B B + B ( B B ) `e laccelerazione di trascinamento (ovvero laccelerazione assoluta del punto del riferimento mobile
occupato da P allistante considerato), 2 B B `e detta accelerazione di
B `e laccelerazione relativa.
Coriolis, mentre
8.4
Capitolo 9
9.1
9.1.1
Coordinate geografiche
Sistemi di riferimento geocentrici
Si consideri una terna di assi solidali alla Terra, assunta per semplicit`a sferica di raggio R , avente come origine O il centro della sfera terrestre, con
lasse z diretto lungo lasse di rotazione della terra (detto asse polare),
, dunque positivo verso
parallelamente alla velocit`
a angolare di rotazione
83
9.1.2
Coordinate geografiche
Un riferimento avente orientamento fisso rispetto alle stelle con origine nel centro
di massa del Sole approssima ancora meglio un riferimento inerziale, in quanto lorigine
della terna geocentrica equatoriale ruota intorno al Sole e dunque non si muove di moto
rettilineo uniforme. Il riferimento eliocentrico, con assi x e y sul piano delleclittica
(scegliendo lasse x parallelo alla linea dellequinozio vernale) e asse z perpendicolare allo
stesso piano viene utilizzato per descrivere il moto di pianeti e comete e di sonde spaziali
interplanetarie rispetto al Sole. Tuttavia, la velocit`
a angolare del moto di rivoluzione
della Terra `e talmente bassa che le accelerazioni apparenti ad esso dovute sono, in genere,
modeste e, dunque, trascurate.
N
r
meridiano di
Greenwich
piano
equatoriale
yR
xR
S
9.2
9.2.1
a
lungo i versori di FT si ottiene
intorno allasse polare. Proiettando
= ( cos )t1 + ( sin )t3
zR
cos
zT
xT
yT
sin
xT
yR
xR
S
9.2.2
(9.1)
9.2.3
Rcos
g app
xT
9.2.4
x
= 2 y sin
y = 2 (x sin + z cos )
z = g + 2 y cos
` possibile integrare ambo i membri della prima e della terza equazione,
E
ottenendo
x(t)
x 0 = 2 y sin
z(t)
z0 = gt + 2 y cos
Esplicitando i valori di x(t)
e z(t)
= x 0 + 2 y 0 t sin
z(t)
= z0 gt + 2 y 0 t cos
e dunque
x(t) = x 0 t + y 0 t2 sin
1
z(t) = z0 t gt2 + y 0 t2 cos
2
!"
"
!"
!"
"
!
!"
2 2
h
xF = zF = 0 ; yF =
h
cos > 0
3
g
che risulta quindi spostato verso Est rispetto a Q. Questo indica che un grave
cade lungo una traiettoria che presenta una deviazione verso Est, legata al
verso della forza (apparente) di Coriolis per una componente verticale di
velocit`a.
9.2.5
Cicloni e anticicloni
Il risultato sopra riportato `e difficilmente verificabile sperimentalmente, anche a causa delle numerose ipotesi semplificative alla base della soluzione
analitica sopra riportata. Esiste tuttavia in natura un fenomeno che dimostra in modo inconfutabile la rilevanza della forza di Coriolis, ed `e il
verso di rotazione delle correnti cicloniche e anticicloniche nei due emisferi
Nord e Sud.
Si considerino delle zone di alta (A) e bassa (B) pressione,2 disposte come
in Fig. 9.4. In assenza di rotazione terrestre, le masse daria sottoposte a una
pressione maggiore tenderebbero a fluire verso le zone di pressione inferiore
(frecce tratteggiate). Nellemisfero boreale la componente di velocit`a ango lungo la verticale locale t3 risulta positiva (diretta verso
lare terrestre
lalto), e provoca una forza di Coriiolis orientata nel verso indicato dalle
2
Lorigine di tali differenze di pressione non ha interesse per la presente discussione e
viene quindi omessa.
9.3
Pendolo di Foucault
Si consideri un pendolo semplice, di massa m, appeso con un cavo inestensibile di massa trascurabile e lunghezza pari a L sulla verticale di un punto Q
sulla superficie terrestre, alla latitudine . Sia Q lorigine di un riferimento
topocentrico FT = {Q; t1 , t2 , t3 }, con i versori disposti in direzione Sud, Est,
e Zenith.
Le forze agenti sul corpo appeso al cavo sono il suo peso mgt3 e la
tensione T nel cavo, espressa nel riferimento topocentrico dalle componenti
T T = (fx , fy , fz )T = T (x/L, y/L, 1 z/L)T . Le equazioni del moto sono
dunque rappresentate tramite il seguente sistema di equazioni differenziali
ordinarie del secondo ordine
T x
+ 2 y sin + 2 R sin cos
mL
T y
y =
2 (x sin + z cos )
mL
T
z
z =
(1 ) g + 2 y cos + 2 R cos2
m
L
x
=
91
variabile complessa u = x + iy, dove i = 1 `e lunit`a immaginaria. Moltiplicando la seconda equazione per i e sommando si ottiene lequazione differenziale lineare del secondo ordine a coefficienti costanti
u
+ 2aiu + p2 u = 0
(9.2)
1,2 = ia
a2 p2 ia ip,
9.4
93
2
366.25636
=
Si definisce crosstrack lo spostamento in direzione laterale rispetto al piano diametrale che contiene la velocit`
a di rientro allistante dellingresso in atmosfera.
6
Date le fluttuazioni nella durata del giorno solare, il secondo ha ora una definizione
diversa, legata al periodo di oscillazione di un fenomeno su scala atomica. Pi`
u esattamente,
il secondo `e definito come la durata di 9 192 631 770 periodi della radiazione corrispondente
alla transizione tra due livelli iperfini, da (F=4, MF=0) a (F=3, MF=0), dello stato
fondamentale dellatomo di cesio-133, qualunque cosa questo voglia dire...
Parte II
Meccanica Lagrangiana
95
Capitolo 10
primo volume del suo testo fondamentale1 sotto lambizioso titolo Equation
differentielle pour la solution de tous les probl`emes de Dinamique. Esso
si basa sulla definizione di due sole quantit`a scalari per lintero sistema,
lenergia cinetica e il lavoro, e consente di ricavare le equazioni che regolano la dinamica del sistema in modo pressoche automatico, una volta
nota la cinematica del sistema stesso. Le forze vincolari, se richiesto, possono essere ricavate in un secondo momento.
Come ulteriore vantaggio, le equazioni di Lagrange possono essere ricavate non solo in termini di coordinate spaziali, ma anche in funzione delle
cosiddette coordinate generalizzate, ovvero di un sistema di parametri
che tenga conto, nella descrizione delle configurazioni del sistema, della
presenza di vincoli che ne condizionano il moto. In questo modo le equazioni
del moto di un generico sistema meccanico vengono ottenute da propriet`a
globali del sistema stesso, e non gi`a dallanalisi individuale delle grandezze
vettoriali che caratterizzano il moto dei suoi componenti (posizione, velocit`a e accelerazione) e delle interazioni fra di esse o con lambiente (forze),
espresse per componenti in un particolare sistema di riferimento cartesiano
ortogonale.
1
J.L. Lagrange, Mcanique Analytique, seconda ed. riveduta e corretta, Parigi, 1811,
che segu` la prima edizione, stampata a Berlino nel 1788.
97
98
`
CAPITOLO 10. SISTEMI VINCOLATI E GRADI DI LIBERTA
La scrittura delle equazioni effettuata in termini di coordinate generalizzate (ovvero quantit`a diverse dalle coordinate vettoriali fisiche della
meccanica newtoniana) porta cos` due vantaggi significativi: una maggiore
versatilit`
a nella scelta delle descrizione matematica del sistema ed una formulazione indipendente dal particolare sistema di riferimento utilizzato per
descrivere il moto. Tali aspetti risultano particolarmente importanti quando
si prendono in esame sistemi con molti gradi di libert`a e/o con vincoli anche
tempovarianti.
10.1
Configurazioni
Si consideri un sistema di N elementi materiali la cui posizione Pi `e rappresentata dai vettori r i , i = 1, 2, ..., N , dove la posizione della iesimo elemento
`e data da r F i = (xi , yi , zi )T in un opportuno riferimento cartesiano F. Il
moto del sistema `e noto se `e possibile definire istante per istante la posizione di ciascun elemento, ovvero se sono note le 3N funzioni del tempo
xi (t), yi (t), zi (t), i = 1, 2, ..., N . Linsieme delle posizioni di tutti gli elementi
del sistema `e detta configurazione, e pu`o essere rappresentata come un
punto P di uno spazio a 3N dimensioni, detto appunto spazio delle configurazioni S. Levoluzione di P in S rappresenta quindi levoluzione del
sistema, ovvero il suo moto, e la curva P (t) S `e detta curva C, ovvero
C = {P = P (t) S, t R} S.
10.2
Vincoli
10.2. VINCOLI
10.2.1
99
(10.1)
dove r i sono le coordinate dei punti occupati dagli elementi del sistema
allistante t in un certo riferimento cartesiano ortogonale.
Se la funzione che esprime il vincolo non dipende espressamente dal
tempo (ovvero t non compare come parametro indipendente fra gli argomenti
di f in Eq. (10.1)) i vincoli sono detti scleronomi. Se invece la funzione
che esprime il vincolo dipende esplicitamente dal tempo, in quanto la natura
del vincolo `e tempovariante, il vincolo `e detto reonomo.
I vincoli che non dipendono solo dalla posizione degli elementi del sistema, e non sono dunque esprimibili nella forma data da Eq. (10.1), sono
detti anolonomi.
10.2.2
Esempi
100
`
CAPITOLO 10. SISTEMI VINCOLATI E GRADI DI LIBERTA
10.3
Gradi di libert`
a
Se i vincoli sono tutti olonomi, il sistema `e detto olonomo. Basta che un solo
vincolo sia anolonomo, e lintero sistema non `e pi`
u olonomo. In alcuni casi,
`
10.3. GRADI DI LIBERTA
101
102
`
CAPITOLO 10. SISTEMI VINCOLATI E GRADI DI LIBERTA
10.4
Coordinate generalizzate
= zN (q1 , q2 , ..., qn )
103
104
`
CAPITOLO 10. SISTEMI VINCOLATI E GRADI DI LIBERTA
Capitolo 11
Estremo di un Funzionale
I =
F (x, y, y )dx
x1
11.2
Approccio analitico
106
I() =
x1
F dy F dy
+
) dx
y d y d
x1
x2 F
F
=
(
+ ) dx
y
y
x1
x2
x1
x2
F d
F x2
d F
dx
=
[
]
(
) dx
y dx
y x1
dx y
x1
dove il termine fra parentesi quadre `e nullo per le condizioni agli estremi su
. Sostituendo tale risultato nellespressione di dI/d si ottiene lespressione
x2 F
dI()
d F
=
[
(
)] dx.
d
y dx y
x1
Dal momento che (x) `e una funziona arbitraria (fatte salve le condizioni agli estremi), dI/d pu`o essere nulla solo se la funzione y(x) soddisfa
lequazione differenziale di EuleroLagrange
F
d F
(
)=0
y dx y
condizione necessaria affinche I abbia un valore stazionario.
107
11.3
Approccio variazionale
d
(y),
dx
d
y
dy
d
d
dy
)=( )( )=
(
y y) =
(y) = (y) .
dx
dx
dx
dx
dx
I =
x1
F
F
y + y ) dx
y
y
I =
[
x1
F
d F
(
)] ydx
y dx y
108
11.4
Lidentit`
a di Beltrami
y
y
dx x y
Ma, moltiplicando lequazione di EuleroLagrange per y , risulta anche
d F
F
y = y
(
)
y
dx y
Sottrando questultima relazione dalla precedente, risulta
d F
dF F F
y y
(
)=0
dx x y
dx y
Tale equazione pu`
o essere riscritta nella forma
F
d
F
+
(F y ) = 0
x dx
y
F
y =C
y
11.5
Storicamente, il primo problema ad essere stato affrontato e risolto adottando il calcolo delle variazioni `e il problema della brachistocrona, ovvero
la determinazione di una curva y = f (x) in un piano cartesiano tale che
una massa m che percorra la curva senza attrito, sotto la sola azione di un
campo di forza costante diretto lungo lasse y impieghi un tempo minimo a
passare la punto P1 = (x1 , y1 ) al punto P2 = (x2 , y2 ).
t2
T =
dt =
t1
P1
x2 dx
ds
=
v
x
x1
v = 2gy
Ma, dette x = dx/dt e y = dy/dt le componenti del vettore velocit`a nel riferimento
cartesiano scelto, si ottiene anche la seguente catena di uguaglianze
v 2 = x 2 + y 2 = x 2 + (
dy 2
x)
= x 2 (1 + y 2 )
dx
2gy/(1 + y 2 )
T =
(1 + y 2 )/(2gy) dx
x1
Dal momento che lintegrando F (x, y, y ) = (1 + y 2 )/(2gy) non dipende esplicitamente da x, `e possibile applicare lidentit`a di Beltrami F y F /y = C.
Sostituendo lespressione di F e
F
1 + y 2
y
=
=
y y
2gy
2gy(1 + y 2 )
nellidentit`
a di Beltrami, si ottiene la relazione
110
1 + y 2
y 2
1
=
=C
2
2gy
2gy(1 + y )
2gy(1 + y 2 )
dy 2
1
) ]=
= k2 > 0
dx
2gC 2
k2
k2
( sin ) ; y() =
(1 cos )
2
2
Capitolo 12
Principi fondamentali
12.1
R
W = f i
ri
i=1
W = f i
ri
i=1
112
il lavoro fatto dalle forze applicate al sistema per spostamenti virtuali infinitesimi nellintorno di una condizione di equilibrio `e nullo.
Ci`
o permette di conferire alla condizione di equilibrio del sistema una
nuova propriet`
a, in quanto, in presenza di vincoli, gli spostamenti virtuali
W = Qk qk
k=1
12.2
Il Principio di dAlembert
Il principio dei lavori virtuali permette di determinare la condizione di equilibrio statico di un sistema, e pertanto non `e adatto alla descrizione della
dinamica di sistemi il cui stato evolva nel tempo. Tuttavia esso pu`o essere
esteso anche allanalisi dinamica tramite il principio di dAlembert.
Secondo il principio di dAlembert, `e possibile considerare la variazione
della quantit`
a di moto di una massa mi come lopposto di una forza dinerzia
I
= mi v
i `e la quantit`a di moto di mi . Lespressione del sec i , dove p
f i = p
ondo principio della dinamica
V
i
f i + f i = p
pu`
o essere quindi riscritta in modo formalmente identico a quella di un
equilibrio statico portando a primo membro la forza dinerzia, ovvero
V
I
f i + f i + f i = 0
113
i )
ri = 0
(f i p
i=1
12.3
Il Principio di Hamilton
ri = 0
(f i mi r i )
i=1
dove
r i = W
f i
i=1
mi r i
r i = mi
114
che pu`
o essere sostituita nellespressione del principio di dAlembert generalizzato, che viene cos` espresso in una nuova forma
N
T + W = mi
i=1
dove
d
(r i
ri )
dt
N
1
T = mi r i r i
i=1 2
t1
t2
(T + W) dt =
t1
t2
N
d
[ mi (r i
r i )] dt = [ mi (r i
r i )
dt
i=1
i=1
t
(T + W) dt = 0
t1
Ldt = 0
t1
Ldt) = I = 0
t1
115
I =
Ldt
t1
12.4
T =
1 L
w(x, t)
] dx
m(x) [
2 0
t
L
0
dx
cos
1
1
1
1
1 w(x, t)
=
=
1+ [
]
2
2
2
cos
2
x
1 tan
w(x,t)
(1 sin )
1 [ x ]
` dunque possibile scrivere
E
116
V = F L =
L w(x, t)
F
] dx
[
2 0
x
Dal momento che il sistema `e olonomo, `e possibile applicare il Principio di Hamilton: la funzione w(x, t) rappresenta il moto della corda se rende stazionario il
funzionale
L
w 2
w 2
1 t2
[m
(
{
)
F
(
) ] dx} dt
2 t1
t
x
0
t2
I =
(T V) dt =
t1
t1
1
w 2
m(
) dt)
2
t
t2
m
t1
w w
(
) dt
t
t
t2
t2
w
w
) w
(m
) wdt
t
t
t
t
1
t1
t2
2w
= m 2 wdt
t
t1
= [(m
1
w 2
F(
) dx)
2
x
F
0
= [(F
L
w w
w
(
) dx = F
(w)dx
x
x
x x
0
L
L
w
2w
) w F
wdt
x
x2
0
0
t2
t1
(T V) dt) =
t1
(F
0
L
2w
2w
w
m
)
wdt
[(F
)
w
} dt = 0
x2
t2
x
0
L
w
) w = 0
x
0
2w
2w
m
=0
x2
t2
Capitolo 13
Equazioni di Lagrange
13.1
r i dqn
=
+
+ ... +
dt
q1 dt q2 dt
qn dt
n
ri
=
qk ,
i = 1, 2, ..., N
k=1 qk
r i =
1 N
mi r i r i
2 i=1
n n
1 n n
1 N
r i
ri
qj qk = Mj,k qj qk
mi
2 i=1 j=1 k=1 qj qk
2 j=1 k=1
Mj,k = Mk,j = mi
i=1
117
r i
ri
qj qk
118
I =
(T V) dt
t1
I = (
t1
Dal momento che le coordinate generalizzate sono state assunte indipendenti, il problema di annullare I pu`o essere risolto imponendo che sia nulla
la variazione indotta dagli incrementi di ciascuna variabile, utilizzando un
approccio simile a quello presentato nel par. 11.1. Pertanto I ha un valore stazionario se le coordinate generalizzate qk soddisfano le Equazioni di
Lagrange
L
d L
(
)
= 0, k = 1, 2, ..., n
dt qk
qk
Se il sistema `e sottoposto ad azioni il cui lavoro non `e riducibile alla
variazione di una quantit`a scalare, occorre utilizzare il Principio di Hamilton
esteso, nella formulazione
t2
(T + W) dt = 0
t1
In questo caso il lavoro virtuale pu`o essere scisso in due contributi, uno,
V, proveniente dalle forze conservative, laltro, W , dovuto allazione di
forze non conservative. Questultimo, definito come
W = f i
ri
N
i=1
pu`
o essere riscritto in funzione della variazione virtuale delle coordinate
generalizzate, in quanto, posto
n
ri
qk
q
k
k=1
ri =
risulta
n
N
n
N
ri
ri
W = f i (
qk ) = ( f i
) qk
qk
i=1
k=1 qk
k=1 i=1
119
W = Qk qk
k=1
le quantit`
a
N
ri
) , k = 1, 2, ..., n
Qk = ( f i
qk
i=1
t1
t2
(T V) dt +
t1
Qk qk dt = 0
Da tale relazione viene ricavato, applicando nuovamente il calcolo delle variazioni per ciascuna variazione indipendente qk , il sistema di Equazioni di
Lagrange nella forma
d L
L
(
)
= Qk , k = 1, 2, ..., n
dt qk
qk
= L1 sin 1
y1
= L1 (1 cos 1 )
x2
= L1 sin 1 + L2 sin 2
y1
= L1 (1 cos 1 ) + L2 (1 cos 2 )
=
=
1
1
m1 (x 21 + y 12 ) + m2 (x 22 + y 22 )
2
2
1
1
2 2
m1 L1 1 + m2 (L21 12 + L22 22 + 2L1 L2 1 2 cos(1 2 ))
2
2
120
= m1 gy1 + m2 gy2
= m1 gL1 (1 cos 1 ) + m2 g (L1 (1 cos 1 ) + L2 (1 cos 2 ))
Nel sistema sono presenti solo forze conservative (la forza peso) il cui potenziale dipende dalla sola posizione e pertanto lequazione di Lagrange assume la
formulazione
T V
d T
(
)
+
=0
dt q
q q
dove
T
V
d T
(
)
+
dt 1
1 1
d T
T
V
(
)
+
dt 2
2 2
= m2 L22 2 + m2 L1 L2 1 cos(1 2 ) +
m2 L1 L2 12 sin(1 2 ) + m2 gL2 sin 2 = 0
= 0
= 0
Il sistema di due equazioni differenziali ordinarie lineari del secondo ordine cos`
ottenuto pu`
o essere scritto in forma matriciale,
[
(m1 + m2 )L21
m2 L1 L2
13.2
m2 L1 L2
(m1 + m2 )gL1
] ( 1 )+[
2
m2 L2
0
2
0
]( 1 ) = (
)
m2 gL2
2
0
121
assume in questo caso che la dissipazione per lelemento iesimo del sistema sia fuzione solo ed esclusivamente del suo stato di moto, e non dipende
quindi dal moto degli altri elementi. Sotto tale ipotesi di disaccoppiamento,
il lavoro virtuale prodotto dalla forze dissipative `e dato da
N
W d = f i
r i = (cxi x i xi + cyi y i yi + czi zi zi )
i=1
i=1
ri
ri =
qk
q
k
k=1
mentre la velocit`
a delliesimo elemento r i risulta pari a
n
ri
r i =
qk
k=1 qk
ri =
1
(cxi x 2i + cyi y i2 + czi zi2 )] qk
i=1 2 qk
= [
k=1
1
(cxi x 2i + cyi y i2 + czi zi2 )
i=1 2
F =
`e possibile scrivere
n
F
qk = Qdk qk
k=1
k=1 qk
W d =
F
.
qk
122
1 n n
1 T
Cj,k qj qk = q C q
2 j=1 k=1
2
Pertanto, in presenza di forze dissipative, le equazioni di Lagrange possono essere riscritte nella forma
d L
L F
(
)
+
= Qk , k = 1, 2, ..., n
dt qk
qk qk
dove il termine relativo alle forze generalizzate Qk include solo le forze non
conservative e nondissipative.
= 0
m2 x
2 c2 x 1 + c2 x 2 k2 x1 + k2 x2
= 0
m1 + m2
0
0
x
c +c
] ( 1 )+[ 1 2
m2
x
2
c2
c2
x
k + k2
] ( 1 )+[ 1
c2
x 2
k2
123
k2
x
0
]( 1 ) = ( )
k2
x2
0
124
Parte III
125
Capitolo 14
128
14.1
Gradi di libert`
a di un corpo rigido
14.2
14.2.1
Terna solidale
14.2.2
14.3
Moti rigidi
14.3.1
Dal momento che il moto del corpo rigido `e rappresentato dal moto della
terna FB ad esso solidale, possiamo definire la velocit`a di ogni punto del
corpo rigido a partire da quanto abbiamo stabilito per i moti relativi di una
terna mobile FB rispetto a una fissa FI . La posizione P in FI di un
elemento di massa m del corpo rigido `e data da
.
r = r O +
1
130
14.3.2
Moto traslatorio
Il moto di un corpo rigido `e detto traslatorio se la velocit`a di tutti i suoi
punti `e uguale.
Teorema 2 Un moto rigido `e traslatorio se e solo se la velocit`
a angolare
del corpo `e costantemente nulla.
Dimostrazione: se la velocit`a angolare `e nulla, i versori della terna solidale conservano la loro direzione in quanto la loro derivata espressa dalla formula di Poisson si
e tutti i punti del corpo rigido hanno
annulla. Pertanto si annulla il termine
una velocit`
a pari a quella dellorigine.
Viceversa, se tutti i punti posseggono la stessa velocit`a, la derivata di tutti i
versori della terna solidale deve essere nulla, il che implica che sia nulla la velocit`a
angolare, QED.
Moto rototraslatorio
solidale
Il moto di un corpo rigido `e detto rototraslatorio se esiste un asse a
al corpo che mantiene un orientamento fisso.
Teorema 3 Un moto rigido `e rototraslatorio se e solo se la velocit`
a angolare
del corpo `e costantemente parallela allasse dei punti che si muovono con la
stessa velocit`
a.
(ci`o `e sempre
Dimostrazione: Si scelga lorigine del riferimento solidale O sullasse a
possibile data larbitrariet`a nella scelta del riferimento solidale). Se la velocit`a
`e diretta come a
=
= k
, essa `e esprimibile come
angolare
a. I punti
a hanno
velocit`
a
=v
O + (
O
v
a) (k
a) = v
non cambia direzione in quanto i suoi punti traslano nello spazio
e dunque lasse a
con velocit`
a pari a quella dellorigine.
131
Moto elicoidale
Il moto rototraslatorio di un corpo rigido `e detto elicoidale se la velocit`a
`e parallela allasse a
stesso.
dei punti dellasse a
Moto rotatorio
Il moto rototraslatorio di un corpo rigido `e detto rotatorio se la velocit`a dei
`e nulla. In questo caso lasse a
`e detto asse di rotazione.
punti dellasse a
Moto polare
Il moto di un corpo rigido `e detto polare se esiste un punto solidale al corpo
rigido che mantiene velocit`
a nulla.
Moto piano
Un moto rigido `e detto piano se esiste un piano solidale al corpo rigido
mantiene costanti la sua giacitura e la sua distanza d da un piano ad esso
parallelo nel riferimento fisso.
14.4
132
P e
Q , con
P
Q = k
P = v
O +
P =
individuata dai vettori
a, risulta v
O +
(
O +
Q = v
Q.
v
Q + k
a) = v
A questo punto `e possibile dimostrare che
Teorema 5 Un atto di moto `e rigido se e solo se `e rototraslatorio.
Q fra
P
Dimostrazione: Un atto di moto `e rigido se la distanza d = P Q =
due elementi di massa generici del corpo rimane costante. Questo `e possibile se la
`e identica, dove g
`e il versore parallelo
componente della velocit`
a di P e Q lungo g
P
Q = d
a
g . Se latto di moto `e rototraslatorio
P
v
g
P ) = g
O + g
(
(
v
= g
vO +
Q + d
g) =
O + g
Q + g
(d
Q) = g
Q
v
(
v
= g
g) = g
vO +
dP dQ
d
g
d
=d
=
= 0,
(P Q ) g
)g
g
dt
dt
dt
dt
il che `e possibile, per una scelta arbitraria dei punti P e Q, e dunque del versore
, solo se la derivata di g
`e espressa dalla formula di Poisson, il che comporta che
g
P = v
Q +
(d
v
g ), ovvero se latto di moto `e rototraslatorio, QED.
14.4.1
Teoremi di Mozzi
P = 12 (
v
O )+
Si consideri lasse del moto, definito dai punti P tali che
con k R uno scalare arbitrario. Lasse del moto `e pertanto parallelo
k ,
alla velocit`
a angolare. Il teorema di Mozzi afferma che
Teorema 7 Un atto di moto rigido `e elicoidale e lasse del moto `e lasse i
cui punti traslano con velocit`
a parallela alla direzione dello stesso asse.
e pertanto i suoi punti hanno la
Dimostrazione: Lasse del moto `e parallelo a
P v
Q = k ,
risulta
stessa velocit`
a in quanto, presi due punti P e Q, tali che v
P = v
Q +
(k )
=v
Q . Inoltre
v
P
v
O +
P = v
O +
[
= v
1
v
O ) + k ]
=
(
2
1
(
v
O )] + k
[
2
1
I
O + 2 [(
v
O )
(
)
vO ] = 2
= v
O +
= v
133
14.4.2
134
Capitolo 15
dh
= m
(15.2)
dt
dove si assume che il polo utilizzato per il calcolo dei momenti coincida
sempre con il centro di massa del sistema.
Nella prima equazione, la quantit`a di moto del corpo `e data semplicemente da
= m
p
vC
C = d
dove v
r C /dt `e la velocit`
a del centro di massa del corpo (v. Eq. (4.2)).
per i sistemi a massa costante, fornisce direttamente la derivata della velocit`a, d
v /dt = f /m, dove m `e la massa totale del corpo rigido. Integrandola
si ottiene levoluzione della velocit`a che, integrata a sua volta, fornisce la
posizione del centro di massa del sistema. Per poter sfruttare la seconda
equazione occorre invece definire in modo appropriato il valore del momento
angolare in funzione dellatto di moto del corpo rigido.
135
136
15.1
O =m
O +
O f
m
O = f (m
O +
O f ) = f m
O = IO .
IO = f m
15.2
A tale proposito si ricorda che, nella visione assiomatica della Meccanica del continuo
proposta da Walter Noll e Clifford Truesdell, uno degli assiomi fondamentali stabilisce che
il lavoro svolto da un sistema di forze su un corpo o sistema di corpi o elementi material
non dipende dal sistema di riferimento considerato. Pertanto le componenti delle forze
agenti sul sistema subiscono i soli effetti cinematici dellassetto di una terna rispetto a
unaltra per quanto riguarda la proiezione delle forze sugli assi coordinati, ma il valore
vettoriale della risultante non dipende dal punto considerato come origine del riferimento.
15.3
137
1
r m
m B
v
C = d
= r r C la posizione
dove v
r C /dt `e la velocit`
a del centro di massa e
dellelemento di massa m rispetto al centro di massa C.
Il momento angolare rispetto a C `e dunque dato da
C =
(m
(
)m
h
v) =
B
(m
m) v
C = 0
v C ) = (
B
relazione in cui le componenti della velocit`a angolare sono le stesse per tutto
il corpo rigido, e dunque sono delle costanti durante loperazione di integrazione su B. Il momento angolare pu`o quindi essere espresso nella forma
2 = I B .
Ixy = xym ,
B
Ixz = xzm ,
B
Iyz = yzm
B
138
15.4
Energia cinetica
1 2
1
) (
) m
v m = (
vC +
vC +
2 B
B 2
1
1
C ) m + v
C (
m) + (
) (
) m
(
vC v
2
2 B
B
B
=
=
In questultima espressione B m = m `e la massa totale del corpo, men m = 0, ragion per cui
tre, dalla definizione di centro di massa, risulta B
il secondo termine nella somma `e nullo. Il terzo addendo pu`o essere trasfor, v
e w,
mato grazie alle propriet`a del prodotto misto. Presi tre vettori u
infatti, si ha che
(
=v
(w
u
).
u
v w)
Grazie a questa relazione, il terzo termine della somma che fornisce T
pu`
o essere riscritto nella forma
1
1
1
(
) m =
) (
) m =
(I )
(
2 B
2
2
B
Si ottiene cos` un caso particolare del Teorema di Koenig (v. Par. 4.5)
sulla decomposizione dellenergia cinetica per il corpo rigido,
1
1
(I )
T = mv 2C +
2
2
in cui il primo termine `e lenergia cinetica associata alla traslazione del baricentro, mentre il secondo contributo `e quello associato al moto di rotazione
dello stesso corpo.
15.5
139
Dalla sua definizione si evince chiaramente cha la matrice di inerzia `e simmetrica. Si pu`
o dimostrare che una matrice simmetrica `e diagonalizzabile,
ovvero ammette una base di autovettori mutuamente ortogonali corrispondenti ad autovalori reali, che costituiscono gli elementi diagonali della matrice di inerzia, quando `e espressa nella base di autovettori, base in cui
si annullano tutti i termini misti extradiagonali.2 Esiste quindi una terna
solidale FP nella quale la matrice di inerzia ha lespressione
Jx 0 0
I = 0 Jy 0
0 0 Jz
dove gli autovalori di I, indicati come Jx , Jy e Jz , sono detti momenti principali di inerzia. Gli autovettori corrispondenti sono detti assi principali
di inerzia.
La matrice I `e anche definita positiva, stante il fatto che la forma
quadratica
TB I B = Trot 0
rappresenta lenergia cinetica rotazionale, quantit`a strettamente positiva, a
meno che la velocit`
a angolare non sia nulla. Di conseguenza tutti i momenti
principali di inerzia sono strettamente positivi.
In una terna di assi principali di inerzia le espressioni delle componenti
del momento angolare e dellenergia cinetica rotazionale diventano particolarmente semplici, in quanto
hP = (Jx 1 , Jy 2 , Jz 3 )T
e
1
Trot = (Jx 12 + Jy 22 + Jz 32 ).
2
Dividendo questultima espressione per Trot si ottiene lespressione dellellissoide
di inerzia (o ellissoide di Poinsot),
12
22
32
+
+
=1
2Trot /Jx 2Trot /Jy 2Trot /Jz
Quando due dei momenti principali di inerzia sono uguali, lellissoide
diviene assialsimmetrico, mentre se sono tutti e tre uguali esso `e una sfera.
Si osservi per`
o che le propriet`
a dellellissoide di inerzia sono s` legate alla
distribuzione di massa del corpo rigido considerato, ma non alla sua forma
geometrica. Ad esempio, un parallelepipedo a base triangolare ha un ellissoide di inerzia assialsimmetrico, pur non essendo un solido assialsimmetrico,
2
140
15.6
15.6.1
Simmetrie
Ci`
P Q parallelo a .
o implica che la densit`a nei punti P e Q sia la stessa.
`e perpendicolare al piano diametrale , `e detto piano
Se la direzione
di simmetria per la distribuzione di massa %(P ).
tale che ogni piano
Se esiste una retta r parallela alla direzione a
contenente r `e un piano di simmetria, allora r `e detta asse di simmetria
per la distribuzione di massa %(P ).
Se esiste un punto O tale che ogni retta rpassante per O `e un asse di
simmetria per la distribuzione di massa %(P ), la distribuzione `e a simmetria
polare rispetto a O.
15.6.2
Centro di massa
1
r m
m B
pu`
o essere semplificata da alcune semplici considerazioni, soprattutto in presenzza di simmetrie per il sistema (o corpo) considerato.
Il centro di massa di due corpi `e collocato sulla congiungente i due corpi,
se i corpi sono puntiformi, o i loro centri di massa, se sono corpi estesi (e.g.
corpi rigidi).
Se un sistema `e piano (o ammette un piano diametrale), il centro di
massa giace sul piano (o sul piano diametrale).
141
15.6.3
Momenti di inerzia
142
Osservando che
T
T
[( )I3 ] m = I
B
`e la matrice di inerzia per una terna baricentrica, e che tutti gli integrali del
tipo
T
T
[(C )] m = C ( m) = 0
B
e
T
T
[(C )] m = C ( m) = 0
B
si ha che
I = I + md2
ea
.
passante per il centro di massa e dC `e la distanza fra i due assi a
, il minimo momento di inerzia
Si osservi che, presa una certa direzione a
passa
lo si ottiene quando il corrispondente asse di rotazione parallelo ad a
per il centro di massa, ragion per cui il corpo presenta momenti di inerzia
minimi per assi di rotazione passanti per il suo centro di massa.
143
`e evidente che
Izz = (x2 + y 2 )m = x2 m + y 2 m = Ixx + Iyy .
B
Propriet`
a 3: La terna di assi principali di inerzia `e ruotata rispetto a una
generica terna di assi avente con 3 perpendicolare al piano xy di un angolo
che soddisfa lequazione tan(2) = 2Ixy /(Ixx Iyy ).
Dimostrazione: Come scritto precedentemente per una generica figura piana, la
posizione di un elemento di massa nella terna di assi principali di inerzia `e data
da
P = R3 () = (x cos + y sin , y cos x sin , 0)T .
Affinche la terna x y z sia principale di inerzia `e necessario che il termine misto
Ix y ()
1
= Ixy cos(2) (Ixx Iyy ) sin(2)
2
si annulli, il che accade per tan(2) = 2Ixy /(Ixx Iyy ).
Se lo spessore del corpo `e finito, come accade nei casi pratici, ma la distribuzione di
massa omogenea lungo z, ci si riconduce a un caso piano prendendo il piano mediano
fra le due superfici che delimitano il corpo, che `e piano di simmetria.
Appendice: Diagonalizzabilit`
a di una matrice reale
simmetrica
Definizione: si dice atuoaggiunta o Hermitiana una matrice A tale che
uT Av = (Au)T v
144
xH Ax
xH x
xH x
xH x
xH x
il che `e
ovvero lautovalore deve essere uguale al suo coniugato ,
possibile solo se la parte immaginaria `e nulla, ovvero se `e reale.
2. Se la matrice A `e reale e simmetrica i suoi autovettori sono mutuamente ortogonali.
Si considerino due autovalori distinti i j e i loro autovettori, xi e
xj . Prese le relazioni
Axi = i xi
Axj
= j xj
si moltiplichi la prima per xTj e la seconda per xTi , dove gli autovalori
sono certamente reali, per quanto detto sopra, e dunque sono reali
anche gli autovettori. Sottraendo le due relazioni si ottiene che
xTj Axi (xTi Axj ) = i xTj xi (j xTi xj )
3
La propriet`
a `e facilmente verificabile effettuando il prodotto scalare di due vettore per
componenti e osservando il primo termine del prodotto scalare nei due casi
145
ovvero
xTj Axi (Axj )T xi = (i j )xTj xi
Ma per la simmetria della matrice A il membro di sinistra `e nullo
e, dunque, deve essere nullo il prodotto xTj xi , il che comporta che
autovettori distinti sono perpendicolari.4
146
Capitolo 16
Equazioni di Eulero
dh
) = h B + B hB = mB
dt B
(16.1)
(16.2)
Jz 3 (Jx Jy )1 2 = m3
16.2
Per poter descrivere il moto del corpo rigido occorre definire la relazione fra
le componenti della velocit`
a angolare e il rateo di variazione degli angoli di
Eulero. Solo in questo modo `e possibile ricavare levoluzione della posizione
angolare della terna solidale al corpo rigido.
Se per descrivere lassetto del corpo rigido (ovvero lassetto della terna
ad esso solidale) si utilizzano gli angoli di Eulero nella sequenza 313 (pre147
148
= e
1 `e la linea dei nodi, data dalla posizione dellasse x lungo il
dove n
3 .
1 dopo la prima rotazione intorno a e
versore e
Proiettando tutti i versori in assi corpo si ha che
0
cos
sin sin
1
3
0
cos
In forma matriciale si ha che
cos
0
1
cos / cos
sin / sin
=
cos
sin
sin / tan cos / tan
Si osservi che per = 0 (condizione di singolarit`a per gli angoli di Eulero), la derivata degli angoli di precessione e rotazione propria diviene indeterminata, in quanto si hanno a denominatore delle quantit`a che tendono a
0. Questo `e il maggiore problema posto dalla rappresentazione degli assetti
tramite rotazioni elementari.
16.3
Se sul corpo rigido non agiscono forze, il suo centro di massa C si muove di
C , per cui un riferimento con origine in
moto rettilineo uniforme a velocit`a v
C, orientato come il riferimento inerziale, `e anchesso inerziale.
Se poi sul corpo rigido non agiscono neppure coppie esterne, i principi di
conservazione del momento angolare e dellenergia cinetica stabiliscono che
`e costante in un riferimento inerziale, per cui nel rifer il vettore h
ma
imento di assi solidali al corpo avr`a modulo costante, h = h,
componenti variabili nel tempo, nel caso pi`
u generale;
I ,
`e costante
lenergia cinetica associata alla rotazione, Trot = 21
(propriet`
a indipendente dal riferimento in cui vengono proiettate le
variabili del moto).
1
La cosa pu`
o essere fatta per via analitica, ma viene omessa e si fornisce solo il risultato
finale.
16.3.1
149
Poloide
= 1
12
22
32
+
+
2Trot /Jx 2Trot /Jy 2Trot /Jz
= 1
16.3.2
Erpoloide
h
T =
2
Dividendo ambo i membri per h e moltiplicandoli per 2 si ha che
h =
h
T
=2
h
h
`e fisso nel riferimento inerziale, il che significa che il vetIl versore h/h
ha una componente h lungo la direzione (costante) del momento
tore
visto
angolare che `e anch essa costante. Dunque lestremo del vettore ,
150
16.3.3
Rotazione pura
16.4
16.4.1
151
chiusa delle equazioni di Eulero per il caso triinerziale, in cui tutti e tre i
momenti principali di inerzia sono diversi fra loro. Tale soluzione, tuttavia,
`e estremamente complessa e viene tralasciata. Ma se il corpo rigido ha una
simmetria assiale, per cui due dei momenti principali di inerzia sono uguali,
le equazioni di Eulero ammettono una soluzione in forma chiusa. Si assuma
3 , per cui Jx = Jy = Jt `e il momento di
che lasse di simmetria inerziale sia
3 , mentre
inerzia trasversale, rispetto a un qualunque asse perpendicolare a
3 .
Jz = Ja `e il momento di inerzia intorno a
In assenza di coppie, le equazioni di Eulero in questo caso assumono la
formulazione
Jt 1 (Jt Ja )2 3 = 0
Jt 2 (Ja Jt )3 1 = 0
(16.3)
Ja 3 = 0
Dalla terza equazione si evince immediatamente che 3 (t) = , dove
`e costante. Sostituendo tale risultato nelle prime due equazioni si ha che
Jt 1 (Jt Ja )2 = 0
Jt 2 (Ja Jt )1 = 0
(16.4)
Dividendo entrambe le equazioni per Jt e ponendo (Jt Ja )/Jt = si ha
1 2 = 0
2 + 1 = 0
Derivando la prima equazione rispetto a t (
1 2 = 0) e sostituendo
lespressione di 2 ricavata dalla seconda ( 2 = 1 ) si ha che levoluzione
temporale di 1 `e governata dallequazione differenziale lineare omogenea
del secondo ordine
1 + 2 1 = 0
che ammette come soluzione
1 (t) = A sin(t) + B cos(t) = A0 sin(t + )
Osservando che 2 = (1/) 1 si ha che
2 (t) = A0 cos(t + )
Le costanti A0 e possono essere determinate dalle condizioni iniziali,
1 (0) = 10 e 2 (0) = 20 .
Si osservi che 12 + 22 = 12 = A20 `e costante. Pertanto le prime due
12 = 1
1 + 2
2
componenti della velocit`
a angolare descrivono un vettore
3 con pulsazione .
che gira intorno a
152
a)
b)
12
12
Fig. 16.1: Coni di Poinsot per corpi oblati (Caso a, Ja > Jt ) e prolati
(Caso b, Ja < Jt ).
3 ,
e il vettore momento angolare
= Jx 1
12 + Ja
1 + J y 2
2 + Jz
h
3 = Jt
3 ,
e
sono entrambi combigiacciono tutti sullo stesso piano , in quanto h
12 e
12 nelle sue rotazioni
3 . Il piano segue
nazioni lineari dei vettori
3 .
intorno a
e
3 `e tale
In un riferimento di assi solidali la distanza angolare fra
12 / = A0 / ed `e dunque costante, ragion per cui
descrive
che tan =
3 (Fig. 16.1). Tale cono
un cono di ampiezza intorno allasse di simmetria
`e fisso rispetto al corpo. Al tempo stesso anche la distanza angolare fra il
e
12 /Ja = (Jt /Ja ) tan .
3 `e costante, e pari a , con tan = Jt
vettore h
Per corpi oblati (Ja > Jt , come un disco o una moneta) si avr`a < , e
giace allinterno dellangolo formato
in questo caso il momento angolare h
e asse di simmetria
3 . Per corpi prolati (Ja < Jt ,
da velocit`
a angolare
come una matita o una mazza da baseball) risulta invece > , per cui `e
e
`e costante in un riferimento
a trovarsi fra h
3 . Dal momento che h
e il vettore velocit`a
inerziale, in tale riferimento il piano ruota intorno ad h
angolare descrive un cono di apertura , fisso nel riferimento inerziale.
Per corpi oblati il cono fisso nel riferimento inerziale si trova allinterno
del cono solidale al corpo (Fig. 16.1.a), mentre per corpi prolati i due coni
sono uno esterno allaltro (Fig. 16.1.b). Dal momento che lasse lungo
153
16.4.2
3 del riferiSi pu`o assumere senza alcuna perdita di generalit`a che lasse e
che, in assenza di coppie
mento inerziale sia parallelo al momento angolare h,
esterne, risulta inerzialmente costante. In questo caso, langolo di nutazione
risulta pari a ed `e costante. Sostituendo tale dato nelle equazioni che
regolano levoluzione degli angoli di Eulero in funzione dellevoluzione (nota)
delle componenti della velocit`
a angolare `e possibile definire in forma analitica landamento nel tempo dei parametri che descrivono lassetto del corpo
rigido.
= 0 e che sono note le evoluzioni temporali di
Ricordando che risulta
1 , 2 e 3 = , si ha che
sin sin
12 sin(t + ) =
sin cos
12 cos(t + ) =
cos +
(16.5)
(16.6)
per cui la derivata temporale dellangolo di precessione `e costante. Dividendo le prime due equazioni fra loro si ha
= = (Jt Ja )/Jt
tan(t + ) = tan = t + + k/2
Quindi anche la derivata dellangolo di rotazione propria `e costante.
Dalla terza equazione si ha che
= = Ja
cos cos Jt
ovvero il rateo di precessione `e proporzionale alla componente (costante) di
velocit`a angolare lungo lasse di simmetria ed `e maggiore di essa per corpi
oblati (Ja > Jt ) e inferiore per quelli prolati. Esplicitando dalla definizione
si ha anche che
di = ,
=
Ja
cos Jt Ja
154
Capitolo 17
17.1
Vincoli
17.1.1
Appoggio e scorrimento
156
forze attive o inerziali che ne provocano il distacco. Si rammenta che il vincolo di appoggio `e anolonomo e rappresenta lesempio pi`
u tipico di vincolo
unilatero.
Se il vincolo risulta invece bilaterale, ovvero il corpo si deve muovere
sulla superficie, si ha un vincolo olonomo, in quanto il vincolo riduce di
uno il numero di gradi di libert`a del corpo (o sistema di corpi) in esame.
Se il corpo `e vincolato a muoversi lungo una data curva, i gradi di libert`a
eliminati diventano due e si parla di vincolo di scorrimento, in quanto il
corpo pu`
o scorrere lungo una guida.
Ideale
Nel caso ideale la reazione vincolare f v `e normale alla superficie (o alla
curva). In presenza della sola forza gravitazionale sono possibili equilibri
solo in punti in cui la superficie sia (almeno localmente) orizzontale. In
generale, la reazione vincolare `e pari alla risultante componente totale delle
forze di massa (forze gravitazionali e inerziali) in direzione normale alla
superficie di appoggio o di scorrimento.
Con attrito
In presenza di attrito si genera una forza in verso opposta o a quello della
del corpo. Nel caso
direzione di avanzamento, identificato dalla velocit`a v
di attrico secco si ha una forza proporzionale alla reazione vincolare f v in
direzione normale, indipendente dallestensione della superficie di contatto.
Si hanno due casi in presenza o in assenza di moto relativo fra corpo e
superficie.
In presenza di moto, lattrito `e dato da
f ad = d f v (
v /v),
dove d `e un coefficiente di attrito dinamico, dipendente dalla scabrosit`a
delle superfici fra le quali si realizza il contatto.
In assenza di moto relativo, lattrito statico mantiene il sistema in quiete
fino a quando la risultante delle forze in direzione tangenziale alla superficie
non diventano superiori alla massima forza di attrito che pu`o essere sviluppata dal contatto fra le due superfici e che va sotto il nome di attrito di
primo distacco. Quindi il corpo rimane in quiete se `e soddisfatta la relazione
f t s f v
dove s `e un coefficiente di attrito statico. Per superfici non lubrificicate si
ha che s > d .
Dal punto di vista pratico, un equilibrio `e possibile se la direzione della
risultante delle forze agenti sul corpo giace allinterno di un cono avente per
17.1. VINCOLI
157
17.1.2
1
v
f aero = SCD v 2 ( )
2
v
dove `e la densit`
a del mezzo fluido, S `e una superficie di riferimento e CD
un coefficiente di resistenza.
17.1.3
Cerniera sferica
158
17.1.4
Cerniera cilindrica
17.1. VINCOLI
159
17.1.5
Cassetto
17.1.6
Incastro
Qualora il corpo rigido sia vincolato a un punto in modo tale che il vincolo
blocchi tutti i gradi di libert`
a di traslazione del punto vincolato e di rotazione
intorno ad assi passanti per esso, il moto del corpo risulta impedito. Un
singolo vincolo di incastro `e dunque in grado di portare a 0 il numero di
gradi di libert`
a residui del corpo e, per fare questo, deve essere in grado di
realizzare 3 componenti di forza e 3 di coppia, nel caso tridimensionale.
17.1.7
Carrello
17.1.8
160
17.2
17.3
Nel caso di un sistema olonomo, tutti gli spostamenti virtuali sono reversibili. Se infatti si rappresenta la configurazione del sistema in termini di
n coordinate generalizzate qk indipendenti, ciascuna coordinata pu`o essere
variata in modo arbitrario e, in particolare, il segno della variazione qk pu`o
essere positivo o negativo. In questo caso si ricorda che vale il principio dei
lavori virtuali secondo il quale il lavoro virtuale delle forze generalizzate si
annulla in una configurazione di equilibrio, il che equivale alla richiesta che,
in una condizione di equilibrio tutte le forze generalizzate Qk si annullino.
Se invece il sistema `e soggetto a vincoli unilateri esso diviene anolonomo.
Un vincolo unilatero pu`
o essere spesso rappresentato da una o due disequazioni del tipo `infk qk `supk , dove `infk e `supk rappresentano i limiti inferiore e superiori alle variazioni ammesse per qk . Questo `e possibile
quando il vincolo opera su una sola delle coordinate generalizzate, del sistema e non sulle altre.2 .
Le configurazioni possibili per un sistema soggetto a vincoli unilateri si
classificano in due categorie: configurazioni ordinarie, tali che nessun
vincolo unilatero `e attivo, e configurazioni di confine, in cui uno o pi`
u
vincoli unilateri risultano attivi. Nel primo caso gli spostamenti virtuali
ammissibili sono reversibili, in quanto non soggetti ad alcun vincolo. Nel
secondo, invece, si potranno avere solo spostamenti virtuali positivi (qk 0),
se `e attivo il vincolo inferiore (`infk = qk ), o spostamenti virtuali negativi
(qk 0), se `e attivo quello superiore (qk = `supk ). Se lo spostamento virtuale
deve avere un segno predeterminato, allora esso diviene irreversibile.
A questo punto la condizione di equilibrio di un sistema soggetto a vincoli unilateri non attivi `e indistinguibile da quella di un sistema olonomo,
in quanto in assenza di vincoli attivi tutti gli spostamenti virtuali sono reversibili e, dunque, arbitrari e la condizione di equilibrio `e sempre espressa
dalla nullit`
a delle forze generalizzate, ovvero dalla nullit`a del lavoro virtuale.
Se viceversa un vincolo unilatero `e attivo, gli spostamenti virtuali di una
o pi`
u coordinate generalizzate devono avere un segno predeterminato. In
questo caso la configurazione `e di equilibrio solo se il lavoro delle componenti
generalizzate delle forze attive risulta negativo, ovvero in verso opposto allo
spostamento virtuale ammissibile. Questo implica che le forze attive tendono a spingere il sistema considerato oltre il vincolo unilatero, e solo la
presenza del vincolo, che esplica una reazione sul sistema stesso, impedisce
tale passaggio. Pertanto per avere un equilibrio la forza generalizzata non
deve pi`
u annullarsi, ma solo avere un segno predeterminato.
2
Pi`
u in generale, un vincolo unilatero sar`
a rappresentato da una disequazione nella
forma f (q1 , q2 , ..., qn ) 0.
162
Capitolo 18
Rotazioni a un grado di
libert`
a
Quando i vincoli presenti nel sistema non bloccano tutti i gradi di libert`a
disponibili, il corpo o i corpi che lo costituiscono possono muoversi e avere
una loro dinamica, compatibile con i vincoli imposti. Si consideri a mo
desempio unasta rigida pesante di lunghezza ` i cui estremi A e B sono
vincolati a scorrere lungo due direzioni prefissate, una orizzontale lungo
lasse x per A e una verticale lungo y per B. Il sistema `e piano (3 gradi di
libert`a) ma la direzione di scorrimento per due punti dellasta `e vincolata,
e questo sottrae due gradi di liberti`a. Ne resta dunque uno solo e si pu`o
scegliere come coordinata generalizzata che rappresenta la posizione del sistema una coordinata fra xA , yB e (dove `e langolo compreso fra lasta
e la verticale locale). Scelta la coordinata generalizzata, `e possibile determinare le altre due in funzione di essa e avere la posizione dellasta descritta
in modo completo.
In generale `e spesso possibile determinare le reazioni vincolari dalle
equazioni della dinamica, in modo tale che forze e momenti sviluppati dai
vincoli inducano la variazione osservata per la posizione e lassetto del corpo
considerato. In questo caso il sistema `e detto labile, in quanto i vincoli cui `e
sottoposto consentono un atto di moto. La condizione di labilit`a `e ottenuta
anche quando latto di moto ammissibile `e solo virtuale.
18.1
Fra tutti i possibili vincoli cui pu`o essere sottoposto un corpo rigido, `e di
particolare interesse nelle applicazioni pratiche il vincolo che impone al corpo
di ruotare intorno a un asse fisso. La tipologia del vincolo `e la stessa sia
per casi in cui il moto di rotazione viene mantenuto a lungo in una data
direzione (le ruote di un treno, il mandrino di un trapano, lingranaggio
di un orologio), sia per rotazioni finite di ampiezza limitata che avvengono
163
164
`
CAPITOLO 18. ROTAZIONI A UN GRADO DI LIBERTA
alternativamente in versi opposti (ad esempio una finestra o una porta, una
biella, etc.).
Dal punto di vista cinematico il problema `e estremamente semplice in
quanto `e presente un solo grado di libert`a, , dove `e langolo compreso
, uno solidale
fra le normali a due semipiani uscenti dallasse di rotazione a
al corpo e laltro fisso rispetto a terra. La velocit`a angolare del corpo sara
a.
=
quindi data da
dh
a+m
v) a
= (m
a
dt
si ha che
Ia = ma
aa
, ma = m
`e la
dove Ia `e il momento di inerzia del corpo intorno ad a
delle coppie attive agenti sul corpo rotante e i vincoli
componente lungo a
v a
, sicche m
= 0.
non producono coppie intorno ad a
Bilanciamento statico
A questo punto `e possibile determinare le forze e i momenti sviluppati dai
vincoli. Consideriamo innanzi tutto il moto del centro di massa C. Se esso
, avremo che
si trova a distanza d da a
m
d
vC
= f a + f v
dt
d
vC
e2 + (d2 )
e2 (d2 )
2 = (d)
= (d)
e3 e
e1
dt
Dunque la forza vincolare risulta pari a
e2 (d2 )
f v = m [(d)
e1 ] f a
165
` evidente che
massa, quando questo si trova fuori dallasse di rotazione. E
1 e si traduce quindi in una sollecitazione
la reazione vincolare ruota con e
ciclica con pulsazione che si scarica sui cuscinetti che garantiscono la
, sollecitazione che si annulla se il centro di massa
rotazione intorno ad a
.
giace su a
`e detto
Definizione 1 Un copro rigido libero di ruotare intorno a un asse a
bilanciato staticamente se il suo centro di massa C giace sullasse di
.
rotazione a
Bilanciamento dinamico
= I ,
dove I `e il tensore di inerzia calcolato
Ricordando che, in generale, h
rispetto a una data terna solidale al corpo (non necessariamente baricen abbia una sola componente diversa da 0 (quella
trica!), nel caso in cui
e
3 nellesempio in esame) si ha che
lungo a
= I13
e1 + I23
e2 + Ia
e3
h
dove I13 e I23 sono i momenti misti. In questo caso, dunque, si ha che
dh
dt
1 + I23 e
2 + Ia e
3 ) + 2 e
3 (I13 e
1 + I23 e
2 + Ia e
3 ) =
= (I13 e
1 + I23 e
2 + Ia e
3 ) + 2 (I13 e
2 I23 e
1 )
= (I13 e
dh
a=
m
dt
1 + (I23 + I13 2 ) e
2
= (I13 I23 2 ) e
v =
m
166
`
CAPITOLO 18. ROTAZIONI A UN GRADO DI LIBERTA
`e detto
Definizione 2 Un copro rigido libero di ruotare intorno a un asse a
`e un asse princibilanciato dinamicamente se il suo asse di rotazione a
e velocit`
pale di inerzia, nel qual caso i vettori momento angolare h
a angolare
risultano entrambi paralleli ad a
stesso.
18.2
167
168
`
CAPITOLO 18. ROTAZIONI A UN GRADO DI LIBERTA
Parte IV
Stabilit`
a
169
Capitolo 19
Stabilit`
a: prime definizioni
Il concetto di stabilit`
a `e familiare e intuitivo anche in contesti non specialistici, dove un sistema `e detto stabile se le perturbazioni cui esso `e inevitabilmente sottoposto non ne alterano in modo considerato critico le condizioni.
Il concetto di stabilit`
a `e fondamentale in moltissime branche delle discipline
tecnicoscientifiche. In questi contesti, per`o, diviene necessario fornire una
formulazione matematica che permetta di definire in termini quantitativi
la condizione di stabilit`
a o instabilit`a di un sistema, e in questo caso le cose
diventano molto meno banali.
Peter Hughes cita ben 35 diverse definizioni matematiche di stabilit`a di
applicabilit`
a generale e un numero indefinito, se la definizione `e riferita ad
ambiti disciplinari specifici. Fortunatamente se si rimane nellambito della
Meccanica, `e sufficiente utilizzare un numero molto inferiore, anche se non
piccolissimo, di definizioni. Uno dei primi ambiti in cui si `e formulato un
requisito matematico per la nozione di stabilit`a `e quello della Meccanica
Celeste, in cui Lagrange defin` il moto di un astro rispetto a un altro come
stabile se la distanza fra i due corpi rimane finita. Si pu`o dunque dire che
un sistema meccanico `e stabile secondo Lagrange se il moto degli elementi
che lo compone rimane limitato.
` evidente che una nozione cos` blanda di stabilit`a `e di scarsa utilit`a
E
in ambiti diversi dalla Meccanica Celeste. Se il moto di un sistema rimane s` limitato, ma una perturbazione causa delle deviazioni molto grandi
dallevoluzione che si avrebbe in assenza di essa, `e chiaro che difficilmente si
` dunque necessario
definirebbe tale situazione come stabile. E
1. fornire un limite superiore allampiezza del moto e
2. fare in modo che tale limite divenga inferiore se la perturbazione diviene meno intensa.
Poste queste premesse, la definizione pi`
u utilizzata della nozione di stabilit`a `e di gran lunga quella proposta dal matematico russo Lyapunov.
171
` PRIME DEFINIZIONI
CAPITOLO 19. STABILITA:
172
19.1
19.1.1
Definizioni preliminari
19.1.2
173
19.1.3
Altre definizioni
19.2
Disturbi e incertezze
174
` PRIME DEFINIZIONI
CAPITOLO 19. STABILITA:
ipotesi che rappresentano inevitabilmente una idealizzazione e/o una semplificazione della realt`
a fisica delloggetto o sistema studiato. Inoltre i coefficienti numerici che ne rappresentano le caratteristiche (e.g. i momenti di
inerzia di un corpo rigido e la posizione del suo centro di massa) sono spesso
noti con un certo grado di approssimazione e quindi il comportamento del
sistema reale zar`
a influenzato dalla differenze fra il valore effettivo di tali
parametri e quelli assunti in sede di modellizzazione.
19.2.1
Stabilit`
a totale
19.2.2
Stabilit`
a strutturale
Leffetto delle discrepanze fra i valori nominali dei parametri e quelli effettivi pu`
o essere studiata esplicitando tale dipendenza nella formulazione del
modello matematico del sistema,
x = f (t, x; p)
` DI MOTI PERIODICI
19.3. STABILITA
175
19.3
Stabilit`
a di moti periodici
176
` PRIME DEFINIZIONI
CAPITOLO 19. STABILITA:
della soluzione nominale e due quella perturbata sono anche solo poco diversi, le due soluzioni tenderanno scostarsi una dallaltra. Tuttavia, in molti
casi pratici, pu`
o essere sufficiente richiedere che le due traiettorie rimangano
vicine nello spazio dello stato, senza preoccuparsi del fatto che lo sfasamento
fra di esse aumenti nel tempo a causa della diversa frequenza di oscillazione.
Il sistema `e certamente stabile dal punto di vista strutturale (la topologia delle soluzioni rimane infatti la stessa), ma le soluzioni dei due sistemi
non rimangono arbitrariamente vicine nel dominio del tempo, ma solo nello
spazio dello stato. Poincare ha introdotto in questo ambito il concetto di
stabilit`
a orbitale.
Definizione 6 Una soluzione periodica xT (t) `e detto orbitalmente stabile se, per ogni > 0 arbitrariamente scelto, esistono un numero reale > 0
e un istante t tali che risulti
(t; t0 , xT (t0 ) + x) xT (t) < , t > t0
se x < .
Definizione 7 Una soluzione periodica xT (t) `e asintoticamente orbitalmente stabile se `e orbitalmente stabile e inoltre risulta
(t; t0 , xT (t0 ) + x) xT (t) 0 per t .
Capitolo 20
Stabilit`
a di un equilibrio
20.1
f
A=
=
x e
f1
x1 e
f1
x2 e
f2
x1 e
f2
x2 e
fn
x1 e
fn
x2 e
f2
xn e
fn
xn e
f1
xn e
` DI UN EQUILIBRIO
CAPITOLO 20. STABILITA
178
20.2
t2 2 t3 3
A + A + ...
2
3!
si ha infatti che
t2
d
([etA ] x0 ) = (A + tA2 + A3 + . . . ) x0
dt
2!
= A (I + tA+
20.3
t2 2
A + . . . ) x0 = A [etA ] x0 = Ax(t)
2!
Se la matrice Jacobiana A calcolata nel punto di equilibrio ammette n autovalori distinti (o almeno n autovettori indipendenti q j organizzati per
colonne nella matrice Q) si pu`o riscrivere la soluzione nella forma
x(t) = Q [et ] Q1 x0
La matrice Q1 proietta la condizione iniziale x0 sugli autovettori che sono
linearmente indipendenti. Detta q j la jesima riga della matrice Q1 , il
coefficiente j = q j x0 `e detto ampiezza modale della condizione iniziale
rispetto al modo j, dove il jesimo modo `e rappresentato dalla coppia autovalore e autovettore, (j , q j ). Ciascun coefficiente j viene poi moltiplicato
per ej t e successivamente per la jesima colonna della matrice Q, ovvero
per il jesimo autovetture. La soluzione `e dunque ottenuta come sovrapposizione dellevoluzione delle ampiezze modali, che vanno a moltiplicare
ciascuna lautovettore corrispondente, ovvero
x(t) = q 1 1 e1 t + q 2 2 e2 t + + q n n en t
Dal momento che la matrice di stato A `e reale, il polinomio caratteristico
ha coefficienti reali e radici che possono essere o reali o a coppie di soluzioni
1
179
complesse coniugate. Ad ogni autovalore reale corrisponde un modo esponenziale. Se lautovalore `e reale e negativo, lampiezza modale corrispondente converge a 0, viceversa se `e reale positivo, essa diverge. Questo implica che il sistema `e instabile, in quanto lo stato si allontana dallorigine,
lungo la direzione dellautovettore corrispondente allautovalore positivo.
A una coppia di autovalori complessi coniugati, cc = a ib, corrispon cc e, di conseguenza, sono
dono due autovettori complessi coniugati, q cc e q
complesse coniugate anche le ampiezze modali iniziali, . Utilizzando la formula di Eulero, ei = cos +i sin `e per`o possibile considerare solo coefficienti
reali. Infatti si ha che
cc
q cc e(a+ib)t + q
e(aib)t = eat [q c cos(bt) + q s sin(bt)]
dove q c e q s sono combinazioni lineari delle parti reale e immaginaria dei
due autovettori coniugati. Dunque il moto associato a una coppia di autovalori e autovettori complessi coniugati corrisponde a un moto oscillatorio di
ampiezza modulata dal termine eat . Se la parte reale dellautovalore complesso `e negativa, lampiezza delle oscillazioni si riduce nel tempo. Viceversa,
se a > 0 lampiezza delle oscillazioni aumenta e lo stato del sistema diverge,
ovvero il sistema `e instabile.
Quindi, se gli autovalori sono tutti a parte reale negativa, tutte le componenti modali della condizione iniziale, = Q1 x0 , decrescono nel tempo
in modo monotono, se lautovalore `e reale negativo, in modo oscillatorio,
con oscillazioni di ampiezza decrescente, per coppie di autovalori complessi
coniugati. Dunque in questo caso la perturbazione iniziale decade nel tempo
e converge asintoticamente a zero. Pertanto `e possibile stabilire che
Teorema 11 Condizione sufficiente affinche un punto di equilibrio sia asintoticamente stabile `e che tutti gli autovalori della matrice Jacobiana A =
f /xe calcolata nel punto di equilibrio xe siano a parte reale negativa.
Teorema 12 Se la parte reale di uno o pi`
u autovalori risulta positiva, il
sistema `e instabile.
Teorema 13 Il sistema lineare risulta Lstabile in presenza di un autovalore a parte reale nulla o un autovalore reale nullo solo se gli autovalori sono
tutti distinti.
In presenza di autovalori a parte reale nulla, la stabilit`a del sistema non
pu`o essere determinata dallo studio degli autovalori se gli autovalori non
sono distinti.
20.4
Stabilit`
a per sistemi nonlineari
180
` DI UN EQUILIBRIO
CAPITOLO 20. STABILITA
20.5
181
20.5.1
20.5.2
Teorema 15 (Lstabilit`
a): lequilibrio nellorigine per il sistema x = f (x)
tale che f (0) = 0 `e Lstabile se esiste una funzione definita positiva V (x)
182
` DI UN EQUILIBRIO
CAPITOLO 20. STABILITA
20.6
Sistemi conservativi
183
[eM ]
= e M .
Tale propriet`
a si ottiene osservando che eM eM = I. Infatti, sviluppando
il prodotto
1
1
1
1
eM eM = (I + M + M 2 + M 3 + . . . ) (I M + M 2 M 3 + . . . )
2
3!
2
3!
si ha che
I
+M
1
M4 +...
+ 21 M 2 + 16 M 3 + 24
M 2
12 M 3 16 M 4
+ 12 M 2
+ 12 M 3
+ 41 M 4
16 M 3
16 M 4
+...
+...
+...
1
+ 24
M4 + = I
0
m2,2
0
M =
0
0
mn,n
` DI UN EQUILIBRIO
CAPITOLO 20. STABILITA
184
allora
em1,1
0
M
e =
em2,2
emn,n
1
1
= I + M + M2 + M3 + ...
2
3!
1
1
= QQ1 + QQ1 + (QQ1 )2 + (QQ1 )3 + . . .
2
3!
1 2 1 3
1
= Q (I + + + + . . . ) Q = Q (e ) Q1
2
3!
Detto
p() = det(I M ) = p0 + p1 + p2 2 + + pn1 n1 + n
Capitolo 21
Stabilit`
a di sistemi meccanici
lineari
In molti casi il risultato della modellizzazione di un sistema meccanico porta
alla definizione di un sistema di equazioni differenziali ordinarie del secondo
ordine scritto nella forma
+ (D + G)q + (K + A)q = Q
Mq
dove
M = M T `e la matrice di massa simmetrica definita positiva;
D = D T `e la matrice di smorzamento simmetrica semidefinita positiva;
G = GT `e la matrice giroscopica antisimmetrica;
K = K T `e la matrice di rigidezza simmetrica;
A = AT `e la matrice vincolare antisimmetrica;
Q `e il vettore di forze generalizzate.
In generale `e possibile convertire il sistema sopra scritto nella forma
x = Ax, ponendo xT = (q T , q T ) e partizionando opportunamente la matrice
A, ma a questo punto si pu`
o arrivare facilmente a sistemi di ordine molto
elevato.
21.1
Sistema massivo
186
21.2
Sistema giroscopico
21.3
Richiamo elastico
21.4
Smorzamento
187
21.5
Frequenze proprie
Data la linearit`
a del sistema, la soluzione per q(t) pu`o essere posta nella
forma
2n
q(t) = ai ei t
i=1
2n
q(t)
= i ai ei t
i=1
2n
(t) = 2i ai ei t .
q
i=1
Data la linearit`
a del sistema, vale il principio di sovrapposizione degli effetti.
188
Capitolo 22
Esercizi
22.1
Dinamica Newtoniana
1.
Un elemento materiale si trova in quiete sulla sommit`a di una sfera liscia
di raggio R. Lelemento comincia a scivolare lungo la sfera sotto lazione
della forza di gravit`
a. Determinare il punto da cui si stacca e la velocit`a
in quel punto. Assumere che gli attriti siano trascurabili e che il campo
gravitazionale sia uniforme e verticale.
Risposta: cos = 2/3.
2.
Determinare la gittata R di un proiettile lanciato con velocit`a u e alzo su
un piano inclinato di un angolo rispetto al piano orizzontale. Assumendo
un campo gravitazionale uniforme e verticale, determinare langolo che
massimizza la gittata.
Risposta: sin(2) = cos , ovvero = /4 + /2.
3.
Un uomo di massa m si trova a poppa di unimbarcazione di lunghezza L e
massa M . La prua tocca la banchina, ma il natante non `e ormeggiato. Cosa
accade quando luomo si sposta da poppa verso prua? Trascurare gli attriti.
Risposta: Limbarcazione si allontana dalla banchina di una distanza
d = L(m/M ).
4.
Tracciare il diagramma dellenergia associata a un carrello di massa m, collegato a due molle di rigidezza k, a riposo quando x = 0, tenendo conto che
189
190
per x > L/8 il carrello entra in contatto con unaltra molla, di rigidezza sempre pari a k. Descrivere il moto che si instaura partendo da una posizione
x = L/4 e determinare landamento della velocit`a in funzione di x.
5.
Due masse mA e mB sono collocate agli estremi di un bilancere inestensibile
di lunghezza L. Determinare il moto dovuto allazione di una forza impulsiva
F che agisce sul bilancere allistante t0 in direzione perpendicolare ad esso.
6.
Due corpi di massa mA e mB orbitano intorno al comune centro di massa C
sotto lazione della mutua attrazione di gravit`a mantenendo una distanza R
costante (orbita relativa circolare) con velocit`a angolare
c = R3 / G(mA + mB ).
Un terzo corpo, di massa m trascurabile, si trova lungo la congiungente. Determinare le sue condizioni di equilibrio lungo la direzione AB, discutendone
qualitativamente la stabilit`a.
Determinare analiticamente la posizione di equilibrio, supponendo m
mB mA , ed effettuando il calcolo per il sistema Terra-Luna ( = mB /(mA +
mB ) = 0.01215 e R = 384 750 km).
22.2
191
Moti relativi
1.
Calcolare la deviazione della verticale apparente rispetto a quella geometrica
alla latitudine di Lecce (40 N circa).
2.
Calcolare la deviazione verso Est di un grave lasciato cadere dalla terrazza
panoramica della Tour Eiffel (h = 310 m, a una latitudine approssimativa
di 49 N).
3.
Il rotore di un elicottero ruota con velocit`a angolare costante pari a 260 rpm.
Quando lelicottero viaggia a una velocit`a di avanzamento di 100 km/h le
pale sono soggette a un moto di flappeggio () fuori dal piano di rotazione
intorno a una cerniera C posta a 0.4 m dallasse del mozzo, con andamento
(in prima approssimazione) cosinusoidale, se langolo di rotazione della pala
= t `e contato a partire dalla direzione che punta verso poppa.
Sapendo che la pala ha un raggio R = 8.2 m e che lampiezza del moto di
flappeggio `e pari a -3 , determinare le componenti della velocit`a e dellaccelerazione
assolute del punto allestremit`
a della pala in funzione di e calcolare il loro
modulo in un riferimento fisso rispetto allelicottero per le 4 posizioni = 0,
90 , 180 e 270 , dove gli assi solidali alla fusoliera sono diretti nel seguente
modo: x verso prua (concorde con la velocit`a di avanzamento), y verso la
destra del pilota e z verso il basso.
Nota: `e possibile assumere langolo di flappeggio come piccolo, tale che
cos 1 e sin .
4.
Ripetere il calcolo proposto nellesercizio precedente in un riferimento solidale con il rotore e (per la sola velocit`a) in un riferimento solidale con la
pala.
192
22.3
Equazioni di Lagrange
1.
Una massa m `e appesa a due molle di uguali caratteristiche, agganciate a due
uniforme. In
punti di sospensione a distanza d in un campo gravitazionale g
condizioni di equilibrio, quando la tensione nelle molle bilancia esattamente
il peso, le molle formano un angolo di 45 con la verticale locale. Derivare
le equazioni di Lagrange per il sistema, ipotizzando che le molle di rigidezza
k abbiano massa trascurabile.
2.
Riscrivere le equazioni, ipotizzando che gli spostamente di m siano piccoli
(x, y d).
3.
Osservando le equazioni `e possibile identificare delle coordinate lungo le
quali le equazioni relative ai due gradi di libert`a si disaccoppiano. Scrivere
` possibile disaccoppiare
le equazioni linearizzate in forma disaccoppiata. E
le equazioni, se langolo formato dalle molle rispetto alla verticale locale
dovesse essere diverso da 45 ?
4.
Ripetere gli esercizi 1. e 2., ipotizzando che la massa sia collegata a terra
anche da due smorzatori, sempre di massa trascurabile, posti sul lato inferiore, anchessi inclinati a 45 rispetto alla verticale locale, che generano una
d
forza f = c`.
5.
Ripetere lesercizio precedente, ipotizzando che uno degli smorzatori sia
rotto. Le equazioni del moto sono ancora smorzate?
6.
Ricavare le equazioni del moto del pendolo elastico nel caso ideale, utilizzando come coordinate generalizzate lo spostamento angolare del pendolo
dalla verticale locale e la sua lunghezza `.
7.
Ricavare le equazioni del moto del pendolo elastico in presenza di attrito
viscoso nel sistema di sospensione e con il mezzo fluido.
22.4
193
1.
Determinare la posizione del baricentro e i momenti di inerzia baricentrici
di un parallelepipedo omogeneo di lati 2a, 2b e 2c.
2.
Determinare i momenti di inerzia di una barra omogenea di lunghezza ` e
sezione trascurabile rispetto alla sua dimensione lineare. Effettuare il calcolo
a) utilizzando una distribuzione lineare di massa per unit`a di lunghezza
uniforme e b) utilizzando le formule del punto precedente, assumendo che
sia c = `/2 a, b. Valutare infine il momento di inerzia rispetto a un estremo
della barra.
3.
Determinare i momenti di inerzia baricentrici di un disco di raggio R.
4.
In una lastra piana di lati 2a e 2b (con a > b) e densit`a uniforme viene
effettuato un foro centrale di forma circolare di raggio R < b). Valutare i
momenti di inerzia della piastra forata.
5.
Per un errore di lavorazione, il foro viene effettuato in una posizione errata,
dove a a e b b sono gli errori nella posizione del centro del foro rispetto
al centro della lastra, misurati rispettivamente in direzione parallela ai lati
2a e 2b. Valutare i momenti di inerzia della piastra forata.
6.
Determinare la posizione del baricentro di una lastra triangolare con un
angolo retto fra i cateti di lunghezza a e b.
7.
Determinare la posizione del centro di massa di un cono di altezza h e raggio
di base R.
Risposta:
194
8.
Determinare i momenti di inerzia baricentrici del cono e quelli rispetto a un
sistema di assi principali centrati sul centro della base.
Risposta:
` E PICCOLE OSCILLAZIONI
22.5. STABILITA
22.5
195
Stabilit`
a e piccole oscillazioni
1.
Dimostrare che una sola forma di dissipazione nellesercizio 7. del Par. 22.3
`e insufficiente per rendere asintoticamente stabile la posizione di equilibrio
con = 0 del pendolo elastico. Quale?
2.
Determinare il periodo delle piccole oscillazioni di un pendolo composto
formato da una barra omogenea sospesa a un suo estremo.
3.
Determinare le equazioni del moto di un doppio pendolo composto (formato
da un pendolo composto sospeso per un suo estremo al punto O fisso a terra
tramite un cavo inestensibile di massa trascurabile) assumento M per la
massa del pendolo, L la sua lunghezza, ` la lunghezza del cavo e trascurabili
gli effetti di dissipazione. Linearizzare le equazioni del moto e analizzare le
frequenze proprie di oscillazione.
4.
Un disco omogeneo di massa M e raggio r rotola senza strisciare sul fondo
di una guida circolare di raggio R. Determinare il periodo delle piccole
oscillazioni.
5.
Lo stesso disco `e sbilanciato al suo estremo da una massa m posta a distanza
r dal suo centro. Determinare le condizioni di equilibrio in funzione del
rapporto = m/M e della posizione angolare che la massa m ha rispetto
alla verticale locale quando il punto di contatto fra disco e guida `e nel
punto pi`
u basso della guida. Studiare la stabilit`a e la frequenza delle piccole
oscillazioni.
6.
Il disco sbilanciato `e ora posizionato su un piano inclinato con un angolo .
Verificare lesistenza o meno di posizioni di equilibrio al variare di e .
7.
Un disco scorre lungo una guida rettilinea di massa trascurabile incernierata
a un suo estremo e sospesa con una molla di costante elastica k e lunghezza
196
8.
Ripetere lesercizio precedente, assumendo che le due molle abbiano costanti
elastiche diverse. Verificare quale impatto ci`o pu`o avere sulla stabilit`a degli
equilibri.
9.
Un carrello di massa M scorre su una guida orizzontale tramite due ruote di
massa m e raggio r, trattenuto da una molla di costante elastica k. Studiare
la frequenza delle piccole oscillazini del carrello intorno al punto di equilibrio.
10.
Sul carrello dellesercizio precedente `e sospeso un pendolo di lunghezza `
e massa M . Determinare le equazioni del moto e le frequenze proprie
nellintorno della posizione di equilibrio.
22.6
197
Moti rigidi
1.
Un robot `e composto da due bracci snodati AB e BC, di lunghezza `1 ed `2 .
Il braccio AB `e collegato a terra nel punto A con uno snodo cilindrico ad
asse orizzontale, montato su un supporto rotante intorno a un asse verticale,
snodi che realizzano un giunto sferico con due gradi di libert`a (azimuth
ed elevazione ), mentre in B i due bracci sono articolati fra loro con una
semplice cerniera cilindrica ad asse orizzontale e formano fra loro un angolo
. In C `e collocato un utensile che deve eseguire una lavorazione con velocit`a
di avanzamento v costante lungo una linea retta che passa a una distanza
d da A. Calcolare levoluzione temporale degli angoli nelle tre cerniere (due
in A e una in B).
198
Bibliografia
[1] Paolo Biscari, Tommaso Ruggeri, GIuseppe Saccomandi e Maurizio
Vianello, Meccanica Razionale per lIngegneria, 2 ed., Monduzzi Editoriale, Bologna, 2007.
[2] Piero G. Bordoni, Lezioni di Meccanica Razionale, Veschi Ed., Roma,
1972.
[3] Bruno Finzi, Meccanica Razionale, Zanichelli, Bologna, 19??.
[4] Herbert Goldstein, Charles Poole, and John Safko, Classical Mechanics,
3rd Edition, Addison Wesley, ???, 19??.
[5] Peter C. Hughes, Spacecraft Attitude Dynamics, Revised Edition,
Dover, Mineola (NY), 2004.
[6] Lev D. Landau e Evgenij M. Lifsits, Fisica teorica: 1 - Meccanica,
Editori Riuniti, Roma, 2009.
[7] Leonard Meirovitch, Methods of Analytical Dynamics, 2nd Edition,
Dover, Mineola (NY), 2003.
199