Visualizza wikitesto di Calcolo telai piani col metodo delle rigidezze

{{C|il tono e le parole usate fanno pensare a un manuale più che a una voce di enciclopedia|ingegneria|febbraio 2014}}

Di seguito si mostra il '''calcolo di telai piani''' ossia si sviluppa una procedura che permette di determinare le caratteristiche di sollecitazione interne di un telaio piano. Per telaio si intende un insieme di travi collegate tra loro con vincoli interni.

In particolare si utilizza il [[metodo degli spostamenti]] che passa attraverso il calcolo della matrice di [[rigidezza]] della struttura. Questo metodo è anche definito [[Metodo delle rigidezze]].

== La matrice di rigidezza di una struttura ==
Su una struttura si evidenziano gli spostamenti ([[Traslazione (geometria)|traslazioni]] e [[Rotazione (matematica)|rotazioni]]) di alcuni punti, definibili attraverso valori scalari reali. 
[[File:Calcolo telai piani-Convenzioni spostamenti.png|center|600px]]
Negli stessi punti si applicano i carichi esterni omologhi a tali spostamenti (traslazioni associate a forze con la stessa direzioni, rotazioni a momenti con stesso piano di applicazione e verso di rotazione).
[[File:Calcolo telai piani-Convenzioni forze.png|center|600px]]
Riportando gli spostamenti nodali su un vettore verticale <math> \begin{Bmatrix} X \end{Bmatrix} </math> e gli sforzi nodali applicati sui medesimi nodi su un vettore verticale <math> \begin{Bmatrix} F \end{Bmatrix} </math>, la relazione che lega le due grandezze è del tipo:
<div style="text-align:center;"><math> \begin{Bmatrix} F \end{Bmatrix} =  \begin{bmatrix} H \end{bmatrix} \begin{Bmatrix} X \end{Bmatrix} </math>.</div>
in cui <math> \begin{bmatrix} H \end{bmatrix} </math> è una matrice associata alla struttura definita matrice di rigidezza della struttura. Per il [[teorema di reciprocità]] o di Betti, la matrice <math> \begin{bmatrix} H \end{bmatrix} </math> è simmetrica.

== Matrice di rigidezza trave ==
Per determinare la matrice di rigidezza dell'intera struttura si parte dal calcolo della matrice di rigidezza di una singola trave.

In tal caso, indicando con <math> \begin{Bmatrix} \eta \end{Bmatrix} </math> gli spostamenti dei punti estremi della trave e con <math> \begin{Bmatrix} S \end{Bmatrix} </math> i carichi applicati nei due punti, le due grandezze saranno legate da una relazione del tipo: 
<div style="text-align:center;"><math> \begin{Bmatrix} S \end{Bmatrix} =  \begin{bmatrix} K \end{bmatrix} \begin{Bmatrix} \eta \end{Bmatrix} </math>.</div>
in cui <math>  \begin{Bmatrix} K \end{Bmatrix} </math> è la matrice di rigidezza della trave che è del tipo:
<div style="text-align:center;"><math>
\begin{bmatrix}K\end{bmatrix} = 
\begin{bmatrix}
k_{1,1} & k_{1,2} & k_{1,3} & k_{1,4} & k_{1,5} & k_{1,6} \\
k_{2,1} & k_{2,2} & k_{2,3} & k_{2,4} & k_{2,5} & k_{2,6} \\
k_{3,1} & k_{3,2} & k_{3,3} & k_{3,4} & k_{3,5} & k_{3,6} \\
k_{4,1} & k_{4,2} & k_{4,3} & k_{4,4} & k_{4,5} & k_{4,6} \\
k_{5,1} & k_{5,2} & k_{5,3} & k_{5,4} & k_{5,5} & k_{5,6} \\
k_{6,1} & k_{6,2} & k_{6,3} & k_{6,4} & k_{6,5} & k_{6,6} \\
\end{bmatrix} 
</math></div>

=== Calcolo matrice di rigidezza locale singola trave ===
La relazione scritta sopra è di carattere generale e si applica a qualsiasi trave. Di seguito si calcola la matrice di rigidezza di una trave a sezione costante, costituita da un materiale omogeneo ed isotropo, seguendo le ipotesi del [[De Saint-Venant|de Saint Venant]]. 

In particolare si suppone che la sezione della trave sia costante e di area <math>A</math>. Il [[momento di inerzia]] calcolato lungo un asse perpendicolare al piano del telaio e passante per il baricentro della sezione è <math>J</math>. Inoltre tale asse è uno degli assi principali di inerzia della sezione, di modo che i momenti applicati nel piano del telaio diano luogo a flessione retta.

Per determinare i valori dei singoli coefficienti della matrice <math> \begin{Bmatrix} K \end{Bmatrix} </math> si considerano gli schemi statici.

