Effetti non lineari

Da "Fisica, onde Musica": un sito web su fisica delle onde e del suono, acustica degli strumenti musicali, scale musicali, armonia e musica.

In quasi tutti gli esempi di oscillazioni ed onde trattati in questo sito ci siamo limitati ad analizzare i fenomeni in regime lineare, anche se alcuni effetti non lineari, o anarmonici sono stati quasi sempre citati tra gli "approfondimenti" delle varie sezioni.

Il motivo principale per cui non abbiamo affrontato una vera trattazione sistematica, consiste nel fatto che gli effetti della non linearità sono molteplici e diversissimi, e la fisica dei sistemi non lineari include fenomeni completamente nuovi, rispetto al semplice moto ondulatorio. Ne sono esempi:

  • l'emergenza di dinamiche caotiche, o di attrattori strani
  • la presenza di nuove eccitazioni elementari del sistema che non sono onde
  • l'interazione tra onde (non vale più il principio di sovrapposizione)
  • la modificazione del comportamento dei modi normali del sistema: essi non sono più indipendenti tra di loro, ma sono accoppiati, e quindi non corrispondono più ad onde stazionarie.
  • la modificazione delle risonanze e delle frequenze proprie del sistema
  • la presenza di fenomeni di isteresi o saturazione

ecc.

Nel seguito cerchiamo di chiarire i termini e portare un solo, ma importante esempio.

Altri esempi, più complessi, ma più interessanti, si trovano nelle sezioni catena di Fermi-Pasta-Ulam e Instabilità modulazionale.

Regime lineare e non lineare

Chiariamo i termini

  • Un singolo oscillatore si trova in regime lineare quando è un oscillatore armonico, e cioè quando la forza di richiamo è direttamente proporzionale alla distanza dell'oscillatore da un punto, che ne costituisce il punto di equilibrio
𝐅=−k(𝐱−𝐱0).

Questa forza quindi ammette un potenziale direttamente proporzionale al quadrato della distanza dal punto di equilibrio

V=12k(𝐱−𝐱0)2.
  • Un sistema di N oscillatori accoppiati si trova in regime lineare quando ciascun oscillatore è soggetto a un insieme di forze armoniche. Per esempio, se ogni oscillatore è legato alla propria posizione di equilibrio da una molla ideale e gli oscillatori sono collegati tra loro da altre molle ideali, l'i-esimo oscillatore sarà soggetto ad una forza complessiva pari a
𝐅i=−kii(𝐱i−𝐱i0)−∑j≠ikij[(𝐱i−𝐱i0)−(𝐱j−𝐱j0)]

dove 𝐱i0 è la posizione di equilibrio dell'oscillatore i. Le forze dipendono quindi linearmente dagli spostamenti di tutti gli oscillatori. Il potenziale dell'intero sistema è

V=12∑ikii(𝐱i−𝐱i0)2+12∑i<jkij[(𝐱i−𝐱i0)−(𝐱j−𝐱j0)]2

Le costanti kij formano una matrice. Gli elementi diagonali kii descrivono la rigidità della molla propria di ciascun oscillatore, che lo richiama verso la sua posizione di equilibrio; gli elementi non diagonali kij=kji (con i≠j) descrivono l'accoppiamento reciproco tra gli oscillatori i e j. Se gli elementi non diagonali sono nulli, gli oscillatori diventano indipendenti l'uno dall'altro: ciascuno oscilla per conto proprio, con la frequenza fissata dalla propria molla e dalla propria massa.

  • Regime non lineare si riferisce a tutte le situazioni che si discostano dal caso precedente.

Cosa accade nella realtà

  • Dal punto di vista della realtà sperimentale, e senza ulteriori ipotesi, i sistemi fisici reali sono da considerarsi tutti non lineari.
    • Ad esempio, le forze elastiche reali non seguono la legge di Hooke per qualunque dislocazione dalla posizione di equilibrio. Ce ne si rende conto facilmente notando che un sistema elastico, al crescere della deformazione, prima risponde elasticamente, ma poi si snerva, rimane deformato in modo perenne, e poi si rompe, quando viene sollecitato oltre un certo limite. Tutti fenomeni evidentemente non lineari.
    • Dal punto di vista fisico-matematico la presenza di termini non lineari si traduce nel fatto che le forze non sono più semplicemente direttamente proporzionali agli spostamenti, ma dipendono anche dal loro quadrato, cubo, ecc. Un esempio tipico è costituito da un oscillatore il cui potenziale è
V=12k1(𝐱−𝐱0)2+13k2(𝐱−𝐱0)3.
  • Tuttavia, la maggior parte dei sistemi fisici possono essere considerati approssimativamente lineari quando le oscillazioni sono abbastanza piccole.
    • Questa proprietà è equivalente ad assumere che le forze (o i potenziali) siano tutte sviluppabili in serie di potenze con un termine al primo ordine. In parole povere è come dire "la terra è piatta a tutti gli effetti, finché ne vediamo una porzione sufficientemente piccola".

