Fisica 1
Fisica 1
MECCANICA CLASSICA
Roberto Bonciani∗
Dipartimento di Fisica, Università di Firenze
e INFN Sezione di Firenze,
via Sansone 1,
50019 Sesto Fiorentino
∗
Email: [Link]@[Link]
Indice
1 INTRODUZIONE 6
1.1 Metodo Scientifico . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.1.1 Il Modello . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.1.2 Grandezze fisiche e loro misura . . . . . . . . . . . . . . . . . . . . . . . 7
1.2 Richiami di calcolo vettoriale . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
1.3 Rappresentazione Cartesiana . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
2 APPUNTI DI CINEMATICA 12
2.1 Sistema di riferimento . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
2.2 Punto Materiale . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
2.3 Parametrizzazione naturale . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
2.4 Moto rettilineo di un punto . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
2.4.1 Moto rettilineo uniforme . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
2.4.2 Moto uniformemente accelerato . . . . . . . . . . . . . . . . . . . . . . . 18
2.4.3 Esempio: caduta dei gravi . . . . . . . . . . . . . . . . . . . . . . . . . . 19
2.5 Moto piano generico di un punto . . . . . . . . . . . . . . . . . . . . . . . . . . 22
2.5.1 Derivata del vettore posizione . . . . . . . . . . . . . . . . . . . . . . . . 23
2.6 Parametrizzazione naturale . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
2.7 Velocità e accelerazione . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
2.7.1 Moto circolare . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
2.8 Moto piano in coordinate cartesiane . . . . . . . . . . . . . . . . . . . . . . . . . 28
2.9 Moto piano in coordinate polari . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
2.10 Moto circolare in coordinate polari . . . . . . . . . . . . . . . . . . . . . . . . . 29
2.10.1 Grandezze angolari . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
2.11 Moto armonico . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
2.12 Esempio di moto piano: balistica . . . . . . . . . . . . . . . . . . . . . . . . . . 33
2.13 Cinematica relativa . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
2.13.1 Esempio: treno-binario . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
2.13.2 Esempio: giostra . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
2.13.3 Esempio: deviazione verso oriente . . . . . . . . . . . . . . . . . . . . . . 41
1
3.4 Vincoli . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48
3.4.1 Vincoli lisci . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48
3.4.2 Vincoli con attrito . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
3.4.3 Esempio . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
2
6 APPUNTI SULLA GRAVITAZIONE 96
6.1 Leggi di Keplero . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 96
6.1.1 Velocità areolare . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 96
6.1.2 Ellisse in coordinate cartesiane . . . . . . . . . . . . . . . . . . . . . . . 97
6.1.3 Ellisse in coordinate polari . . . . . . . . . . . . . . . . . . . . . . . . . . 98
6.2 Da Keplero a Newton . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 99
6.2.1 Le orbite sono piane ... . . . . . . . . . . . . . . . . . . . . . . . . . . . . 99
6.2.2 La velocità areolare è costante ... . . . . . . . . . . . . . . . . . . . . . . 100
6.2.3 Andamento 1/r 2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 100
6.2.4 Dalla Terza Legge ... . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 101
6.3 Da Newton a Keplero . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 102
6.3.1 Traiettoria e terza legge di Keplero . . . . . . . . . . . . . . . . . . . . . 106
6.4 Teorema di Gauss . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 107
6.4.1 Distribuzione di Massa a simmetria sferica . . . . . . . . . . . . . . . . . 109
6.4.2 Misura della costante di gravitazione universale . . . . . . . . . . . . . . 110
6.4.3 Forza peso . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 110
6.5 Esercizi . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 112
6.5.1 Potenziale centrale elastico . . . . . . . . . . . . . . . . . . . . . . . . . . 112
6.5.2 Potenziale centrale lineare in r . . . . . . . . . . . . . . . . . . . . . . . . 112
3
9.7.3 Pendolo di torsione . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 140
9.7.4 Rotolamento puro . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 140
9.8 Energia cinetica del corpo rigido . . . . . . . . . . . . . . . . . . . . . . . . . . . 144
9.8.1 Asse fisso . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 145
9.8.2 Moto piano . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 145
9.9 Lavoro e Teorema delle forze vive . . . . . . . . . . . . . . . . . . . . . . . . . . 145
9.10 Espressione generale del momento angolare. Tensore d’inerzia . . . . . . . . . . 146
9.11 Equazioni di Eulero . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 147
9.12 Espressione generale dell’energia cinetica . . . . . . . . . . . . . . . . . . . . . . 148
4
Questi appunti non vogliono essere una
trattazione esaustiva della meccanica,
ma soltanto una schematica presen-
tazione di ciò che è stato trattato a
lezione.
5
Capitolo 1
INTRODUZIONE
4. Esperimento. Una volta fatte delle predizioni teoriche sulla base del modello, tali predi-
zioni devono ripassare al vaglio dell’Esperimento, il quale potrà confermare il modello per
altri processi (corroborandone la struttura) oppure falsificarlo. Se l’Esperimento mette in
evidenza un processo per il quale il Modello non dà una descrizione consistente, il Modello
deve essere raffinato e cambiato per includere nella descrizione il nuovo fenomeno. Ciò
permette al Modello di raggiungere un livello di generalità più elevato.
1.1.1 Il Modello
Il Modello (o la Teoria) si esplica in “Principi”, affermazioni semplici di caratere generale che pos-
sano essere formalizzate con delle “Leggi Fisiche”, cioè delle relazioni quanitative fra “grandezze
fisiche”. Il Modello è una schematizzazione della realtà, che descrive in forma semplificata.
6
1.1.2 Grandezze fisiche e loro misura
Le “grandezze fisiche” sono le grandezze che servono per la descrizione oggettiva del fenome-
no sotto studio. Per esempio, nello studio del movimento di un corpo la distanza del corpo
stesso dall’osservatore oppure il tempo che il corpo impiega a raggiungere un punto sulla sua
traiettoria, sono entrambe grandezze fisiche.
Misura
Nella descrizione quantitativa di un fenomeno si pone il problema di come associare un numero
ad una grandezza fisica. Ciò viene definito come procedimento di “misura”, che deve essere il
più possibile univoco (non ambiguo) e deve essere riproducibile.
Praticamente, il procedimento di misura si quantifica con il confronto della grandezza fisica
con un “campione”, che chiameremo “unità di misura”. Al campione verrà fatto corrispondere
il numero 1 in quella unità di misura.
Il rapporto fra grandezza e campione, ovvero il numero di volte in cui il campione “entra”
nella grandezza è la misura di quest’ultima. Questo concetto si rifà direttamente alla misura di
una lunghezza, nella quale per esempio contiamo quante volte entra il campione “metro” nella
distanza che devo misurare. Per il tempo, le cose sono un po’ più complicate. La misura di un
intervallo di tempo si deve fare per confronto con due eventi consecutivi. È chiara la necessità
di avere eventi che facciano parte di un processo periodico, come per esempio un pendolo: Se
diamo una certa durata alla oscillazione completa del pendolo, potremo misurare un intervallo
di tempo fra due eventi, contando il numero di oscillazioni pendolari che li separa.
Si parla di misurazione diretta o indiretta. Nel primo caso la grandezza viene misurata
con un confronto diretto col campione (esempio: lunghezza). Il secondo caso invece è rilevante
quando la grandezza in questione è una funzione di grandezze “più semplici”, che diremo fonda-
mentali, e quindi la misura di tale grandezza deriva dalla funzione delle misure delle grandezze
che la definiscono. Per esempio, la velocità, che è data dal rapporto di una lunghezza e di un
tempo, sarà misurata tramite tale rapporto una volta che vengano indipendentemente misurate
le grandezze lunghezza e tempo pertinenti al processo in esame.
Grandezze Fondamentali
Le grandezze fondamentali sono le grandezze in termini delle quali possono essere espresse tutte
le altre. In meccanica storicamente si usano come grandezze fondamentali lunghezza, tempo e
massa e si indicano le rispettive dimensioni con i simboli [l], [t] e [m].
Le altre grandezze si esprimono in funzione delle fondamentali. Per esempio la velocità sarà
espressa dal rapporto di una lunghezza e di un intervallo di tempo e quindi avrà dimensione
[v] = [l] [t]−1 .
Grandezze prive di dimensioni vengono dette adimensionali. Per esempio il rapporto di due
grandezze omogenee, ovvero con la stessa dimensione, è adimensionale. Un angolo, definito in
radianti come rapporto di due lunghezze, è una quantità adimensionale: [α] = [l] [l]−1 = 0.
La scelta delle grandezze fondamentali è legata alla riproducibilità e alla precisione di mi-
sura delle grandezze stesse. Storicamente si usano lunghezza, tempo e massa come grandezze
fondamentali perché è immediato pensare a come misurare una lunghezza col metro campione,
il tempo come confronto con un fenomeno periodico (esempio il pendolo) e la massa con la
bilancia a due bracci. Queste definizioni però sono cambiate nel tempo e la necessità di avere
7
un campione sempre più preciso e riproducibile ha fatto evolvere il concetto di metro, secondo
e chilogrammo massa nel corso dei secoli.
Unità di misura
Il sistema di unità di misura è l’insieme dei campioni scelti per le grandezze fisiche fondamen-
tali. Il Sistema Internazionale, per esempio, utilizza come unità di misura della lunghezza il
metro, che ha cambiato nei secoli diverse definizioni operative, come unità di misuta del tempo
il secondo e come unità di massa il kg massa. Il metro, inizialmente definito come la lunghezza
di un pendolo il cui semiperiodo sia pari ad un secondo, venne poi legato alla fine del ’700 alla
misura del meridiano terrestre (1m=1/10000000 del semi-meridiano terrestre, dal Polo Nord
all’Equatore) e nel 1889 alla distanza fra due tacche incise su una barra di platino-iridio con-
servata in condizioni standard al Bureau International des Poids et Mesures di Sèvres, vicino
a Parigi. La precisione di tale definizione è limitata, per esempio, dallo spessore delle tacche.
Con l’avvento dei laser, negli anni ’60 si definisce il metro come 1650763.73 volte la lunghezza
d’onda nel vuoto della radiazione conseguente alla transizione 2p10 → 5d5 del kripton-86. Suc-
cessivamente, il metro diviene effettivamente una grandezza derivata, essendo definito dal 1983
come la distanza che la luce percorre in 1/299792458 di secondo. Tale definizione permette
una precisione dell’ordine di ∼ 10−10 . Anche per il secondo si passa dalla sua definizione come
1/86400 del giorno solare medio, poi specificato come “giorno solare medio dell’anno 1900”, al-
l’introduzione dell’orologio atomico nel 1967, con cui il secondo viene definito come 9192631770
periodi della radiazione corrispondente alla transizione fra due livelli iperfini dello stato fon-
damentale dell’atomo di cesio-133. Questa definizione permette una precisione di ∼ 10−12 . Il
chilogrammo è definito come la massa di un campione in lega di platino-iridio, conservato al
Bureau International des Poids et Mesures di Sèvres ed ha una precisione intrinseca di ∼ 10−8 .
Per facilitare le misure col campione, se ne definiscono multipli e sottomultipli.
La branca della fisica che studia la corretta definizione delle grandezze indipendenti e le
corrispondenti unità di misura è detta metrologia.
• vettori: sono grandezze caratterizzate da un’intensità (il modulo del vettore), da una
direzione e da un verso. Per esempio la velocità: per la caratterizzazione del moto non
basta dire che procedo a 50 km/h ma occorre precisare anche la mia direzione (Roma-
Milano) e il mio verso (da Roma a Milano oppure viceversa da Milano a Roma).
8
Definizione 1.2.1 Dati due vettori a e b, si dice che a = b se kak = kbk e se hanno la stessa
direzione e lo stesso verso. Quindi a = −b se kak = kbk e se hanno stessa direzione e verso
opposto.
Definizione 1.2.2 Si definisce come somma fra due vettori a e b il vettore c = a + b che si
trova con la regola del parallelogramma. Analogamente a = c − b.
a= c−b
a c
c= a+b
b −b
z
d
c
a
b
• Proprietà commutativa: a + b = b + a.
• Proprietà associativa: a + (b + c) = (a + b) + c.
• (λ + µ)v = λv + µv.
• λ(v + w) = λv + λw
• λv = 0 implica λ = 0 o v = 0.
v
Definizione 1.2.4 Si dice versore un vettore di modulo unitario: v̂ = kvk
, kv̂k = 1.
Definizione 1.2.5 Si dice prodotto scalare fra i due vettori v e w, appartenenti ad uno spazio
vettoriale sul campo reale, una funzione che associa ai due vettori un numero del campo così
definito:
P (v, w) = v · w = kvk kwk cos α , (1.1)
dove α è l’angolo fra i due vettori.
• v · v = kvk2 ,
• v·w=w·v ,
9
• v · (w + z) = v · w + v · z ,
• v · (λw) = λv · w .
Definizione 1.2.6 Si dice prodotto vettoriale fra i due vettori v e w, appartenenti ad uno
spazio vettoriale sul campo reale, una funzione che associa ai due vettori un altro vettore dello
stesso spazio vettoriale così definito:
v ∧w = z, (1.2)
tale che kzk = kvk kwk sin α, la direzione è quella perpendicolare al piano individuato dai
vettori v e w, e il verso è quello della regola della mano destra, ovvero il vettore z punta nel
verso da cui si veda ruotare v su w in senso antiorario coprendo l’angolo acuto.
• v∧v=0 ,
• v ∧ (−w) = −v ∧ w ,
• v ∧ w = −w ∧ v ,
• v ∧ (w + z) = v ∧ w + v ∧ z ,
• v ∧ (λw) = λv ∧ w .
• a ∧ (b ∧ c) = b(a · c) − c(a · b) (“bac men cab”).
.. nota che a ∧ (b ∧ c) 6= (a ∧ b) ∧ c.
v
k̂ î · ĵ = ĵ · k̂ = k̂ · î = 0 ,
y î ∧ ĵ = k̂ , ĵ ∧ k̂ = î , k̂ ∧ î = ĵ .
0 ĵ
î
10
Allora si trova che:
• v · w = vx wx + vy wy + vz wz ,
î ĵ k̂
• v ∧ w = vx vy vz = (vy wz − vz wy )î + (vz wx − vx wz )ĵ + (vx wy − vy wx )k̂.
wx wy wz
11
Capitolo 2
APPUNTI DI CINEMATICA
In queto capitolo tratteremo la Cinematica, ovvero lo studio del movimento di un sistema fisico
in un certo sitema di riferimento, indipendentemente dalle cause che generano il moto.
Corpo
z
SdR
k̂
y
0 ĵ
î
Dovendo misurare delle distanze e degli intervalli di tempo, saranno necessarie alcune
precisazioni sulla struttura dello spazio in cui opereremo e sul tempo stesso.
• Lo spazio in cui intendiamo ambientare il nostro studio è uno spazio vettoriale Euclideo,
ovvero valgono gli assiomi della geometria Euclidea. Che questo sia vero, almeno nel
grado di approssimazione che ci interessa, può essere verificato sperimentalmente (Es:
esperimento di Gauss sulla somma degli angoli interni di un triangolo). Inoltre lo spazio
è omogeneo (cioè non ci sono punti o zone “privilegiate” rispetto ad altre) e isotropo (cioè
non ci sono direzioni privilegiate rispetto ad altre).
12
• In meccanica Newtoniana il tempo costituisce un parametro esterno ed assoluto, che se-
gna il susseguirsi degli eventi. Per assoluto, si intende il fatto che in ogni punto del nostro
spazio, il tempo scorre nello stesso identico modo. Due orologi che vengano sincroniz-
zati in un punto del nostro SdR, in un centro istante, rimangono sincronizzati sempre e
ovunque. Questa in realtà è un’approssimazione che vale finché le velocità in gioco sono
piccole rispetto alla velocità della luce e finché i campi gravitazionali considerati sono
sufficientemente poco intensi, ma per lo studio della meccanica Newtoniana questa sarà
un’assunzione esatta.
Definizione 2.2.1 Il punto material è un oggetto privo di dimensioni (punto geometrico) che
eventualmente ha massa m.
Definizione 2.2.2 I gradi di libertà di un sistema sono i parametri indipendenti necessari per
definirne univocamente la configurazione rispetto al SdR scelto.
Il punto materiale P nel sistema di riferimento SdR sarà individuato dal suo raggio vettore
r(t) (un vettore in R3 , quindi appunto 3 parametri indipendenti):
z
SdR
P
r(t)
k̂ γ
y
0 ĵ
î
r = r(t) , (2.1)
13
detta equazione vettoriale del moto. In componenti, rispetto al SdR ortonormale si ha
Allo scorrere del tempo (parametro t) l’Eq. (2.2) descrive la traiettoria γ, ovvero il luogo
geometrico dei punti toccati da P nel suo moto.
0
s P
Tale lunghezza cambia nel tempo, dipendendo dal moto di P su γ. Quindi sarà una fun-
zione del tempo. Supponiamo che sia una funzione monotona e derivabile. Allora, possiamo
parametrizzare il moto di P su γ in termini di s e l’equazione vettoriale del moto potrà essere
scritta come segue:
r = r(s) , (2.3)
s = s(t) . (2.4)
In questo modo abbiamo diviso la parte puramente geometrica del moto da quella cinematica.
La prima è descritta dall’Eq. (2.3) e riguarda la traiettoria come luogo geometrico, parametriz-
zata dalla lunghezza della curva s. La seconda, in Eq. (2.4), mi dice come si muove il punto su
tale traiettoria ed è detta equazione oraria.
P
î
0 x(t) x
14
Il moto del punto P è dato dal cambiamento della sua coordinata, x, rispetto al tempo:
x = x(t) . (2.5)
Nota: in tale sistema di riferimento non useremo la notazione vettoriale x(t) = x(t)î, ma
utilizzeremo semplicemente la funzione x(t) col suo segno. Ciò è sufficiente per la descrizione
completa del moto. Infatti la direzione è fissata a priori, rimane dunque soltanto l’arbitrarietà
del verso, che si traduce nel segno di x.
La funzione x(t) determina il “diagramma orario” del punto P :
x(t)
x0
0 t
v = v(t) (2.8)
15
per cui conviene, allo stesso modo di come è stato fatto per la coordinata x(t), introdurre la
variazione di v col tempo, che si chiamerà accelerazione. Avremo un’accelerazione media data
da
∆v
ā = (2.9)
∆t
e un’accelerazione istantanea, in cui si prende il limite in cui i ∆ti diventano sempre più piccoli
(e più numerosi):
∆v
a = lim . (2.10)
∆t→0 ∆t
In realtà anche l’accelerazione a sua volta è una funzione di t e quindi potremmo introdurre
la sua variazione rispetto al tempo, ma la meccanica si ferma alle accelerazioni ...
Dal punto di vista matematico si ha:
16
La costante x(t0 ) è la costante d’integrazione che necessariamente appare quando risolvo l’equa-
zione differenziale (2.15). Per trovare una soluzione particolare, ovvero per determinare proprio
il moto che sto studiando, devo imporre una condizione iniziale. Per esempio sia
x(t0 ) = x0 . (2.17)
x(t) = x0 + v0 t . (2.18)
Possiamo riportare su due grafici il moto del punto P disegnando la velocità e lo spazio come
funzione del tempo. Prendendo una velocità v0 positiva si ha:
x(t) v(t)
v0
x0
O t O t
Esercizio
Supponiamo che due punti si muovano di moto rettilineo uniforme sullo stesso asse x. A t = 0
il primo punto, P1 , si trova in x = x1 ed ha una velocità costante v = v1 , mentre il secondo
punto, P2 , si trova in x = x2 ed ha velocità costante v = v2 . Discutiamo le situazioni in cui i
due punti ad un certo istante si urtano e determiniamo dove e quando ciò avverrà.
Il moto del punto P1 sarà dato dalla seguente legge oraria:
x1 (t) = x1 + v1 t , (2.19)
17
Per trovare dove si urteranno, basta sostituire l’espressione di t0 in uno dei due moti (quello
di P1 oppure equivalentemente quello di P2 ):
x1 − x2 x1 v2 − x2 v1
x0 = x(t0 ) = x1 + v1 = . (2.23)
v2 − v1 v2 − v1
Si hanno quindi tre casi (prendiamo per semplicità e senza perdita di generalità x1 > 0 e
x2 > 0):
1. Primo caso. Si ha x2 > x1 e prendendo entrambe le velocità positive, v1 > 0 e v2 > 0,
si deve avere v1 > v2 . I due punti vanno entrambi nella stessa direzione (delle x positive).
Il punto P1 rincorre il punto P2 e lo raggiunge in x0 , che quindi sarà oltre x2 .
2. Secondo caso. Si ha x2 > x1 e prendendo v1 > 0 e v2 < 0, si avrà che i due punti vanno
l’uno verso l’altro e si incontreranno in x0 tale che x1 < x0 < x2 .
x1 x1 x1
x0
0 t0 t 0 t0 t 0 t0 t
Primo Caso Secondo Caso Terzo Caso
18
1
= x(0) + v(0)t + a0 t2 . (2.27)
2
La (2.27) è la soluzione generale dell’equazione del moto (2.25). Abbiamo adesso dovuto ri-
solvere un’equazione differenziale del secondo ordine e quindi abbiamo due costanti arbitrarie
d’integrazione, x(0) e v(0). Per determinare la soluzione particolare, dobbiamo imporre le
condizioni iniziali. Sia per esempio
x(0) = x0 , (2.28)
v(0) = v0 . (2.29)
Si ha quindi:
1
x(t) = x0 + v0 t + a0 t2 , (2.30)
2
v(t) = v0 + a0 t . (2.31)
Come nel caso del moto rettilineo uniforme possiamo riportare in grafico la posizione, la
velocità e l’accelerazione come funzioni del tempo. Se prendiamo la velovità v0 negativa e
l’accelerazione a0 positiva, ottenendo i grafici seguenti.
v0
a0
x0
O t O t O t
−
→
g
O x
si ha
g = −g ĵ . (2.33)
19
L’equazione del moto per un punto materiale che sia soggetto alla sola forza di gravità è
quindi un moto uniformemente accelerato lungo l’asse delle y
ÿ = −g . (2.34)
Integriamo la 2.34. Si ottiene per la velocità
ẏ(t) = ẏ(0) − gt , (2.35)
e infine l’equazione oraria:
1
y(t) = y(0) + ẏ(0)t − gt2 . (2.36)
2
Vediamo adesso due applicazioni della (2.36).
Caso 1
Supponiamo di lasciar cadere a t = 0 il grave da fermo da un’altezza h. Ciò si traduce nelle
condizioni iniziali
ẏ(0) = 0 , (2.37)
y(0) = h . (2.38)
La soluzione particolare della (2.36) è quindi
1
y(t) = h − gt2 . (2.39)
2
Il grave arriverà a terra ad un istante t1 tale che
1
y(t1 ) = h − gt21 = 0 , (2.40)
2
ovvero s
2h
t1 = ± . (2.41)
g
La soluzione negativa non ha significato fisico (dall’istante t = 0 ammettiamo che il tempo
scorra per valori positivi) quindi abbiamo
s
2h
t1 = . (2.42)
g
La velocità del grave quando tocca terra sarà
p
ẏ(t1 ) = −gt1 = − 2hg , (2.43)
negativa, come deve essere nel SdR che abbiamo scelto.
y(t) ẏ(t)
t1
h
O t
√
− 2hg
O t1 t
20
Caso 2
Supponiamo adesso di lanciare il grave verso l’alto da y = 0 con una velocità iniziale pari a
v0 = v0 ĵ. Le condizioni iniziali sono quindi date da
ẏ(0) = v0 , (2.44)
y(0) = 0 . (2.45)
La soluzione particolare della (2.36) è quindi
1
y(t) = v0 t − gt2 . (2.46)
2
Quindi il grave torna al suolo al tempo t1 tale da annullare y(t1):
1
y(t1) = v0 t1 − gt21 = 0 , (2.47)
2
ovvero
1
t1 (v0 − gt1 ) = 0 , (2.48)
2
che ha soluzioni
2v0
t1 = 0 , .
t1 = (2.49)
g
La prima soluzione corrisponde al valore iniziale. Infatti abbiamo imposto che per t = 0 il
grave sia al suolo. La seconda soluzione, invece, corrisponde a momento in cui il grave ritocca
il suolo.
La velocità con cui torna a terra è
v(t1 ) = v0 − gt1 = −v0 , (2.50)
che in modulo è pari alla velocità iniziale, ma ovviamente di verso opposto.
Siccome ci sono due istanti successivi per i quali v(0) > 0 e v(t1 ) < 0 e la funzione v(t) è
continua (non stiamo considerando urti), si dovrà avere un istante t∗ , con 0 < t∗ < t1 , tale che
v(t∗ ) = 0. Questo è dato dalla
v(t∗ ) = v0 − gt∗ = 0 , (2.51)
cioè
v0
t∗ = . (2.52)
g
Per t∗ = vg0 , il grave raggiunge l’altezza massima h, pari a
1 v0 1 v 2 v2
h = y(t∗ ) = v0 t∗ − g(t∗)2 = v0 − g 02 = 0 . (2.53)
2 g 2 g 2g
21
2.5 Moto piano generico di un punto
Supponiamo adesso che il punto materiale P si muova nel piano. Il nostro SdR sarà dato da
una coppia di assi cartesiani con origine O:
y
SdR
P
r(t)
ĵ γ
O î x
k∆rk
kv(t)k = , (2.56)
∆t
e come direzione e verso quelle di ∆r (come mostrato in figura).
y
SdR ∆r(t)
P
r(t)
ĵ r(t + ∆t) γ
O î x
22
2.5.1 Derivata del vettore posizione
Per definire la velovità istantanea, dovremo fare il limite ∆t → 0. La velocità istantanea sarà
quindi definita come:
r(t + ∆t) − r(t)
v = ṙ = lim (2.57)
∆t→0 ∆t
È facile vedere come al tendere di ∆t → 0, il vettore v(t) tenda ad un vettore tangente alla
traiettoria nel punto r(t). Infatti, il vettore ∆r ha il modulo che tende a zero come ∆t, in
maniera tale che il loro rapporto tenda a kvk; la sua direzione invece tende alla tangente alla
traiettoria nel punto t (come si vede per costruzione in figura) e il verso è quello del moto.
∆r(t)
P
r(t)
γ
O r(t + ∆t)
Se t̂ è il versore tangente alla traiettoria γ nel punto r(t), si potrà esprimere la velocità
come
ṙ = kvkt̂ . (2.58)
Per capire meglio il significato di kvk, è conveniente utilizzare la parametrizzazione naturale
della curva γ.
Esempio
Per esempio consideriamo una circonferenza di raggio r. Un’equazione parametrica per questa
curva è data dalla
r(t) = (r cos (t), r sin (t)) , (2.60)
da cui
ṙ(t) = (−r sin (t), r cos (t)) (2.61)
e quindi
kṙ(t)k = r . (2.62)
23
La lunghezza della curva è quindi
Z 2π Z 2π
s(t = 2π) = kṙ(t)k dt = r dt = 2πr . (2.63)
0 0
Proprietà 2.6.1 La lunghezza di una curva NON dipende dalla scelta della sua parametrizza-
zione.
Infatti, sia τ un nuovo parametro tale che t = t(τ ) sia una funzione almeno C1 nell’intervallo
d’interesse (cioè sia invertibile). Allora r(t(τ )) = r(τ ) con t ∈ (a, b) e τ ∈ (a′ , b′ ) tali che
t(a′ ) = a e t(b′ ) = b. Allora
b b′ b′
dr dt dr
Z Z Z
s(t) = kṙ(t)kdt = dτ = dτ . (2.64)
a a′ dt dτ a′ dτ
Proprietà 2.6.2 Qualunque curva differenziabile e non singolare ammette una parametrizza-
zione “naturale” tramite la lunghezza d’arco.
Infatti, siccome Z t
s(t) = kṙ(t)kdt , (2.65)
t0
si avrà che:
ds
= ṙ(t) ≥ 0 . (2.66)
dt
In particolare ci sarà un intervallo intorno a t in cui
ds
= ṙ(t) > 0 , (2.67)
dt
ciaoè in cui s(t) è monotona crescente e quindi invertibile. L’inversa è
dt 1
= > 0. (2.68)
ds kṙ(t)k
che è tangente alla traiettoria per costruzione (ragionamento identico a quanto fatto per la
derivata rispetto a t) e ha banalmente modulo 1.
24
2.7 Velocità e accelerazione
Tramite la parametrizzazione naturale della traiettoria è possibile esprimere in maniera più
chiara la velocità v = ṙ(t) e anche l’accelerazione.
Parametrizzando la curva con s, invece che con t, si trova:
dr dr ds
v= = = ṡ t̂ . (2.71)
dt ds dt
La ṡ è detta velocità scalare e ci dice effettivamente che kvk è la velocità di spostamento lungo
la curva γ, ovvero “quanta lunghezza di curva P percorre nell’unità di tempo”. Se potessimo
“stirare” la curva γ ad una traiettoria rettilinea, la ṡ sarebbe proprio la ẋ definita nel paragrafo
del moto rettilineo.
Per quanto riguarda l’accelerazione, la situazione in due dimensioni invece è più complicata.
Abbiamo:
dv dt̂ dt̂ ds dt̂
a= = s̈ t̂ + ṡ = s̈ t̂ + ṡ = s̈ t̂ + ṡ2 . (2.72)
dt dt ds dt ds
Dobbiamo vedere a cosa corrisponde la derivata rispetto ad s del versore tangente t̂. Facendo
lo stesso ragionamento di quello fatto per t, si vede che ddst̂ è un vettore perpendicolare a t̂ che
punta verso la concavità della curva γ.
t̂(s0 )
t̂(s) t̂(s0 )
∆t̂
t̂(s)
γ
dt̂
= k(s)n̂(s) , (2.74)
ds
dove il fattore di proporzionalità k(s) è detto curvatura della curva γ in s ed è l’inverso del
raggio del cerchio osculatore di γ in s:
1
k(s) = (2.75)
ρ
ρ è detto raggio di curvatura.
25
t̂(s)
n̂(s)
ρ
γ
Infine si ha:
v = ṡ t̂ , (2.76)
2
ṡ
a = s̈ t̂ + n̂ . (2.77)
ρ
Si vede quindi che, mentre la velocità è tangente alla traiettoria, l’accelerazione ha due com-
ponenti: una componente tangente ed una normale, diretta verso la concavità della γ. La
componente tangenziale è quella che mi dice come varia la velocità lungo la curva. Se si pren-
desse il limite ρ → ∞, cioè se la traiettoria diventasse rettilinea (curvatura nulla), il termine s̈
sarebbe esattamente l’ẍ del moto retilineo che abbiamo trattato nei paragrafi precedenti. Il ter-
mine di curvatura, ci dice che per cambiare direzione nel piano c’è bisogno di una componente
“centripeta” dell’accelerazione.
L’Eq. (2.74) è la prima parte del teorema di Frenet.
Teorema 2.7.1 Per una curva differenziabile non singolare di parametrizzazione naturale s si
ha:
dt̂
= k(s)n̂(s) , (2.78)
ds
dn̂
= −k(s)t̂(s) . (2.79)
ds
Per dimostrare la seconda relazione, Eq. (2.79), basta notare che da n̂ · n̂ = 1 segue che
dn̂/ds ∝ t̂. Inoltre da n̂ · t̂ = 0 segue che
dn̂ dt̂
· t̂ = −n̂ · = −k(s) , (2.80)
ds ds
da cui la (2.79).
26
con r = cost e θ = θ(t). Il moto circolare ha quindi un grado di libertà, ovvero è completamente
determinato una volta che conosciamo l’equazione oraria θ = θ(t) (r è fissato).
La velocità di P sarà data da
di modulo q
kvk = r 2 θ̇2 sin2 (θ) + r 2 θ̇2 cos2 (θ) = r θ̇ . (2.83)
θ+
x
Quindi
dt̂ dt̂ dt̂
= (−θ̇ cos (θ), −θ̇ sin (θ)) = ṡ = r θ̇ . (2.88)
dt ds ds
Da questa relazione si deduce che
dt̂ 1 1 1
= − cos (θ), − − sin (θ) = n̂ . (2.89)
ds r r r
Quindi si ha
1
,
k(s) = (2.90)
r
ovvero, per il moto circolare su una circonferenza di raggio r si ha:
ρ = r. (2.91)
27
2.8 Moto piano in coordinate cartesiane
In coordinate cartesiane, come abbiamo già detto, il vettore posizione r(t) è individuato dalle
sue componenti (x(t), y(t))>
y
SdR
P
r(t)
ĵ γ
O î x
dr d
v(t) = = (x(t)î + y(t)ĵ) , (2.92)
dt dt
= ẋ(t)î + ẏ(t)ĵ , (2.93)
dv d
a(t) = = (ẋ(t)î + ẏ(t)ĵ) , (2.94)
dt dt
= ẍ(t)î + ÿ(t)ĵ . (2.95)
Possiamo quindi scomporre il moto del punto P lungo î e ĵ. I due moti sono “indipendenti” e
trattabili separatamente (vedi esempi di balistica).
y
ûθ
ûr
r(t)
P γ
θ
O x
28
Per individuare la parte vettoriale, introdurremo una coppia di versori mobili, solidali con P,
ûr e ûθ tali che ûr sia il versore prolungamento del raggio vettore e ûθ il versore perpendicolare
ad ûr ottenuto ruotando ûr di π2 in senso antiorario:
Come si vede, si ha
dûr
= −θ̇ sin θî + θ̇ cos θĵ = θ̇ûθ
dt
dûθ
= −θ̇ cos θî − θ̇ sin θĵ = −θ̇ûr (2.98)
dt
Infine si ha:
dr d
v = = (rûr ) = ṙûr + r θ̇ûθ , (2.99)
dt dt
dv
a = = r̈ûr + 2ṙθ̇ûθ + r θ̈ûθ − r θ̇2 ûr = (r̈ − r θ̇2 )ûr + (2ṙ θ̇ + r θ̈)ûθ . (2.100)
dt
θ+
x
v = r θ̇ ûθ ,
a = −r θ̇2 ûr + r θ̈ ûθ , (2.102)
che sono le (2.99,2.100) in cui abbiamo posto ṙ = 0. La velocità, quindi, è tangente alla
traiettoria e di modulo r θ̇, mentre l’accelerazione continua ad avere una parte tangente ed una
normale.
29
Moto circolare uniforme
Se il moto è uniforme, cioè se
θ̇ = ω = cost , (2.103)
che quindi, integrando direttamente, dà la seguente legge oraria
θ(t) = θ0 + ωt , (2.104)
v = r θ̇ ûθ = rω ûθ ,
a = −r θ̇2 ûr = −rω 2 ûr , (2.105)
con la velocità costante in modulo. L’accelerazione adesso è solo normale e diretta verso il
centro del cerchio. Per i moduli si ha:
r 2 θ̇2 kvk2
kak = r θ̇2 = = . (2.106)
r r
Moto periodico
Il moto circolare uniforme è un primo esempio di moto periodico, ovvero di un moto che dopo
un certo lasso di tempo, che chiameremo periodo e indicheremo con T , si riproduce uguale a sé
stesso.
Sia n ∈ N. Allora un moto periodico è caratterizzato dalla seguente equazione
Nel nostro caso particolare (moto circolare uniforme) ad ogni periodo T l’angolo θ ha un
incremento di 2π:
θ(t + T ) = θ(t) + 2π . (2.108)
Si definisce frequenza del moto periodico, l’inverso del periodo
1
ν= (2.109)
T
e pulsazione
2π
ω = 2πν = . (2.110)
T
Sia la frequenza che la pulsazione hanno le dimensioni dell’inverso di un tempo:
30
Quindi la velocità sarà data da:
dr dûr
v= =r = r θ̇ ûθ . (2.113)
dt dt
Il versore ûθ è perpendicolare al versore ûr , tale che ûr ruoti su ûθ di π/2 in direzione antioraria.
Ricordando le proprietà del prodotto vettoriale, si può introdurre un vettore che abbia modulo
pari a θ̇, direzione perpendicolare al piano {ûr , ûθ } e verso tale da vedere ûr andare su ûθ con
la regola della mano destra:
→
−ω = θ̇ n̂ , (2.114)
dove n̂ è il versore perpendicolare al piano {û , û }. Usando −
r θ
→
ω si può scrivere
dûr
= θ̇ ûθ = −
→
ω ∧ ûr . (2.115)
dt
Ovvero
dr
v= = r−→ω ∧ ûr = r θ̇n̂ ∧ ûr . (2.116)
dt
Per trovare →
−
ω basta moltiplicare vettorialmente la (2.115) per ûr :
dûr
ûr ∧ = ûr ∧ (−→ω ∧ ûr ) , (2.117)
dt
= | utilizzando la regola del doppio prodotto vettoriale, Cap. 1 |
= (ûr · ûr ) −
→
ω − (ûr · −
→
ω ) ûr = −
→
ω. (2.118)
−
→ d−
→ω dω
α = = −̇
→
ω = n̂ = θ̈ n̂ . (2.120)
dt dt
Per il moto circolare (r = cost) abbiamo quindi
dv d − d−
→ω dr
a = = (→
ω ∧ r) = ∧r+−
→
ω ∧ , (2.121)
dt dt dt dt
= −̇
→
ω ∧r+− →
ω ∧ (−
→
ω ∧ r) , (2.122)
→
− →
− →
−
= α ∧ r + ω ∧ ( ω ∧ r) . (2.123)
31
dove abbiamo posto ω = θ̇. Proiettiamo l’Eq. (2.125) sugli assi cartesiani. Si ottiene
ẍ + ω 2 x = 0 , (2.126)
ÿ + ω 2 y = 0 . (2.127)
Su entrambi gli assi cartesiani si compie lo stesso moto. L’equazione del secondo ordine
f¨ + ω 2f = 0 , (2.128)
si dice equazione del moto armonico.
È facile verificare che una soluzione generale della (2.128) è data dalla
f (t) = A cos (ω t + φ) . (2.129)
Infatti
f˙ = −A ω sin (ω t + φ) , (2.130)
f¨ = −A ω 2 cos (ω t + φ) = −ω 2 f . (2.131)
Le due costanti A e φ sono le costanti arbitrarie d’integrazione e dovranno essere fissate
imponendo le condizioni iniziali del moto. A è detta ampiezza e φ la fase. Supponiamo di avere
le seguenti condizioni iniziali: (
f (0) = f0 ,
(2.132)
f˙(0) = f˙0 .
Allora, considerando la soluzione generale (2.129), avremo
A cos φ = f0 , (2.133)
−Aω sin φ = f˙0 (2.134)
e, risolvendo il sistema, si trova !
f˙0
φ = − arctan (2.135)
ω f0
e considerando che dalle (2.134) si ottiene
f˙02
A2 sin2 φ + A2 cos2 φ = A2 = f02 + , (2.136)
ω2
si trova anche l’ampiezza s
f˙02
A= f02 + . (2.137)
ω2
La soluzione particolare del moto è data quindi da
s !!
f˙2 f˙0
f (t) = f02 + 02 cos ωt − arctan (2.138)
ω ω f0
La funzione f (t) è mostrata in figura. Il moto è periodico di periodo
2π
T = (2.139)
ω
e frequenza
1 ω
ν= = . (2.140)
T 2π
32
f (t)
T
A
0 t
−A
v0
ĵ θ
O î x
ẍ(t) = 0 , (2.144)
ÿ(t) = −g . (2.145)
x(0) = 0 , (2.146)
ẋ(0) = v0x = v0 cos θ , (2.147)
y(0) = 0 , (2.148)
33
ẏ(0) = v0y = v0 sin θ . (2.149)
dove abbiamo introdotto le prime due costanti arbitrarie ẋ(0) e ẏ(0), che possono essere fissate
usando le (2.147,2.149).
Integrando un’altra volta otteniamo
Z t
x(t) = x(0) + ẋ(0)dt = x(0) + ẋ(0)t , (2.152)
0
Z t
1
y(t) = y(0) + (ẏ(0) − gt)dt = y(0) + ẏ(0)t − gt2 . (2.153)
0 2
Imponendo quindi le condizioni iniziali, si hanno le seguenti soluzioni particolari delle
equazioni del moto:
Traiettoria
Le (2.154,2.155) costutiscono una forma parametrica della traiettoria del punto. Possiamo
trovare una forma “esplicita” come y = y(x) ricavando il parametro da una delle due equazioni
e sostituendolo nell’altra:
x
t = , (2.156)
v0 cos θ
g
y = − 2 x2 + tan (θ) x . (2.157)
2v0 cos2 θ
ĵ
h
O î x0 x1 x
34
Il punto in cui P tocca di nuovo il suolo è dato dal’equazione y(x) = 0:
g
x − 2 x + tan (θ) = 0 , (2.158)
2v0 cos2 θ
che ha per soluzioni (0, 0), punto da cui P parte, e (x1 , 0) con
2v02 v2
x1 = sin θ cos θ = 0 sin (2θ) . (2.159)
g g
La quota massima raggiunta nel moto può essere trovata estremando la y = y(x). Si ha un
massimo per
dy g
=0=− 2 x + tan (θ) , (2.160)
dx v0 cos2 θ
da cui
v2
x = x0 = 0 sin θ cos θ , (2.161)
g
che sostituita nella (2.157) dà
v02
h = y(x0 ) = ... = sin2 θ . (2.162)
2g
Gittata
Un’altra possibile domanda che ci possiamo porre è quale sia la gittata (cioè l’x1 ) massima, a
parità di v0 , in funzione di θ.
Dalla (2.159), estremando rispetto a θ si può cercare il massimo:
dx1 v2
= 2 0 cos (2θ) = 0 , (2.163)
dθ g
che dà la seguente soluzione (ci possiamo limitare, per la geometria del problema, a cercare una
soluzione nell’intervallo 0 < θ < π/2):
π
θ= . (2.164)
4
Ciò vuol dire che se v0 6= 0, la gittata massima si ottiene sempre sparando a 45◦ , indipenden-
temente dalla velocità iniziale del proiettile.
Alzata
Chiediamoci adesso quale sia l’alzata, ovvero l’angolo θ, che mi permetta di colpire il punto di
coordinate (l, 0).
Imponendo il passaggio della traiettoria da questo punto, troviamo
g 2 l gl
0=− 2 l + tan (θ) l = 2 sin θ cos θ − 2 . (2.165)
2v0 cos2 θ 2 cos2 θ v0
Se θ 6= π/2, cosa che è ben verificata, visto che se spariamo il proiettile in verticale non
abbiamo nessuna possibilità di prendere il bersaglio a (l, 0), allora cos θ 6= 0 e quindi la (2.165)
è soddisfatta se
gl
sin (2θ) = 2 . (2.166)
v0
35
La (2.166) ammette soluzioni se
gl
≤ 1, (2.167)
v02
ovvero, se, data una certa distanza l alla quale
√ vogliamo colpire il bersaglio, la velocità iniziale
del proiettile√è “sufficientemente alta (v0 ≥ lg).
Se v0 = lg, la√(2.166) ha un’unica soluzione, per θ = π/4.
Se invece v0 > lg, si hanno due soluzioni:
1 gl
θ1 = arcsin , (2.168)
2 v02
π
θ2 = − θ1 . (2.169)
2
Il fatto che ci siano due soluzioni, ci dice che se la velocità iniziale del proiettile è sufficien-
temente alta, possiamo colpire il bersaglio sparando con un’alzata inferiore ai 45◦ (in tal caso
la traiettoria si dirà diretta) oppure con un’alzata superiore a tale angolo limite (traiettoria
inversa).
Tempi
Quanto tempo impiega il proiettile a raggiungere di nuovo terra?
Utilizzando l’Eq. (2.155) e imponendo y(t) = 0, si trova l’istante desiderato:
1
t − gt + v0 sin θ = 0 , (2.170)
2
Avremmo potuto anche considerare il moto uniforme lungo le x, Eq. (2.154), e imporre
x = x1 trovando
x1 2v0
t= = sin θ , (2.172)
v0 cos θ g
dove abbiamo usato la (2.159).
36
rO (t), poi, individua l’origine del SdR mobile nel SdR′ fisso. Il tempo (assoluto) scorre allo
stesso modo per i due sistemi di riferimento.
z SdR
′
z
P y
SdR’ r(t) k̂
r′ (t)
O ĵ
k̂′ î
rO (t)
O ′ ĵ′ y′
′
î x
x′
Avremo:
r′ (t) = rO (t) + r(t) , (2.173)
dove quindi
rO (t) = xO (t)î′ + yO (t)ĵ′ + z0 (t)k̂′ , (2.174)
r(t) = x(t)î + y(t)ĵ + z(t)k̂ . (2.175)
Ricordiamo che i versori solidali al SdR mobile sono funzioni del tempo e cambiano con
esso la loro posizione, rispetto a SdR′ .
Cerchiamo una relazione che leghi fra loro la velocità che misurerebbe l’osservatore nel SdR′
fisso e quella che misurerebbe l’osservatore nel SdR′ mobile. Derivando rispetto al tempo la
relazione (2.173) si ottiene
d dî
(î · î) = 2 · î = 0 (2.178)
dt dt
e questo vale anche per ĵ e k̂. Cioè le derivate temporali dei tre versori sono rispettivamente
ortogonali ai versori stessi. Possiamo, allora, introdurre tre vettori, −
→
ω 1, −
→
ω 2 ed −
→
ω 3 , tali che:
37
Teorema 2.13.1 (di Poisson) Esiste uno e un solo →
−
ω tale che →
−
ω1 =−
→
ω2 = −
→
ω3 = −
→
ω.
Per dimostrare il teorema dobbiamo mostrare che si può trovare tale → −
ω e che questo è unico.
Notiamo che
î · ĵ = 0 , ĵ · k̂ = 0 , k̂ · î = 0 , (2.180)
che, derivate, danno le
ω 1 ∧ î · ĵ = −
−
→ →
ω 2 ∧ î · ĵ , −
→
ω 2 ∧ ĵ · k̂ = −
→
ω 3 ∧ ĵ · k̂ , −
→
ω 3 ∧ k̂ · î = −
→
ω 1 ∧ k̂ · î . (2.182)
Da cui si ottiene
ω1,3 = ω2,3 = ω3 . (2.186)
Dalla seconda delle (2.182), analogamente, si ottiene la
e dalla terza
ω3,2 = ω1,2 = ω2 . (2.188)
Quindi i tre vettori si possono scrivere come segue:
−
→
ω 1 = ω1,1 î + ω2 ĵ + ω3 k̂ , (2.189)
→
−
ω 2 = ω1 î + ω2,2 ĵ + ω3 k̂ , (2.190)
→
−
ω = ω î + ω ĵ + ω k̂ . (2.191)
3 1 2 3,3
Rimangono da determinare le tre componenti ω1,1 , ω2,2 e ω3,3 , che non sono vincolate da nessuna
relazione. Quindi potremo scegliere
trovando così →
−
ω1 = −
→
ω2 =−
→
ω3 = −
→
ω , dove
−
→
ω = ω1 î + ω2 ĵ + ω3 k̂ . (2.193)
Ponendo un solo →
−
ω nelle (2.179), si ottengono le così dette formule di Poisson:
dî
= →
−
ω ∧ î , (2.194)
dt
38
dĵ
= −
→
ω ∧ ĵ , (2.195)
dt
dk̂
= →
−
ω ∧ k̂ . (2.196)
dt
Il vettore −
→
ω si dice velocità angolare di SdR rispetto ad SdR′ . Per ogni versore ŵ solidale con
SdR si ha che
dŵ −
=→ω ∧ ŵ . (2.197)
dt
Se −→
ω non fosse unico, si potrebbe definire un −
→ω ′ tale che valga anche
dŵ −
=→
ω ′ ∧ ŵ . (2.198)
dt
Ma allora
(−
→
ω −−
→
ω ′ ) ∧ ŵ = 0 , (2.199)
e siccome ŵ è arbitrario, ne segue che ciò è verificato solo se (−
→
ω −−
→
ω ′ ) = 0, ovvero se −
→
ω =− →
ω ′.
Ritornando alla formula per la velocità, allora, si ha:
v = vO + vR + −
→
ω ∧r (2.200)
Si usa porre
vt = vO + −
→
ω ∧r. (2.201)
vt viene detta velocità di trascinamento di P. In tal modo si può porre
v = vt + vR . (2.202)
Deriviamo la formula per le velocità e otteniamo una relazione per le accelerazioni. Si tro-
verà:
a = aO + aR + 2−
→
ω ∧ vR + −̇
→
ω ∧r+−
→
ω ∧ (−
→
ω ∧ r) . (2.203)
Si può scrivere, allora
a = aR + aC + at , (2.204)
dove abbiamo posto
aC = 2−
→
ω ∧ vR , (2.205)
detta accelerazione di Coriolis e
at = ao + −̇
→
ω ∧r+−
→
ω ∧ (ω ∧ r) , (2.206)
39
2.13.1 Esempio: treno-binario
Due sistemi di riferimento traslano l’uno rispetto all’altro. Si chiede come descrivere il moto
del punto P , soggetto all’accelerazione di gravità, nei due SdR.
Supponiamo che in SdR′ (“fisso”) il punto P si muova con la seguente equazione del moto:
(
ÿ ′ = −g ,
(2.207)
ẍ′ = 0 .
Essendo −→
ω = 0, in questo caso, la relazione che lega il moto di P in sDR′ al moto in SdR è
a′ = a0 + aR . (2.208)
Allora nel caso in cui si abbia semplicemente a0 = a0 ß̂ (traslazione lungo la direzione x), si
avrà
aR = a′ − a0 . (2.209)
Quindi in SdR si ha (
ÿ = −g ,
(2.210)
ẍ = −a0 .
Supponiamo ora che a0 = 0, ma v0 6= 0. Il moto ne risulterebbe immutato:
(
ÿ = −g ,
(2.211)
ẍ = 0 .
(Questo porta alla Legge d’Inerzia, di Galileo).
z′ z −
→
ω
R h y
0 P
x
x′
Supponiamo che un SdR’, con origine nel centro della Terra e asse z ′ verso il Polo Nord e
altri due assi x′ e y ′ sul piano equatoriale, sia inerziale. Prendiamo un SdR solidale con un
punto O dell’equatore (quindi in rotazione con la Terra) con asse z parallelo a z ′ , x tangente
all’equatore nel punto O e y perpendicolare alla superficie e quindi lungo il raggio terrestre.
Consideriamo la velocità angolare della terra costante (quindi −̇→
ω = 0), tale che
−
→
ω = ωk , (2.220)
con
2π
ω= = 7.3 10−5s−1 . (2.221)
86400s
In SdR’, P è soggetto alla forza di gravitazione
F = −mGr̂ , (2.222)
dove
MT
G=G , (2.223)
(R + y)2
con MT massa della Terra, R il suo raggio, G costante di gravitazione universale e y altezza del
punto materiale sulla superficie terrestre.
Scriviamo le equazioni del moto nel SdR in rotazione. Si avrà
dove
aC = 2−
→
ω ∧ vR , (2.225)
= 2ω k̂ ∧ (ẋi + ẏj + żk) , (2.226)
41
= 2ω ẋj − 2ω ẏi , (2.227)
at = a0 + − →
ω ∧ (−→ω ∧ r) , (2.228)
2 2
= −ω Rj + ω k ∧ (k ∧ (xi + yj + zk)) , (2.229)
= −ω 2 (R + y) j − ω 2 x i . (2.230)
Le equazioni del moto in SdR sono quindi date dal seguente sistema:
2
ẍ = 2ω ẏ + ω x ,
ÿ = −G − 2ω ẋ + ω 2(R + y) , (2.231)
z̈ = 0 ,
2. Inoltre, abbiamo che g = 9.8ms−2 , mentre il termine ω ẋ è dell’ordine di 10−4 ms−2 . Quindi
ω ẋ ≪ g e possiamo ulteriormente semplificare l’equazione differenziale:
ÿ = −g , (2.236)
con soluzione
1
y(t) = − gt2 + h . (2.237)
2
3. Consideriamo l’equazione per ẍ. Si ha anche in questo caso una gerarchia fra i termini di
destra: ω 2 x ≪ ωgt. Infatti, considerando spostamenti finali dell’ordine dei centimetri,
ω 2x ≃ 10−11 ms−2 , (2.238)
mentre per tempi dell’ordine dei 10s,
ωgt ≃ 10−2 ms−2 . (2.239)
Quindi l’equazione da risolvere è
ẍ = −2ωgt , (2.240)
con soluzione
1
x(t) = − ωgt3 . (2.241)
3
42
Se tF è il tempo di caduta del grave, si ha
s
1 2h
y(tF ) = 0 = − gt2F + h → tF = (2.242)
2 g
e quindi
32
1 2h
x(tF ) = − ωg , (2.243)
3 g
che, per h = 100m, dà x(tF ) = 2.5 · 10−2 m.
43
Capitolo 3
In questo capitolo definiamo e studiamo alcune grandezze fondamentali che poi ritroveremo
nello studio della dimanica vera e propria, ovvero nello studio del moto dei corpi sottoposti
ad interazione con l’ambiente o con altri corpi. L’interazione sarà la causa del moto. Data
l’interazione si può prevedere il moto del sistema fisico in esame. Vice versa dato il moto si
può cercare di risalire alla forma dell’interazione. L’interazione si esplica quantitativamente
attraverso l’introduzione di una quantità fisica detta “forza”, che studieremo in dettaglio in
questo capitolo.
3.1 Forze
Il concetto di forza si può empiricamente far risalire a quello dello sforzo muscolare che compia-
mo nel cercare di mutare lo stato (configurazione) di un corpo. Da questo dobbiamo astrarre
ciò che ci interessa da un punto di vista fisico, cioè:
Graficamente tutto viene riassunto dalla freccetta in figura, che avrà lunghezza proporzionale
al modulo, retta di applicazione individuata dal prolungamento tratteggiato, verso indicato
dalla punta della freccia e punto di applicazione P.
44
L’altro fatto che evidenzia il comportamento vettoriale di una forza è la composizione di
più forze agenti sullo stesso punto di applicazione. Infatti: due forze si sommano con la regola
del parallelogramma, cioè se su P agiscono le due forze F1 ed F2 , la loro azione è equivalente
all’azione di una forza F = F1 + F2 che sia la composizione delle due secondo la regola del
parallelogramma.
F1 F
P
F2
Questo si può verificare sperimentalmente! (Esperienza dei cavalli o dei tre dinamometri).
Definizione 3.1.1 Si dice che il punto P è fermo (nel SdR considerato), se la sua configura-
zione non muta nel tempo.
Se sul punto P agiscono due forze lungo la stessa direzione e con verso opposto, tali che P sia
fermo, allora si dice che le due forze si equilibrano, hanno la stessa intensità (proprietà intuitiva,
che deriva dall’esperienza).
Il secondo concetto riguarda una caratteristica importante di alcuni corpi quando sono
soggetti ad una forza; una caratteristica che mette in evidenza l’azione della forza stessa e
che quindi può essere utilizzata per misurarla. È esperieza quotidiana il fatto che se una forza
agisce su un corpo reale, può deformarlo. Si pensi per esempio all’azione di una forza che agisca
sull’estremità di una sbarretta metallica fissata, dall’altra estremità, ad una morsa. Spingendo
in direzione ortigonale alla lunghezza della sbarretta, ci accorgiamo che la barretta si flette. La
stessa cosa avviene se consideriamo una molla fissata ad una sua estremità. Agendo con una
forza all’altra estremità, possiamo allungare (o comprimere) la molla. Alcuni corpi hanno la
proprietà importante di tornare alla loro configurazione iniziale una volta che la sollecitazione
(forza applicata) cessi (molti corpi reali hanno un comportamento elastico purché l’intensità
della forza applicata sia modesta e non deformi il corpo stesso). Questi corpi si dicono elastici.
Questa proprietà, insieme al concetto di equilibrio introdotto sopra, può essere utilizzata per
misurare la forza.
Consideriamo, infatti, una molla, attaccata al muro da una estremità, e al punto P dall’altra
estremità:
P F
O x0 x1 x
45
Una forza F agisce sul punto P , tirandolo nella direzione positiva delle x. Se F = 0, la
molla non sarà sollecitata e il punto P sarà fermo all’ascissa x = x0 . Applichiamo adesso F 6= 0
(forza di piccola intensità, tale da garantire il comportamento elastico della molla). La molla
si allungherà fino a quando non si sarà trovata una nuova situazione di equilibrio in x = x1 . La
forza F continua ad agire su P , ma un’altra forza di uguale intensità, stessa direzione e verso
opposto, f, agisce su P in maniera tale da farlo stare fermo in x = x1 , equilibrando F.
L’elongazione della molla, (x1 − x0 ), dà una misura della forza F applicata in P . Infatti,
agendo su P con una forza più intensa, F1 , mi accorgerei che il punto di equilibrio trovato
nella nuova situazione sarebbe x = x2 , tale che (x2 − x0 ) > (x1 − x0 ). Ponendo di nuovo
F = 0 riotterrei l’equilibrio in x = x0 . Allora, ad elongazione nulla corrisponde una forza di
intensità nulla, mentre ad elongazione non nulla corrisponde una forza di intensità dipendente
dall’elongazione stessa. Si potrebbe in prima approssimazione pensare ad una dipendenza
lineare dall’elongazione: kFk ∝ (x − x0 ).
Per andare oltre nel nostro raionamento abbiamo bisogno di introdurre la forza peso. Li-
mitiamoci al campo di gravità. È un fatto sperimentale che un corpo vicino alla superficie
terrestre (nel laboratorio) sia soggetto ad una forza che lo attrae verso il basso. Per un corpo
omogeneo, questa forza dipende dal volume del corpo stesso. Per corpi di natura chimica di-
versa si vede che dipende anche dalla natura stessa del corpo (se prendiamo lo stesso volume di
due corpi diversi la forza peso è diversa). Sappiamo inoltre che tutti i corpi cadono nel campo
di gravità con la stessa accelerazione g = 9.81 m/s2 . Questi fatti sperimentali si combinano
nella legge matematica che individua la struttura della forza peso. La forza peso agente su
un corpo risulta essere proporzionale alla accelerazione di gravità tramite una costante che si
chiama massa gravitazionale del corpo stesso:
P = mg g . (3.1)
Sapendo che per un corpo omogeneo il peso è proporzionale al volume, possiamo agevolmente
definire multipli e sottomultipli di una dato peso P,
che dovrà equilibrare la forza di richiamo della molla,
proporzionale all’elongazione della stessa.
Su questi principi, possiamo costruire uno strumento
di misura della forza, detto dinamometro, mostrato in molla interna
figura. Si definisce un’unità di forza, data per esem-
pio dalla forza che equilibra la forza peso di un certo
campione, e si procede a definire una scala lineare di
misura, tarata con i multipli e sottomultipli dell’unità.
L’unità di misura della forza è il Newton (N), che
però nel Sistema Internazionale (SI) si definisce in ma-
scala graduata
niera dinamica attraverso il Secondo Principio, come lineare
grandezza derivata (cioè in dipendenza dalle grandezze
fondamentali lunghezza, massa e tempo).
Col dinamometro è possibile provare la legge di
composizione delle forze (parallelogramma) in manie- F
ra quantitativa. Infatti, collegando tre dinamometri
allo stesso punto di applicazione possiamo cambiare P
intensità e direzione/verso alle tre forze
Figura 3.1: Dinamometro.
46
3.2 Macchine
Per poter fare esperimenti con le forze, cambiare punto
di applicazione, direzione, verso e anche intensità usiamo delle macchine.
Le due macchine semplici sono la carrucola e la leva.
F2 b = F1 a . (3.2) F2
Nel caso per esempio che b > a, il modulo F1
della forza F2 viene “moltiplicato” per il rap-
porto b/a. Quindi con la leva è possibile fare
equilibrio ad una forza più intensa.
cioè che la risultante delle forze applicate sul punto sia nulla.
47
La (3.3) è detta Prima Cardinale della statica del punto materiale.
3.4 Vincoli
Un elemento importante nello studio del comportamento meccanico di un sistema è il concetto
di vincolo.
Definizione 3.4.1 Si dice vincolo una qualunque limitazione imposta alle configurazioni che
può assumere il mio sistema fisico nello spazio.
Ciò che ci interesa dei vincoli è, ovviamente, il loro comportamento fisico: il vincolo è capace
di applicare una forza sul mio sistema, che si dirà reazione vincolare (reazione, perché viene
in conseguenza all’azione del mio sistema sul vincolo, ma dobbiamo considerare sempre solo il
sistema che isoliamo e quindi ci interessa la reazione vincolare). La forza impressa dal vincolo
sul sistema è tale da impedire al sistema le configurazioni proibite dal vincolo.
La reazione vincolare R può essere scomposta in una reazione normale, Rn , ed una tangen-
ziale, Rt .
Rn R
P Rt
• Se, in ogni caso, si ha Rt = 0, il vincolo è detto liscio (è capace solo di opporre resistenza
normale alla superficie vincolare).
P N = −P . (3.4)
48
3.4.2 Vincoli con attrito
Attrito statico
Consideriamo, prima di tutto, l’attrito radente, cioè l’attrito che si ha per sfregamento contro
una superficie scabra. Il vincolo con attrito sarà un vincolo con Rt 6= 0. Vediamo, allora, come si
comporta questa forza, sempre nel caso statico. L’attrito è stato studiato fenomenologicamente
da Coulomb. Si vede che:
• Se cerchiamo di mettere in moto il punto facendolo strisciare sul piano scabro, la forza di
attrito si oppone alla possibile direzione del moto, ovvero alla forza che imprimiamo sul
punto: ha la stessa direzione, verso opposto e stessa intensità.
• Variando l’intensità della forza applicata al punto, ci accorgiamo che c’e’ un intervallo di
tale intensità per il quale la forza di attrito riesce a contrastare la forza applicata e si ha
equilibrio. Superato un certo valore, il corpo si mette in moto.
• Se per esempio consideriamo il caso di un punto materiale (quindi soggetto alla forza di
gravità) su un ripiano scabro e cerchiamo di metterlo in moto, ci accorgiamo che anche il
peso del corpo gioca un ruolo in questo processo. Se il corpo è più pesante, la forza che
dobbiamo applicare per metterlo in moto ha un’intensità maggiore.
Allora, possiamo caratterizzare la forza d’attrito statico, dicendo che tale forza ha stessa di-
rezione e verso contrario alla forza che agisce sul punto. Inoltre finché c’è equilibrio la sua
intensità è pari a quella della forza applicata, fino ad un valore limite. Questo comportamento
è descritto dalla seguente relazione:
Attrito dinamico
Nel caso dinamico, abbiamo che la forza di attrito si oppone al moto con una forza Rt tale che
abbia la direzione della velocità, verso opposto ad essa e modulo pari a
3.4.3 Esempio
Per illustrare meglio la natura dell’attrito statico consideriamo il seguente esempio. Su un
piano scabro ad inclinazione variabile α poniamo un punto materiale di massa m e studiamone
l’equilibrio al variare di α.
49
y Se α = 0, l’unica sollecitazione che riceve il vincolo
Rn è lungo la normale alla superficie stessa, quindi non ci
sono sforzi tangenziali e si ha semplicemente
O x Rn = mg . (3.7)
P
Cominciamo adesso ad aumentare l’inclinazione del
piano rispetto all’orizzontale.
Possiamo scomporre le forze nel nostro SdR, lungo
y i due assi cartesiani. Imponendo l’equilibrio si ha:
Rn
P sin θ − Rt = 0 , (3.8)
Rt
O P cos θ − Rn = 0 . (3.9)
α Si può avere equilibrio finché Rt ≤ µs Rn , ovvero,
x siccome dalla (3.8) si ha Rt = P sin θ, troviamo la
P
condizione
P sin θ ≤ µs P cos θ , (3.10)
ovvero
tan θ ≤ µs =⇒ α ≤ arctan (µs ) . (3.11)
Se α supera il valore limite α0 = arctan (µs ), il punto si mette in movimento. Per esempio se
µs ∼ 0.2, risulta α0 = 11◦ .
50
Capitolo 4
In questo capitolo ci occuperemo delle cause del moto del punto materiale, cioè dello studio della
relazione che intercorre fra forze e funzioni cinematiche. La dinamica del punto materiale, come
la dinamica di un qualsiasi sistema, è basata sui tre Principi di Newton, che ora brevemente
introdurremo e commenteremo.
Principio 4.1.1 Ogni corpo persevera nel suo stato di quiete o di moto rettilineo
uniforme, se non è soggetto a forze che mutino il suo stato.
Prima di Newton, il Principio d’Inerzia era già stato studiato da Galileo, in netto contrasto
con le “credenze” scientifiche del suo tempo, che volevano il movimento di un oggetto legato
alla “vis insita movendi”, una volta esaurita la quale, il corpo si ferma (a Galileo il merito di
aver riconosciuto il carattere accidentale degli attriti). Il principio d’inerzia è importante anche
perché permette di introdurre un problema teorico di notevole interesse, quello relativo ad una
classe di sistemi di riferimento, nei quali si compiono le misure, che vengono detti inerziali. Sono
definiti come quella classe di SdR in cui vale esattamente il principio d’inerzia e sono individuati
dal fatto che, come vedremo meglio dopo l’enunciato del Secondo Principio di Newton, facendo
una trasformazione che somma alla velocità del punto materiale una velocità costante, non
cambia la fisica del problema (rimane immutata la seconda legge della dinamica).
51
In formule:
F = ma , (4.1)
con la quale si può risolvere un qualsiasi problema di meccanica!
La (4.1) è un’equazione differenziale del secondo ordine, integrata la quale, con le condi-
zioni iniziali s(0) = so e v(0) = vo , si trovano traiettoria ed equazione oraria (cioè si risolve
completamente il problema dello studio del moto).
Per integrare la (4.1) abbiamo bisogno di conoscere la struttura della forza, ovvero la
funzione
F = F(x, ẋ, t) , (4.2)
ovvero come la forza dipende dalla posizione, o dalla velocità (per esempio la forza di attrito
viscoso), o anche dal tempo, in maniera esplicita.
Viceversa, conoscendo il moto (cioè la traiettoria e la legge oraria) è possibile risalire alla
foma dell’interazione.
Nell’enunciato del Secondo Principio compare la costante m, detta massa inerziale del punto
a cui è applicata la forza F. La massa inerziale è una grandezza estensiva, cioè per un corpo
omogeneo si vede che è proporzionale al volume del corpo. Inoltre, la massa di un corpo è
da considerarsi una costante (cioè non dipende dal tempo) qualunque sia lo stato di moto del
corpo.
Principio 4.3.1 Per ogni azione esiste sempre una reazione uguale e contraria;
cosicché le mutue azioni di due corpi sono sempre uguali e di verso opposto.
Si può dire meglio come segue: se un corpo A agisce su un corpo B per mezzo di una forza
F, contemporaneamente il corpo B agisce su A con -F, tale che il sistema (F, -F) costituisca
una coppia di braccio nullo e risultante nulla.
Definizione 4.4.1 Si dice inerziale un SdR in cui valgano esattamente le leggi di Newton,
nella forma in cui le abbiamo esposte. In particolare, si dice inerziale un SdR in cui valga il
Principio d’Inerzia.
Per fare un esempio, consideriamo un punto materiale nello spazio, sufficientemente lontano
da altri corpi in modo tale che si possa considerare isolato (F = 0). Noi osservatori saremo
in un SdR inerziale se il punto, che non è soggetto a forze, resterà fermo o si muoverà di
moto rettilineo uniforme. Supponiamo di aver individuato un tale SdR. Consideriamo, adesso,
un’altro SdR che ruoti rispetto a quello inerziale con velocità angolare ω. Se nel SdR inerziale
il punto materiale era fermo, adesso, nel nuovo SdR, descriverà un’orbita circolare alla velocità
angolare −ω. In questo SdR il punto materiale, non soggetto a forze, si muove di moto circolare
52
uniforme. È quindi soggetto ad un’accelerazione centripeta che non ha una causa “reale”, poiché
F = 0. Il SdR così individuato non è inerziale, poiché non vi valgono le leggi di Newton, il
principio d’inerzia in testa.
Esistono infiniti SdR inerziali. In particolare, una volta definito un SdR inerziale, qualunque
altro sistema di riferimento che si muova di moto traslatorio a velocità costante v rispetto al
primo, sarà anch’esso inerziale.
Prendiamo due sistemi di riferimento, come in figura, e supponiamo che il secondo trasli di
velocità costante v rispetto al primo. Prendiamo per semplicità gli assi del secondo SdR paralleli
a quelli del primo (al massimo dovremo fare una trasformazione ortogonale per riportarci a
questo caso, ma una tale trasformazione in ogni caso è indipendente dal tempo e non inficia il
seguente ragionamento).
z
SdR
z′ P
SdR’ r(t)
r′ (t)
O y
rO (t) = vt
O′ y′ x
x′
Per l’osservatore in SdR’, il punto materiale P è individuato dal raggio vettore r′ (t). Per
l’osservatore in SdR, invece, da r(t). La relazione che lega r′ (t) a r(t) è:
r′ = r + vt , (4.3)
con v = cost.
Se il punto materiale P si muove di velocità costante in SdR’
′ d 2 r′ d
a = 2 = (ṙ + v) = r̈ = a . (4.6)
dt dt
Quindi, nei due sistemi di riferimento l’accelerazione è la stessa. Se ne deduce che in SdR e
SdR’ vale il Secondo Principio nella sua forma semplice enunciata da Newton, F = ma. La
fisica nei due SdR è la stessa!
53
Il ragionamento ora fatto si può ripetere infinite volte, per cui potremo affermare che: Sia
S un SdR inerziale: allora è inerziale un qualunque altro SdR che trasli di moto rettilineo
uniforme rispetto ad S. Studiare uno stesso fenomeno meccanico in due SdR inerziali, dà gli
stessi identici risultati (a meno delle condizioni iniziali).
La discussione fatta finora può essere riespressa nel modo seguente. Definiamo la classe di
trasformazioni (4.3):
Definizione 4.4.2 Le trasformazioni r → r′ = r + vt con v = cost sono dette Trasformazioni
di Galileo.
Allora, la teoria, che è rappresentata dalla legge F = ma, è invariante sotto trasformazioni di
Galileo.
Questa affermazione, che riassume tutto il discorso sui SdR inerziali, è detta Principio di
Relatività Galileiana.
Misura dinamica
Per misurare la massa inerziale, possiamo utilizzare lo stesso Secondo Principio. Infatti, sia
data la forza costante F. Sottoponiamo due corpi A e B alla stessa forza F e supponiamo di
farli muovere nello stesso modo, ovvero per esempio lungo la stessa direzione e verso della F.
Misuriamo i moduli delle accelerazioni aA e aB . Siccome la forza è la stessa, si avrà
mA aB
= . (4.7)
mB aA
Possiamo allora definire l’unità di massa mB e musirare in unità mB la massa mA
aB
mA = mB . (4.8)
aA
Abbiamo definito una procedura dinamica di misura della massa inerziale. Inoltre, potremo
studiare anche le proprietà della massa inerziale, come per esempio il fatto che la massa è uno
scalare. Infatti, supponiamo che applicando F alla massa mA si misuri l’accelerazione aA tale
che:
|F|
mA = , (4.9)
aA
e l’analoga relazione per mB . Allora, se uniamo i due corpi a formare un corpo unico e vi
applichiamo ancora una volta F, misureremo un’accelerazione aC tale che
|F| |F| |F|
mC = = + = mA + mB . (4.10)
aC aA aB
Benché in linea di principio una misura dinamica della massa sia possibile, in pratica non
è ben realizzabile. In particolare, la precisione che può essere raggiunta non è molto grande
e anche la riproducibilità della misura vacilla. Per questo si può cercare di pensare ad una
procedura di misura statica.
54
Misura statica
È un fatto sperimentale che tutti i corpi cadono nel campo di gravità con la stessa accelerazione
g. Allora, per il Secondo Principio si avrà:
P = mg , (4.11)
dove P è la forza peso che agisce sull’oggetto di massa m, attraendolo verso il suolo.
Utilizzando una bilancia a due bracci (vedi sezione 3.2), possiamo pensare di misurare in
maniera statica la massa. Infatti, poniamo la massa da misurare m su uno dei due bracci e
poniamo dei campioni di massa sull’altro, in maniera tale da raggiungere l’equilibrio. Sappiamo
allora che all’equilibrio
P a = P1 b , (4.12)
dove a e b sono i bracci della bilancia. Se costruiamo una bilancia con a = b, avremo:
P = P1 , ovvero mg = m1 g , (4.13)
m = m1 . (4.14)
55
(g)
che non dipende dal corpo B. Il Secondo Principio ci dice che, indicando con mB la massa
(I)
gravitazionale e con mB la massa inerziale di B
(I) (g)
mB a = −mB g , (4.18)
nel nostro SdR con asse z positiva verso l’alto. Siccome sperimentalmente a = g per tutti i
corpi, ne deriva che per tutti i corpi si ha anche
Questa relazione deriva quindi dall’esperimento. Attualmente è verificata con una precisione
molto spinta, di una parte su 1012 .
Quindi, si può misurare la massa gravitazionale nelle stesse unità di misura della massa
inerziale.
L’uguaglianza fra massa gravitazionale e inerziale è detto Principio di Equivalenza ed è
carico di conseguenze fisiche (relatività di Einstein).
Oppure, come faremo nel capitolo della Gravitazione Universale, dato il moto di un pianeta
intorno al Sole potremmo risalire alla forma della forza ed ottenere esattamente la (4.15).
1
ẍ = Fx (x, ẋ, t) , (4.20)
m
1
ÿ = Fy (x, ẋ, t) , (4.21)
m
1
z̈ = Fz (x, ẋ, t) , (4.22)
m
con le condizioni iniziali
56
y(0) = y0 ẏ(0) = ẏ0 , (4.24)
z(0) = z0 ż(0) = ż0 . (4.25)
Dalla teoria delle equazioni differenziali ordinarie, sappiamo che un tale sistema ammette, date
le condizioni iniziali (4.23,4.24,4.25), un’unica soluzione. In questo senso la fisica classica è
completamente deterministica: dati i valori della posizione e della velocità ad un certo istante,
che noi abbiamo fatto coincidere con l’origine dei tempi (t = 0), sapendo la natura (cioè
struttura matematica) delle forze che agiscono sul sistema, c’è un solo moto che soddisfa F =
ma.
F = ṗ , (4.27)
con F forza applicata al punto materiale. Se il mio sistema (il punto materiale in questo caso)
è isolato, si ha F = 0 e quindi
ṗ = 0 , (4.28)
che ha come soluzione
p = cost . (4.29)
Per sistemi isolati la quantità di moto si conserva (nel caso del punto materiale è una banale
conseguenza del fatto che se si annulla la forza, il punto si muove di moto rettilineo uniforme
v = cost, quindi, siccome m è costante, si ha anche p = mv = cost).
57
4.5.2 Momento angolare
Dato un arbitrario centro di riduzione O ′ in un SdR, si dice momento angolare (o momento
della quantità di moto) rispetto ad O ′ il vettore
LO′ = (P − O ′ ) ∧ p . (4.30)
4.5.3 Impulso
Si parla di impulso di una forza, quando questa agisce per un certo lasso di tempo, breve, sul
sistema. L’impulso sarà allora definito come
Z t
I= Fdt . (4.35)
to
Siccome F = ma, si ha
Z t Z t
I= ma dt = m dv = mv(t) − mv(to ) = ∆p , (4.36)
to to
cioè
Teorema 4.5.1 L’impulso di F agente sul punto P nel lasso di tempo ∆t = t − t0 è pari alla
variazione della quantità di moto fra l’istante in cui comincia e l’istante in cui cessa di agire.
Gli impulsi entrano in gioco, soprattutto, negli urti, dove si hanno forze intense, applicate
però per intervalli di tempo molto brevi.
58
4.6 Dinamica relativa
In un SdR inerziale sappiamo che è valida la seconda legge di Newton F = ma, ma se vogliamo
applicarla in un SdR non inerziale, S, dobbiamo tener conto del così detto teorema di Coriolis,
secondo il quale
a = aR + at + aC . (4.37)
Quindi, in S
maR = ma − mat − maC = F + Ft + FC , (4.38)
dove F = ma è la forza realmente e direttamente agente sul punto P, che si registrerebbe nel
SdR inerziale; Ft = −mat e FC = −maC sono la forza di trascinamento e la forza di Coriolis,
rispettivamente. Queste vengono dette forze apparenti, in quanto non si può trovare il corpo
che le applica sul punto P. In pratica sono dei termini correttivi, che servono all’osservatore del
SdR S, non inerziale, per poter spiegare il moto di P e poter ancora scrivere la seconda legge
nella forma
Φ = maR , (4.39)
dove adesso però Φ = F + Ft + FC .
In S si potrà ancora parlare di principio d’inerzia, basterà considerare l’annullarsi di Φ e
non di F. In poche parole, se io applico a P una forza F tale che F = maC + mat , nel SdR
S, non inerziale, vedrò P fermo, o muoversi di moto rettilineo ed uniforme (cosa non vera per
l’osservatore nel SdR inerziale).
4.7 Esempi
4.7.1 Piano inclinato liscio
Consideriamo un punto materiale di massa m che si
y muova su un piano inclinato di un angolo θ rispetto al-
Rn l’orizzontale sotto l’azione della forza di gravità. Sup-
poniamo che il vincolo (il piano) sia liscio. Ciò vuol di-
O re che l’unica reazione vincolare cui è soggetto il punto
materiale è normale alla superficie del vincolo.
θ Con riferimento alla figura, scegliamo un sistema di
x riferimento Cartesiano il cui asse x sia lungo la penden-
P
za del piano e l’asse y ovviamente perpendicolare alle
x. Il verso positivo dell’asse x sia verso valle e quello
dell’asse y verso l’alto. Inoltre supponiamo di far partire il punto da fermo (ẋ(0) = 0) in
x(0) = x0 .
Studiamo il moto del punto materiale. Scriviamo l’equazione del moto, F = ma, nel nostro
sistema di riferimento. Proiettando l’equazione vettoriale dui due assi cartesiani, otterremo le
due equazioni differenziali
0 = Rn − P cos θ , (4.40)
mẍ = R sin θ , (4.41)
dove abbiamo considerato il fatto che il moto si svolge solo lungo l’asse delle x (quindi l’acce-
lerazione lungo le y è nulla) e dove
P = mg . (4.42)
59
L’equazione (4.40) ci dice come si comporta il vincolo:
Ovvero, il vincolo reagisce alla componente normale del peso, P cos θ, con una reazione vincolare
in modulo pari a P cos θ, con la stessa direzione ma con verso opposto, necessaria per far sì che
il punto materiale di muova sul vincolo e non sprofondi in esso.
L’equazione (4.41), invece, ci dà la vera equazione del moto
ẍ = g sin θ . (4.44)
Il moto quindi risulta essere un moto uniformemente accelerato con accelerazione costante pari
a g sin theta.
Integrando si ottiene la velocità
x(0) = x0 , (4.47)
ẋ(0) = 0 , (4.48)
0 = Rn − P cos θ , (4.51)
mẍ = R sin θ − Rt , (4.52)
60
L’equazione (4.51) ci dà la reazione normale del vincolo
che è ancora un moto uniformemente accelerato con accelerazione costante pari a g(sin θ −
µd cos θ).
Rn = P = mg , (4.59)
ẍ = −µd g . (4.61)
x(0) = 0 , (4.63)
ẋ = v0 , (4.64)
si ottiene
1
x(t) = v0 t − µd gt2 . (4.65)
2
61
Se vogliamo trovare lo spazio di frenata sf , ovvero lo spazio che percorre il punto materiale
finché non si ferma, si può trovare l’istante t0 in cui ciò avviene imponendo l’annullarsi della
velocità
ẋ(t0 ) = v0 − µd gt0 = 0 , (4.66)
che dà
v0
t0 = (4.67)
µd g
e calcolare lo spazio raggiunto a t = t0
1 v2 v2 v2
sf = x(t0 ) = v0 t0 − µd gt20 = 0 − 0 = 0 . (4.68)
2 µd g 2µd g 2µd g
Un altro modo per trovare direttamente sf consiste nell’utilizzare (sostanzialmente) il teore-
ma delle forze vive (che vedremo nel prossimo capitolo). Sappiamo che il punto ha accelerazione
costante a = −µd g. Si può scrivere
dv dx dv
a dx = dt = v dt = v dv . (4.69)
dt dt dt
Integrando la (5.16) fra l’istante iniziale e quello finale si ha
Z sf Z 0
a dx = v dv , (4.70)
0 v0
m2 ÿ = m2 g − T . (4.72)
62
Consideriamo adesso il moto di m1 . Il punto sarà soggetto alla forza peso, P1 e alla reazione
vincolare normale che dovrà annullare l’azione del peso per far sì che m1 si muova sul piano
senza sprofondare. Inoltre, m1 sarà soggetto alla tensione T della fune (che non fa altro che
comunicare la tensione da un capo all’altro della sua lunghezza). Prendiamo l’asse delle x
diretta positivamente verso destra (nel senso del moto). Le equazioni del moto saranno date da
0 = Rn − m1 g , (4.73)
m1 ẍ = T . (4.74)
Adesso bisogna notare che la variazione di x è esattamente data dalla variazione di y. La fune è
insestensibile, quindi di quando si abbassa la quota di m2 di tanto aumenta la x di m1 . Potremo
allora scrivere che
ẍ = ÿ . (4.75)
La soluzione del sistema (4.73,4.74), allora, è possibile con la sostituzione del modulo della
tensione T da (4.73) a (4.74). Si ottiene
(m1 + m2 )ẍ = m2 g (4.76)
e quindi
m2
ẍ = g. (4.77)
(m1 + m2 )
Il moto è uniformemente accelerato, con accelerazione pari alla (4.77):
1 m2
x(t) = x0 + ẋ0 t + gt2 , (4.78)
2 (m1 + m2 )
dove le costanti x0 e ẋ0 andranno trovate imponendo le condizioni iniziali.
La tensione della fune è data da
m1 m2
T = m1 ẍ = g, (4.79)
(m1 + m2 )
ed è quindi costante lungo tutto il moto.
63
Il secondo punto, m2 , è soggetto anche lui al peso, P2 e alla tensione della fune T. Questa
volta la coordinata y sarà positiva verso l’alto. Siccome la fune è inestensibile, se la massa m1
si abbassa, la massa m2 si alzerà, e vice versa. Allora
m2 ÿ = T − m2 g . (4.81)
T = m2 ÿ + m2 g , (4.82)
(m1 + m2 )ÿ = (m1 − m2 )g . (4.83)
O ı̂ l0 x x
64
Quindi k sarà data in N m−1 o in kg s−2 . l0 è la lunghezza di riposo della molla.
Scriviamo l’equazione del moto. Il punto materiale sarà soggetto alla forza peso, che sarà
compensata dalla reazione normale del vincolo e lungo le x alla forza di richiamo elastico:
0 = Rn − mg , (4.89)
mẍ = −k(x − l0 ) . (4.90)
Quindi l’equazione del moto è
k
(x − l0 ) = 0 .
ẍ + (4.91)
m
Poniamo r
k
ω= , (4.92)
m
e faccaimo un cambio di variabile per riportare l’Eq. (4.91) in forma “canonica”:
ξ = (x − l0 ) , (4.93)
tale che ξ˙ = ẋ e ξ¨ = ẍ. Allora la (4.91) diventa
ξ¨ + ω 2ξ = 0 , (4.94)
equazione del moto armonico.
Possiamo cercare la soluzione generale nella forma
ξ(t) = A cos (ωt + φ) , (4.95)
con l’ampiezza A e lo sfasamento φ che dovranno essere fissate imponendo le condizioni iniziali
del moto.
Tornando alla variabile x si ha dunque
x(t) = l0 + A cos (ωt + φ) . (4.96)
Supponiamo di lasciare il punto materiale da fermo nell’origine degli assi. Allora
x(0) = 0 , (4.97)
ẋ(0) = 0 . (4.98)
Le condizioni iniziali si traducono in un sistema per la determinazione di A e φ:
x(0) = lo + A cos φ = 0 , (4.99)
ẋ(0) = −Aω sin φ = 0 , (4.100)
ovvero
A = −l0 , (4.101)
φ = 0. (4.102)
Il moto allora sarà dato dalla
x(t) = l0 (1 − cos (ωt)) , (4.103)
moto armonico di periodo r
2π m
T = = 2π . (4.104)
ω k
65
2π
ω
2x0 ωx0 ω 2 x0
3π π 3π
2ω ω 2ω
x0 π π π 2π
O 2ω ω t O 2ω ω t
π π 3π 2π
−ωx0 −ω 2 x0
O 2ω ω 2ω ω t
−
→
g
y0
O
Prendiamo un sistema di riferimento in cui l’asse delle y è diretta verso l’alto. Se y0 è il
punto di riposo della molla, si avrà per il Secondo Principio
66
Quindi: il moto è sempre un moto armonico, con la stessa pulsazione e periodo. La forza
di gravità è una forza costante, quindi non cambia queste caratteristiche del moto. Interviene
però nello spostamento del punto di equilibrio. Infatti, la soluzione della (4.108) sarà
x
O x0 x
Il punto è soggetto, quindi, oltre alla forza di richiamo elastico, anche alla forza d’attrito.
Siccome il punto è massivo, subisce la forza attrattiva di gravità, in conseguenza della quale il
vincolo reagisce con una reazione normale
Rn = mg . (4.112)
Ft = µd Rn = µd mg . (4.113)
Il verso della forza d’attrito si oppone al verso della velocità del punto. Allora, potremo scrivere
l’equazione di moto nel modo seguente:
ẋ
mẍ = −k(x − x0 ) − µd mg . (4.114)
|ẋ|
Il termine ẋ/|ẋ| è un segno, ±1, che quindi sta ad indicare il fatto che quando ẋ > 0, cioè
quando il punto si muove nerso le x positive, la forza d’attrito si oppone al moto ed è quindi
diretta verso l’origine. Vice versa, quando ẋ < 0, la forza d’attrito è diretta verso le x positive.
67
La soluzione della (4.114) si ottiene spezzettando il moto in tratti in cui il termine d’attrito
ha un segno definito.
Supponiamo di partire con le condizioni iniziali
x(0) = x1 = 0 , (4.115)
ẋ(0) = ẋ1 = 0 . (4.116)
ẍ + ω 2 (x − x0 ) + µd g = 0 , (4.117)
dove abbiamo diviso la (4.114) per m e abbiamo considerato il fatto che ẋ/|ẋ| = 1.
Ridefinendo la variabile come
µd g
ξ = x − x0 + 2 , (4.118)
ω
si ha un’equazione del moto armonico in forma “canonica” che ha come soluzione
ẍ + ω 2(x − x0 ) − µd g = 0 , (4.125)
69
x(t)
2x0
2x0 − 2λ
x0
4λ
π 2π
O ω ω t
4.9 Pendolo
Consideriamo un punto materiale di massa m appeso ad una corda ideale di lunghezza r,
soggetto alla forza di gravità (vedi la figura). Il problema ha un grado di libertà. Infatti, data
la lunghezza della corda, che è costante, il punto sarà costretto a muoversi sulla circonfnza di
raggio r. Siccome il problema ha simmetria circolare, è conveniente passare alla coordinate
polari.
−
→
g
θ+
T
ûθ
h ûr V =0
P
v = r θ̇ûθ , (4.141)
a = −r θ̇2 ûr + r θ̈ûθ . (4.142)
70
L’equazione lungo ûr dà la tensione della fune in funzione di θ
La seconda equazione invece dà l’equazione del moto vera e propria, dalla quale possiamo
trovare θ = θ(t):
g
θ̈ + sin θ = 0 . (4.146)
r
Piccole Oscillazioni
La soluzione esatta della (4.146) esula dallo scopo di questo corso. Però, se consideriamo delle
oscillazioni per le quali
θ ≪ 1, (4.147)
dette piccole oscillazioni, si può ottenere una soluzione per la (4.146) in approssimazione. Se
vale la (4.147), si ha anche (arrestando l’espansione di Taylor al primo ordine)
sin θ ≃ θ . (4.148)
Pendolo di Newton
Notiamo che siamo giunti alla (4.149) utilizzando l’uguaglianza fra massa gravitazionale e massa
inerziale. Se avessimo tenuto separate le due, avremmo trovato piuttosto
g mg
θ̈ + θ = 0, (4.151)
r mI
e di conseguenza
r mI
r
T = 2π . (4.152)
g mg
La misura del periodo in Eq. (4.152) può quindi essere utilizzata per mettere in evidenza
eventuali discrepanze dal risultato sperimentale
mI
= 1. (4.153)
mg
71
del corpo in movimento. È esperienza comune che se il corpo ha una superficie che si estende
perpendicolarmente al moto, più è estesa tale superficie più è intensa la resistenza che il corpo
incontra nel suo moto. Inoltre, tale resistenza dipende dalla natura del fluido che il corpo at-
traversa. Per esempio se consideriamo il moto in aria, la resistenza incontrata sarà inferiore a
quella di un moto in acqua.
Queste considerazioni si riassumono nella seguente forma fenomenologica per la forza che
agisce sul corpo che si muove con velocità v nel fluido:
F = −βv , (4.154)
dove β è detto coefficiente d’attrito viscoso, ha le dimensioni di una forza su una velocità
[F ] [m][L][t]−2 [m]
[β] = = = , (4.155)
[v] [L][t]−1 [t]
e quindi si può dare per esempio in kg s−1 .
Il coefficiente d’attrito viscoso, dipendendo sia dalla geometria del corpo che si muove nel
fluido sia dalla natura del fluido, si può esprimere attraverso la relazione seguente:
β = γη. (4.156)
Il coefficiente γ tiene in conto la geometria del corpo, mentre il coefficiente η è detto viscosità
del fluido e tiene in conto la natura di quest’ultimo. Per esempio, la viscosità dell’aria in con-
dizioni standard è η = 1.7 · 10−6 kg m−1 s−1 . In acqua invece η ∼ 10−3 kg m−1 s−1 .
y F
−
→
g
P
O x
72
NB. Il segno della forza peso è chiaro (lo abbiamo già discusso precedentemente); il segno
della forza d’attrito deve essere anch’esso negativo. Infatti, la forza si deve opporre al moto;
se ẏ > 0, ovvero se nel nostro SdR il punto si sta muovendo verso l’alto, la forza d’attrito deve
puntare verso il basso: da qui il segno meno. Allo stesso tempo, se ẏ < 0, ovvero se nel nostro
SdR il punto si sta muovendo verso il basso, la forza deve puntare verso l’alto. Infatti, −ẏ sarà
positivo e la forza, scritta così, punterà verso l’alto (come deve). La stessa cosa avverrebbe se
per esempio prendessimo un SdR con le y positive verso il basso (cambierebbe di segno soltanto
la forza peso, che è costante ...).
Dividendo per la massa, si ha
β
ÿ + ẏ + g = 0 . (4.160)
m
Per risolvere l’equazione (4.160), notiamo che il termine y(t) è assente. La (4.160) è un’e-
quazione differenziale del secondo ordine, ma ha una forma semplice e in realtà è un’equazione
del primo ordine nella velocità. Allora si ha:
ẏ = v , (4.161)
v̇ + γv + g = 0 , (4.162)
β
dove abbiamo posto γ = m . La (4.162) è un’equazione del primo ordine non omogenea.
Facciamo il seguente cambio di variabile
g
ξ=v+ , (4.163)
γ
e quindi
g
v(t) = ξ0 e−γt −. (4.166)
γ
La (4.166) dipende da una costante arbitraria, la ξ0 , che bisognerà trovare imponendo le
condizioni iniziali del moto.
Supponiamo di lasciare il grave da una certa altezza h da fermo. Si ha quindi
y(0) = h , (4.167)
ẏ(0) = 0 . (4.168)
73
La (4.171) ci descrive il cambiamento della velocità col tempo. A t = 0 si ha v = 0. come
deve essere (lo abbiamo imposto noi). Via via che il tempo scorre, la velocità diventa negativa
(il grave si muove verso il basso nel nostro SdR), il termine esponenziale diventa sempre più
piccolo e quindi per t → ∞ la velocità tende ad un valore limite
g
v∞ = − . (4.172)
γ
v(t)
x0
O t
Fisicamente, tale limite si può spiegare così: quando il punto è fermo subisce solo la forza
di gravità; si mette in moto e, siccome la velocità diventa 6= 0, comincia a farsi sentire anche
la forza d’attrito. Via via che il corpo cade la velocità aumenta (in modulo) e quindi aumenta
anche la forza d’attrito, che va con v. La differenza fra la forza peso e la forza d’attrito, che
divisa per la massa dà l’accelerazione del corpo, diventa sempre più piccola e, per t → ∞, tende
a zero. Quindi ad un certo punto (a t = ∞) il moto non è più soggetto a forze e quindi diventa
un moto rettilineo uniforme, che ha quindi velocità costante.
Possiamo integrare la (4.171) per ottenere y(t):
Z t
g 1 −γt 1
y(t) = y(0) + v(t)dt = y(0) + − e + −t . (4.173)
0 γ γ γ
mẍ = −β ẋ − kx . (4.175)
Ponendo β/m = 2λ e k/m = ω 2 , ci riduciamo alla ricerca di una soluzione della seguente
equazione differenziale ordinaria del secondo ordine omogenea:
74
Sostituendo la (4.177) nella (4.176), e imponendo che la (4.177) sia soluzione, si ottiene un’e-
quazione del secondo ordine per α (detta equazione indiciale)
La soluzione generale (4.180) descrive moti che sono molto diversi a seconda del valore
dell’espressione sotto radice. Si hanno i seguenti tre casi: λ2 −ω 2 > 0, λ2 −ω 2 = 0 e λ2 −ω 2 < 0.
• Oscillatore sovra-smorzato
Supponiamo quindi che λ2 − ω 2 > 0. Da un punto di vista fisico, ciò vuol dire che la parte
d’attrito viscono “vince” sul richiamo elastico della molla.
Se λ2 − ω 2 > 0, le soluzioni (4.180) sono reali.
√ √
λ2 −ω 2 t λ2 −ω 2 t
x(t) = e−λt Ae−
+ Be , (4.181)
x(0) = x0 , (4.182)
ẋ(0) = ẋ0 . (4.183)
x0 = A + B , (4.184)
√
ẋ0 = −λ(A + B) + λ2 − ω 2 (B − A) , (4.185)
75
e
λe−λt x0 h √ √
2 2
√ √
2 2
i
ẋ(t) = − √ −(λ − λ2 − ω 2 )e− λ −ω t + (λ + λ2 − ω 2 )e λ −ω t
2 λ2 − ω 2
e−λt x0 h √ √
− λ2 −ω 2 t
√ √
λ2 −ω 2 t
i
+ 2
(λ − λ − ω )e 2 2 2
+ (λ + λ − ω )e , (4.190)
2 √ √
e−λt x0 (λ + λ2 − ω 2 )(λ − λ2 − ω 2) −√λ2 −ω2 t
= √ e
2 λ2 − ω 2
√ √
e−λt x0 (λ − λ2 − ω 2 )(λ + λ2 − ω 2 ) √λ2 −ω2 t
− √ e , (4.191)
2 λ2 − ω 2
e−λt x0 ω 2 −√λ2 −ω2 t √
2 2
= √ e − e λ −ω t , (4.192)
2 λ2 − ω 2
che è sempre negativa e non si annulla mai.
Se invece anche ẋ(0) 6= 0, si può avere il caso in cui la velocità si annulla una volta, ovvero
il moto fa una oscillazione e poi converge a zero, come nel caso precedente.
Il moto si dice sovra-smorzato.
• Smorzamento critico
Supponiamo adesso che λ2 − ω 2 = 0. In questo caso l’equazione indiciale è degenere, cioè
ha due soluzioni coincidenti α± = −λ. In questo caso la soluzione generale è
x(0) = x0 , (4.194)
ẋ(0) = ẋ0 , (4.195)
si ha
A = x0 , (4.196)
ẋ0
B = −λ (4.197)
x0
e quindi
−λt ẋ0
x(t) = x0 e 1+ −λ t . (4.198)
x0
Se supponiamo per esempio di lasciare da fermo il punto in x0 con velocità iniziale nulla,
ẋ0 = 0, si ha
x(t) = x0 e−λt [1 − λ t] , (4.199)
che quindi fa un’oscillazione prima di convergere a zero per t → ∞. Infatti
che si annulla
0 = ẋ(t) = −λx0 e−λt (2 − λt) , (4.201)
76
se
2
t= . (4.202)
λ
Se anche ẋ0 6= 0, si ha
−λt ẋ0 ẋ0
ẋ(t) = x0 e −2λ + −λ −λ t , (4.203)
x0 x0
• Moto smorzato
Supponiamo che λ2 − ω 2 < 0. Da un punto di vista fisico, ciò vuol dire che la partedi
richiamo elastico (la molla) vince sull’attrito viscoso.
Ponendo β/m = 2λ, k/m = ω 2 e Ã/m = A, ci riduciamo alla ricerca di una soluzione della
seguente equazione differenziale ordinaria del secondo ordine non omogenea:
La soluzione dell’Eq. (4.206) sarà data dalla somma della soluzione generale dell’omogenea
associata, ẍ + 2λẋ + ω 2 x = 0, e di una soluzione particolare. La soluzione dell’omogenea è:
√ √
2 2 2 2
x(t) = e−λt A1 e− λ −ω t + A2 e λ −ω t , (4.207)
77
dove adesso b ∈ C. Si ha
ẋ = ibeiΩt ẍ = −bΩ2 eiΩt . (4.210)
Sostituendo le (4.210) nella (4.208) e imponendo che (4.209) sia soluzione, otteniamo la seguente
equazione per b:
−bΩ2 + i2λΩb + ω 2 b = Aeip , (4.211)
dove abbiamo semplificato a destra e a sinistra dell’uguale il termine eiΩt . Si trova quindi:
Aeip
b= . (4.212)
ω 2 − Ω2 + i2λΩ
Riscriviamo adesso b come
b = BeiC (4.213)
con B, C ∈ R. Avremo
Akeip k A
B = kbk = 2 2
=p , (4.214)
kω − Ω + i2λΩk (ω 2 − Ω2 )2 + 4λ2 Ω2
2λΩ
C = p + arctan , (4.215)
ω 2 − Ω2
dove abbiamo usato il fatto che se a, b ∈ C allora kabk = kak kbk e ka/bk = kak/kbk.
Quindi, la soluzione dell’Eq. (4.206) è:
√ √ A
2 2 2 2
x(t) = e−λt A1 e− λ −ω t + A2 e λ −ω t + p cos (Ωt + C) (4.216)
(ω − Ω2 )2 + 4λ2 Ω2
2
dove √
ω′ = ω 2 − λ2 . (4.218)
Se aspettiamo un tempo sufficientemente lungo, t ≫ 1, il termine esponenziale diventa piccolo
e abbiamo "accesso" alla parte risonante:
t≫1 A
x(t) = p cos (Ωt + C) . (4.219)
(ω 2 − Ω2 )2 + 4λ2 Ω2
A
f (Ω) = p , (4.220)
(ω 2 − Ω2 )2 + 4λ2 Ω2
78
Figura 4.1: Curve di risonanza.
df (ω 2 − Ω2 − 2λ2 )Ω
= 2A p = 0, (4.222)
dΩ (ω 2 − Ω2 )2 + 4λ2 Ω2
e infine la (4.221).
Studiamo la concavità della f (Ω) al picco. Si ha:
d2 f 4λ2 Ω2R
=− 3 (2λ2 + 3Ω2R ) < 0 . (4.224)
dΩ2 Ω=ΩR (4λ2 (λ2 + Ω2R )) 2
Un’altro importante parametro della curva (4.220) è la larghezza a mezza altezza. Per
trovarla cerchiamo le soluzioni del sistema
A
f (Ω) = √ 2 ,
2λ ω − λ2 (4.226)
A
f (Ω) = p .
(ω 2 − Ω2 )2 + 4λ2 Ω2
Si ha:
(ω 2 − Ω2 )2 + 4λ2 Ω2 − 16λ2 (ω 2 − λ2 ) = 0 , (4.227)
che puo’ essere scritta nella forma
e quindi √ √ √
Ω2± = ω 2 − 2λ2 ± 2 3λ ω 2 − λ2 = Ω2R ± 2 3λω ′ . (4.230)
Per trovare Ω dobbiamo estrarre la radice, trovando quindi 4 soluzioni possibili. Imponendo
però che Ω sia positivo, le due soluzioni negative vanno scartate e infine si ottiene:
√ √
q
Ω± = ω 2 − 2λ2 ± 2 3λ ω 2 − λ2 . (4.231)
79
Per cui:
(∆Ω)2 = (Ω+ − Ω− )2
p
= 2(ω 2 − 2λ2 ) − 2 (ω 2 − 2λ2 )2 − 12λ2 (ω 2 − λ2 )
p
= 2(ω 2 − 2λ2 ) − 2 (ω 2 − 4λ2 )2 − 8λ2 ω 2
q
2 2 2 2 2 2
= 2 ΩR − (ΩR − 2λ ) − 8λ (ΩR + 2λ ) . 2 (4.232)
e quindi √
∆Ω ≃ 2 3λ . (4.234)
80
Capitolo 5
P Z B
dr L= F · dr . (5.1)
r(t)
B Aγ
F
γ
O x
Il lavoro della forza F dipende quindi dal cammino γ e dai punti iniziale e finale, A e B.
Quindi, in generale, il lavoro infinitesimo F · dr non è un differenziale esatto. Per sottolineare
ciò, si indica con il seguente simbolo:
6 dL = F · dr . (5.4)
Noto il moto del punto, (x(t), y(t), z(t)) e nota la F(r, ṙ, t), posso ricondurre l’integrale
curvilineo ad un comune integrale definito, tramite la parametrizzazione della traiettoria:
Z B Z t1
LAB = F · dr = F · vdt (5.5)
Aγ to
81
Il termine
W =F·v (5.6)
è detto potenza della forza F, si indica con W ed indica, quindi, il lavoro compiuto da F
nell’unità di tempo
6 dL = W dt . (5.7)
Da un punto di vista dimensionale, il lavoro si esprime come segue:
[L] = [F ] [l] = [m] [l] [t]−2 [l] = [m] [l]2 [t]−2 , (5.8)
Definizione 5.1.2 Il joule è il lavoro fatto da una forza che sposta il proprio punto di appli-
cazione (parallelamente a sé stessa) per un metro
1 J = 1 Nm . (5.9)
Per la potenza si ha
[W ] = [L] [t]−1 = [m] [l]2 [t]−3 . (5.10)
L’unità di misura che si utilizza per la potenza nel SI è il watt.
Definizione 5.1.3 Il watt è la potenza di una forza capace di compiere il lavoro di un joule in
un secondo
1 W = 1 J/s . (5.11)
È bene sottolineare un’altra volta quanto già asserito implicitamente nelle formule: il lavoro
di una forza dipende dal cammino su cui è calcolato. Per esempio, pensiamo al trascinamento di
un masso pesante per un certo tratto, in linea retta, sull’asfalto. La forza che noi imprimiamo
al masso deve vincere la forza di attrito che c’è fra le due superfici in contatto. Se riusciamo a
spostare il masso per un certo numero di metri, avremo fatto in certo lavoro (o meglio, la forza da
noi applicata avrà fatto un certo lavoro). Se adesso ripetiamo l’impresa, raddoppiando il numero
di metri del trascinamento, il lavoro, fermo restando le condizioni di attrito, sarà raddoppiato
anch’esso. Non solo. Potremo ripetere l’esperimento tenendo fissi i punti di partenza e di arrivo
del tragitto e variando quest’ultimo, per lunghezza e forma. In tutti questi casi il lavoro sarà
diverso, a seconda del cammino su cui ci siamo mossi.
F = ma . (5.12)
82
e ricordando che 6 dL = F · vdt, si ha
2
v
6 dL = d m . (5.14)
2
Il termine
1
T = mv 2 (5.15)
2
è detto energia cinetica del punto materiale P di massa m e velocità v.
Se integriamo lungo la traiettoria γ del punto, otteniamo:
Z B
m B m
Z
LAB = 6 dL = d(v 2 ) = [v 2 (B) − v 2 (A)] = T (B) − T (A) , (5.16)
Aγ 2 A 2
cioè il lavoro fatto dalla forza F, agente sul punto materiale P, per spostarlo dal punto A al
punto B della traiettoria γ, è pari alla variazione di energia cinetica del punto materiale stesso
fra il punto A e il punto B. L’espressione (5.16) è il così detto teorema delle forze vive (l’energia
cinetica veniva chiamata forza viva).
Dal punto di vista della risoluzione del moto di P, il teorema delle forze vive non ci dice
assolutamente nulla. Infatti, per trovare il lavoro che fa F per andare da A a B, bisogna
conoscere già il moto di P. È però un importante risultato, poiché permette di identificare
l’omogeneità fra energia cinetica e lavoro e in alcuni casi particolari (vedi forze conservative),
dà un integrale primo del moto.
γ1
A γn
B
γ2
Inoltre si ha anche la
LAγB = −LBγA , (5.18)
83
poiché Z B Z A Z A
F · dr = F · (−dr) = − F · dr (5.19)
Aγ Bγ Bγ
• 1) F · dr = 0.
H
Infatti, sia F conservativa. Presi γ1 e γ2 due cammini qualsiasi che vanno dal punto A al
punto B (come nella figura precedente), si ha che
I
LAγ1 B = LAγ2 B ⇒ LAγ1 B + LBγ2 A = F · dr = 0 . (5.20)
Vice versa, se l’integrale ciclico lungo un cammino chiuso è nullo, posso spezzare questo
cammino in due, γ1 e γ2 ed avere quindi
Infatti, sia F conservativa. Fissiamo un punto qualsiasi dello spazio, O (come in figura).
A
B
Ponendo
Z
LOA = F · dr = U(A) , (5.23)
Z OA
LOB = F · dr = U(B) , (5.24)
OB
si ha che
Z Z Z
LAB = F · dr = F · dr + F · dr
AB AO OB
84
Z Z
= − F · dr + F · dr =
OA OB
= U(B) − U(A) . (5.25)
Il vice versa è banale. Infatti se è valida la (5.25), vuol dire che il lavoro fatto da F
dipende soltanto dai punti A e B e non dal cammino fatto per arrivare dall’uno all’altro.
La funzione U(x, y, z) è detta Potenziale della forza F.
• 3) Esiste una funzione del posto U(x, y, z), con U ∈ C 2 (Ω) cioè continua con derivate
continue sul dominio Ω, tale che
∂U ∂U ∂U
F = ∇U = i+ j+ k. (5.26)
∂x ∂y ∂z
Infatti, sia F conservativa. Consideriamo uno spostamento infinitesimo lungo l’asse x,
dx. Allora
6 dL = F · dx i = Fx dx , (5.27)
ma si ha anche, per il punto 2), che
Quindi
U(x + dx, y, z) − U(x, y, z) ∂U
Fx = = (5.29)
dx ∂x
e così per le altre componenti.
Vice versa, se F = ∇U, si ha che
Z B Z B B B
∂U ∂U ∂U
Z Z
F · dr = ∇U · dr = ( dx + dy + dz) = dU ,
Aγ Aγ Aγ ∂x ∂y ∂z Aγ
= U(B) − U(A) , (5.30)
∂ ∂ ∂U ∂2U
Fz = =
∂y ∂y ∂z ∂y∂z
= per il teorema di Schwartz
∂2U ∂
= = Fy . (5.32)
∂z∂y ∂z
85
Allora la componente x del rotore si annulla. Così per tutte le componenti del rotore e
quindi
rotF = 0 . (5.33)
Forza di gravità
Un primo esempio di forza conservativa è la forza peso. Il campo corrispondente è uniforme.
86
z Se consideriamo il SdR in figura, in ogni
punto di R3 il campo è definito come
g = −g k̂ . (5.41)
−
→
g
Se il punto materiale P ha massa m, su di lui
agisce una forza
F = mg = −mg k̂ , (5.42)
y
O
quindi le componenti di F sono
Fx = Fy = 0 , Fz = −mg . (5.43)
Forza elastica
Un altro esempio di forza conservativa è dato dalla forza elastica. Prendiamo il caso generale
tridimensionale (è una forza centrale):
F = −kr , (5.50)
cioè Fx = −kx, Fy = −ky, Fz = −kz. Le derivate miste sono tutte nulle, quindi F è
irrotazionale e quindi conservativa.
L’energia potenziale è data da
1 1
V (x, y, z) = kr 2 = k(x2 + y 2 + z 2 ) , (5.51)
2 2
tale che
F = −∇V . (5.52)
87
Forza di gravitazione
La forza di gravitazione universale è una forza centrale definita come segue:
Mm
F = −γ ûr . (5.53)
r2
Quindi, in coordinate cartesiane si ha
x
Fx = −γMm 3 , (5.54)
(x2 + y 2 + z 2 ) 2
y
Fy = −γMm 3 , (5.55)
(x2 + y 2 + z 2 ) 2
z
Fz = −γMm 3 . (5.56)
(x2 + y 2 + z 2 ) 2
dûr
⊥ ûr . (5.59)
dt
Siccome la forza è lungo ûr , come si vede in Eq. (6.58), si ha che
Mm Mm
F · ṙ dt = −γ 2
ṙ dt = −γ 2 dr . (5.60)
r r
Quindi
rB
Mm 1 1
Z
LAγB = − γ 2 dr = γMm − , (5.61)
rA r rB rA
= U(B) − U(A) = V (A) − V (B) . (5.62)
88
5.5 Conservazione dell’energia meccanica
Noto il risultato del teorema delle forze vive, LAB = T (B) − T (A), dove LAB è il lavoro della
forza F applicata al punto materiale P, di massa m, consideriamo F conservativa.
Per quanto detto due paragrafi fa, si avrà
LAB = U(B) − U(A) (5.64)
e quindi sostituendo nella formula del teorema delle forze vive:
T (B) − U(B) = T (A) − U(A) . (5.65)
Questo importante risultato ci dice che se la forza che agisce su P è conservativa, la quantità
(T − U), che chiameremo energia meccanica totale del punto, si conserva durante il moto.
Conviene introdurre l’energia potenziale di P, che non è altro che U cambiata di segno:
V (x, y, z) = −U(x, y, z) , (5.66)
con la quale, fra l’altro, F = −∇V .
Se conosco l’energia meccanica totale in un certo punto O del moto, E = T (O) + V (O),
l’equazione (5.65) mi dice che in ogni altro punto, la somma dell’energia cinetica e dell’energia
potenziale uguaglierà E:
1 2
mv + V = E . (5.67)
2
L’Equazione (5.67) è un’equazione differenziale del primo ordine (ci sono le velocità), detta
equazione della conservazione dell’energia, che descrive il moto di P. È quindi un’equazione
che può essere utilizzata nella soluzione di questo. È un’equazione scalare, quindi rimuove
un grado di libertà. Per sistemi ad un grado di libertà, quindi, dà la soluzione del moto.
Essendo un’equazione del primo ordine, è, in generale, più semplice da integrare. Di qui la
grande importanza nella ricerca degli integrali primi del moto, ovvero delle quantità che, sotto
particolari condizioni, si conservano (cioè sono costanti del moto).
Integrando la (5.70), con le condizioni iniziali che determinano E, si trova la funzione x(t)
ovvero l’equazione oraria del punto P.
89
−
→
g
θ+
T
ûθ
h ûr V =0
P
5.6 Il pendolo
Consideriamo il pendolo in figura 5.1. Il problema è ad un grado di libertà. Inoltre l’unica
forza agente sul punto che fa lavoro è la forza peso, che è conservativa (la tensione del filo è
sempre perpendicolare allo spostamento e quindi fa lavoro nullo). Quindi, possiamo utilizzare
l’equazioni di conservazione dell’energia
E =T +V , (5.71)
dove
1
T = m|ṙ|2 , (5.72)
2
V = mgh = mgr(1 − cos θ) . (5.73)
In coordinate polari si ha
r = rûr , (5.74)
ṙ = ṙûr + r θ̇ûθ = r θ̇ûθ , (5.75)
visto che ṙ = 0, per cui |ṙ|2 = r 2 θ̇2 . Allora l’equazione della conservazione dell’energia si riscrive
come segue
1 2 2
mr θ̇ + mgr(1 − cos θ) = E . (5.76)
2
L’energia meccanica E è da determinarsi sapendo le condizioni iniziali. Supponiamo di lasciare
il pendolo a t = 0 in θ = θ0 da fermo, θ̇(0) = 0. L’energia meccanica iniziale è semplice da
calcolare. Infatti, siccome il pendolo viene lasciato da fermo, l’energia cinetica è nulla. Si avrà
quindi soltanto energia potenziale
90
Quindi la (5.76) diventa
1 2 2
mr θ̇ + mgr(cos θ0 − cos θ) = 0 , (5.78)
2
ovvero
2g
θ̇2 +
(cos θ0 − cos θ) = 0 . (5.79)
r
Per provare che la (5.79) è un integrale primo del moto, deriviamola rispetto al tempo e
otteniamo l’equazione del moto del pendolo. Si ottiene
2g
0 = 2θ̇θ̈ + θ̇ sin θ ,
rg
= 2θ̇ θ̈ + sin θ , (5.80)
r
e quindi per θ̇ 6= 0,
g
θ̈ +
sin θ = 0 . (5.81)
r
La soluzione generale della (5.79) è la seguente:
r
dθ 2g
=± (cos θ − cos θ0 ) , (5.82)
dt r
ovvero Z θ Z t
dθ
q = dt = t . (5.83)
θ0 ± 2g
r
(cos θ − cos θ0 ) 0
Piccole oscillazioni
Se θ ≪ 1 si può espandere la (5.79) in serie di Taylor. Tenendo il secondo ordine (cos θ ∼
1 − θ2 /2) si ha:
g
θ̇2 + (θ2 − θ02 ) = 0 , (5.84)
r
che integrata dà la Z θ
dθ
q = t. (5.85)
θ0 ± g 2 2)
(θ
r 0
− θ
p
Moltiplicando ambo i membri per ω = g/r, si ottiene
Z θ
dθ
p
2
= ωt . (5.86)
2
θ0 ± (θ0 − θ )
Siccome
θ θ/θ0
dθ dξ θ
Z Z
p =± p = ± arcsin , (5.87)
θ0 ± (θ02 − θ2 ) 1 ± (1 − ξ 2 ) θ0
si ha
θ
arcsin = ±ωt , (5.88)
θ0
e quindi
θ(t) = θ0 sin (ωt) , (5.89)
dove abbiamo preso la soluzione positiva in quanto per t = 0 abbiamo posto θ = θ0 e non
θ = −θ0 .
91
5.7 Oscillatore armonico
Consideriamo adesso il punto materiale soggetto alla forza elastica di una molla (moto unidi-
mensionale). L’equazione di conservazione dell’energia si esprime come segue:
1 1
E = mẋ2 + kx2 , (5.90)
2 2
da cui r
2E k
ẋ = ± − x2 . (5.91)
m m
Quindi la soluzione generale è
x t
dx
Z Z
q = dt = t − t0 . (5.92)
x0 ± 2E
m
− k 2
m
x t0
92
5.8 Teorema delle forze vive forze conservative e forze dis-
sipative
Abbiamo visto che
LAγB = TB − TA . (5.101)
D’altra parte Z B
LAγB = F · dr . (5.102)
Aγ
e Z B
Fcons · dr = VA − VB . (5.104)
Aγ
e quindi Z B
TB − TA = VA − VB + Fd · dr , (5.106)
Aγ
o anche Z B
(TB + VB ) − (TA + VA ) = Fd · dr . (5.107)
Aγ
NB. Scriversi sempre il teorema delle forze vive per esteso, come in (5.106) e poi svolgere
l’integrale stando attenti al segno del prodotto scalare (se la forza è dissipativa si oppone al
moto e quindi verrà fuori un segno negativo). Solo così si evitano errori di segno nell’applicazione
della (5.106).
Supponiamo, adesso, che tutte le forze Fi che appaiono nella (5.108) siano conservative. Allora
sarà possibile trovare un’energia potenziale, eventualmente somma delle varie energie potenziali
delle singole forze, se queste non hanno origine comune, V (x, y, z) tale che
X
Fi = F = −∇V . (5.109)
i
93
Adesso, F = 0 se e solo se ∇V = 0. Quindi, le condizioni di equilibrio vanno ricercate tra
quelle configurazioni, compatibili con i vincoli, che rendono stazionaria l’energia potenziale.
Questo risultato è molto importante, poiché dà uno strumento potente nello studio dell’e-
quilibrio dei sistemi vincolati. Infatti, quando scriviamo l’energia potenziale, automaticamente
teniamo conto dei vincoli introducendo l’opportuno numero di parametri indipendenti, atti ad
individuare la configurazione del sistema. La differenziazione sarà fatta, allora, rispetto a questi
parametri.
Il punto stazionario può essere di massimo, di minimo o di flesso. Nei tre casi cambiano
le caratteristiche dell’equilibrio. Si dice di equilibrio stabile una configurazione per cui una
variazione infinitesima dei gradi di libertà, fa sì che il sistema rimanga in un intorno della confi-
gurazione di equilibrio. Per intenderci meglio, consideriamo un punto materiale in equilibrio nel
punto più basso di una valle, sotto l’azione della forza peso. Questa è una posizione di equilibrio
stabile, perché, se proviamo ad allontanare un po’ il punto da questa, non se ne discosterà mai
indefinitamente. Se invece mettiamo il punto sulla cima di una collina (il punto di massimo
è di equilibrio), questa sarà una configurazione di equilibrio instabile. Infatti allontanatolo un
po’, se ne discosterà indefinitamente. Anche il flesso è instabile, poiché da una parte si ripete
la condizione appena analizzata. La stabilità si ottiene se la derivata seconda dell’energia po-
tenziale è positiva, cioè in corrispondenza dei punti di minimo relativo. Si può avere, infine,
anche una situazione di equilibrio indefinito. Ciò avviene in configurazioni in cui non solo la
derivata prima si annulla, ma anche le derivate di ordine superiore (per esempio un tratto in
cui il potenziale è piatto).
Prendiamo per semplicità un caso unidimensionale V = V (x) e visualizziamo i punti di
equilibrio.
V (x)
x1 x2 x3 x4 x
94
il punto materiale verso x2 . Analogamente, se il dx è negativo (si va verso x1 ) la derivata
dell’energia potenziale è negativa e quindi la forza è positiva e richiama il punto verso x2 .
95
Capitolo 6
La Legge di Grazitazione Universale è uno dei più grandi successi di Newton, che riesce a
sintetizzare in una formula matematica le osservazioni sperimentali fatte dai grandi astronomi
che lo precedono.
In particolare, si deve sottolineare l’importanza di Copernico che nella prima metà del Cin-
quecento cambia la prospettiva nella concezione dell’Universo, e le raffinate osservazioni di
Tycho Brahe, negli ultimi 20 anni del Cinquecento. Keplero, allievo di Brache, sulla base di
queste osservazioni enuncia, nei primi 20 anni del 1600, le sue tre Leggi. Sono Leggi “descritti-
ve”, “fenomenologiche”, non sono formulate come “leggi fisiche”, non costituiscono un modello.
Ciò non di meno, sulla base di queste leggi empiriche e del secondo Principio, Newton dà la
spiegazione dinamica del movimento dei pianeti intorno al Sole, facendo anche di più, ovve-
ro capendo la portata del suo modello che assurge a Legge Universale, cioè valida anche per
fenomeni che non si pensava potessero essere correlati.
Seconda Legge. Il raggio vettore (che unisce il Sole al pianeta durante il suo moto)
spazza aree uguali in tempi uguali (la velocità areolare è costante.
Terza Legge. Il rapporto fra il quadrato del periodo di rivoluzione e il cubo del
semiasse maggiore dell’ellisse è una costante uguale per tutti i pianeti.
Introduciamo adesso alcuni concetti di cui abbiamo bisogno in seguito e andiamo poi a
derivare la Legge di Gravitazione Universale partendo dalle Leggi di Keplero.
96
e quindi
1 1
I I
S= (x dy − y dx) = k r ∧ dr k , (6.2)
2 2
dove r = xß̂ + y æ̂ and dr = dxß̂ + dy æ̂. Essendo dr = dr
dt
dt = vdt, si ha anche
1
I
S= k r ∧ v k dt , (6.3)
2
da cui
1
Ṡ = k r∧v k . (6.4)
2
Esprimendo r in coordinate polari, si ha
Definizione 6.1.1 Si dice Ellisse il luogo dei punti P tali che P F1 + P F2 = 2a con a ∈ ℜ+
tale che a > c.
F1 e F2 sono detti fuochi dell’ellisse e c è la distanza focale.
In coordinate cartesiane quindi si ha
p p
(x + c)2 + y 2 + (x − c)2 + y 2 = 2a , (6.8)
o anche
(a2 − c2 ) x2 + a2 y 2 = a2 (a2 − c2 ) . (6.11)
Poniamo
a2 − c2 = b2 (6.12)
che è positivo visto che a > c. Allora
b2 x2 + a2 y 2 = a2 b2 , (6.13)
97
Figura 6.1: Ellisse in forma canonica relativa agli assi
1 + cos (2θ)
= cos2 θ = ,
2
Z π
ab
= (1 + cos x) dx = πab . (6.17)
2 −π
Il rapporto
c
e= (6.18)
a
è detto eccentricità dell’ellisse. Siccome a > c si ha 0 ≤ e ≤ 1. Per e = 0 si ha c = 0 e quindi
l’ellisse degenera in una circonferenza. Per e → 1, invece, c → a e quindi l’ellisse è sempre più
schiacciata e al limite degenera nel segmento che unisce i fuochi.
Ponendo c = ae, si ha
b2 = a2 − c2 = a2 (1 − e2 ) (6.19)
e quindi
b √
= 1 − e2 . (6.20)
a
θ
F1 H F2
Abbiamo
2
P F1 = r2 , (6.21)
98
2
P F2 = r 2 + (2c)2 − 2r(2c) cos θ , (6.22)
= r 2 + 4c2 − 4rc cos θ , (6.23)
per cui √
r 2 + 4c2 − 4rc cos θ + r = 2a . (6.24)
Portando r a dx dell’uguale e quadrando, si ha
• Prendiamo un SdR centrato sul Sole e facciamo l’ipotesi che sia “inerziale”. In questo
sistema il Sole è fisso e vale il Secondo Principio nella sua forma più semplice, F = ma.
• Il pianeta ed il Sole sono entrambi approssimabili come punti materiali (ottima approssi-
mazione).
• Il moto di un pianeta, sul quale ci concentriamo, non viene perturbato dal moto degli altri
pianeti che che girano intorno al Sole. Praticamente considero un sistema Sole-pianeta
isolato.
Fatte queste approssimazioni iniziali (che rispecchiano in realtà quello che Brache e Keplero
vedevano, con l’accuratezza dell’epoca, passiamo ad analizzare le tre Leggi e vediamo cosa ne
possiamo ricavare.
Siccome r e v sono sempre sul piano del moto, LO è un vettore con direzione costante.
99
6.2.2 La velocità areolare è costante ...
Abbiamo visto che
1
Ṡ = r 2 θ̇ = cost . (6.28)
2
Allora LO è costante anche in verso e modulo. Quindi è un vettore conservato durante il moto
L̇O = 0 . (6.29)
Siccome L̇O = M0 , dove M0 è il momento della forza che agisce su P rispetto all’origine O, si
ha automaticamente che
L̇O = 0 =⇒ MO = 0 . (6.30)
Essendo
MO = r ∧ F , (6.31)
con r, F 6= 0, se ne deduce che r k F. Cioè il vettore forza gravitazionale è lungo il raggio
vettore che unisce il Sole al pianeta.
Un altro modo per vedere la stessa cosa è il seguente. Scriviamo l’accelerazione generica in
coordinate polari. Avremo
d
a = ṙûr + r θ̇ûθ , (6.32)
dt
= (r̈ − r θ̇2 )ûr + (2ṙθ̇ + r θ̈)ûθ , (6.33)
1d 2
= (r̈ − r θ̇2 )ûr + (r θ̇)ûθ . (6.34)
r dt
Siccome r 2 θ̇ = 2Ṡ = cost, il pezzo lungo ûθ si annulla, e quindi
e
F = ma ∝ ûr . (6.36)
d2 1 d2
e e 1 1
= − cos θ = cos θ = − . (6.44)
dθ2 r dθ2 p p p r
Quindi, in totale
L2 L2 1 L2 1
2 1 1
a = r̈ − r θ̇ = − 2 2 − − = − . (6.45)
mr p r m2 r 2 r m2 p r 2
T2
= k. (6.49)
a3
101
Ma b2 /a = p, per cui
4π 2 m2
p = k, (6.53)
L2
da cui infine
L2 4π 2
= m. (6.54)
mp k
La (6.54) ci dice che
m
F = −k ′ ûr . (6.55)
r2
Per considerazioni di simmetria si trova che
k ′ = γM , (6.56)
dove M è la massa del Sole. Per cui infine si trova la Legge di Gravitazione Universale
Mm
F = −γ ûr , (6.57)
r2
dove γ è la costante di Gravitazione Universale e vale γ = 6.670 · 10−11 N m2 kg−2 .
Mm
V (r) = −γ , (6.60)
r
tale che
F = −∇V . (6.61)
102
Siccome la forza è centrale, il suo momento rispetto all’origine del sistema di riferimento
(cioè rispetto alla posizione del Sole) è nullo
MO = r ∧ F = 0 , (6.62)
r = r(t) , (6.66)
θ = θ(t) . (6.67)
Come equazioni differenziali abbiamo a disposizione due integrali primi del moto, che de-
rivano dalle due leggi di conservazione: la conservazione dell’energia meccanica, visto che F è
conservativa e la conservazione del momento angolare, visto che F è anche centrale.
Troviamo l’energia cinetica del pianeta. Avremo
1 1
T = mkṙk2 = m(ṙ 2 + r 2 θ̇2 ) . (6.68)
2 2
L’energia potenziale dipende solo da r ed è data da
Mm
V = −G . (6.69)
r
Quindi, l’equazione di conservazione dell’energia meccanica è
1 Mm
m(ṙ 2 + r 2 θ̇2 ) − G =E. (6.70)
2 r
Consideriamo adesso il momento angolare del pianeta rispetto al sole, fisso in O:
L = kLO k , (6.72)
1 2 L2 Mm
mṙ + 2
−G =E. (6.75)
2 2mr r
Risolvendo questa equazione, si potrà ricavare r(t) e sostituirlo nella (6.73), per poi integrare
e trovare anche θ(t).
Studiamo l’Eq. (6.75). Si tratta di un’equazione per un moto unidimensionale, parametriz-
zato da r(t), dove l’energia cinetica
1
T = mṙ 2 (6.76)
2
si somma ad un’energia potenziale “efficace”
L2 Mm
Vef f = 2
−G , (6.77)
2mr r
che dipende solo da r. Vef f è la somma di due pezzi: l’energia potenziale della forza gravitazio-
nale e l’energia potenziale della forza centrifuga che agisce sul pianeta nel SdR che gira con θ(t)
intorno al sole insieme al pianeta. In questo SdR non inerziale, noi vediamo il pianeta muoversi
solo lungo il raggio (moto unidimensionale, appunto).
Se rappresentiamo Vef f in funzione di r, otteniamo il grafico in figura
Vef f ≤ E . (6.78)
104
Se, dato sempre L 6= 0, si ha E = Emin , la traiettoria è circolare (moto circolare uniforme)
con r = r0 . Il valore di r0 si trova estremando la funzione Vef f (r):
dVef f L2 Mm
= − 3 + G 2 = 0, (6.79)
dr mr r
da cui si ricava (r 6= 0)
L2
r0 = (6.80)
GMm2
e di conseguenza
M 2 m3
Vef f (r0 ) = ... = −G 2 . (6.81)
2L2
Se E = 0 la traiettoria è parabolica, con v → 0 per r → ∞.
Se E > 0 la traiettoria è iperbolica e v → const per r → ∞.
NB in tutti questi casi (L 6= 0) il pianeta NON cade mai sul Sole.
Allora, se il pianeta parte ad una certa distanza r0 dal Sole ed ha una velocità tale che
E < 0, si allontanerà al massimo fino a rmax e poi tornerà a cadere sul Sole. Se v è tale che
E = 0, riuscirà ad allontanarsi indefinitamente arrivando con v = 0 all’infinito (r → ∞). Tale
velocità è detta “velocità di fuga” vf e si ricava da
1 2 Mm
mvf − G = 0, (6.83)
2 rin
ovvero r
M
vf = 2G . (6.84)
rin
Si vede che la velocità necessaria per sfuggire all’attrazione gravitazionale dipende dall’inverso
della radice quadrata di rin . Quindi, più si parte distanti dal Sole, più piccola è la velocità di
fuga.
105
6.3.1 Traiettoria e terza legge di Keplero
Vediamo adesso come risolvere il sistema di equazioni del moto ottenendo direttamente la
traiettoria r = r(θ). Dall’Eq. (6.75) si ha
r
dr 2E L2 M
=± − 2 2 + 2G . (6.85)
dt m mr r
Questa equazione insieme all’equazione di θ̇ dà la forma parametrica della traiettoria. Noi però
vogliamo la forma esplicita r = r(θ). Abbiamo
r
dr dr 1 mr 2 dr mr 2 2E L2 M
= dθ
= = ± − 2 2
+ 2G , (6.86)
dθ dt dt L dt L m mr r
da cui si trova
L dr
dθ = ± 2
q . (6.87)
mr 2E 2
− mL2 r2 + 2G M
m r
La radice è reale se
2E L2 M
> 2 2 − 2G . (6.88)
m mr r
Integrando termine a termine, si ottiene:
θ(r) r
1 dr
Z Z
dθ = θ(r) − θ(r0 ) = ±L 2
q . (6.89)
θ(r0 ) r0 r 2m E + 2G Mm
− L2
r r2
Poniamo r
L2 2EL2
p= , e = 1 + . (6.90)
GMm2 G 2 M 2 m3
Allora l’Eq. (6.89) si trasforma in
Z r Z r
dr 1 dr p
θ(r) − θ(r0 ) = ± 2
q = ± 2
q . (6.91)
r0 r e2 −1
+ p2 1r − r12 r0 r e 1 − 1e − p
p2 er
Infine:
− pr + 1 − rp0 + 1
θ(r) − θ(r0 ) = ∓ arccos ± arccos . (6.92)
e e
Poniamo r0 = a + c e θ(r0 ) = 0, cioè facciamo una rotazione degli assi per avere la conica in
forma canonica. Allora ...
Si ottiene
p
= (1 − e cos (θ)) , (6.93)
r
che è l’equazione di una conica in coordinate polari.
L’eccentricità è data dall’equazione
r
2EL2
e= 1+ 2 2 3 , (6.94)
G M m
106
quindi a seconda del valore di E possiamo avere i seguenti casi:
E<0 → e<1 → ellisse , (6.95)
E=0 → e=1 → parabola , (6.96)
E>0 → e>1 → iperbole . (6.97)
Se
G 2 M 2 m3
E=− (6.98)
2L2
si ha e = 0 e quindi una circonferenza.
v
n̂
θ
dS
La normale alla superficie ha verso che sarà definito in modo convenzionale. Per esempio si
può decidere di definire la normale “esterna” (ovvero che punti dalla parte della convessità della
107
superficie stessa). Se la superficie ’è chiusa, per esempio, la normale in ciascun punto quindi
rimane sempre esterna alla superficie stessa.
Il campo gravitazionale
M
G = −G ûr , (6.105)
r2
definisce un campo vettoriale in R3 . Troviamo il flusso di G attraverso una superficie chiusa
contenente la massa M che genera il campo:
cos (θ)dS
Z Z Z
φ = G · dS = G cos (θ)dS = −GM
S S S r2
dSn
Z
= −GM 2
, (6.106)
S r
r1
n̂1 n̂2
dSn1 dSn2
M
dΩ
r2
Quindi, abbiamo Z
φ = −GM dΩ = −4π GM , (6.108)
S
indipendente dalla posizione che M occupa dentro S.
Nel caso in cui la superficie chiusa attraverso la quale vogliamo valutare il flusso non contenga
la massa M, il flusso si annulla. Infatti per ogni elemento di superficie con dSn positivo, ce ne
sarà uno con dSn negativo tali che il contributo dei corrispondenti dΩ all’integrale si annullino.
Abbiamo dimostrato il
108
Teorema 6.4.2 (di Gauss). Il flusso del campo gravitazionale attraverso una superficie chiusa
S, che contenga la massa M che genera il campo, è proporzionale ad M e indipendente dalla
posizione che M occupa dentro S. Se M è esterna alla superficie chiusa S, si ha che il flusso
è nullo.
NB1 . Il fatto che il campo vada come 1/r 2 è cruciale per la ricostruzione dell’angolo solido,
che quindi non dipende da r.
NB2 . Il teorema di Gauss vale per qualunque massa interna alla superficie S. Se abbiamo
più masse, M1 , M2 , ... , Mn interne ad S, si avrà che il flusso totale sarà dato da
X X
Φ= φi = −4π G Mi . (6.109)
i i
109
Andamento del campo gravitazionale interno alla sfera
Ci possiamo chiedere quale sia l’andamento del campo gravitazionale all’interno della sfera P ,
ovvero per un R tale che R < r. Per capire questo andamento, consideriamo una sfera S
concentrica a P ma di raggio R < r e consideriamo il flusso del campo gravitazionale generato
da M attraverso questa sfera.
Il teorema di Gauss ci dice che il flusso attraverso tale superficie dipenderà soltanto dalla
massa interna alla superficie stessa e quindi
φ = −G 4π M(R) , (6.117)
dove adesso
4 R3
M(R) = ρ πR3 = M 3 . (6.118)
3 r
Troveremo quindi
M(R) MR
G(R) = −G 2
ûr = −G 3 ûr , (6.119)
R r
con R < r. In totale (
−G Mr3R ûr R<r
G(R) = M
(6.120)
−G R2 ûr R ≥ r.
Quindi, l’andamento in R del campo gravitazionale è lineare all’interno della sfera e all’esterno
ha l’usuale andamento con l’inverso del quadrato di R.
|G(R)|
GM
r2
r R
110
Consideriamo la Terra come una distribuzione di massa a simmetria sferica. Il campo
gravitazionale a distanza h dalla superficie terrestre è quindi determinato dal teorema di Gauss
ed è dato da
MT MT MT 2h
G(RT + h) = −G ûr = −G 2 ûr ≃ −G 2 1 − ûr , (6.121)
(RT + h)2 RT RT
RT2 1 + RhT
dove abbiamo sviluppato al primo ordine in h/RT ∼ 10−5 . Si capisce quindi come, in questa
approssimazione di distribuzione di massa perfettamente sferica, G sia sempre diretto verso il
centro della Terra e il suo modulo cambi al cambiare dell’altezza relativa alla superficie terrestre.
Il cambiamento però è molto piccolo sulle scale tipiche degli esperimenti in laboratorio. Infatti,
il termine 2h/RT è dell’ordine di 10−5 su una variazione di altezza dell’ordine di 100m.
Possiamo riscrivere la (6.121) come segue
2h
G(RT + h) = −G0 1 − ûr , (6.122)
RT
dove
MT 6.66 · 10−11 Nm2 kg−2 · 5.97 · 1024 kg
G0 = G 2 = ≃ 9.8 ms−2 . (6.123)
RT (6.37 · 10 ) m
6 2 2
FP = Fg + Fc , (6.124)
dove
Quindi
111
Si vede dalla (6.128) come, tranne all’Equatore o al Polo Nord, la forza peso NON sia diretta
esattamente verso il centro della Terra (neanche nell’approssimazione molto simmetrice che
abbiamo fatto). Inoltre, in modulo la forza peso è maggiore al Polo Nord (dove si annulla la
forza centrifuga) e minima all’Equatore, dove la forza centrifuga è massima. Questa variazione
è comunque piccola (dell’ordine del per mille).
6.5 Esercizi
6.5.1 Potenziale centrale elastico
6.5.2 Potenziale centrale lineare in r
112
Capitolo 7
Fi = mi ai , (7.1)
dove i sta ad indicare l’i-esimo punto ed Fi è la risultante di tutte le forze che agiscono su
questo. Cerchiamo di ricavare qualche informazione utile, sul moto del sistema, da queste n
equazioni. Sommiamo su i:
X X d X
R= Fi = mi ai = mi vi = Q̇ , (7.2)
i i
dt i
dove Q è la quantità di moto totale del sistema e R la risultante delle forze agenti sul sistema.
. Conviene, a questo punto, fare una distinzione fra le forze che agiscono sull’i-esimo punto
materiale. Ci sono forze che agiscono sul sistema dall’esterno, che chiameremo forze esterne e
sono quelle le cui reazioni, nel senso del terzo principio, non sono applicate a punti del sistema
stesso. Ci sono poi forze, che agiscono sui singoli punti, dovute alle interazioni con gli altri
punti del sistema: queste forze verranno chiamate forze interne ed avranno reazioni applicate
a punti del sistema stesso. Si potrà scrivere, allora:
X X X
R= Fi = Fint
i + Fext
i = Rint + Rext . (7.3)
i i i
in quanto le forze interne sono a due a due uguali in modulo e direzione, hanno la stessa retta
di applicazione e verso opposto, come in figura 7.1.
Si arriva così alla Prima Cardinale della dinamica dei sistemi:
X
Rext = Fext
i = Q̇ . (7.5)
i
113
j
fkj
fij fji = −fij
fki = −fik
fjk = −fkj
...
fji
fjk
i fki fik k
È chiaro che la prima cardinale è un’equazione che dà delle informazioni parziali sulla dinamica
del sistema studiato. Infatti siamo partiti da n equazioni vettoriali e siamo arrivati ad una sola.
In più si sono perse completamente tutte le informazioni relative alle forze interne. L’utilità
di una trattazione sì fatta risiede soprattutto nello studio del corpo rigido, per il quale si
dimostra che le due cardinali (la seconda sarà derivata nel prossimo paragrafo) sono sufficienti
a determinarne univocamente il moto.
Troviamo adesso una forma un po’ più comoda della (7.5).
Introduciamo un punto particolare, il centro di massa del sistema, individuato dall’equazione
P P
im i (Pi − O) mi (Pi − O)
(G − O) = P = i , (7.6)
i mi M
cioè la media, pesata sulle masse, dei raggi vettore dei singoli punti del sistema. Il centro di
massa ha una notevole proprietà: derivando rispetto al tempo la (7.6) si trova:
X
MvG = mi vi = Q (7.7)
i
ovvero, il centro di massa si muove come un punto che abbia la massa totale del sistema e la sua
quantità di moto. Derivando ancora rispetto al tempo, si trova la forma della prima cardinale
che stavamo cercando: X
MaG = Q̇ = Fext
i = Rext . (7.8)
i
Il centro di massa del sistema è un punto geometrico che riassume le caratteristiche del
sistema stesso. Non è necessariamente un punto del sistema.
Se il nostro sistema è un sistema continuo, per trovare il centro di massa, invece di una
somma bisognerà fare un integrale
1 1
Z Z
(G − O) = (P − O) dm = ρ(P − O) dV . (7.9)
M ρV V
Supponendo la densità ρ indipendente dal volume (oggetto omogeneo):
1
Z
(G − O) = (P − O) dV . (7.10)
V V
114
7.2 Conservazione della quantità di moto
Supponiamo adesso che il sistema di punti sia isolato (ovvero che non agiscano forze esterne su
di esso). Allora X
Fext
i =0 (7.11)
i
e quindi
Q̇ = 0 =⇒ Q = cost . (7.12)
La quantità di moto del sistema si conserva, ovvero rimangono costanti nel tempo il suo modulo,
direzione e verso.
Siccome
Q = MvG , (7.13)
ciò vuol dire che il centro di massa del sistema si muove di moto rettilineo uniforme. In
particolare, se a t = 0 si ha vG = 0, allora vG = 0 sempre, e i punti del sistema si muovono a
causa delle interazioni interne in maniera tale da lasciare il loro centro di massa fermo.
Esempio: il rinculo del cannone (o come uscire dal centro di un lago ghiacciato)
Un cannone di massa M spara un proiettile di
massa m in direzione inclinata di θ sull’oriz-
y zontale. Il cannone è inizialmente fermo su un
piano orizzontale liscio. Sapendo che la velo-
cità con la quale viene espulso il proiettile è in
m modulo v0 , determinare il moto del cannone
v0
dopo lo sparo.
θ Supponiamo che lo sparo avvenga nell’i-
M stante t = 0. Prima di questo istante il can-
O x none e il proiettile sono in quiete sulla super-
ficie liscia. Le unichge forze che agiscono sul
sistema sono: la forza di gravità e la reazione vincolare, che equilibra la forza di gravità, poiché
il cannone col proiettile non sprofondano al di sotto di tale superficie.
Siccome la forza peso è
y P = −(M + m)g ĵ , (7.14)
N = (M + m)g . (7.15)
M +m
Lungo la direzione x, invece, le forze che
x
agiscono sul sistema sono nulle, per cui si ha
P
Q̇x = 0 =⇒ Qx = (m + M)vcm = cost . (7.16)
Qx = 0 . (7.17)
115
Dopo lo sparo le forze esterne agenti sul sistema in direzione x rimangono nulle. L’esplosione
che permette di sparare il proiettile genera delle forze interne. Quindi la prima cardinale ci dice
comunque che
Qx = 0 . (7.18)
Ma adesso il sistema è costituito dai due pezzetti (cannone e proiettile), che si muovono con
velocità incognita Vx il primo e con v0 cos (θ) (velocità lungo le x) il secondo. Allora vremo
ovvero
m
Vx = − v0 cos (θ) . (7.20)
M
Quindi se v0 > 0, cioè il proiettile viene sparato verso destra, il cannone rincula verso sinistra
con velocità Vx . Ovviamente più grande è la massa del cannone, rispetto a quella del proiettile,
più la velocità di rinculo del cannone è piccola rispetto alla velovità con cui parte il proiettile.
Nell’istante subito dopo lo sparo, si nota una variazione della quantità di moto lungo le y.
Infatti il proiettile, essendo sparato con un angolo θ rispetto all’orizzontale, ha una componente
di quantità di moto pari a mv0 sin (θ)ĵ lungo le y. Questa variazione della quantità di moto
dovrà essere assorbita dal vincolo.
Fi = mi ai , (7.21)
Con alcuni accorgimenti ne abbiamo ricavato la prima cardinale, muovendo dal concetto di
risultante. Ora partiamo dal concetto di momento risultante MO′ rispetto ad un centro di
riduzione O ′.
z
Fi
Pi
ri
r′i
y
O
rO ′ O′
116
Definiamo il momento risultante come il vettore libero dato dalla somma vettoriale di tutti gli
Mi O′ . Allora si ha X X X
MO′ = Mi O′ = r′i ∧ Fi = r′i ∧ mi ai . (7.23)
i i i
Analogamente, definiamo il momento risultante della quantità di moto rispetto ad O ′ come il
vettore libero somma vettoriale dei singoli momenti angolari:
X
LO′ = r′i ∧ mi vi . (7.24)
i
Fext
i
Pi
r′i
Fint
ij
O′
Fint
ji
r′j Pj
Fext
j
La risultante delle forze agenti su Pi è stata scomposta nella somma della risultante delle forze
esterne agenti su Pi , Fext
i e la forza interna che agisce da Pj su Pi , Fint
ij . Lo stesso è stato fatto
per Pj . Per il terzo principio si ha
Fint int
ij = −Fji . (7.32)
117
Il momento risultante rispetto ad O ′ lo possiamo scomporre nella somma vettoriale di Mext
O′
e MO′ come segue
int
dove
Mext
O′ = r′i ∧ Fext
i + r′j ∧ Fext
j , (7.36)
Mint
O′ = r′i ∧ Fint ′ int
ij + rj ∧ Fji . (7.37)
O ′ . Utilizzando la (7.32) si ha
Prendiamo in considerazione Mint
Mext
O ′ = L̇O ′ + vO ′ ∧ Q . (7.39)
Mext
O ′ = L̇O ′ . (7.41)
Pi = −mi g k̂ . (7.42)
118
La risultante delle forze peso è un vettore che agisce su un punto di massa M = mi pari
P
i
alla massa totale del sistema. Siccome per la prima cardinale dei sistemi si ha
ne deriva che
−Mg k̂ = MaG , =⇒ aG = −g k̂ , (7.45)
cioè il centro di massa del sistema si muove come un punto materiale di massa M soggetto alla
forza peso.
Vediamo adesso cosa possiamo dire della seconda cardinale. Prendiamo i momenti rispetto
al centro del SdR inerziale O, che è quindi fermo. Possiamo scrivere
X X
MO = ri ∧ P i = − ri ∧ mi g k̂ (7.46)
i i
X
= − mi ri ∧ g k̂ = −MrG ∧ g k̂ (7.47)
i
Pi
ri
r′i
y
O
rG G
x
Si ha:
ri = rG + r′i (7.50)
119
e quindi
vi = vG + vi′ (7.51)
Troviamo l’espressione del momento angolare rispetto ad O:
X X
LO = ri ∧ mi vi = (rG + r′i ) ∧ mi (vG + vi′ ) ,
i i
X X X X
= rG ∧ mi vG + rG ∧ mi vi′ + r′i ∧ mi vG + r′i ∧ mi vi′ . (7.52)
i i i i
LO = LG + L′ , (7.56)
dove
LG = rG ∧ MvG (7.57)
e X
L′ = r′i ∧ mi vi′ . (7.58)
i
120
7.6 Equazioni Cardinali e sistemi isolati
Le equazioni cardinali della dinamica sono quindi
Rext = Q̇ , (7.59)
Mext
O = L̇O , (7.60)
dove abbiamo preso come centro di riduzione per i momenti O il centro di massa del sistema
oppure un punto fermo nel SdR inerziale considerato.
Supponiamo di avere un sistema isolato, ovvero un sistema di punti che non sia soggetto a
forze esterne. Allora Rext = 0 e Mext
O = 0. Ne deriva che
( (
Q̇ = 0 Q = cost
=⇒ (7.61)
L̇O = 0 LO = cost .
Se il sistema è isolato, si conservano la sua quantità di moto e il suo momento rispetto al centro
di massa del sistemad O.
Si può leggere questa legge di conservazione anche al contrario e utilizzarla per dare una
formulazione alternativa al terzo principio della dinamica.
Partiamo dall’osservazione sperimentale che in un SdR inerziale per un sistema isolato (cioè
non soggetto a forze esterne, posto suficientemente lontano da qualsiasi altro corpo che possa
interagire con i punti del sistema) si conservano i vettori Q e LO . Allora, possiamo riottiene il
principio di azione e reazione. Infatti, prendiamo per semplicità un sistema costituito da due
punti (ma il ragionamento si può estendere ad un numero qualsivoglia di punti). Si avrà
Q = q1 + q2 = cost , (7.62)
da cui segue che
Q̇ = q̇1 + q̇2 = 0 . (7.63)
Siccome q̇i = Fi , ovvero la derivata della quantità di moto del punto i-esimo è pari alla risultante
delle forze interne agenti sul punto (trattandosi solo di due punti, è pari alla forza con cui Pj
agisce su Pi ), si avrà anche
F1 + F2 = 0 =⇒ F1 = −F2 . (7.64)
Cioè la forza che agisce da Pj su Pi deve avere modulo uguale, direzione uguale e verso opposto
alla forza che agisce da Pi a Pj .
D’altra parte si avrà anche
L̇O = 0 = MO = m1 O + m2 O = r1 ∧ F1 + r2 ∧ F2 . (7.65)
Usando la (7.64), si ottiene
(r1 − r2 ) ∧ F1 = 0 oppure (r2 − r1 ) ∧ F2 = 0 , (7.66)
cioè F1 e F2 sono lungo la congiungente di Pi con Pj (costituiscono una coppia di braccio
nullo!).
Abbiamo ritrovato il principio di azione e reazione. Quindi un modo equivalente di enunciare
il Terzo Principio della Dinamica è
Principio 7.6.1 In un sistema di riferimento inerziale, per un sistema isolato la quantità di
moto e il momento angolare si conservano.
121
7.7 Teorema delle forze vive
Partiamo dalle solite
Fi = mi ai (7.67)
e moltiplichiamo scalarmente a destra per vi e sommiamo su i:
X X
Fi · vi = mi ai · vi , (7.68)
i i
dove per T intendiamo l’energia cinetica totale del sistema. Allora, integrando
dove LAB è il lavoro totale, fatto sugli n punti materiali, dalla configurazione A alla configu-
razione B. Per mettere in evidenza il diverso contributo delle Fi nel lavoro LAB , operiamo la
stessa distinzione, fatta nel paragrafo precedente, fra forze esterne ed interne:
XZ B
T (B) − T (A) = (Fext
i + Fint ext int
i ) · dri = LAB + LAB (7.71)
i A
In generale, pero’, il lavoro fatto dalle forze interne è diverso da zero (le distanze relative fra i
punti cambiano). Quando faremo il corpo rigido, vedremo che in quel caso particolare Lint AB = 0,
in quanto le distanze relative fra i punti del sistema rigido rimangono costanti nel tempo.
dove vi è la velocità dell’i-esimo punto pateriale, ovvero la derivata del suo raggio vettore ri
rispetto al tempo. Consideriamo la figura che abbiamo utilizzato per il teorema di König per
il moemnto angolare e riproponiamo la scomposizione del raggio vettore che individua il punto
i-esimo come somma del raggio vettore del centro di massa e del raggio vettore relativo al centro
di massa
ri = rG + r′i . (7.73)
Per le velocità si ha
vi = vG + vi′ (7.74)
e quindi, sostituendo nella (7.72)
X1
T = mi (vG + vi′ ) · (vG + vi′ ) , (7.75)
i
2
122
1X 2 1X X
= mi vG + mi vi′ 2 + mi vG · vi′ , (7.76)
2 i 2 i i
1 2 1 X X
= MvG + mi vi′ 2 + vG · mi vi′ (7.77)
2 2 i i
L’ultimo termine è nullo per le considerazioni già fatte in precedenza: sarebbe la derivata
temporale del raggio vettore del centro di massa nel sistema di riferimento del centro di massa
stesso e quindi il vettore nullo.
Quindi, come risultato, si ottiene
T = TCM + T ′ . (7.78)
Per quanto riguarda il lavoro delle forze interne, prendiamo in considerazione due punti ma-
teriali, Pi e Pj . Siamo nell’ipotesi che Fij e Fji siano conservative. Per il terzo principio si
avrà:
Fij = −∇i Vij = ∇j Vij = −Fji , (7.80)
con
Vij = V (k ri − rj k) (7.81)
e ∇i gradiente rispetto alla variabile ri .
Il lavoro totale delle forze interne sarà dato dalla somma
XZ B XZ B XZ B XZ B
int
2 LAB = Fij · dri + Fji · drj = − ∇i Vij · dri − ∇j Vji · drj (7.82)
i A j A i A j A
e ponendo ∇ij Vij = ∇i Vij = −∇j Vij e drij = dri − drj , si ottiene:
B
1X 1X
Z
Lint
AB =− ∇ij Vij · drij = − (Vij (B) − Vij (A)) (7.83)
2 ij A 2 ij
123
mediante la quale, una volta definita l’energia cinetica totale come
1X
T = mi vi2 (7.85)
2 i
7.10 Il problema dei due corpi nel SdR centrato sul Sole
Consideriamo il problema Sole-Pianeta (S-P), studiando il moto relativo del pianeta rispetto al
Sole.
z
z′
S
M y
F1 F2
r
P
m
O′ x y′
x′
Supponiamo che il sistema di riferimento (O ′, x′ , y ′, z ′ ) sia inerziale e consideriamo un sistema
di riferimento non inerziale centrato sul Sole, (O ≡ S, x, y, z), con assi paralleli agli assi di
(O ′ , x′ , y ′, z ′ ). Il SdR centrato sul Sole non è inerziale, ma ha velocità angolare nulla. Ciò vuol
dire che l’accelerazione del pianeta nel SdR inerziale è data dalla somma della sua accelerazione
in quel SdR, aR , e dell’accelerazione di trascinamento, che sarà soltanto l’accelerazione del Sole
nel SdR inerziale.
Per il terzo principio si ha, che, se F1 è la forza con cui il pianeta attrae il Sole ed F2 quella
con cui il Sole attrae il pianeta
F1 + F2 = 0 . (7.87)
Scriviamo l’equazione del moto per il pianeta P nel SdR centrato sul Sole. Avremo
maR = F2 + Ft = F2 + maS . (7.88)
Nel SdR inerziale, si ha
MaS = F1 = −F2 , (7.89)
da cui
F2
aS = − . (7.90)
M
Sostituendo la (7.90) nella (7.88) si ottiene
(m + M)
maR = F2 , (7.91)
M
124
ovvero
mM
F2 = aR = µ aR , (7.92)
(m + M)
dove
mM
µ= (7.93)
(m + M)
è detta massa ridotta del sistema S-P. Nel sistema di riferimento centrato sul Sole quindi il
moto si riconduce a quello di una massa µ sottoposta alla forza gravitazionale
Mm
F2 = −G ûr . (7.94)
r2
Se M ≫ m, allora µ ≃ m. Nel caso Sole-Terra, abbiamo
e quindi
m
= 3 · 10−6 . (7.96)
M
Allora
mM m
µ= ≃m 1− ≃m (7.97)
(m + M) M
a meno di qualche parte su un milione.
7.11 Il problema dei due corpi nel SdR del centro di massa
Consideriamo il sistema Sole-Pianeta nel SdR
z centrato nel centro di massa. Abbiamo che,
se il sistema è isolato
S Q̇ = 0 (7.98)
M
F1 e quindi la quantità di moto totale
cm F2 y
P Q = (M + n)vcm (7.99)
m
è costante. Ne deriva il fatto che vcm è costan-
x te e quindi il SdR centrato sul centro di massa
è anch’esso un SdR inerziale. Individuiamo il
Sole col raggio vettore r1 e il pianeta col raggio
vettore r2 . Allora, in questo sistema di riferimento, le equazioni del moto sono
(
Mr̈1 = −G Mr2m ûr1 ,
(7.100)
mr̈2 = −G Mr2m ûr2 ,
125
da cui derivano le due equazioni per i moduli:
m
r1 = r2 , (7.102)
M
M
r2 = r1 , (7.103)
m
ovvero le due possibilità
m+M
r = r2 , (7.104)
M
m+M
r = r1 . (7.105)
m
Sostituendo le (7.105) nelle equazioni del moto (7.100), si ottengono due equazioni disaccoppiate
nei due raggi vettore:
m2
( 1
Mr̈1 = −G Mm (m+M û ,
)2 r12 r1
M2 1
(7.106)
mr̈2 = −G Mm (m+M )2 r2 ûr2 .
2
Quindi entrambi i corpi percorrono orbite ellittiche, che sono ovviamente diverse in ragione
della massa.
Notiamo che, siccome la quantità di moto si conserva,
ovvero il Sole è praticamente fermo nel centro del nostro SdR (il centro di massa coincide col
Sole, infinitamente massivo) e il pianeta si muove in un campo centrale determinato dalla massa
del Sole.
126
Capitolo 8
127
Capitolo 9
z S
′
z
P y
r(t) k̂
r′ (t)
O ĵ
k̂′ î
rO (t)
O ′ ĵ′ y′
′
î x
x′
La configurazione del corpo rigido sarà allora determinata dalla configurazione di S, fissato sul
128
corpo, rispetto al SdR “fisso” in cui si vuole studiare il corpo rigido. Tale configurazione è
fissata da 6 parametri: tre parametri per individuare l’origine del sistema e tre angoli di Eulero
per la rotazione rispetto al sistema fisso.
Per studiare la cinematica di un sistema rigido, basta, quindi, rifarsi a quanto detto nel
paragrafo della cinematica relativa, considerando il vincolo di rigidità. Infatti se fissiamo il
sistema di riferimento S sul corpo, questo parteciperà del moto del corpo stesso.
Il punto O, origine del sistema di riferimento S, sarà un punto generico del corpo rigido. La
velocità di un altro punto P appartenente al corpo rigido sarà data da
vP = vO + −
→
ω ∧ (P − O) , (9.3)
in quanto, per la stessa definizione di corpo rigido vR = 0. La precedente è detta formula del
Mozzi. La − →
ω , che chiaramente è funzione del tempo, è la velocità angolare propria del corpo
rigido e dipende esclusivamente dal suo moto. Infatti, supponiamo di prendere un altro punto
di riferimento del corpo rigido O ′′ , invece di O. Rispetto ad O ′′ , si avrà
vP = vO′′ + −
→
ω ′′ ∧ (P − O ′′ ) , (9.4)
ma dovrà essere −
→
ω ′′ = →
−
ω poiché O ′′ è pur sempre un punto del corpo, per cui anche per lui
vale la
vO′′ = vo + −→ω ∧ (O ′′ − O) . (9.5)
Sottraendo la (9.5) dalla (9.3), si ha
vP − vO′′ = →
−
ω ∧ (P − O) − − →
ω ∧ (O ′′ − O)
= −
→
ω ∧ (P − O ′′ ) , (9.6)
aP = ao + −̇
→
ω ∧ (P − O) + −
→
ω ∧ [−
→
ω ∧ (P − O)] (9.7)
9.2 Considerazioni su ω e vo
Facciamo due esempi banali di moti rigidi, per capire il ruolo della velocità angolare ω e di vo .
• Se −
→ω = 0, avremo che, presi in qualsiasi modo P ed O appartenenti al corpo rigido,
vP = vo ; cioè tutti i punti si muoveranno con la stessa velocità. Il moto che ne deriva si
dice traslatorio.
• Se vo = 0, avremo vP = − →ω ∧ (P − O).
Con −→ω costante, i punti della retta passante per il centro O e parallela ad − →ω hanno
velocità nulla: il corpo rigido ruota attorno a questa retta (asse di rotazione); i punti del
corpo rigido formano traiettorie circolari, che giacciono su piani perpendicolari all’asse
di rotazione, con centro il punto di intersezione del piano con l’asse e velocità angolare
k ω k.
Con −→
ω non costante abbiamo un moto di precessione, in cui, istante per istante, c’è un’asse
di rotazione parallela ad −→ω , attorno alla quale si compie un atto di moto rotatorio.
129
9.3 Prima Equazione Cardinale
Siccome il corpo rigido è pur sempre un sistema di punti materiali, possiamo ripetere ciò che
abbiamo detto in occasione dello studio della Prima Cardinale dei sistemi di punti materiali.
La risultante delle forze sarà anche in questo caso data dalla sola risultante delle forze esterne
e avremo
Rext = Q̇ = Macm , (9.8)
dove per un sistema di punti discreto
P
mi ri
rcm = Pi , (9.9)
i mi
ri ∧ mi (vO + −
→
X
= ω ∧ ri ) (, 9.14)
i
mi ri ∧ (−
→
X
= ω ∧ ri ) . (9.15)
i
LO = LO k + LO ⊥ , (9.16)
130
dove
−
→
" ! #
ω →
− mi ri ∧ (−
→
X
LO k = LO · − ω = ω ∧ ri ) · k̂ k̂ , (9.17)
k→ωk i
" ! #
X
= mi ri ∧ ωri sin (θi ) ûθi · k̂ k̂ , (9.18)
i
!
X
= mi (ri sin (θ))2 ω k̂ , (9.19)
i
mi h2i −
→
X
= ω. (9.20)
i
La quantità X
I= mi h2i , (9.21)
i
si chiama momento d’inerzia, rispetto all’asse z passante per O, del corpo rigido e mi “descrive”
l’inerzia del corpo a mettersi in rotazione una volta sollecitato da un momento risultante non
nullo (vedremo più in dettaglio il significato nei paragrafi successivi).
La componente perpendicolare, invece, è
" ! #
m r ∧ (−→
X
L O⊥= (L · û ) û =
O r r ω ∧ r ) · û û ,
i i i r(9.22)
r
i
X
= − mi ri2 sin (θi ) cos (θi )ω ûr . (9.23)
i
Nel caso di corpo rigido simmetrico rispetto all’asse di rotazione, il contributo di ogni
punto materiale ad LO ⊥ viene cancellato, nella somma vettoriale, da quello del suo simmetri-
co. Quindi, in questo caso, rimane solamente la componente parallela all’asse z del momento
angolare
LO = LO k = I −
→ω. (9.24)
In generale, invece
LO = LO k + LO ⊥ = I −
→
ω + LO ⊥ . (9.25)
Siccome il sistema è rigido, tutti i vettori
131
9.5 Momento d’Inerzia e Teorema di Huygens-Steiner
Il momento d’inerzia gioca nelle rotazioni un ruolo analogo a quello della massa inerziale nelle
traslazioni. Ovvero, rappresenta l’“inerzia” a mettersi in moto (rotatorio) sotto la sollecitazione
di un momento.
Per un sistema rigido discreto abbiamo che
X
I= mi h2i , (9.28)
i
dove hi è la distanza dall’asse di rotazione del punto i-esimo. Se il sistema è continuo, invece,
si avrà un integrale: Z
I= h2 dm . (9.29)
V
Se infine il corpo continuo è omogeneo con densità costante ρ, si può scrivere
Z
I= ρ h2 dV . (9.30)
V
Pi
r′i ri
rO O
cm ≡ O ′
y ≡ y′
x′ x
Supponiamo che il corpo rigido sia formato da un insieme discreto di punti (se si tratta di
un continuo dovremo semplicemente passare agli integrali, invece che utilizzare le somme). Il
132
contributo al momento angolare del punto i-esimo, di massa mi , sarà dato, nel sistema SdR’ da
Ii = mi (x′2 + y ′2 ) (9.32)
e quindi X
IA1 = Icm = mi (x′i 2 + yi′ 2 ) , (9.33)
i
Siccome r′i = rO + ri , la relazione fra le componenti di Pi nei due sistemi di riferimento è data
dalle seguenti equazioni
′
xi = xi ,
yi′ = yi + d , (9.35)
′
zi = zi .
Quindi si ha
X X
IA 2 = mi (x2i + yi2 ) = mi [x′i 2 + (yi′ − d)2 ] , (9.36)
i i
X X
= mi (x′i 2 + yi′ 2 ) + Md2 − 2d mi x′i , (9.37)
i i
X
= mi (x′i 2 + yi′ 2 ) + Md , 2
(9.38)
i
= Icm + Md2 , (9.39)
poiché i mi x′i = 0, essendo la “x” del centro di massa nel sistema di riferimento centrato sul
P
centro di massa stesso (quindi dove il cm coincide con l’origine).
Corollario 9.5.2 Se A1 è un asse non passante per il centro di massa e distante d1 da un asse
parallelo passante per cm e A2 è un asse, parallelo ad A1 e distante d2 dall’asse passante per
cm, si ha
IA1 = IA2 + M(d21 − d22 ) . (9.40)
e analogamente
IA2 = Icm + Md22 . (9.42)
Ricavando quindi Icm dalla (9.42) e sistituendolo nella (9.41) si ottiene la .
133
9.5.1 Una formula importante per i sistemi piani
Supponiamo di avere un sistema rigido piano (cioè bidimensionale). Se l’asse di rotazione è
perpendicolare al piano, potremo identificarla come l’asse z di un sistema cartesiano con gli
assi x e y sul piano del corpo rigido. Allora
X
Iz = mi (x2i + yi2) . (9.43)
i
Se ora consideriamo la rotazione del corpo rispetto all’asse x o risspetto all’asse y, però, abbiamo
X
Ix = mi yi2 (9.44)
i
e X
Iy = mi x2i . (9.45)
i
Quindi
Iz = Ix + Iy . (9.46)
−
→
ω
z
B
l
α
l x
134
Si ha
l l
dx
Z Z
2
I = x dm = ρ x2 dl = abbiamo che dl =
−l −l cos (α)
xB l cos (α)
ρ ρ l2
Z Z
2
= x dx = x dx = m cos2 (α) .
2
(9.48)
cos (α) xA cos (α) −l cos (α) 3
−
→
ω
z
Si ha:
Z Z
2
I = r dm = ρ r 2 ds = | abbiamo che ds = r dθ |
γ γ
Z 2π
= ρ r 3 θ = 2πρ r 3 = mr 2 . (9.49)
0
−
→
ω
z
Si ha:
R 2π R
R4 R2
Z Z Z Z Z
2 2 2
I= r dm = r ρ dS = ρ r ρr dr dθ = 2πρ r 3 dr = 2πρ =M . (9.50)
S S 0 0 0 4 2
135
Banderuola triangolare di densità superficiale costante
Cilindro di densità volumetrica costante
Sfera di densità volumetrica costante
LO = LO k + LO ⊥ = I −
→
ω + LO ⊥ . (9.51)
Mext
0 = L̇O . (9.52)
ext d−
→ω dLO ⊥
M0 = I + , (9.55)
dt dt
dove
!
dLO ⊥ d X
= − mi zi ωhi ûr , (9.56)
dt dt i
X X dûr
= − mi zi hi ω̇ ûr − − mi zi ωhi , (9.57)
i i
dt
ω̇
= LO ⊥ + −
→
ω ∧ LO ⊥ . (9.58)
ω
In totale
d−
→ω ω̇
Mext
0 = I + LO ⊥ + −→
ω ∧ LO ⊥ . (9.59)
dt ω
Se il corpo rigido è simmetrico rispetto all’asse di rotazione (cioè questa è un asse di
simmetria) e Mext
0 = 0, allora si conserva il momento angolare rispetto a questo asse:
ovvero
I θ̇ = cost . (9.61)
Il sistema continua a ruotare a θ̇ = cost intorno all’asse di simmetria indefinitamente (asse
libero di rotazione).
136
9.7 Moti rigidi piani
Il moto descritto nella prima parte del secondo punto del paragrafo precedente, cioè in cui
• →
−
ω = costante
• vo · −
→
ω =0
si dice moto rigido piano. I punti del sistema hanno velocità perpendicolare alla direzione
di ω. Anche se non si tratta di un corpo rigido piatto (che non esiste), il moto del sistema
resta descritto da quello dei punti di un piano perpendicolare ad ω, che viene detto piano
rappresentativo.
9.7.1 Volano
Consideriamo il sistema in figura, costituito da un manicotto di massa m (raggio r), un disco di
massa M (raggio R) e un altro manicotto di massa m (raggio r) coassiali, rigidamente collegati,
che ruotino attorno all’asse di simmetria del sistema (asse z). Su uno dei due manicotti, viene
avvolta una fune ideale collegata ad una massa m′ , soggetta alla forza di gravità. Lasciando
il sistema da fermo, il volano (l’insieme dei due manicotti col disco) si metterà a ruotare,
sollecitato dalla massa m′ che tira la corda.
m m
z M
−
→
g
M R
r
m′
m
m′
m′ ÿ = m′ g − T . (9.63)
137
T
P = m′ g
y
Rn
T P
m′ (g − r θ̈)r = I0 θ̈ , (9.67)
ovvero
m′ gr 2m′ gr
θ̈ = = Ω. (9.68)
I0 + m′ r 2 MR2 + 2(m + m′ )r 2
Il moto è quindi uniformemente accelerato.
Con la prima equazione cardinale possiamo trovare quali sono le reazioni del vincolo (l’asse
che tiene sospeso il volano).
138
9.7.2 Pendolo fisico
Consideriamo un sistema rigido bidimensionale (o tridimensionale con un asse fisso di rotazione)
che ruoti attorno ad un asse fisso passante per O (perpendicolare al piano della figura). Sia il
centro di massa non coincidente con O. Sia inoltre il sistema soggetto alla forza di gravità.
O
−
→
g
θ
c.m.
Mg
Per trovare l’equazione del moto, applichiamo la sseconda cardinale nel punto O. Il mo-
mento risultante delle forze esterne rispetto ad O risulta essere soltanto quello dovuto alla forza
peso, che possiamo schematizzare come agente nel centrop di massa (c.m.) del corpo rigido.
Prendiamo come verso di rotazione positivo quello antiorario. Quindi:
dove per k̂ si intende il versore perpendicolare al foglio che punta nell’occhio dell’osservatore.
La seconda cardinale sarà allora
ovvero
Mgrcm
θ̈ + sin (θ) = 0 . (9.71)
IO
La (9.71) corrisponde al moto di un pendolo semplice con “lunghezza equivalente” pari a
IO
leq = . (9.72)
Mrcm
Nel limite delle piccole oscillazioni, avremo che la pulsazione sarà data da
r
Mgrcm
ω= (9.73)
IO
e quindi il periodo da s
2π IO
T = = 2π . (9.74)
ω Mgrcm
139
9.7.3 Pendolo di torsione
Il pendolo di torsione è costituito da un disco appeso ad un filo, che reagissca in maniera lineare
alla torsione (sottile filo d’acciaio o filo di quarzo).
−
→
ω
z
Ruotando il disco attorno al suo asse (coincidente col filo), si sposta dalla situazione di
equilibrio, tocendo il filo. Se θ è l’angolo di tale rotazione (rispetto ad una direzione fissa nel
piano del disco), il filo reagisce con un momento proporzionale a θ (Legge di Hooke):
Mz = −kθ , (9.75)
che si oppone a tale torsione. L’equazione del moto allora sarà un’equazione del moto armonico
in θ
Iz θ̈ = −kθ , (9.76)
dove Iz = MR2 /2, ovvero
2k
θ̈ + θ = 0. (9.77)
MR2
θ̇+
r
O
x
Per ricavare una formula cinematica che caratterizzi il moto (velocità istantanea di ogni
punto P del disco), possiamo applicare la formula del Mozzi identificando O col punto di
contatto fra disco e vincolo (che è istantaneamente fermo). Si ha, allora
v =−
P
→ω ∧ (P − O) , (9.78)
140
che determina la distribuzione di velocità istantanea in figura
θ̇+
C vC
x
O
In particolare la velocità del centro del cerchio è legata alla velocità angolare θ̇ dalla
vC = ẋC = r θ̇ , (9.79)
che risponde all’intuizione che ci dice che, se il disco non può strisciare ma solo rotolare, lo
spostamento che fa il centro è esattamente l’arco di circonferenza r∆θ di cui ha ruotato il
disco.
y
θ̇+
−
→
g
C
Mg
Rn
Rt
α x
Abbiamo posto la reazione vincolare tangente al piano nel verso negativo delle x. Ciò
deriva dal fatto che effettivamente la forza d’attrito, per evitare che il disco strusci sul piano,
deve essere diretta lungo il piano e avere quel verso. Comunque, è da notare che potremmo
anche posizionare la forza nel verso opposto e scrivere le equazioni di conseguenza. Alla fine,
troveremmo un Rt negativo, che ci dice che il verso è effettivamente quello della figura.
141
Applichiamo le cardinali della dinamica per descrivere il moto del disco. La prima cardinale,
scomposta lungo le direzioni x e y degli assi coordinati, dà le seguenti due equazioni:
Rt r = IC θ̈ . (9.82)
ẋC = r θ̇ (9.83)
e quindi
IC
Mg sin (θ) ≤ µS Mg cos (θ) . (9.88)
Mr 2
+ IC
Da questa relazione possiamo per esempio trovare il valore minimo per µS , dato l’angolo θ, per
garantire il rotolamento puro:
IC
µS ≥ tan (θ) . (9.89)
Mr 2 + IC
142
Rotolamento puro con momento (motore)
Consideriamo adesso il caso in cui si abbia una ruota che rotola su un piano sotto l’azione di
un momento esterno τ , che considereremo agente in C (schematizzazione della ruota soggetta
al momento dato da un motore, per esempio) e tale da far rotolare il disco nel senso positivo
dell’angolo θ. Le forze che agiscono sul disco sono schematizzate in figura. Questa volta la
reazione vincolare tangente al piano è ha verso positivo, perchè si deve opporre ad un possibile
slittamento del disco sotto l’azione del momento τ .
θ̇+
C
Mg
Rn Rt
x
La seconda cardinale in C dà
τ − Rt r = IC θ̈ , (9.90)
mentre la prima cardinale scomposta lungo le due direzioni x e y dà
Rn = Mg , (9.91)
Rt = M ẍC = Mr θ̈ , (9.92)
143
Attrito volvente
Si parla di attrito volvente quando una ruota invece di toccare il piano in un punto solo, lo
tocca lungo una zona estesa. Esagerando, si può schematizzare la situazione come in figura.
θ̇+
C
Mg Rn
δ
L’attrito volvente sarà dato da un momento M che si oppone al moto, pari in modulo a δRn ,
dove δ è un parametro caratteristico (una lunghezza) che schematizza il fatto che la reazione
normale non avviene più lungo il raggio che collega C perpendicolarmente al piano:
M = −δRn . (9.97)
L’attrito volvente, schematizzato in questo modo, incide solo sulla seconda cardinale.
144
9.8.1 Asse fisso
Vediamo cosa succede nel caso semplificato in cui il corpo sia in rotazione attorno ad un asse
fisso. Prendiamo il punto 0 su tale asse e quindi v0 = 0.
→
−
ω
z
LO,i
θi Pi
hi
θi
ri
O y
Allora
1X 1X
T = mi (−
→
ω ∧ r i ) · (−
→ω ∧ ri ) = mi (ωri sin (θi ))2 , (9.101)
2 i 2 i
!
1 X 1X 1
= mi ri2 sin2 (θi ) ω 2 = mi h2i ω 2 = Iω 2 , (9.102)
2 i
2 i 2
dove
fiext · (vcm + −
→
X X
6 dL = fiext · vi dt = ω ∧ ri )dt , (9.104)
i i
fiext · −
→
X X
= ( fiext ) · vcm dt + ω ∧ ri dt , (9.105)
i i
X X
= ( fiext ) · vcm dt + ·(ri ∧ fiext )dt , (9.106)
i i
145
→
−
= Rext · drcm + Mext
cm · d θ . (9.107)
I due contributi sono dunque un lavoro dovuto al movimento di traslazione del centro di massa
e uno dovuto alla rotazione attorno a − →ω.
Nel caso di rotazione attorno ad un asse fisso, quindi, prendiamo 0 sull’asse (quindi v0 = 0)
e avremo che il lavoro fatto per passare dalla configurazione A alla configurazione B sarà dato
da Z B Z B Z B
ext →
−
LAB = 6 dL = M0 · d θ = Mkext dθ . (9.108)
A A A
ri ∧ mi (vcm + −→
X X
Lcm = ri ∧ mi vi = ω ∧ ri ) , (9.109)
i i
r i ∧ (→
−
X
= ω ∧ ri ) , (9.110)
i
poiché i ri ∧ mi vcm = 0.
P
Svolgendo il doppio prodotto vettoriale, si trova
mi ri2 −
→ mi (ri · −
→
X X
Lcm = ω − ω )ri . (9.111)
i i
ri = xi î + yi ĵ + zi k̂ , (9.112)
−
→
ω = ωx î + ωy ĵ + ωz k̂ . (9.113)
Analogamente si trova
X X X
Lcm,y = mi (x2i + zi2 )ωy − mi yi zi ωz − mi xi yi ωx , (9.116)
i i i
X X X
Lcm,z = mi (x2i + yi2 )ωz − mi zi yi ωy − mi xi zi ωx . (9.117)
i i i
146
Individuando la matrice
Ix −Ixy −Ixz
σ = −Ixy Iy −Iyz , (9.118)
−Ixz Iyz −Iz
dove Ix , Iy e Iz sono rispettivamente i momenti d’inerzia rispetto all’asse x, y e z:
X
Ix = mi (yi2 + zi2 ) , (9.119)
i
X
Iy = mi (x2i + zi2 ) , (9.120)
i
X
Iz = mi (x2i + yi2 ) , (9.121)
i
e dove gli altri sono i momenti misti:
X
Ixy = mi xi yi , (9.122)
i
X
Ixz = mi xi zi , (9.123)
i
X
Iyz = mi yi zi , (9.124)
i
si trova
Lcm = σ −
→
ω). (9.125)
La matrice σ viene detta “Tensore d’inerzia”. è una matrice 3×3 simmetrica e quindi diagona-
lizzabile.
Rispetto ad un SdR che abbia come terna gli autovettori di σ, normalizzati (versori), v̂1 ,
v̂2 , v̂3 , potremo porre la (9.125) in forma diagonale
Lcm = I1 ω1 v̂1 + I2 ω2 v̂2 + I3 ω3 v̂3 . (9.126)
I tre assi individuati dagli autovettori di σ sono detti “assi principali d’inerzia”.
147
9.12 Espressione generale dell’energia cinetica
Consideriamo adesso l’energia cinetica in Eq. (9.100) e consideriamo il caso in cui 0 ≡ c.m..
Avremo in generale
1 1X
T = Mvcm 2
+ mi (−
→
ω ∧ r i ) · (−
→
ω ∧ ri ) . (9.134)
2 2 i
Consideriamo il secondo termine e scriviamo la decomposizione di ri e →
−
ω rispetto ad un SdR
generico centrato nel centro di massa del corpo rigido:
ri = xi î + yi ĵ + zi k̂ , (9.135)
−
→
ω = ωx î + ωy ĵ + ωz k̂ . (9.136)
Allora
î î k̂
−
→
ω ∧ ri = ωx ωy ωz , (9.137)
xi yi zi
= (zi ωy − yi ωz )î + (xi ωz − zi ωx )î + (yi ωx − xi ωy )k̂ . (9.138)
Quindi
(−
→
ω ∧ ri )2 = (zi ωy − yi ωz )2 + (xi ωz − zi ωx )2 + (yi ωx − xi ωy )2 , (9.139)
2 2 2 2 2 2 2 2 2
= (yi + zi )ωx + (xi + zi )ωy + (xi + yi )ωz − 2xi yi ωx ωy − 2xi zi ωx ωz − 2yi zi(9.140)
ωy ωz ,
che è una forma quadratica (detta Ellissoide d’inerzia) che si può scrivere nella seguente forma
matriciale:
1 1 →t −
T = Mvcm 2
+ − ω σ→ω. (9.141)
2 2
La matrice σ è il tensore d’inerzia (9.118).
148