Partendo dal calcolo dei coefficienti che riguardano i carichi nodali <math>S_1</math> e <math>S_4</math>, si considera una trave vincolata sul secondo estremo, libera sul primo e sottoposta su quest'ultimo estremo ad un carico assiale <math>S_1</math>. Supponendo di trascurare gli effetti del second'ordine (non linearità geometriche), si è in presenza di uno spostamento <math>\eta_{1} \ne 0 </math> e <math> \eta_{2...6} = 0</math>. La relazione che lega <math>S_1</math> a <math>\eta_1</math> è:
<div style="text-align:center;"> 
<math>S_{1} = E A \varepsilon = \frac {E A}{l} \eta_1</math> 
</div>
in cui <math>\varepsilon </math> è la deformazione unitaria dell'asta.

La reazione vincolare assiale sul secondo estremo è data da:
<div style="text-align:center;"> 
<math>S_{4} = -S_1 = -\frac {E A}{l} \eta_1</math>
</div>
Procedendo analogamente per <math>\eta_4</math>, si può ricavare:
<div style="text-align:center;"> 
<math>
\begin{matrix}
k_{1,1} = \frac{EA}{l} & k_{1,2} = 0 & k_{1,3} = 0 & k_{1,4} = -\frac{EA}{l} & k_{1,5} = 0 & k_{1,6} = 0 \\
k_{2,1} = 0 & k_{2,4} = 0 & k_{3,1} = 0 & k_{3,4} = 0 \\
k_{4,1} = -\frac{EA}{l} & k_{4,2} = 0 & k_{4,3} = 0 & k_{4,4} = \frac{EA}{l} & k_{4,5} = 0 & k_{4,6} = 0 \\
k_{5,1} = 0 & k_{5,4} = 0 & k_{6,1} = 0 & k_{6,4} = 0  \\
\end{matrix}
</math>
</div>

Si considera ora una trave semplicemente appoggiata, sottoposta nei due estremi a due momenti <math> S_3 </math> e <math> S_6 </math>.

Con alcuni calcoli si ottiene che la relazione che lega gli spostamenti nodali <math> \eta_3 </math> e <math> \eta_{6} </math> agli sforzi nodali <math> S_3 </math> e <math> S_6 \,\! </math> è:
<div style="text-align:center;"><math> \begin{Bmatrix} \eta_3 \\ \eta_6 \end{Bmatrix} =  \begin{bmatrix} \frac{l}{3EJ} & -\frac{l}{6EJ} \\ -\frac{l}{6EJ} & \frac{l}{3EJ} \end{bmatrix} \begin{Bmatrix} S_3 \\ S_6 \end{Bmatrix} </math></div>

Invertendo la relazione suddetta si ha:
<div style="text-align:center;"><math>
\begin{Bmatrix} S_3 \\ S_6 \end{Bmatrix} =  \begin{bmatrix} \frac{4EJ}{l} & \frac{2EJ}{l} \\ \frac{2EJ}{l} & \frac{4EJ}{l} \end{bmatrix} \begin{Bmatrix} \eta_3 \\ \eta_6 \end{Bmatrix}
</math> </div>
che permette di trovare <math> k_{3,3}</math> <math> k_{3,6}</math> <math> k_{6,3}</math> e <math> k_{6,6}</math>.

Calcolando le reazioni vincolari verticali sui due estremi, si trova: 
<div style="text-align:center;"><math>
\begin{matrix}
S_2 =  \frac{S_6+S_3}{l} =  \frac{1}{l} (\frac{6 E J}{l} \eta_3 + \frac{6 E J}{l} \eta_6 ) \\
S_5 = -\frac{S_6+S_3}{l} = -\frac{1}{l} (\frac{6 E J}{l} \eta_3 + \frac{6 E J}{l} \eta_6 ) \\
\end{matrix}
</math> </div>
da cui:
<div style="text-align:center;"><math>
\begin{matrix}
k_{2,3} = k_{2,6} =  \frac{6 E J}{l^2} \\
k_{5,3} = k_{5,6} = -\frac{6 E J}{l^2} \\
\end{matrix}
</math> </div>

Per trovare i restanti coefficienti si analizza una trave orizzontale incastrata sul primo estremo e vincolata sul secondo con una guida che consenta la sola traslazione verticale. Si applica sul secondo vertice una forza verticale <math>S_5 \,\! </math>. In tal caso lo spostamento verticale <math> \eta_5</math> del secondo estremo è dato dalla relazione:
<div style="text-align:center;"> 
<math>\eta_5 = \frac{S_5 l^3}{12 E J}</math>
</div>
da cui:
<div style="text-align:center;"> 
<math>k_{5,5} = \frac{12 E J }{l^3}</math>
</div>
Sempre facendo riferimento alla configurazione sopra descritta, la reazione vincolare verticale nel primo estremo, è pari a <math>-S_5</math>. Pertanto: 
<div style="text-align:center;"> 
<math>k_{2,5} = -\frac{12 E J }{l^3}</math>
</div>

Per simmetria si trova <math>k_{2,2}</math> e <math>k_{5,2}</math> pari rispettivamente a <math>\frac{12 E J }{l^3}</math> e <math>-\frac{12 E J }{l^3}</math>.

