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 (traslazioni e rotazioni) di alcuni punti, definibili attraverso valori scalari reali.
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).
Riportando gli spostamenti nodali su un vettore verticale e gli sforzi nodali applicati sui medesimi nodi su un vettore verticale , la relazione che lega le due grandezze è del tipo: . in cui è una matrice associata alla struttura definita matrice di rigidezza della struttura. Per il teorema di reciprocità o di Betti, la matrice è 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 gli spostamenti dei punti estremi della trave e con i carichi applicati nei due punti, le due grandezze saranno legate da una relazione del tipo: . in cui è la matrice di rigidezza della trave che è del tipo:
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.
In particolare si suppone che la sezione della trave sia costante e di area . Il momento di inerzia calcolato lungo un asse perpendicolare al piano del telaio e passante per il baricentro della sezione è . 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 si considerano gli schemi statici.
Partendo dal calcolo dei coefficienti che riguardano i carichi nodali e , si considera una trave vincolata sul secondo estremo, libera sul primo e sottoposta su quest'ultimo estremo ad un carico assiale . Supponendo di trascurare gli effetti del second'ordine (non linearità geometriche), si è in presenza di uno spostamento e . La relazione che lega a è:
in cui è la deformazione unitaria dell'asta.
La reazione vincolare assiale sul secondo estremo è data da:
Procedendo analogamente per , si può ricavare:
Si considera ora una trave semplicemente appoggiata, sottoposta nei due estremi a due momenti e .
Con alcuni calcoli si ottiene che la relazione che lega gli spostamenti nodali e agli sforzi nodali e è:
Invertendo la relazione suddetta si ha:
che permette di trovare e .
Calcolando le reazioni vincolari verticali sui due estremi, si trova:
da cui:
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 . In tal caso lo spostamento verticale del secondo estremo è dato dalla relazione:
da cui:
Sempre facendo riferimento alla configurazione sopra descritta, la reazione vincolare verticale nel primo estremo, è pari a . Pertanto:
Per simmetria si trova e pari rispettivamente a e .
Ricapitolando quanto visto finora si ha la seguente matrice di rigidezza:
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 e rispettivamente spostamenti e sollecitazioni nodali nel sistema di riferimento globale, i due dovranno essere legati tra loro dalla relazione:
.
Per fare questo si inizia caratterizzando la trave con un versore bidimensionale orientato come la stessa asta ed avente per definizione modulo unitario. Supponendo che gli estremi della trave siano i due punti e , le componenti di sono date da:
Si noti che assumendo pari ad l'inclinazione della trave rispetto all'asse delle ascisse, le componenti diventano:
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 e gli spostamenti nel sistema di riferimento globale è:
in cui la matrice è data da:
è una matrice simmetrica con una particolarità: la sua inversa è uguale alla sua trasposta .
Analogamente la relazione che lega le sollecitazioni nodali nel sistema locale a quelle nel sistema globale è:
Invertendo le suddette relazioni, si ottiene:
che sostituite in:
forniscono:
da cui si ottiene la matrice di rigidezza cercata:
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 travi, si uniscono i vettori spostamenti nodali delle singole travi ottenendo un vettore degli spostamenti nodali del telaio disarticolato così definito:
Lo scopo è passare dal vettore (spostamenti nodali di ogni singola trave) al vettore degli spostamenti nodali del telaio. Perciò si ricerca un operatore che leghi le due grandezze secondo:
è definita matrice di incidenza ed è una matrice , in cui è il numero degli spostamenti nodali del telaio, è il numero complessivo di spostamenti nodali per ogni trave (4 traslazioni e 2 rotazioni) e è il numero di aste. Il generico elemento sarà pari a 1 se allo spostamento interno della trave corrisponde lo spostamento del telaio, altrimenti sarà pari a 0.
Passando ad analizzare i carichi nodali, si possono unire i carichi nodali in un unico vettore così definito:
La relazione che lega il vettore (carichi nodali di ogni singola trave) al vettore (carichi nodali sull'intero telaio) è:
Ricompare così la matrice di incidenza.
Assemblando opportunamente le matrici di rigidezza delle singole travi, si può notare che anche e sono legati tra loro dalla relazione:
in cui:
Invertendo la relazione che lega ed e sostituendovi le relazioni precedenti, si ha:
da cui si evince che la matrice di rigidezza della struttura è data dalla relazione:
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:
Il vettore dipende dai carichi concentrati e distribuiti applicati sulla struttura.
Per capire il significato fisico di questo nuovo vettore, si considera il caso:
In tal caso:
Quindi il vettore si definisce come il valore degli sforzi nodali in presenza di spostamenti nulli. Il vettore 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 ; si trasformano i risultati così ottenuti dal sistema di riferimento locale a quello globale utilizzando la matrice ; si è così arrivati ad avere un vettore per ogni trave della struttura; per ottenere un unico vettore si dovrà utilizzare nuovamente la matrice di incidenza .
Segue il calcolo del vettore . 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 ed una curvatura totale .
Si calcolano le reazioni vincolari nel vincolo incastrato, indicate con , ed , tali valori vengono inseriti nel vettore , così definito:
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:
Le forze nodali possono essere calcolate servendosi della matrice di rigidezza secondo la relazione:
Per il principio di sovrapposizione degli effetti il vettore cercato sarà dato dalla relazione:
Il vettore così ottenuto si riferisce al sistema di riferimento locale della trave analizzata. Per il passaggio al sistema di riferimento globale utilizziamo la matrice sopra definita, secondo la relazione:
Rimane infine da assemblare il vettore 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:
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.