Meccanica Razionale per Allievi Ingegneri “Aerospaziali” Esercitazioni numeriche 2006

1. Meccanica Razionale per Allievi Ingegneri “Aerospaziali”Esercitazioni numeriche 2006 Realizzazione dei “Ritratti in Fase” per alcuni problemi classici della Meccanica

2. Obiettivi del Corso Scopo del presente corso è fornire agli allievi strumenti concettuali ed un primo training per lo studio di problemi della Meccanica Razionale con i metodi del calcolo computerizzato. In particolare si fornisce una introduzione alla dinamica di sistemi non lineari e di sistemi dotati di cicli limite e si propone agli allievi l’uso critico di integratori numerici per tracciare i ritratti in fase di alcuni dispositivi classici particolarmente rappresentativi.

3. Programmazione “ben fatta” • Un programma di calcolo “ben fatto” ha un punto di inizio e una condizione di termine ben definiti • L’ obbiettivo del programma è definito ed enunciato (espressamente) • Sono assegnati i dati di ingresso e sono ben definiti i risultati attesi

4. Primi programmi • Programma uno: stampare un titolo con nome e matricola dell’ allievo • Programma due: scrivere su un file le coordinate di nove punti equidistanti su una circonferenza ( usare le funzioni sin e cos) • Programma tre : come due, ma ordinando i punti in modo che siano i vertici di una stella, una volta congiunti successivamente con corde opportune.P.es. mediante l’uso del software “gnuplot” in ambiente Linux o mediante applicazioni Excel in ambiente Window .

5. Equazioni Differenziali Ordinarie • Metodi numerici per ”Problemi di Caucy” (vedi il corso di Analisi2). • Vogliamo considerare equazioni differenziali: • - lineari: particella libera, moto per gravità, crescita Malthusiana, decadimento radioattivo, oscillatore e repulsore armonico; • - non lineari: pendolo, crescita di popolazioni con risorse limitate, moto kepleriano, bipendolo; • -dotate di ciclo limite: orologio meccanico, problema di Lotka-Volterra, problema di Van Der Pol • Con esiti caotici: pendolo forzato, mappa standard

6. O.D.E Lineari. Part.Libera • Particella libera: d2x/dt2=0 ; Per la soluzione introduciamo la variabile v, velocità, definita da v= dx/dt: • l’equazione del secondo ordine (d2x/dt2=0) diviene una coppia di equazioni del primo ordine: dv/dt=0 ; dx/dt=v ; con le condizioni iniziali x(0)=x0 , v(0)=v0 • Vedi più avanti le traiettorie rettilinee parallele nel piano delle fasi.

7. Moto per gravità • d2x/dt2 = g. Come sopra introduciamo la variabile v, velocità, definita da v= dx/dt: • l’equazione del secondo ordine diviene una coppia di equazioni del primo ordine: dv/dt=g; dx/dt=v ; con le condizioni iniziali x(0)=x0, v(0)=v0. • Più avanti sarà richiesto di tracciare anche per questo caso le traiettorie nel piano delle fasi

8. Oscillatore e Repulsore Armonico L’equazione d2x/dt2 = - k /m x , detta dell’oscillatore armonico, descrive il comportamento di un punto materiale soggetto al richiamo di una forza elastica. L’equazione è rappresentativa del comportamento di sistemi fisici in prossimità della posizione di equilibrio stabile Se il coefficiente del termine lineare fosse invece positivo, si avrebbe il “repulsore armonico”, rappresentativo del comportamento di sistemi fisici in prossimità degli equilibri instabili. Ne è esempio il moto di un punto in un campo centrifugo, p.es. il moto relativo di un punto materiale su una giostra.