Ricapitolando quanto visto finora si ha la seguente matrice di rigidezza:
<div style="text-align:center;"><math>
\begin{bmatrix}K\end{bmatrix} = 
\begin{bmatrix}
\frac{EA}{l} & 0 & 0 & -\frac{EA}{l} & 0 & 0 \\
0 & \frac{12EJ}{l^3} & \frac{6EJ}{l^2} & 0 & -\frac{12EJ}{l^3} & \frac{6EJ}{l^2} \\
0 & \frac{6EJ}{l^2} & \frac{4EJ}{l} & 0 & -\frac{6EJ}{l^2} & \frac{2EJ}{l} \\
-\frac{EA}{l} & 0 & 0 & \frac{EA}{l} & 0 & 0 \\
0 & -\frac{12EJ}{l^3} & -\frac{6EJ}{l^2} & 0 & \frac{12EJ}{l^3} & -\frac{6EJ}{l^2} \\
0 & \frac{6EJ}{l^2} & \frac{2EJ}{l} & 0 & -\frac{6EJ}{l^2} & \frac{4EJ}{l} \\
\end{bmatrix} 
</math></div>

=== Sistema di riferimento globale ===
Quanto appena descritto si riferisce ad un [[sistema di riferimento]] locale della trave avente l'asse delle ascisse coincidente con l'asse della trave. Dovendo successivamente assemblare la matrice di rigidezza globale della struttura è però necessario coordinare tra loro i sistemi di riferimenti delle varie travi. Per fare questo si riportano i sistemi di riferimento delle singole travi ad un sistema di riferimento comune. Questo vorrà dire trasportare spostamenti e sollecitazioni dal sistema di riferimento locale della trave al sistema di riferimento globale.

Indicando con <math>\begin{Bmatrix}\eta_G\end{Bmatrix} </math> e <math>\begin{Bmatrix}S_G\end{Bmatrix} </math> rispettivamente spostamenti e sollecitazioni nodali nel sistema di riferimento globale, i due dovranno essere legati tra loro dalla relazione:
<div style="text-align:center;"><math>
\begin{Bmatrix}S_G\end{Bmatrix} = \begin{bmatrix}K_G\end{bmatrix} \begin{Bmatrix}\eta_G\end{Bmatrix}
</math></div>

.{{Chiarire|Nel seguito si vuole individuare i coefficienti di|incompleto}} 

Per fare questo si inizia caratterizzando la trave con un versore bidimensionale <math>\vec t</math> orientato come la stessa asta ed avente per definizione modulo unitario. Supponendo che gli estremi della trave siano i due punti <math>P_1 (x_1, y_1)</math> e <math>P_2 (x_2,y_2)</math>, le componenti di <math>\vec t</math> sono date da: 
<div style="text-align:center;"><math>
\begin{matrix}
t_x = \frac {x_2 - x_1}{\sqrt{ (x_2 - x_1)^2 + (y_2 - y_1)^2}} \\
t_y = \frac {y_2 - y_1}{\sqrt{ (x_2 - x_1)^2 + (y_2 - y_1)^2}} \\
\end{matrix} 
</math></div>

Si noti che assumendo pari ad <math>\alpha</math> l'inclinazione della trave rispetto all'asse delle ascisse, le componenti diventano:
<div style="text-align:center;"><math>
\begin{matrix}
t_x = \cos \alpha \\
t_y = \sin \alpha \\
\end{matrix} 
</math></div>

In tal caso, considerando che le rotazioni rimangono immutate con le variazioni di sistema di riferimento, il legame tra gli spostamenti nel sistema di riferimento locale <math>\begin{Bmatrix}\eta \end{Bmatrix} </math> e gli spostamenti nel sistema di riferimento globale <math>\begin{Bmatrix}\eta_G \end{Bmatrix} </math> è:
<div style="text-align:center;"><math>
\begin{Bmatrix} \eta_G \end{Bmatrix} = \begin{bmatrix} G \end{bmatrix} \begin{Bmatrix} \eta \end{Bmatrix}
</math></div>
in cui la matrice <math>\begin{bmatrix} G \end{bmatrix}</math> è data da:
<div style="text-align:center;"><math>
\begin{bmatrix} G \end{bmatrix} = 
\begin{bmatrix} 
t_x & -t_y & 0 & 0 & 0 & 0 \\
t_y & t_x & 0 & 0 & 0 & 0 \\
0 & 0 & 1 & 0 & 0 & 0 \\
0 & 0 & 0 & t_x & -t_y & 0 \\
0 & 0 & 0 & t_y & t_x & 0 \\
0 & 0 & 0 & 0 & 0 & 1 \\
\end{bmatrix} = 
\begin{bmatrix} 
\cos \alpha & -\sin \alpha & 0 & 0 & 0 & 0 \\
\sin \alpha & \cos \alpha & 0 & 0 & 0 & 0 \\
0 & 0 & 1 & 0 & 0 & 0 \\
0 & 0 & 0 & \cos \alpha & -\sin \alpha & 0 \\
0 & 0 & 0 & \sin \alpha & \cos \alpha & 0 \\
0 & 0 & 0 & 0 & 0 & 1 \\
\end{bmatrix}
</math></div>

<math>\begin{bmatrix} G \end{bmatrix} </math> è una [[matrice simmetrica]] con una particolarità: la sua inversa è uguale alla sua trasposta <math>(G^{-1} = G^{T})</math>.

