Lo studio della stabilità dei sistemi meccanici analizza il comportamento del moto nelle vicinanze delle [[Stabilità dell'equilibrio|configurazioni di equilibrio]]. Attraverso la linearizzazione delle equazioni del moto, è possibile determinare non solo se un equilibrio è stabile o instabile, ma anche caratterizzare le frequenze di oscillazione e i modi normali del sistema. ### Sistemi con un grado di libertà Per un sistema olonomo a vincoli fissi con un solo grado di libertà (coordinata $q$) come l'[[Oscillatore meccanico]], le configurazioni di equilibrio $q^\circ$ corrispondono ai punti stazionari del potenziale: $U'(q^\circ) = 0$. Per studiare i moti "piccoli" attorno all'equilibrio, si introduce una perturbazione $\eta(t)$ tale che $q(t) = q^\circ + \varepsilon \eta(t)$. Sostituendo questa espressione nell'equazione di Lagrange e sviluppando in serie di Taylor fino al primo ordine in $\varepsilon$, si ottiene l'**equazione di moto linearizzata**: $ \ddot{\eta}(t) - \frac{U''(q^\circ)}{a(q^\circ)} \eta(t) = 0 $ dove $a(q^\circ) > 0$ è il coefficiente della matrice di massa valutato all'equilibrio. Il comportamento del sistema dipende esclusivamente dal segno della derivata seconda del potenziale $U''(q^\circ)$: 1. **$U''(q^\circ) < 0$ (Massimo del potenziale)**: L'equazione descrive un moto armonico. L'equilibrio è **stabile** e il sistema oscilla con una **frequenza delle piccole oscillazioni** pari a $\omega = \sqrt{-\frac{U''(q^\circ)}{a(q^\circ)}}$. 2. **$U''(q^\circ) > 0$ (Minimo del potenziale)**: L'equazione ammette soluzioni esponenziali (iperboliche). L'equilibrio è **instabile**. 3. **$U''(q^\circ) = 0$ (Flesso o caso degenere)**: L'approssimazione lineare prevede un moto uniforme, ma non è sufficiente per concludere. È necessario analizzare le derivate di ordine superiore. Tuttavia, si dimostra che se esiste un intorno in cui $U(q) \geq U(q^\circ)$, la posizione è instabile. ### Sistemi con più gradi di libertà Per sistemi con $N$ gradi di libertà, la linearizzazione porta a un sistema di equazioni differenziali accoppiate. Sviluppando l'energia cinetica e il potenziale attorno alla configurazione di equilibrio $\mathbf{q}^\circ$, si ottiene l'equazione matriciale: $ \mathbf{A}^\circ \ddot{\boldsymbol{\eta}} = \mathbf{B}^\circ \boldsymbol{\eta} $ dove $\mathbf{A}^\circ$ è la [[Matrice di massa]] (definita positiva) e $\mathbf{B}^\circ$ è la matrice Hessiana del potenziale, entrambe valutate in $\mathbf{q}^\circ$. Cercando soluzioni del tipo $\boldsymbol{\eta}(t) = \boldsymbol{\eta}_0 e^{\lambda t}$, il problema si riduce alla ricerca degli autovalori $\lambda^2$ e degli autovettori $\boldsymbol{\eta}_0$ generalizzati: $ \det(\mathbf{B}^\circ - \lambda^2 \mathbf{A}^\circ) = 0 $ Le radici $\lambda^2$ determinano la natura dei **modi normali** del sistema: - **$\lambda^2 < 0$**: Modo normale **oscillatorio**. Tutte le coordinate oscillano in fase con la stessa autofrequenza $\omega = \sqrt{-\lambda^2}$. - **$\lambda^2 > 0$**: Modo normale **iperbolico**. Il moto diverge esponenzialmente, indicando instabilità. - **$\lambda^2 = 0$**: Modo normale **lineare**. Indica un caso degenere che richiede analisi non lineari. #### Stabilità e contributo dei modi normali Se *tutte* le radici $\lambda^2$ sono negative, l'equilibrio è stabile e il moto generale è una sovrapposizione di oscillazioni armoniche pure (i modi normali). Se esiste *almeno una* radice positiva, l'equilibrio è instabile. Gli autovettori $\boldsymbol{\eta}_0$ formano una base ortogonale rispetto al prodotto scalare definito dalla matrice di massa ($\mathbf{u} \cdot \mathbf{A}^\circ \mathbf{v}$). Scegliendo opportunamente le condizioni iniziali (posizioni e velocità proporzionali a un singolo autovettore), è possibile "eccitare" un solo modo normale alla volta, facendo sì che il sistema oscilli con una singola frequenza pura. ### Esempi ed esercizi Immagina una corda di chitarra. Quando la pizzichi, non vibra in modo casuale. Il suo movimento complesso è in realtà la somma di movimenti molto semplici chiamati "modi normali". - **Il primo modo (fondamentale)**: Tutta la corda va su e giù insieme, come una corda per saltare. Ha una frequenza specifica ($\omega_1$). - **Il secondo modo (armonica)**: La metà sinistra va su mentre la metà destra va giù. Il centro sta fermo. Ha una frequenza più alta ($\omega_2$). La matematica della linearizzazione fa esattamente questo per qualsiasi sistema meccanico complesso (come un'auto con le sue sospensioni): scompone il movimento caotico in una somma di "vibrazioni fondamentali" indipendenti, ognuna con la sua frequenza. Se anche solo una di queste "vibrazioni" tende a crescere all'infinito invece di oscillare, l'intero sistema è instabile. ##### Domande di teoria - Perché nello sviluppo in serie di Taylor dell'equazione di Lagrange, la matrice di massa viene valutata solo all'ordine zero (costante), mentre il potenziale viene sviluppato fino al secondo ordine? - Qual è il significato fisico di un autovettore $\boldsymbol{\eta}_0$ associato a un autovalore $\lambda^2 < 0$? - Se l'Hessiano del potenziale ha un autovalore nullo, perché l'approssimazione lineare fallisce nel determinare la stabilità? ##### Esercizi - Esempi ed esercizi | Lezione 59 + 64 Turzi - Esempi capitoli 13.6-13.7 | Biscari - - **Esercizio 1**: Un punto materiale di massa $m$ è vincolato a muoversi sull'asse $x$ sotto l'azione del potenziale $U(x) = -kx^2 + \alpha x^4$ (con $k, \alpha > 0$). Trovare le posizioni di equilibrio stabile e calcolare la frequenza delle piccole oscillazioni attorno ad esse. - **Esercizio 2**: Due masse $m$ sono collegate tra loro e a due pareti fisse da tre molle identiche di costante elastica $k$, potendosi muovere solo lungo una retta orizzontale. Scrivere le matrici $\mathbf{A}^\circ$ e $\mathbf{B}^\circ$, calcolare le autofrequenze del sistema e descrivere fisicamente i due modi normali di oscillazione. ### Collegamenti --- > [!info]- Risorse > ![[!Meccanica Razionale#Risorse#Bibliografia]] > ![[!Meccanica Razionale#Risorse#Approfondimenti]]