Il pendolo

È straordinariamente facile ottenere comportamenti complessi anche in sistemi apparentemente molto semplici. Il comune pendolo ne è un esempio.

L'equazione del moto del pendolo determina la legge oraria del moto dell'angolo di oscillazione. Tuttavia la forza netta di richiamo (dovuta alla composizione della forza di gravità e della tensione della corda) dipende dal seno dell'angolo di oscillazione. Perciò l'equazione esatta del moto del pendolo (ferme restando le altre approssimazioni: corda ideale senza massa, massa oscillante puntiforme, moto in un piano) è pertanto

α+glsin⁡θ=0, (1)

dove α indica l'accelerazione angolare.

La forma comune, che riduce il pendolo ad un oscillatore armonico, invece, è in realtà la forma approssimata al prim'ordine nello sviluppo in serie della funzione seno

α+glθ=0. (2)

Utilizzando la (2) anziché la (1) si ottiene la nota espressione per il periodo del pendolo

T=2πlg, (3)

Questa legge del moto è in accordo con l'osservazione di Galileo che il periodo non dipende dall'ampiezza dell'oscillazione, ma solo l'equazione (1) descrive correttamente anche le oscillazioni di grande ampiezza. Per queste oscillazioni il periodo è dato da

T=2πlg(1+14sin2θmax2+964sin4θmax2+⋯). (4)

Dall'esame delle seguenti figure si nota che le oscillazioni di maggiore ampiezza hanno periodo maggiore rispetto a quelle di minore ampiezza. Inoltre la loro forma non è più sinusoidale.

Oscillazioni di 15° e 150° a confronto. Il periodo dell'oscillazione più ampia è molto maggiore di quello dell'oscillazione più piccola. Diagramma velocità angolare-angolo (spazio delle fasi) per un'oscillazione di 15°. L'oscillazione sinusoidale genera un'ellisse (figura di Lissajous). Diagramma velocità angolare-angolo (spazio delle fasi) per un'oscillazione di 150°. L'oscillazione non è sinusoidale, e non genera pertanto una figura di Lissajous.

L'oscillatore di Van der Pol

Un tipo di oscillatore non lineare modello che trova diverse applicazioni pratiche è l'oscillatore di Van der Pol. Esso soddisfa l'equazione

ω˙+θ+μ(θ2−1)ω=0

dove ω=θ˙ e α=ω˙ indicano rispettivamente velocità ed accelerazione (angolari) dell'oscillatore.

L'equazione fu introdotta da Balthasar van der Pol nel 1926 per descrivere i circuiti oscillanti a valvole termoioniche[1]. È il prototipo degli oscillatori autosostenuti, che continuano a oscillare grazie a una sorgente di energia costante, con un'ampiezza fissata dalla non linearità e non dalle condizioni iniziali. Si comportano in modo simile, per esempio, la corda strofinata dall'arco e l'ancia degli strumenti a fiato, alimentate dal moto continuo dell'arco o dal fiato del suonatore.

I primi due termini descrivono un oscillatore armonico di periodo 2π, poiché la costante di proporzionalità tra accelerazione angolare ed angolo ha modulo uguale a 1; l'ultimo termine è responsabile del comportamento non lineare. Vediamo secondo quali modalità.

Analisi grafica

Per farci un'idea delle proprietà di questo oscillatore abbiamo risolto numericamente l'equazione di Van der Pol, e abbiamo rappresentato il moto in due modi nelle figure che seguono. Nella prima colonna si trova la solita legge del moto, in cui l'angolo di oscillazione è disegnato in funzione del tempo. Nella seconda colonna abbiamo invece rappresentato la velocità angolare in funzione dell'angolo. Questa rappresentazione (detta spazio delle fasi) ha diversi vantaggi rispetto a quella tradizionale. In particolare permette di accorgersi rapidamente se esistono dei cicli limite, o attrattori, verso i quali il moto del sistema è appunto attratto a partire da diverse condizioni iniziali.

Per avere un riferimento abbiamo tracciato la traiettoria dell'oscillatore di Van der Pol in nero, e quella dell'oscillatore armonico con le stesse condizioni iniziali in rosso.

Esaminiamo ora i tracciati. Un'ulteriore breve analisi si trova nel paragrafo seguente.

  • Per la condizione iniziale scelta si nota che l'ampiezza delle oscillazioni aumenta sempre gradualmente fino a raggiungere un ciclo limite
  • La forma del ciclo limite non dipende dalle condizioni iniziali, ma solo dal valore del parametro μ, e si discosta tanto più dal ciclo armonico (circonferenza) quanto maggiore è il valore di μ.


Legge del moto Spazio delle fasi
μ=0.1
μ=0.5
μ=1
μ=4

La seguente osservazione suggerisce che l'oscillatore sia robusto rispetto alle condizioni iniziali. Infatti, vediamo che un ciclo limite esiste sia che iniziamo il moto all'interno o all'esterno di esso. In altre parole, se le oscillazioni iniziali sono troppo piccole esse verranno amplificate, mentre se sono troppo grandi, verranno smorzate fino ad un valore stabile. Il ciclo limite raggiunto è lo stesso nei due casi.