9. Crescita Malthusiana e decadimento nucleare • I fenomeni sono descritti dalla equazione differenziale ordinaria lineare del primo ordine dy/dt = A y , con A positivo nel caso della crescita di popolazioni viventi. A rappresenta il numero di nuovi nati nella unità di tempo (vedi anche il modello di crescita del capitale bancario). • Per A negativo il modello può rappresentare il decadimento radioattivo o lo svuotamento di una spugna (o di uno strato di terreno) impregnati di acqua (o l’evoluzione di una popolazione senile in cui A è il numero di decessi nella unità di tempo) • Gli esiti numerici dei due problemi possono essere sorprendenti in termini di precisione di calcolo e di dominanza degli errori numerici sulla vera soluzione

10. Sviluppo di Taylor per il problema di Cauchy: crescita Malthusiana. • Gli allievi hanno studiato ampiamente la teoria delle serie e delle successioni nei corsi precedenti. Il Calcolo Numerico richiede effettive utilizzazioni di tale teoria. Invitiamo a rileggere nei testi di analisi la teoria dello sviluppo in serie di Taylor-McLaurin per le funzioni più comuni • Come primo esempio si consideri l’equazione dy/dt = A y , rappresentativa per numerosi problemi concreti, dagli interessi bancari fino alla crescita delle popolazioni biologiche. • Nel nostro esempio, semplificato e standardizzato, si ponga la condizione iniziale y(0) = 1 ed anche il coefficiente A sia costante (ed uguale ad 1). Vogliamo capire cosa significa il termine “ordine di precisione del solutore numerico”.

11. Il grafico dell’ esponenziale • Si cerca y(T). • Cerco uno sviluppo di Taylor, quindi ho bisogno delle successive derivate di Y, calcolate nel punto iniziale. • A è costante, da dy/dt = A y, si ha: d2y/dt2 = A dy/dt = A2 y; d3y/dt3 = A d2y/dt2 = A2 dy/dt = A3 y . Quindi:y(T)= y(0) ( 1+ AT + ½ (AT)2+1/6 (AT)3+ 1/24…). Con le consuete osservazioni per la convergenza della serie (prodotto AT < 1) e per la valutazione del resto.