Analogamente la relazione che lega le sollecitazioni nodali nel sistema locale <math>\begin{Bmatrix}S\end{Bmatrix}</math> a quelle nel sistema globale <math>\begin{Bmatrix}S_G\end{Bmatrix}</math> è:
<div style="text-align:center;"><math>
\begin{Bmatrix} S_G \end{Bmatrix} = \begin{bmatrix} G \end{bmatrix} \begin{Bmatrix} S \end{Bmatrix}
</math></div>

Invertendo le suddette relazioni, si ottiene:
<div style="text-align:center;"> <math>
\begin{matrix}
\begin{Bmatrix} \eta \end{Bmatrix} = \begin{bmatrix} G^T \end{bmatrix} \begin{Bmatrix} \eta_G \end{Bmatrix}\\
\begin{Bmatrix} S \end{Bmatrix} = \begin{bmatrix} G^T \end{bmatrix} \begin{Bmatrix} S_G \end{Bmatrix}\\
\end{matrix}
</math></div>
che sostituite in:
<div style="text-align:center;"> <math>
\begin{Bmatrix} S \end{Bmatrix} = \begin{bmatrix} K \end{bmatrix} \begin{Bmatrix} \eta \end{Bmatrix}
</math></div>
forniscono:
<div style="text-align:center;"> <math>
\begin{bmatrix} G^T \end{bmatrix} \begin{Bmatrix} S_G \end{Bmatrix} =
\begin{bmatrix} K \end{bmatrix}
\begin{bmatrix} G^T \end{bmatrix} \begin{Bmatrix} \eta_G \end{Bmatrix}
</math></div>
da cui si ottiene la matrice di rigidezza <math> \begin{bmatrix} K_G \end{bmatrix} </math> cercata:
<div style="text-align:center;"> <math>
\begin{bmatrix} K_G \end{bmatrix} =
\begin{bmatrix} G \end{bmatrix} \begin{bmatrix} K \end{bmatrix} \begin{bmatrix} G^T \end{bmatrix}
</math></div>

== Assemblaggio matrice di rigidezza globale ==

Nei paragrafi precedenti si è analizzato le travi singolarmente, ora l'intero telaio. Il telaio si presenta come un insieme di travi che interagiscono tra loro e con l'ambiente esterno. All'interno del formalismo matematico tale interazione si manifesta attraverso i vincoli interni (interazione tra le travi) e quelli esterni (interazione con l'ambiente esterno).

Analiticamente si dovrà passare dalla matrice di rigidezza delle singole travi alla matrice di rigidezza del telaio.

Si suppone che il telaio sia composto da <math>n</math> travi, si uniscono i vettori spostamenti nodali delle singole travi ottenendo un vettore <math>\begin{Bmatrix} \Eta \end{Bmatrix} </math> degli spostamenti nodali del telaio disarticolato così definito:
<div style="text-align:center;"> <math>
\begin{Bmatrix} \Eta \end{Bmatrix} = 
\begin{Bmatrix}
  \begin{Bmatrix} \eta_G \end{Bmatrix}_1 \\
  \begin{Bmatrix} \eta_G \end{Bmatrix}_2 \\
     ... \\
  \begin{Bmatrix} \eta_G \end{Bmatrix}_n \\
\end{Bmatrix}
</math></div>

Lo scopo è passare dal vettore <math>\Eta </math> (spostamenti nodali di ogni singola trave) al vettore <math> X </math> degli spostamenti nodali del telaio. Perciò si ricerca un operatore <math>\begin{bmatrix} I \end{bmatrix}</math> che leghi le due grandezze secondo:
<div style="text-align:center;"> <math>
\begin{Bmatrix} \Eta \end{Bmatrix} = \begin{bmatrix} I \end{bmatrix} \begin{Bmatrix} X \end{Bmatrix}
</math></div>
<math>\begin{bmatrix} I \end{bmatrix}</math> è definita matrice di incidenza ed è una matrice <math>m \times 6n</math>, in cui <math>m</math> è il numero degli spostamenti nodali del telaio, <math>6</math> è il numero complessivo di spostamenti nodali per ogni trave (4 traslazioni e 2 rotazioni) e <math>n</math> è il numero di aste. Il generico elemento <math>i_{i, 6(j-1)+k} \quad (i=1...m, j=1... n, k=1...6) </math> sarà pari a 1 se allo spostamento <math>\eta_k</math> interno della trave <math>j</math> corrisponde lo spostamento <math>X_i</math> del telaio, altrimenti sarà pari a 0.

Passando ad analizzare i carichi nodali, si possono unire i carichi nodali in un unico vettore <math>\begin{Bmatrix} \Sigma \end{Bmatrix}</math> così definito:
<div style="text-align:center;"> <math>
\begin{Bmatrix} \Sigma \end{Bmatrix} =
\begin{Bmatrix} 
     \begin{Bmatrix} S \end{Bmatrix}_n \\
  \begin{Bmatrix} S \end{Bmatrix}_2 \\
  ...\\
  \begin{Bmatrix} S \end{Bmatrix}_n \\