μ=1.5 Legge del moto Spazio delle fasi
θ(0)=30∘;ω(0)=1
θ(0)=−80∘;ω(0)=2

Approfondimenti

Vediamo come possiamo giustificare le osservazioni precedenti a partire dall'analisi dell'equazione del moto. Ricordiamo che il solo termine non lineare è quello proporzionale a μ. Dall'esame dell'equazione

ω˙+θ+μ(θ2−1)ω=0

vediamo che

  • se μ=0 l'oscillatore si riduce evidentemente ad un oscillatore armonico. Il suo moto sarà sinusoidale puro, e il diagramma angolo-velocità seguirà una traiettoria circolare. Nelle immagini di sopra riportiamo in rosso queste linee come riferimento.
  • se μ>0 il termine aggiunto ha due effetti che si compensano:
    • per oscillazioni al di sopra di un certo angolo limite (in questo caso |θ|>1) il termine è positivo, e corrisponde ad una sottrazione di energia dall'oscillatore. In sostanza equivale ad un termine di attrito proporzionale alla velocità del moto. In sostanza, dice che le oscillazioni di grande ampiezza tendono ad essere smorzate
    • per oscillazioni al di sotto dell'angolo limite (in questo caso |θ|<1) il termine è negativo, e corrisponde ad un'aggiunta di energia all'oscillatore. In sostanza esso equivale ad una forza esterna proporzionale alla velocità del moto. L'aggiunta di energia significa che le oscillazioni di piccola ampiezza tendono ad essere aumentate.

Si intuisce che l'effetto netto del termine non lineare, qualunque siano le condizioni iniziali, consiste nello stabilizzare le oscillazioni su un ciclo limite, che costituisce il punto di equilibrio tra oscillazioni troppo grandi, che tendono a smorzarsi, ed oscillazioni piccole, che invece tendono ad aumentare di ampiezza. La forma del ciclo limite si discosterà tanto più dalla circonferenza (oscillatore armonico) quanto maggiore è il valore di μ.

Verifica la comprensione

Problema. Un pendolo lungo 1 m oscilla con un'ampiezza di 60°. Calcola il periodo con la formula (3) delle piccole oscillazioni e poi con i primi termini della (4). (g = 9,81 m/s²) Soluzione. Per piccole oscillazioni T=2πl/g ≈ 2,006 s. Con θmax = 60° si ha sin 30° = 0,5, quindi il fattore tra parentesi vale 1 + 1/4 · 0,25 + 9/64 · 0,0625 ≈ 1,071 e il periodo è circa 2,15 s, il 7% in più (il valore esatto è 2,153 s).

Problema. Un orologio a pendolo è regolato per segnare l'ora esatta quando il pendolo oscilla con un'ampiezza di 5°. Se l'ampiezza sale a 10°, quanto ritarda in un giorno? Soluzione. Con la (4), limitandosi al primo termine correttivo, il periodo cambia del fattore relativo 14(sin25∘−sin22,5∘) ≈ 0,0014, cioè dello 0,14%. In un giorno (86400 s) l'orologio ritarda di circa 0,0014 · 86400 s ≈ 120 s, cioè di circa 2 minuti: per questo negli orologi di precisione l'ampiezza deve restare costante.

Problema. Un oscillatore ha il potenziale V=12k1x2+13k2x3 (con x0=0), con k1 = 100 N/m e k2 = 1000 N/m². Calcola la forza di richiamo per x = +5 cm e per x = −5 cm. Soluzione. La forza è F=−dV/dx=−k1x−k2x2. Per x = +0,05 m: F = −5 N − 2,5 N = −7,5 N. Per x = −0,05 m: F = +5 N − 2,5 N = +2,5 N. La forza riporta sempre verso l'equilibrio, ma non è simmetrica: dalla parte delle x negative la molla è molto più "morbida" e l'oscillazione non è più sinusoidale.

Problema. Nell'oscillatore di Van der Pol con μ = 1, il termine non lineare fornisce o sottrae energia quando l'angolo vale 0,5 rad? E quando vale 2 rad? Soluzione. Il termine è μ(θ2−1)ω. Per θ = 0,5 il coefficiente di ω è 0,25 − 1 = −0,75: negativo, quindi il termine agisce come una forza che spinge nel verso del moto e fornisce energia. Per θ = 2 il coefficiente è 4 − 1 = 3: positivo, il termine agisce come un attrito e sottrae energia. Per questo le oscillazioni piccole crescono e quelle grandi si smorzano fino al ciclo limite.

Note

  1. ↑ B. van der Pol, "On relaxation-oscillations", Philosophical Magazine 2, 978-992 (1926), doi:10.1080/14786442608564127

"Fisica, onde Musica": un sito web su fisica delle onde, acustica degli strumenti musicali, scale musicali, armonia e musica.