12. Accorgimenti per il calcolo Una volta scelto il numero di termini della serie, il calcolo va organizzato evitando di calcolare grandi fattoriali con i loro inversi e potenze successive di AT, ma raccogliendo i fattori comuni così da ridurre l’errore di troncamento che si verifica quando si sommano termini via via più piccoli y=yo (1+AT(1+1/2 AT(1+1/3 AT(1+1/4 AT(1+1/5 AT(…). Come esercitazione richiesta si costruisca una procedura che chieda in input l’ordine massimo e realizzi il prodotto dei termini annidiati a partire dal più interno. In uscita la procedura fornirà, oltre al risultato approssimato voluto, anche una stima del primo termine trascurato.

13. Esempi : passo passo • Per valutare y a tempi t > R (raggio di convergenza), è possibile valutarla per t < R e quindi scegliere il nuovo valore di y come nuovo punto iniziale per lo sviluppo in serie di Taylor. L’espressione dello sviluppo in serie resta la medesima ed è sufficiente aggiornare i valori della funzione e delle derivate • La scelta del valore di y al termine del primo passo come valore iniziale per un passo successivo illustra la filosofia generale dei metodi numerici che intendiamo adottare. • Attenzione: gli errori commessi in ciascun passo vanno sommati in valore assoluto per ottenere la stima dell’errore complessivo

14. DA SVOLGERE PER L’ESAME • Può convenire accorciare il passo elementare, aumentando il numero dei passi, e simultaneamente riducendo il numero dei termini dello sviluppo in serie. La scelta dipende dal problema particolare ed è conveniente quando diviene grande l’onere del computo delle successive derivate. • La semplicità dell’esempio suggerito consente almeno in questi primi casi di valutare la soluzione dell’equazione proposta dapprima con passi piuttosto lunghi, prossimi ad R, e numerosi termini della serie di Taylor, poi con serie di Taylor più brevi e con un maggiore numero di brevi passi di calcolo.Si può scegliere A=1,T=0.5. Sviluppo fino alla 4 o alla 6 potenza, poi, per confronto T =1, T=0.25,T=0.125

15. Esempi : seni e coseni • Consideriamo ora la coppia di equazioni: dx/dt= -y ; dy/dt= x ; con x(0)= 1 ; y(0)=0. • È immediato calcolare le derivate successive delle due espressioni e quindi ottenered2x/dt2= - dy/dt= - x d2y/dt2= dx/dt= - yd3x/dt3= - dx/dt= y d3y/dt3= - dy/dt= -x eccetera. • Svolgendo quindi lo sviluppo in serie di Taylor:x(T)= x(0) – y(0) T – ½ x(0) T2 + 1/6 y(0) T3 …… y(T)= y(0) + x(0) T – ½ y(0) T2 – 1/6 x(0) T3 ……

16. Esempi : seni & coseni, segue • E quindi si possono riorganizzare i termini:x(T)=x(0) (1-1/2 T2+1/(2x3x4)T4–1/(2x3x4x5x6) T6 ….)– y(0)T(1-1/(2x3)T2 +1/(2x3x4x5)T4–1/5040 T6….)y(T)=y(0) (1-1/2 T2 + 1/24 T4 – 1/720 T6 ….) + + x(0)T(1-1/6 T2 +1/120 T4 – 1/5040 T6 ….) • Anche ora è opportuno trasformare le somme di termini rapidamente decrescenti in prodotti come p.esx(T)=x(0) (1-(1/2 T2(1-1/(3x4)T2(1+1/(5x6) T2 (….)– y(0)T(1-1/(2x3)T2(1+(1/(4x5)T2(1-1/(6x7)T2 (….)etc... e quindi valutare se la scelta di compiere numerosi passi di piccola ampiezza può essere complessivamente vantaggiosa in termini di precisione dei risultati.

17. Seni & coseni, integrali primi Risulta particolarmente utile, nello studio dei sistemi dinamici, considerare gli eventuali “Integrali Primi” o “Costanti del moto” ammessi dal sistema. Per la coppia di equazioni dx/dt= -y ; dy/dt= x ; si può notare, moltiplicando la prima per x e la seconda per y e sommando i prodotti, x dx/dt= -x y ; y dy/dt= y x ; d(x2)/dt +d(y2)/dt= 0 da cui si deduce che deve restare costante la somma x2+y2, ed uguale al valore iniziale ( in questo caso x(0)2+y(0)2=1 Valutare nel corso della integrazione numerica tale espressione consente di avere una indicazione dell’errore numerico complessivamente commesso.

18. Metodi numerici generali • Gli argomenti presentati oggi servono per costruire un linguaggio didattico comune • Dovrebbero essere tutti a voi completamente noti e serviranno nelle prossime lezioni per costruire strumenti di lavoro più sofisticati ed efficaci, • Ma dovrebbero altresì essere utilizzabili per verificare e per scegliere metodi risolutivi adeguati ai problemi che si porranno alla vostra attenzione.

19. DA PORTARE ALL’ESAMEPer casa: • Organizzare il calcolo della funzione esponenziale per passi successivi, impiegando sviluppi in serie di Taylor a 4, a 6 ed a 10 termini, • Con passi temporali di ampiezza prossima ad R fino a 4R , poi con passi di R/2 e di R/8 • Illustrare brevemente i risultati in termini di precisione e rapidità di calcolo. • Organizzare similmente il calcolo di seni e coseni fino a raggiungere 4 pigreco e 8 pigreco

20. Oscillatore armonico 1 • L’esempio sviluppato per seni e coseni si trasferisce immediatamente alla descrizione del sistema dinamico chiamato “ oscillatore armonico” • La realizzazione fisica elementare è rappresentata da un punto materiale dotato di massa, soggetto a forza centrale elastica • Inizialmente lo immaginiamo vincolato su una retta. Si potrà constatare che il moto deve essere piano e che opportune condizioni iniziali consentono il moto rettilineo.

21. Oscillatore Armonico 2 • L’equazione del moto può essere scritta : m d2x/dt2 = - k x • Con m massa, k costante elastica, origine delle coordinate nel centro della forza elastica e coordinate iniziali compatibili con il moto rettilineo che si intende descrivere.

22. Velocità; equazioni del primo ordine. • Introduco la variabile v =dx/dt • Il sistema di O.D.E. diviene • mdx/dt=mv • mdv/dt=-kx • Una prima trasformazione usualmente operata è dividere entrambe le equazioni per il fattore m (massa: costante positiva)

23. Normalizzazione • Si ottiene così: dx/dt = v dv/dt = - k/m xAnche k/m è una costante positiva e ciò permette di introdurre il nuovo parametro ω con ω2 = k/m.La seconda equazione diviene: dv/dt= -ω2x. Se cambio l’unità di misura dei tempi tale equazione diviene dv/dt= - x dove in luogo dei secondi utilizzo come unità dei tempi la “pulsazione” ω.

24. Condizioni iniziali • L’ equazione dunque può scriversi : d2x/dt2 = -x • Oppure dx/dt = v dv/dt = -x Alla equazione vanno aggiunte le condizioni iniziali, che consentono di selezionare una soluzione particolare nella famiglia (a due parametri) di tutte le soluzioni possibili.

25. Serie di Taylor per le funzioni circolari • Scelgo come come condizioni iniziali i valori di x e di v=dx/dt al tempo iniziale, t0 e li chiamo x0 e v0 . Posso subito esprimere la soluzione nella forma dello sviluppo in serie di Taylor. • x = x0 + dx/dt (t- t0) + ½ d2x/dt2 (t- t0)2 + + (½ *1/3) d3x/dt3 (t- t0)3 + . . . con le successive derivate espresse dalla equazione originale: d2x/dt2 = dv/dt = -x; d3x/dt3=- dx/dt =-v e calcolate nel tempo t0

26. . • Si ottiene : • x(t) = x0 + v0 (t- t0) - ½ x0 (t- t0)2 + -(1/6) v0(t- t0)3 + 1/24 . x0 (t- t0)4 + . . . cioè • x(t) = x0 ( 1- (t- t0)2 (1/2-1/12 (t- t0)2 (….))) + v0 (t- t0) ( 1- 1/6 (t- t0)2 (…))) • Si riconoscono nelle espressioni che moltiplicano x0 e v0 le espressioni di Taylor per cos(ωt) e sin(ωt), con ω = 1 per la particolare scelta effettuata per l’unità di misura del tempo. • Si noti che abbiamo scelto anziché la forma di somma di termini (rapidamente) decrescenti, la forma di prodotto di fattori sempre più prossimi ad 1, per ragioni di precisione numerica complessiva.

27. DA PORTARE ALL’ESAMEIl programma richiesto • Le espressioni sopra individuate permettono di valutare x(t) e v(t) per successivi valori di t. Si suggerisce di scegliere l’intervallo 0<=t<=10 ed in tale intervallo 100 punti. • Si costruisca un file testo con 101 record di terne t,x(t),v(t) e, mediante il software gnuplot, si rappresentino su un grafico le curve x(t) ed y(t) sovrapposte e su un altro grafico la curva ottenuta unendo tutte le coppie (x,t) di ciascun tempo successivo.

28. Metodo di Eulero • A partire dal punto P(x(0),v(0)) = P, si valuti il vettore ΔP = (dx/dt*Δt, dv/dt*Δt), mediante le espressioni dx/dt=v;dv/dt=-x. Si valuti ora il nuovo punto P1 = P(0+Δt) = P0 + ΔP. • Il vettore ΔP è ortogonale a OP ( il prodotto scalare dei due è x dx+v dv = xv – vx =0 ) ed è tangente alla traiettoria. Ripetere la costruzione ottenendo P2=P(0+2Δt). • Nel programma (programma 5) si registrino sul file uscita le successive coppie di coordinate relative ai successivi valori t+n Δt e si tracci la traiettoria così ottenuta, confron-tandola con quella ottenuta mediante lo sviluppo di Taylor.

29. Eulero con correttore • Nel metodo di Eulero l’arco di traiettoria percorso nell’intervallo t -> t+ Δt viene approssimato con un segmento rettilineo tangente alla traiettoria in t. è una approssimazione del primo ordine e l’errore vale circa derivata seconda per dt2 . • Una approssimazione migliore si otterrebbe sostituendo al segmento di retta un arco di parabola, che oltre a partire da P(t) con tangente dP(t), giungesse a P(t+Δt ) con tangente dP(t+Δt). • Tale risultato potrebbe essere ottenuto se si potesse calcolare la tangente nel punto P(t+1/2Δt)

30. Predictor-corrector • Proviamo ad utilizzare come valore approssimato della tangente nel punto P(t+Δt) il valore dato dalla derivata calcolata in P+ΔP. Si valuti un nuovo valore di ΔP come valor medio di quello calcolato al tempo t e di quello valutato nella approssimazione al primo ordine. L’errore che si commette ora è del terzo ordine e consente di ricalcolare la traiettoria con errore ancora minore. • Il metodo complessivo può essere definito un “predictor-corrector-evaluation-corrector”, esemplare di metodi più elaborati che si studieranno in seguito.

31. Implementazione • STAMPA DI TITOLO ED AUTORE • Float tempo,tmax,delt,x,v,u,w,xa,va • Integer i,j,k,npunti,ndiv,nru • Leggi tmax, npunti,ndiv,nru • delt=0.5*(tmax/npunti)/ndiv • Leggi x0,v0 ;Tempo=0 ; x=x0; v=v0 • (Oscar: ) dx=v;dv=-x • For i=1 to npunti ;Registra tempo,x,v,dx,dv • For j=1 to ndiv ;xa=x+dx*delt;va=v+dv*delt; • For k=1 to nru ;x=xa+dx*delt;v=va+dv*delt • (Oscar:) dx=v;dv=-x; Next k; Next j; Next i; FINE

32. Suggerimenti per l’esecuzione • xa e va sono variabili ausiliarie. Si userà la medesima struttura di programma anche per un altro metodo. • Si scelga inizialmente tmax=10, npunti=100, ndiv =1, 2, 4, 8, 16. • Nru =1 realizza il metodo di Eulero; = 2 Eulero con correttore; =3 con ulteriore correzione • Si chiede di confrontare le traiettorie ottenute con nru=1 e con ndiv diversi e le traiettorie con ndiv=1 ed nru diversi e scegliere e motivare la scelta combinata preferibile.

33. Serie di Taylor (complementi) Si chiede di trovare sull’intervallo [to, to+T] la soluzione dell’equazione differenziale: dy/dt = G ( y , t ) con la condizione iniziale y(to)=yo. G è analitica nel punto (yo,to). Vale lo sviluppo y(t)= (1/i!) y(i) (to) (t-to)i Se [t-to] è più grande del raggio di convergenza della serie , l’errore di troncamento non tende verso zero al crescere del numero dei termini. Ma si può suddividere l’intervallo [to to+T] in una catena di intervalli parziali definiti dai successivi punti tj . Supponendo ora di avere già trovato nel punto tj la y e tutte le derivate necessarie, posso scrivere l’approssimazione di Taylor per t compreso tra tj e tj+1 : y(t)  zj(t) = (1/i!)y(i)(t-to)i .

34. Formule di integrazione del 2do ordine Nel caso si intenda arrestare lo sviluppo di Taylor al termine di grado 1, si ottiene yj+1 =zj(tj+1) = yj + h G(yj,tj). È la formula che descrive il “metodo di Eulero”, impiegato per la dimostrazione costruttiva della esistenza ed unicità della soluzione y(t) al problema di Cauchy proposto. Considero l’uguaglianza y(t+h)=y(t)+G(y,t+) d  Se l’integrale è sostituito dal prodotto h·G(y(t),t), l’errore è O(h2) .Una approssimazione dell’ordine di O(h3) si ottiene valutando l’integrale con la formula dei Trapezi: ½ h (G(y(t),t) + G(y(t+h),t+h) Ma non è ancora noto il valore di y al tempo t+h.Ne conosco però una approssimazione: y(t+h)= y(t)+h G(y,t) +O(h2) L’integrale può essere ottenuto come prodotto di h per G(y,t) con y valutato in un opportuno punto interno all’intervallo h di integrazione.

35. Errore O(h2) ed O(h3) • Una prima approssimazione di y è: y* = yi + h G(yi,ti) = yi+1 + O(h2)tale approssimazione consente di valutare :yi+1= yi + ½ h (G(yi,ti)+G(y*,ti+1)). Che errore si commette utilizzando y* anziché yi+1 ? • L’errore è espresso da:½ h (G(yi+1,ti)-G(y*,ti+1))= ½ hGy(yM,ti)(y*-yi+1))dove la derivata parziale di G va valutata in un opportuno yM , intermedio tra y* e yi+1.. L’errore quindi è O(h3).

36. Altre formule di 2do ordine • Un’altra formula al 2do ordine è:y* = yi + ½ h G(yi,ti)yi+1 =yi + h(Gy*,ti + ½ h) • Nei casi in cui le derivate terze sono piccole, una formula del secondo ordine con errore minimo è data da: • k1=h G(y,t) ; k2=h G(y+2/3 k1,t+2/3 h) ;Δy= ¼ ( k1+3 k2)

37. Formule del 4to ordine • È molto diffusa la formula del 4to ordine:k1=hG(y,t); k2=hG(y+ ½ k1,t+ ½ h);k3=h(G(y+ ½ k2, t+ ½ h); k4=h(G(y+k3, t+h);yi+1= yi + 1/6 (k1+2k2+2k3+k4)

38. Formule parametriche (complementi) • Tale formula del 4to ordine è dedotta per s=1 da una formula nel parametro libero s:k1=hG(y,t); • k2=hG(y+ ½ k1,t+ ½ h);k3=h(G(y+ ½k1+1/(2s)(k2-k1), t+ ½ h); k4=h(G(y+k2 + s (k3-k2), t+h);yi+1= yi + 1/6 (k1+4k2+2s(k3-k2)+k4)

39. Altre formule (complementi) La formula Runge – Kutta – Gill che suggeriamo impiega il valore s=(0.5) e la vedremo esplicitamente più avanti. Un’altra formula al 4to ordine è:k1=hG(y,t);k2=hG(y+ 1/3 k1,t+ 1/3 h);k3=h(G(y+k2-1/3 k1, t+2/3 h); k4=h(G(y+k3+k1-k2, t+h);yi+1= yi + 1/8 (k1+3k2+3k3+k4).

40. Metodi di Runge – Kutta di ordine > 4 (complementi) Negli esempi visti si sono costruiti i vettori ki di incrementi attribuibili ad y, per valutare diverse determinazioni di G(y,t). I metodi visti si generalizzano proponendo i vettor α e p e la matrice β costanti, di dimensione q e q x q-1, mediante i quali si valuterà il vettore k degli incrementi da assegnare ad y, per calcolare G in punti diversi del dominio, in modo simile a come si sono calcolate le k1 e k2 nei metodi del secondo ordine. Posto qm=q-1 k1= h G(y,t); k2= h G(y+b21 k1, t+a2 h); k3= h G(y+b31 k1 + b32 k2 , t+a3 h ); kq= h G(y+bq1 k1+…+ b q qm k qm , t+aq h). y(t+h)  z(h) = y(t)+ i pi ki

41. Metodi di Runge – Kutta espliciti (complementi) La matrice beta sarà triangolare inferiore, giacché si valutano i successivi elementi di k solo mediante gli elementi già calcolati. (è questo il significato dell`aggettivo “espliciti “). Le componenti di alfa, p e beta vanno scelte in modo da realizzare la precisione prevista. A tal fine si definisce le scarto E(h) tra la soluzione cercata y(t+h) ed il valore fornito dall’algoritmo, z(h). Tale E può essere sviluppata in serie di Taylor e, per metodi di ordine s, E deve annullarsi ad h=0 con le sue prime s derivate Per metodi del primo ordine (q=1) si ha z(h)=y + p1 h G(y,t); E(h) = y(t+h) - z(h) = y(t+h)-y(t)- p1h G(y,t). Si può scrivere la serie di Taylor per E in funzione di h e si vede che p1 deve valere 1.

42. Il piano delle fasiLe immagini dei “ritratti in Fase” sono tratte dagli appunti delle lezioni di G.Benettin per gli allievi ingegneri dell’Informazione • Uno strumento prezioso per lo studio di problemi della Fisica Matematica è l’uso del “piano delle fasi”, cioè la rappresentazione delle curve ‘posizione-velocità’ nel piano cartesiano x-y. Nel caso dell’oscillatore armonico tali curve sono ellissi concentriche, o, con scelta opportuna della scala, cerchi, con il raggio determinato dal valore dell’invariante x2+v2 . Le condizioni iniziali con x0 =0 sono rappresentate da punti sulla retta verticale, quelle con v0 = 0 sulla retta orizzontale.

43. Alcuni ritratti in fase

44. Verso di percorrenza dell’orbita Nei sistemi meccanici il verso con cui il punto rappresen-tativo percorre le orbite è orario: un punto in x=0 con v>0 parte dal ramo positivo dell’asse verticale ed invade il primo quadrante man mano che x aumenta (perché v è >0) e v diminuisce (allorché dv/dt è ~-x, quindi < 0). Poi quando v=0, x raggiunge il massimo positivo ed il punto entra nel quarto quadrante, con x >0, decrescente e v che diviene <0. Nella immagine sono illustrate le ellissi dell’oscillatore armonico, così percorse, le iperboli del repulsore armonico, per le quali vale la medesima osservazione, le rette della particella libera, percorse da sinistra a destra per v >0 e da destra verso sinistra per v<0. L’asse delle x (v=0) non è una traiettoria ma un insieme di punti di equilibrio isolati.

45. Moto per gravità • L’equazione è d2x/dt2 = -mg (con g = 9.81) • Le orbite sono parabole ad asse orizzontale. • Presentare come esercizio un esempio di moto per gravità del punto libero, con condizioni iniziali diverse ( p. es. velocità iniziali positive di valori diversi e posizioni iniziali positive con v=0 )

46. Equazioni differenziali ordinarie non lineari • Pendolo, crescita con risorse limitate, moto kepleriano, bipendolo • Equazioni con ciclo limite: Van der Pol, Lotka Volterra • Equazioni con esiti caotici: pendolo forzato, mappa standard

47. Moto del pendolo • Si distinguono tre tipi di moto assai diversi: • Le librazioni, od oscillazioni, rappresentate da ellissi, (con centro in 0,0 ed in 2nπ,0, n intero positivo o negativo), ellissi che assomigliano sempre più a cerchi per oscillazioni di piccola ampiezza. • Una curva limite, rappresentata dal punto isolato v=0, x=±π, in equilibrio instabile, e da due traiettorie che partono e terminano quanto si vuole vicino al punto isolato, passano per x=0,v=±2 e vengono percorse in tempo infinito!!! • Le rotazioni, orarie p.es. passanti per x=0, v >2 ed antiorarie, p.es. passanti per x=0, v <-2

48. Ritratto in fase: oscillatore armonico,buca di potenziale

49. Buca con colle, doppia buca

50. Ritratto in fase per il pendolo,relazione con il potenziale