Meccanica analitica/La lagrangiana e le equazioni di Eulero-Lagrange

Sulla scia del modulo precedente, cerchiamo di dare una forma dipendente dalle variabili lagrangiane all'equazione di Newton. Per prima cosa, diciamo che questo modulo è molto derivante, cioè ci saranno derivate, derivate, derivate. Poi, iniziamo subito a ricordare qualche elemento da tenere in considerazione:

$\begin{array}{cc}\text{eq. simbolica di d'Alembert}& \sum _{i=1}^N(m_i\ddot{x_i}-\overrightarrow{f_i})\delta O{P}_i=\sum _{i=1}^N\overrightarrow{r_i}\cdot \delta O{P}_i=0\\ \text{velocita'}& v_i=\sum _{h=1}^n\frac{\partial O{P}_i}{\partial q_h}\dot{q_h}+\frac{\partial O{P}_i}{\partial t}\\ \text{energia cinetica}& T=\frac{1}{2}\sum _{h,k=1}^n{a}_{hk}\dot{q_h}\dot{q_k}+\sum _{h=1}^n{a}_{0h}\dot{q_h}+\frac{1}{2}{a}_{00}\end{array}$

Premessa immediata: se i vincoli sono indipendenti dal tempo, i valori $a_{00},a_{0h}$ sono nulli, e quindi l'energia cinetica ha un'espressione semplice. Consideriamo proprio questo caso: un sistema i cui vincoli siano indipendenti dal tempo.

Deriviamo rispetto a $q_h$ e alle $\dot{q_h}$ velocità e energia cinetica:

$\begin{array}{cc}& \frac{\partial v_i}{\partial \dot{q_h}}=\frac{\partial O{P}_i}{\partial q_h}\\ & \frac{\partial v_i}{\partial q_h}=\frac{d}{dt}\frac{\partial O{P}_i}{\partial q_h}\\ & \frac{\partial T}{\partial \dot{q_h}}=\frac{1}{2}\sum _{i=1}^N{m}_{i}2v_i\frac{\partial v_i}{\partial \dot{q_h}}=\sum _{i=1}^N{m}_i{v}_i\cdot \frac{\partial O{P}_i}{\partial q_h}\\ & \frac{\partial T}{\partial q_h}=\frac{1}{2}\sum _{i=1}^N{m}_{i}2v_i\frac{\partial v_i}{\partial q_h}=\sum _{i=1}^N{m}_i{v}_i\cdot \frac{d}{dt}\left(\frac{\partial O{P}_i}{\partial q_h}\right)\end{array}$

Nell'ultima espressione, non abbiamo scritto per esteso il valore di $\frac{\partial v_i}{\partial q_h}$ perché risulta più comodo vederla così, in vista del prossimo passo che stiamo per compiere, ovvero derivare totalmente rispetto al tempo l'ultima espressione:

$\frac{d}{dt}\frac{\partial T}{\partial \dot{q_h}}=\sum _{i=1}^N(m_i\ddot{x_i}\cdot \frac{\partial O{P}_i}{\partial q_h}+\underset{\frac{\partial T}{\partial q_h}}{\underbrace{m_i{v}_i\frac{d}{dt}\frac{\partial O{P}_i}{\partial q_h}}})=\sum _{i=1}^N{m}_i\ddot{x_i}\frac{\partial O{P}_i}{\partial q_h}+\frac{\partial T}{\partial q_h}$

Da cui si ottiene che:

$\frac{d}{dt}\frac{\partial T}{\partial \dot{q_h}}-\frac{\partial T}{\partial q_h}=\sum _{i=1}^N{m}_i\ddot{x_i}\frac{\partial O{P}_i}{\partial q_h}$

Riprendiamo ora l'equazione simbolica di d'Alembert e imponiamo al problema di soddisfarla, ovvero di compiere il suo moto naturale; in questa, sostituiamo l'espressione di $\delta O{P}_i$ ottenendo:

$\sum _{i=1}^N(m_i\ddot{x_i}-f_i)\cdot \sum _{h=1}^n\frac{\partial O{P}_i}{\partial q_h}\delta q_h=0$

Poiché il sistema soddisfa l'equazione simbolica di d'Alembert, i vincoli saranno perfetti bilaterali e il lavoro virtuale è quindi nullo. È lecito quindi scrivere l'espressione di sopra $-\delta L=0$:

$\sum _{h=1}^n[\sum _{i=1}^N\left(m_i\ddot{x_i}\frac{\partial O{P}_i}{\partial q_h}\delta q_h\right)-\sum _{i=1}^N{f}_i\frac{\partial O{P}_i}{\partial q_h}\delta q_h]=0$

Ricordando l'espressione scritta poco fa, e ricordando l'espressione del lavoro virtuale che abbiamo dato nel precedente modulo, otteniamo:

$\sum _{h=1}^n[(\frac{d}{dt}\frac{\partial T}{\partial \dot{q_h}}-\frac{\partial T}{\partial q_h})\delta q_h-Q_h\delta q_h]=0$

Adesso, ricordiamo che, per la loro definizione, i $\delta q_h$ sono valori completamente indipendenti tra loro, e possiamo quindi semplificarli nella precedente espressione, ottenendo quelle che vengono chiamate equazioni di Lagrange:

$\frac{d}{dt}\frac{\partial T}{\partial \dot{q_h}}-\frac{\partial T}{\partial q_h}=Q_h$

Queste sono $n$ equazioni differenziali. Il caso interessante è quello in cui le componenti lagrangiane derivano da forze derivanti da potenziali; per queste, è possibile scrivere:

$Q_h=-\frac{\partial U}{\partial q_h}$

Ed è quindi possibile sostituire questa espressione nelle equazioni di Lagrange:

$\frac{d}{dt}\frac{\partial T}{\partial \dot{q_h}}-(\frac{\partial T}{\partial q_h}-\frac{\partial U}{\partial q_h})=0$

Ora ragioniamo sul fatto che, se i potenziali $U(q_1,\cdots ,q_n)$ sono funzioni solo delle $q_h$, la loro derivata rispetto a queste è nulla, poiché le $q_h$ sono indipendenti tra loro. Se questo accade, allora è lecito scrivere:

$\frac{\partial T}{\partial q_h}=\frac{\partial (T-U)}{\partial q_h}$

Dove ovviamente non ci interessa minimamente se si deriva in $q_h$ o in $\dot{q_h}$: è la stessa cosa, è lecito poter scrivere $T-U$ al posto di $T$. Allora, a questo punto, introduciamo la funzione.

$L(q_h,\dot{q_h};t)=T-U$

Che chiameremo lagrangiana. Ripetendo il ragionamento compiuto finora, otteniamo l'espressione:

$\frac{d}{dt}\frac{\partial L}{\partial \dot{q_h}}-\frac{\partial L}{\partial q_h}=0$

Queste vengono chiamate equazioni di Eulero-Lagrange; come per le equazioni di Lagrange, sono <math>n<math> equazioni differenziali. La soluzione di questo sistema differenziale fornisce le equazioni del moto del sistema. Adesso, ragionando al di là del percorso teorico affrontato per poter raggiungere questo strumento, è immediato comprendere come, utilizzando queste equazioni, si possa facilmente e velocemente risolvere un problema di meccanica in cui è presente un sistema vincolato. Attraverso l'opportuna scelta delle variabili lagrangiane, che ricordiamo essere arbitrarie, risulta poi facile determinare le equazioni del moto e, quindi, come evolve il sistema nel tempo. Basta semplicemente scriversi la funzione lagrangiana, derivarla e il lavoro è fatto. Template:Avanzamento