\end{Bmatrix}
</math></div>

La relazione che lega il vettore <math>\begin{Bmatrix} \Sigma \end{Bmatrix}</math> (carichi nodali di ogni singola trave) al vettore <math>\begin{Bmatrix} F \end{Bmatrix}</math> (carichi nodali sull'intero telaio) è:
<div style="text-align:center;"> <math>
\begin{Bmatrix} \Sigma \end{Bmatrix} = \begin{bmatrix} I \end{bmatrix}  \begin{Bmatrix} F\end{Bmatrix}
</math></div>

Ricompare così la matrice di incidenza.

Assemblando opportunamente le matrici di rigidezza delle singole travi, si può notare che anche <math>\begin{Bmatrix}\Sigma \end{Bmatrix}</math> e <math>\begin{Bmatrix}\Eta\end{Bmatrix}</math> sono legati tra loro dalla relazione:
<div style="text-align:center;"> <math>
\begin{Bmatrix} \Sigma \end{Bmatrix} = \begin{bmatrix} \Kappa \end{bmatrix} \begin{Bmatrix} \Eta \end{Bmatrix}
</math></div>
in cui:
<div style="text-align:center;"> <math>
\begin{bmatrix} \Kappa \end{bmatrix} =
\begin{bmatrix}
\begin{bmatrix} \Kappa \end{bmatrix}_1 && 0 && 0 &&... && 0 \\
0 && \begin{bmatrix} \Kappa \end{bmatrix}_2 && 0 &&... && 0 \\
0 && 0 && \begin{bmatrix} \Kappa \end{bmatrix}_3 &&... && 0 \\
... && ... && ... &&... && ... \\
0 && 0 && 0 &&... && \begin{bmatrix} \Kappa \end{bmatrix}_n \\
\end{bmatrix} =
\begin{bmatrix}
\frac{EA}{l} & 0 & 0 & -\frac{EA}{l} & 0 & 0 \\
0 & \frac{12EJ}{l^3} & -\frac{6EJ}{l^2} & 0 & -\frac{12EJ}{l^3} & -\frac{6EJ}{l^2} \\
0 & -\frac{6EJ}{l^2} & \frac{4EJ}{l} & 0 & \frac{6EJ}{l^2} & \frac{2EJ}{l} \\
-\frac{EA}{l} & 0 & 0 & \frac{EA}{l} & 0 & 0 \\
0 & -\frac{12EJ}{l^3} & \frac{6EJ}{l^2} & 0 & \frac{12EJ}{l^3} & \frac{6EJ}{l^2} \\
0 & -\frac{6EJ}{l^2} & \frac{2EJ}{l} & 0 & \frac{6EJ}{l^2} & \frac{4EJ}{l} \\
\end{bmatrix}
</math></div>

Invertendo la relazione che lega <math>\begin{Bmatrix} \Sigma \end{Bmatrix}</math> ed <math>\begin{Bmatrix} F \end{Bmatrix}</math> e sostituendovi le relazioni precedenti, si ha:
<div style="text-align:center;"> <math>
\begin{bmatrix} F \end{bmatrix} = \begin{bmatrix} I \end{bmatrix}^T \begin{bmatrix} \Sigma \end{bmatrix} =
\begin{bmatrix} I \end{bmatrix}^T \begin{bmatrix} \Kappa \end{bmatrix} \begin{Bmatrix} \Eta \end{Bmatrix} =
\begin{bmatrix} I \end{bmatrix}^T \begin{bmatrix} \Kappa \end{bmatrix} \begin{bmatrix} I \end{bmatrix} \begin{Bmatrix} X \end{Bmatrix}=
\begin{bmatrix} H \end{bmatrix} \begin{Bmatrix} X \end{Bmatrix}
</math></div>
da cui si evince che la matrice di rigidezza della struttura è data dalla relazione:
<div style="text-align:center;"> <math>
\begin{bmatrix} H \end{bmatrix} = \begin{bmatrix} I \end{bmatrix}^T \begin{bmatrix} \Kappa \end{bmatrix} \begin{bmatrix} I \end{bmatrix}
</math></div>

== Carichi concentrati e distribuiti applicati lungo la trave (carichi infranodali) ==

Fin qui si è supposto la presenza di soli carichi applicati nei nodi. Di fatto questa ipotesi in presenza di carichi concentrati lungo l'asta sarebbe un limite facilmente superabile. In tal caso basterebbe ridefinire la struttura includendo il punto in questione tra i nodi del telaio.

In presenza di carichi distribuiti il limite sopra indicato non può essere superato. In tal caso è necessario rivedere l'impostazione fin qui vista introducendo nell'equazione che regola il problema un nuovo vettore di modo da avere:
<div style="text-align:center;"><math>
\begin{Bmatrix} F \end{Bmatrix} = \begin{bmatrix} H \end{bmatrix} \begin{Bmatrix} X \end{Bmatrix} + \begin{Bmatrix} F_0 \end{Bmatrix}
</math></div>
Il vettore <math>\begin{Bmatrix} F_0 \end{Bmatrix}</math> dipende dai carichi concentrati e distribuiti applicati sulla struttura.

Per capire il significato fisico di questo nuovo vettore, si considera il caso:
<div style="text-align:center;"><math>
\begin{Bmatrix} X \end{Bmatrix} = \begin{Bmatrix} 0 \end{Bmatrix}
</math></div>

In tal caso:
<div style="text-align:center;"><math>
\begin{Bmatrix} F \end{Bmatrix} = \begin{Bmatrix} F_0 \end{Bmatrix}
</math></div>

Quindi il vettore <math>\begin{Bmatrix}F_0\end{Bmatrix}</math> si definisce come il valore degli sforzi nodali in presenza di spostamenti nulli. Il vettore <math>\begin{Bmatrix}F_0\end{Bmatrix}</math> si presenta dipendente dalle caratteristiche geometriche e fisiche del telaio oltre che dai carichi esterni infranodali applicati.

Per il calcolo di tale vettore si procede in maniera analoga a quanto fatto sopra per il calcolo della matrice di rigidezza del telaio: si parte quindi dall'analisi di una singola trave, calcolando all'interno del suo sistema di riferimento locale i valori dei carichi nodali che annullano gli spostamenti nodali <math>\begin{Bmatrix} S_0 \end{Bmatrix}</math>; si trasformano i risultati così ottenuti dal sistema di riferimento locale a quello globale utilizzando la matrice <math>\begin{bmatrix}G\end{bmatrix}</math>; si è così arrivati ad avere un vettore <math>\begin{Bmatrix} S_0 \end{Bmatrix}</math> per ogni trave della struttura; per ottenere un unico vettore <math>\begin{Bmatrix} F_0 \end{Bmatrix}</math> si dovrà utilizzare nuovamente la matrice di incidenza <math>\begin{bmatrix} I \end{bmatrix}</math>.

Segue il calcolo del vettore <math>\begin{Bmatrix} S_0 \end{Bmatrix}</math>. A tale scopo si suppone per la trave un qualsiasi tipo di vincolo, tale da rendere la trave isostatica. Per semplicità ad esempio si suppone di avere un incastro sul primo vertice e che il secondo vertice sia libero. A seguito dell'applicazione dei carichi esterni infranodali, si ha sulla trave una deformazione assiale totale <math>\lambda(z) = \frac{N(z)}{EA(z)} + \bar{\lambda}(z) </math> ed una curvatura totale <math>\mu (z) = \frac{M(z)}{EJ(z)} + \bar{\mu}(z)</math>.

Si calcolano le reazioni vincolari nel vincolo incastrato, indicate con <math>S_1^* </math>, <math>S_2^* </math> ed <math>S_3^* </math>, tali valori vengono inseriti nel vettore <math>S^* </math>, così definito:
<div style="text-align:center;"><math>
\begin{Bmatrix} S^* \end{Bmatrix} = \begin{Bmatrix} S_1^* \\ S_2^* \\ S_3^* \\ 0 \\ 0 \\ 0 \end{Bmatrix}
</math></div>

Si calcolano ora gli spostamenti nodali nel vertice libero. Per fare questo si utilizza il [[principio dei lavori virtuali]], assumendo come struttura di confronto sempre una trave incastrata sul primo vertice e libera nel secondo. Gli spostamenti della trave saranno: 
<div style="text-align:center;"><math>
\begin{Bmatrix} \eta_1^{**} \\ \eta_2^{**} \\ \eta_3^{**} \\ \eta_4^{**} \\ \eta_5^{**} \\ \eta_6^{**} \end{Bmatrix} = \begin{Bmatrix} 0 \\ 0 \\ 0 \\ \int_0^l \lambda(z)\,dz \\ \int_0^l l \left ( \frac{z}{l} - 1 \right ) \mu(z) \,dz \\ \int_0^l \mu(z) \,dz \end{Bmatrix}
</math></div>

Le forze nodali <math>\begin{Bmatrix} S^{**} \end{Bmatrix}</math> possono essere calcolate servendosi della matrice di rigidezza secondo la relazione:
<div style="text-align:center;"><math>
\begin{Bmatrix} S^{**} \end{Bmatrix} = \begin{bmatrix} K \end{bmatrix} \begin{Bmatrix} \eta^{**} \end{Bmatrix} 
</math></div>

Per il [[principio di sovrapposizione degli effetti]] il vettore <math>\begin{Bmatrix} S_0 \end{Bmatrix} </math> cercato sarà dato dalla relazione:
<div style="text-align:center;"><math>
\begin{Bmatrix} S_0 \end{Bmatrix} = \begin{Bmatrix} S^{*} \end{Bmatrix} + \begin{Bmatrix} S^{**} \end{Bmatrix} 
</math></div>

Il vettore <math>\begin{Bmatrix} S_0 \end{Bmatrix}</math> così ottenuto si riferisce al sistema di riferimento locale della trave analizzata. Per il passaggio al sistema di riferimento globale utilizziamo la matrice <math>\begin{bmatrix} G \end{bmatrix}</math> sopra definita, secondo la relazione:
<div style="text-align:center;"><math>
\begin{Bmatrix} S_{0,G} \end{Bmatrix} = \begin{bmatrix} G \end{bmatrix} \begin{Bmatrix} S_0 \end{Bmatrix} 
</math></div>

Rimane infine da assemblare il vettore <math>\begin{Bmatrix} F_{0} \end{Bmatrix}</math> dell'intera struttura. Per fare questo si dovrà analizzare l'equilibrio di ciascun nodo della struttura. In questo modo si ottengono nuovamente i termini della matrice di incidenza, secondo:
<div style="text-align:center;"><math>
\begin{Bmatrix} F_0 \end{Bmatrix} = \begin{bmatrix} I \end{bmatrix} \begin{Bmatrix} S_{0,G} \end{Bmatrix} 
</math></div>

== Non linearità ==

Finora si è sviluppato il modello assumendo un comportamento elastico lineare del materiale; più nello specifico si è supposto l'esistenza di una relazione lineare tra curvatura, deformazione media assiale e caratteristiche di sollecitazione. 

I materiali utilizzati nella tecnica corrente si discostano molto da tale comportamento teorico. Tale aspetto viene genericamente individuato con l'espressione '''non-linearità meccanica'''.

Un'altra ipotesi fatta riguarda l'influenza delle deformazioni sullo stato finale di equilibrio: si calcolano le caratteristiche di sollecitazione della trave riferendoci alla configurazione indeformata ma se si considera ad esempio il caso di un pilastro snello incastrato alla base e soggetto sull'estremo libero a due carichi, uno verticale ed uno orizzontale, le deformazioni assunte dal pilastro fanno sì che il carico verticale determini la comparsa di sollecitazioni flessionali sulla trave per [[carico di punta]].

Nel seguito si modificherà il modello di modo che riesca a considerare anche queste problematiche.

=== Non linearità meccanica ===
Per analizzare le problematiche legate alla non linearità meccanica si dovrà scendere ad un maggior livello di dettaglio, arrivando ad analizzare il comportamento dei singoli punti costituenti la sezione della trave.

Si suppone di avere una sezione definita tramite un insieme di punti riferiti ad un certo sistema di riferimento. All'interno dell'ipotesi di conservazione delle sezioni piane, la deformazione di un generico punto costituente la sezione può essere ottenuta a partire da quella dell'origine della sezione secondo:

<div style="text-align:center;"><math>
\begin{matrix}
\varepsilon_z (y) = \lambda + \mu \ y \\
\gamma = \eta
\end{matrix}
</math></div>

Assumendo per il materiale costituente la sezione un [[Relazione costitutiva (meccanica)|legame costitutivo]] di tipo elastico-lineare isotropo, si dovrà associare a tali deformazioni delle tensioni ottenute secondo:

<div style="text-align:center;"><math>
\begin{matrix}
\sigma_z (y) = E \left( \lambda + \mu \ y  \right ) \\
\tau = G \ \eta
\end{matrix}
</math></div>

Tali tensioni si traducono in un legame lineare tra le caratteristiche di sollecitazione ed i parametri deformativi.

<div style="text-align:center;"><math>
\begin{matrix}
\begin{Bmatrix}N \\ M \\ T \end{Bmatrix} =
\begin{bmatrix}E A & ES & 0\\ ES & EJ & 0\\ 0 & 0 & tGA\end{bmatrix} \begin{Bmatrix} \lambda \\ \mu \\ \eta \end{Bmatrix}
\end{matrix}
</math></div>

Si trascura il taglio (per travi snelle il taglio ha un contributo ridotto rispetto alle altre sollecitazioni interne) e ci si concentra solo su [[momento flettente]] e sforzo normale (ipotesi più che plausibile nella maggior parte dei casi). La relazione sopra esposta mostra un legame di tipo lineare tra sforzo normale, momento flettente e le caratteristiche di deformazione della sezione.

Sia l'acciaio che il calcestruzzo hanno un comportamento diverso da quello lineare fin qui descritto. Per entrambi i materiali la non-linearità del legame costitutivo del materiale si traduce in una non-linearità del legame tra sforzo normale, momento flettente e caratteristiche deformative. Si dovrà quindi modificare il modello teorico per fare in modo che tenga conto di quest'ultimo aspetto. 

Purtroppo non è disponibile una formula chiusa che permetta di calcolare il telaio sotto queste ipotesi, servirà perciò utilizzare dei metodi iterativi.

Si possono individuare nella trave un certo numero significativo di sezioni e su queste andare a verificare la discrepanza tra modello teorico lineare e comportamento effettivo del materiale.

A questo punto un primo approccio è quello di modificare i parametri di rigidezza della sezione, sostanzialmente area e momento di inerzia, per riportare il comportamento teorico in quello effettivo. 

<div style="text-align:center;"><math>
\begin{matrix}
\begin{Bmatrix}N \\ M \end{Bmatrix} =
\begin{bmatrix}E A^* & E S^* \\ E S^* & E J^* \end{bmatrix} \begin{Bmatrix} \lambda \\ \mu \end{Bmatrix}
\end{matrix}
</math></div>

A questo punto si deve ricalcolare la matrice di rigidezza della trave come se si trattasse di una nuova trave a sezione variabile o analogamente, un insieme di travi più piccole collegate tra loro. Si deve quindi ricalcolare tutte le matrici di rigidezza delle singole travi e successivamente riassemblare la matrice di rigidezza del telaio. Risolvere il sistema lineare così definito e verificare nuovamente quale discrepanza c'è tra il modello adottato e quello effettivo del materiale. Si procede così fino a quando lo scarto tra i due sarà sufficientemente piccolo.

Computazionalmente questo approccio presenta non poche difficoltà. Il difetto principale del metodo appena descritto è nella necessità di ricalcolare più volte la matrice di rigidezza e successivamente di doverla risolvere. Per telai minimamente realistici questo vuol dire avere a disposizione una notevole potenza di calcolo.

Un approccio meno dispendioso sotto questo punto di vista lascia invece invariati i parametri di rigidezza della sezione e trasla la legge costitutiva fino a portarla a coincidere con quella effettiva del materiale. Traducendo il tutto in formule, questo vuol dire introdurre due coefficienti <math> \bar \lambda </math> e <math> \bar \mu </math> tali per cui il legame della sezione diventi:

<div style="text-align:center;"><math>
\begin{Bmatrix}N \\ M \end{Bmatrix} =
\begin{bmatrix}E A & E S \\ E S & E J \end{bmatrix} \begin{Bmatrix} \lambda + \bar \lambda \\ \mu + \bar \mu \end{Bmatrix}
</math></div>

Le grandezze <math> \bar \lambda </math> e <math> \bar \mu </math> si presentano come deformazioni impresse sulla sezione e quindi alla stregua di normali carichi infranodali. Utilizzando l'approccio descritto sopra, si può tradurre quindi tali parametri in vettori <math> \begin{Bmatrix} S_0 \end{Bmatrix} </math> e quindi in un vettore <math> \begin{Bmatrix} F_0 \end{Bmatrix} </math> proprio della struttura.

A questo punto si è evitato il ricalcolo della matrice di rigidezza della struttura però si deve comunque risolvere il sistema:
<div style="text-align:center;"><math>
\begin{Bmatrix} F \end{Bmatrix} = \begin{bmatrix} H  \end{bmatrix} \begin{Bmatrix} X \end{Bmatrix} + \begin{Bmatrix} F_0 \end{Bmatrix}
</math></div>

A questo punto bisogna ricalcolare nuovamente lo scarto tra modello teorico e modello reale e ripetere il procedimento sopra indicato fino ad ottenere un livello sufficiente di precisione.

Volendo risparmiare tempo si può apportare una piccola modifica al procedimento sopra esposto: la prima volta che si risolve il sistema <math> \begin{Bmatrix} F \end{Bmatrix} = \begin{bmatrix} H  \end{bmatrix} \begin{Bmatrix} X \end{Bmatrix} + \begin{Bmatrix} F_0 \end{Bmatrix} </math>, invertire la matrice <math> \begin{bmatrix} H  \end{bmatrix} </math>. Il maggiore onere computazionale necessario tornerà utile nei passaggi successivi, allorché sarà sufficiente calcolare:
<div style="text-align:center;"><math>
\begin{Bmatrix} X \end{Bmatrix} = \begin{bmatrix} H^{-1}  \end{bmatrix} \left( \begin{Bmatrix} F \end{Bmatrix} - \begin{Bmatrix} F_0 \end{Bmatrix} \right)
</math></div>
in cui i termini a destra del segno uguale sono tutti noti.

Non è detto che il procedimento iterativo così descritto porti sempre ad una soluzione; si può ad esempio verificare che il telaio non è in grado di sopportare i carichi imposti.

=== Non linearità geometrica ===
Nella stragrande maggioranza dei casi ingegneristici si possono ritenere trascurabili gli effetti delle deformazioni della struttura sull'entità delle [[Azione esterna|sollecitazioni]] e sui fenomeni di instabilità.

Pertanto fino ad ora l'analisi del telaio si è potuta effettuare con la [[teoria del primo ordine]] ('''analisi del primo ordine''') e cioè imponendo l'equilibrio sulla configurazione iniziale della struttura (struttura indeformata).

In presenza di elementi strutturali snelli gli spostamenti non più trascurabili prodotti dalle azioni applicate determinano l'insorgere di un'eccentricità del carico assiale agente, con conseguente formazione di un momento flettente ([[carico di punta]]) o incremento di quello già presente sull'elemento strutturale (strutture pressoinflesse).

In questo caso di parla di '''effetti del secondo ordine''' e il relativo momento flettente addizionale è detto '''momento del secondo ordine'''.

Tale fenomeno influenza notevolmente sia la deformabilità in esercizio, sia la capacità resistente ultima di una struttura snella.

Le verifiche di stabilità si devono condurre attraverso un''''analisi del secondo ordine''', imponendo l'equilibrio sulla configurazione deformata della struttura.

== Voci correlate ==
*[[Statica delle strutture]]
*[[Metodo delle forze]]
*[[Metodo delle rigidezze]]
{{Portale|ingegneria}}

[[Categoria:Analisi strutturale]]