Relazione Sismica Finita
Relazione Sismica Finita
Relazione Sismica Finita
FACOLTÀ DI INGEGNERIA
Corso di laurea in Ingegneria Civile
Specialistica in Strutture
Dipartimento di Modellistica per l’Ingegneria
Corso di Costruzioni in Zona Sismica
Progetto di:
Analisi sismiche di un telaio spaziale
Docente: Studenti:
prof. ing. Alfonso VULCANO Marco LANZA
Paolo TALARICO
Esercitatore: Giuseppe TOMAINO
ing. Fabio MAZZA
INDICE
1. PREDIMENSIONAMENTO.......................................................................................1
1.1. Introduzione..........................................................................................................1
1.2. Dati geometrici.....................................................................................................1
1.3. Caratteristiche dei materiali..................................................................................3
1.4. Azioni sulla costruzione.......................................................................................3
1.5. Analisi dei carichi.................................................................................................4
1.6. Combinazione dei carichi.....................................................................................5
1.7. Predimensionamento pilastri................................................................................6
1.8. Predimensionamento travi....................................................................................9
2. ANALISI STATICA (SHEAR-TYPE).......................................................................22
2.1. Introduzione........................................................................................................22
2.2. Calcolo baricentri delle masse e delle rigidezze................................................22
2.3. Ripartizione taglianti di piano............................................................................24
2.4. Ripartizione taglianti di piano con eccentricità accidentale...............................27
2.5. Commenti...........................................................................................................31
3. ANALISI STATICA (CONDENSAZIONE STATICA)............................................33
3.1. Introduzione........................................................................................................33
3.2. Condensazione statica........................................................................................33
3.3. Ripartizione taglianti di piano............................................................................35
3.4. Commenti...........................................................................................................38
4. ANALISI DINAMICA (PSEUDO-DINAMICA)......................................................40
4.1. Introduzione........................................................................................................40
4.2. Analisi modale....................................................................................................40
4.3. Ripartizione taglianti di piano............................................................................43
4.4. Commenti...........................................................................................................48
i
Capitolo 1 - Predimensionamento
1. PREDIMENSIONAMENTO
1.1. Introduzione
Tutte le grandezze geometriche del problema sono riassunte in tab. 1.1 e fanno riferimento
alle schematizzazioni riportate in fig. 1.1.
a 5.00 m
b 6.00 m
c 4.00 m
d 5.00 m
h1 4.00 m
h2 3.30 m
h3 3.30 m
Tab. 1.1 – Dati geometrici struttura
1
Capitolo 1 - Predimensionamento
Per quanto riguarda il solaio è stato utilizzato come criterio per il predimensionamento il
seguente:
1
HS ≥ {c , d }
25
Ottenendo:
HS 20.0 cm
s 4.0 cm
hL 16.0 cm
sI 1.5 cm
sM 3.0 cm
b0 10.0 cm
i 50.0 cm
bL 40.0 cm
Tab. 1.2 – Dati geometrici solaio
2
Capitolo 1 - Predimensionamento
3
Capitolo 1 - Predimensionamento
• Carichi variabili
2.00 kN/m2 carico di esercizio
neve zona III, a S =200 m
V N =50 anni
CU =1
V R =V N⋅C U=50 anni
In base a quanto visto finora è stata svolta l'analisi dei carichi, che determina l'influenza che
hanno pesi e azioni relativamente ad un metro lineare.
• Carichi permanenti strutturali
g S=25⋅s⋅1m soletta
g N =25⋅b 0⋅h L⋅2 nervature
g L=8⋅b L⋅h L⋅2 laterizi
• Carichi permanenti non strutturali
g M =15⋅s M⋅1m massetto
g I =15⋅s I⋅1m intonaco
g P =0.4⋅1m pavimento
g TRAM=1.2⋅1m tramezzi
g TAMP=2.7⋅h i tamponature
g IMP=0.3⋅1m impermeabilizzazione
• Carichi variabili
QK1=i⋅q SK⋅cE⋅cT neve
QK2=2.0⋅1m manutenzione
I valori ottenuti sono riportati nella tab. 1.4.
4
Capitolo 1 - Predimensionamento
gS 1.00 kN/m
gN 1.00 kN/m
gL 1.02 kN/m
gM 0.45 kN/m
gI 0.23 kN/m
gP 0.40 kN/m
gTRAM 1.20 kN/m
gTAMP,3 8.91 kN/m2
gTAMP,2 8.91 kN/m2
gIMP 0.30 kN/m
QK1 0.48 kN/m2
QK2 2.00 kN/m2
Tab 1.4 – Riepilogo analisi dei carichi
Facendo riferimento alla sola combinazione per stati limite ultimi (SLU), come prescritto da
normativa si è fatto riferimento alle seguenti tre combinazioni di carichi:
• Combinazione 1 – Soli carichi verticali
n
F d1=G1⋅G1G2⋅G2 Q1⋅QK1∑ Qi⋅ 0i⋅QKi
i=1
• Combinazione 2/3 – Soli carichi verticali in presenza di sisma nelle due direzioni
n
F d2,3 =±EG1G 2∑ 2i⋅QKi
i=1
5
Capitolo 1 - Predimensionamento
Questo metodo ripartisce per ogni impalcato il peso del solaio, il peso delle travi trasversali e
il peso della tamponatura secondo delle aree d'influenza proprie per ognuno dei 9 pilastri.
Ogni pilastro regge mezza luce.
Tali aree non rappresentano la reale ripartizione e ipotizzano trave schematizzata come
6
Capitolo 1 - Predimensionamento
In seguito vengono tabellati i risultati ottenuti finora facendo riferimento alla fig. 1.6.
7
Capitolo 1 - Predimensionamento
Pilastro N1 N2 N3 N4 N5 N6 N7 N8 N9
Liv. 3 73.16 190.70 83.91 188.21 98.82 166.19 85.99 162.84 365.97
Liv. 2 197.17 467.09 226.19 450.58 260.81 415.10 228.99 396.87 780.58
Liv. 1 327.71 754.50 373.64 723.97 429.81 675.31 378.76 641.63 1217.35
Tab 1.6 – Sforzo normale alla base dei pilastri [kN]
Pilastro N1 N2 N3 N4 N5 N6 N7 N8 N9
Liv. 3 0.02 0.03 0.02 0.03 0.02 0.03 0.02 0.03 0.05
Liv. 2 0.05 0.08 0.05 0.08 0.06 0.07 0.05 0.07 0.11
Liv. 1 0.08 0.13 0.09 0.13 0.10 0.12 0.09 0.11 0.17
Tab 1.7 – Aree teoriche pilastri [m2]
8
Capitolo 1 - Predimensionamento
- risega pari a 5.0 cm o 10.0 cm, per avere una graduale ripartizione delle rigidezze e
una migliore risposta sismica;
N
- sforzo normale adimensionalizzato n= ≤0.5 .
b PIL⋅h PIL⋅f c ,d
Per quanto precedentemente detto si è ottenuto:
Pilastro A1 A2 A3 A4 A5 A6 A7 A8 A9
Liv. 3 0.09 0.09 0.09 0.09 0.09 0.09 0.09 0.09 0.09
Liv. 2 0.11 0.12 0.11 0.12 0.11 0.12 0.11 0.12 0.16
Liv. 1 0.12 0.15 0.12 0.15 0.12 0.15 0.12 0.15 0.25
2
Tab 1.8 – Aree effettive pilastri [m ]
Pilastro n1 n2 n3 n4 n5 n6 n7 n8 n9
Liv. 3 0.06 0.15 0.07 0.15 0.08 0.13 0.07 0.13 0.29
Liv. 2 0.13 0.28 0.15 0.27 0.18 0.25 0.15 0.23 0.35
Liv. 1 0.19 0.36 0.22 0.34 0.25 0.32 0.22 0.30 0.35
Tab 1.9 – Sforzo normale adimensionalizzato
9
Capitolo 1 - Predimensionamento
Tali aree non rappresentano la reale ripartizione e ipotizzano trave schematizzata come
semplicemente appoggiata, quindi nella valutazione si tengono conto di coefficienti correttivi
da applicare ad ogni traverso.
k tra, Esterna =1.00 k tra , Interna =1.25
Sono stati ottenuti i carichi ripartiti linearmente sulle varie travi per ogni condizione di carico,
riportati in allegato 1.
In tab. 1.10 è riportato un sommario di tali carichi ripartiti.
10
Capitolo 1 - Predimensionamento
7,00
6,00
5,00
4,00
Sa [m/s2]
elastico
inelastico
3,00
2,00
1,00
0,00
0,00 0,50 1,00 1,50 2,00 2,50 3,00 3,50 4,00
To [s]
11
Capitolo 1 - Predimensionamento
∑ Wi
i=1
F h= ⋅Sa T1
g
dove:
- Wi è il peso sismico di ogni impalcato;
- g è l'accelerazione di gravità;
- =0.85 se n P≥3 e T1≤2 TC ;
- Sa è la pseudo accelerazione spettrale.
Infine è stata ripartita questa forzante sismica totale, applicata in modo statico, ai vari piani in
maniera pesata.
W i⋅z i
F hi =F h nP
∑ W j⋅z j
j=1
W3 674.11 kN
W2 1078.92 kN
W1 1134.51 kN
T1 0.44 s
Sa 1.26 m/s2
Fh 314,61 kN
Fh3 114.93 kN
Fh2 126.68 kN
Fh1 72.99 kN
Tab 1.11 – Riepilogo risultati
Le forze sismiche di impalcato sono state ripartite su tutti i telai del piano considerato
mediante aliquote ottenute per esperienza, come rappresentato in figura.
È stato riportato in allegato 3 una schematizzazione.
12
Capitolo 1 - Predimensionamento
Sisma in x Sisma in y
esterno 38.31 45.97
Liv. 3
interno 38.31 22.99
esterno 42.23 50.67
Liv. 2
interno 42.23 25.34
esterno 24.33 29.20
Liv. 1
interno 24.33 14.60
Tab 1.12 – Forze concentrate sui telai [kN]
Le due modalità di carico sui telai sono state risolte separatamente per la determinazione del
momento flettente sulle travi. La prima, carichi verticali, è stata risolta ipotizzando trave
continua appoggiata; la seconda, carichi orizzontali, è stata risolta ipotizzando telaio con
schema shear-type.
Si sono ottenuti gli andamenti riportati in allegato 2 e 3.
Si riorganizzano i risultati di momento flettente massimo nella seguente tabella in base
all'impalcato, telaio, condizione di carico e direzione del sisma.
13
Capitolo 1 - Predimensionamento
- h=
M
⋅
1
b f cd⋅0.8⋅0.259⋅1−0.4⋅0.259
M 1
- b= ⋅
h f cd⋅0.8⋅0.259⋅1−0.4⋅0.259
2
14
Capitolo 1 - Predimensionamento
Telaio in x Telaio in y
esterno 0.38 0.24
Liv. 3
interno 0.43 0.17
esterno 0.50 0.41
Liv. 2
interno 0.53 0.32
esterno 0.56 0.50
Liv. 1
interno 0.59 0.55
Tab. 1.15 – Dimensioni teoriche travi [m]
Si sono tenute in considerazione, in questa fase, alcune accortezze per quanto riguardano le
travi in spessore:
- b≤b PIL H TRA ;
- b≤2⋅bPIL .
Telaio in x Telaio in y
esterno 0.50 0.50
Liv. 3
interno 0.50 0.50
esterno 0.55 0.55
Liv. 2
interno 0.55 0.55
esterno 0.60 0.60
Liv. 1
interno 0.60 0.60
Tab. 1.16 – Dimensioni effettive travi [m]
Viene riportata di seguito la carpenteria delle travi e dei pilastri in scala 1:100.
15
Capitolo 1 - Predimensionamento
16
Capitolo 1 - Predimensionamento
17
Capitolo 1 - Predimensionamento
18
Capitolo 1 - Predimensionamento
19
Capitolo 1 - Predimensionamento
20
Capitolo 1 - Predimensionamento
21
Capitolo 2 – Analisi statica (shear-type)
2.1. Introduzione
In questo capitolo verrà affrontata una tipologia di analisi sismica che assume forzante
sismica applicata in modo statico sulla struttura. Tale analisi è molto semplificata e ipotizza,
su ogni telaio, un solo grado di libertà per ogni piano. È una forte semplificazione quella di
considerare il numero minimo di cinematismi della struttura che, per alcune tipologie di
strutture, come la presente, non effettua errori grossolani.
Questa analisi permette di ripartire attraverso le rigidezze dei telai la forzante sismica su ogni
telaio e livello.
Per lo sviluppo di questa analisi sismica e delle successive, bisogna determinare la posizione
del baricentro delle masse (G) e del baricentro delle rigidezze (C) per ogni impalcato.
In G si ipotizza concentrata l'intera massa dell'impalcato e quindi la relativa forza sismica.
Vista la difficile determinazione si schematizza coincidente con il baricentro geometrico
dell'impalcato, valutabile come:
ab cd
xG= yG =
2 2
In C si ipotizza concentrata l'intera rigidezza dell'impalcato e quindi intorno a questo punto
esso ruoterà.
22
Capitolo 2 – Analisi statica (shear-type)
nTEL nTEL
∑ k yi⋅x i ∑ k xj⋅y j
i=1
x C= nTEL
yC = j=1n TEL
∑ k yi ∑ k xj
i=1 j=1
dove:
- xj è yi sono le distanze dei singoli telai dal baricentro G;
- kxj e kyi sono le rigidezze laterali dei singoli telai.
9 9
E⋅I Xn E⋅I Yn
k yi =12 ∑ 3
k xj =12 ∑
n =1 h n n=1 h 3n
Si ottiene allora che G coincide per ogni piano, mentre C tende a spostarsi ad ogni piano verso
i telai più rigidi.
In fig. 2.1 viene riportato il caso in esame con le grandezze introdotte finora.
23
Capitolo 2 – Analisi statica (shear-type)
Rispetto G
x1 x2 x3 y1 y2 y3
-5.50 -0.50 5.50 -4.50 0.50 4.50
Rispetto C
x1' x2' x3' y1' y2' y3'
Liv. 3 -5.33 -0.33 5.67 -4.67 0.33 4.33
Liv. 2 -5.32 -0.32 5.68 -4.68 0.32 4.32
Liv. 1 -5.30 -0.30 5.70 -4.70 0.30 4.30
Tab 2.2 – Posizione telai rispetto i baricentri [m]
Per il cemento armato è possibile fare l'ipotesi di telaio che lavora secondo il modello shear-
type. Questo ipotizza:
- inestensibilità per tutti gli elementi ( EA ∞ );
- infinita rigidezza flessionale per le travi ( EJ TRA ∞ );
- infinita rigidezza tagliante per tutti gli elementi ( GA ∞ ).
Ne consegue che esso potrà solo subire una traslazione orizzontale in corrispondenza dei
traversi.
Trasportando le stesse ipotesi sulla struttura 3D e quindi sull'impalcato, ne consegue che esso
è inestensibile e infinitamente rigido flessionalmente e potrà solo subire due traslazioni nel
proprio piano e una rotazione secondo l'asse verticale.
Si conclude che ogni struttura 3D, sotto queste ipotesi, ha 3 gradi di libertà per impalcato.
La cinematica viene schematizzata nella seguente figura.
24
Capitolo 2 – Analisi statica (shear-type)
Scrivendo le equazioni di equilibrio per ogni piano è stato possibile determinare le tre
incognite cinematiche u, v e φ dell'impalcato per sisma orientato lungo x e lungo y.
n TEL n TEL
u ∑ k xj ∑ k xj⋅y j=Fh
j=1 j=1
nTEL nTEL
v ∑ k yi − ∑ k yi⋅x i=0
i=1 i=1
n TEL nTEL nTEL n TEL
Le precedenti equazioni di equilibrio sono state scritte rispetto il baricentro delle masse G. È
possibile semplificare considerevolmente l'intero sistema scrivendolo rispetto il baricentro
delle rigidezze C, annullando così tutti i momenti statici. Si ottiene:
nTEL
u ' ∑ k xj=F h
j=1
nTEL
v ' ∑ k yi =0
i=1
n TEL nTEL
La ripartizione della forzante sismica è stata fatta, per sisma in direzione x e y, come:
F xj =k xj u⋅y j F yi =k yi v−⋅xi
In ultimo i taglianti alla base dei pilastri sono stati valutati come la sommatoria delle forzanti
25
Capitolo 2 – Analisi statica (shear-type)
Sisma lungo x
X1 X2 X3 Y1 Y2 Y3
Liv. 3 39.19 38.25 37.49 -1.01 -0.06 1.07
Liv. 2 41.44 45.83 39.41 -1.20 -0.08 1.29
Liv. 1 22.62 29.04 21.33 -0.76 -0.06 0.82
Sisma lungo y
X1 X2 X3 Y1 Y2 Y3
Liv. 3 0.88 -0.06 -0.82 37.30 38.25 39.38
Liv. 2 1.06 -0.08 -0.98 39.18 45.83 41.67
Liv. 1 0.67 -0.06 -0.62 21.19 29.04 22.76
Tab. 2.4 – Forzante sismica [kN]
Sisma lungo x
X1 X2 X3 Y1 Y2 Y3
Liv. 3 39.19 38.25 37.49 -1.01 -0.06 1.07
Liv. 2 80.64 84.08 76.90 -2.21 -0.15 2.36
Liv. 1 103.26 113.13 98.23 -2.97 -0.20 3.18
Sisma lungo y
X1 X2 X3 Y1 Y2 Y3
Liv. 3 0.88 -0.06 -0.82 37.30 38.25 39.38
Liv. 2 1.94 -0.15 -1.80 76.48 84.08 81.05
Liv. 1 2.62 -0.20 -2.41 97.67 113.13 103.82
Tab. 2.5 – Taglianti alla base [kN]
26
Capitolo 2 – Analisi statica (shear-type)
Considerato che nella precedente analisi eseguita la posizione del baricentro G è approssimata
(v. par. 2.2), la normativa si cautela a riguardo introducendo un'ulteriore ipotesi. Pone un
range di variazione al baricentro G (eccentricità accidentale) in dipendenza della direzione del
sisma e della geometria dell'impalcato.
L'eccentricità accidentale da considerare è pari al ±5% della dimensione dell'edificio valutata
ortogonalmente alla direzione del sisma. È riportata in fig. 2.3 la rappresentazione.
Ne consegue che bisogna ripetere identicamente l'analisi del par. 2.3 ulteriori 4 volte
spostando il baricentro G, contemporaneamente su tutto l'edificio, dell'eccentricità
accidentale.
Per sisma lungo x si è ottenuto:
27
Capitolo 2 – Analisi statica (shear-type)
X1 X2 X3 Y1 Y2 Y3
Liv. 3 36.81 38.42 39.70 1.71 0.11 -1.82
Liv. 2 38.81 46.04 41.83 1.79 0.12 -1.91
Liv. 1 21.10 29.17 22.72 0.96 0.07 -1.03
Tab. 2.10 – Forzante sismica +5% [kN]
28
Capitolo 2 – Analisi statica (shear-type)
X1 X2 X3 Y1 Y2 Y3
Liv. 3 41.58 38.08 35.28 -3.73 -0.23 3.96
Liv. 2 44.08 45.63 36.98 -4.20 -0.29 4.48
Liv. 1 24.14 28.91 19.94 -2.48 -0.19 2.66
Tab. 2.11 – Forzante sismica -5% [kN]
X1 X2 X3 Y1 Y2 Y3
Liv. 3 36.81 38.42 39.70 1.71 0.11 -1.82
Liv. 2 75.62 84.46 81.54 3.50 0.23 -3.73
Liv. 1 96.72 113.63 104.26 4.46 0.30 -4.76
Tab. 2.12 – Taglianti alla base +5% [kN]
X1 X2 X3 Y1 Y2 Y3
Liv. 3 41.58 38.08 35.28 -3.73 -0.23 3.96
Liv. 2 85.65 83.71 72.26 -7.93 -0.52 8.45
Liv. 1 109.80 112.62 92.19 -10.40 -0.71 11.11
Tab. 2.13 – Taglianti alla base -5% [kN]
29
Capitolo 2 – Analisi statica (shear-type)
X1 X2 X3 Y1 Y2 Y3
Liv. 3 3.79 -0.27 -3.52 33.98 38.04 42.92
Liv. 2 4.28 -0.33 -3.95 35.52 45.59 45.58
Liv. 1 2.53 -0.22 -2.32 19.09 28.88 25.02
Tab. 2.18 – Forzante sismica +5% [kN]
X1 X2 X3 Y1 Y2 Y3
Liv. 3 -2.03 0.14 1.88 40.63 38.46 35.85
Liv. 2 -2.16 0.17 1.99 42.84 46.08 37.77
Liv. 1 -1.19 0.10 1.09 23.28 29.20 20.51
Tab. 2.19 – Forzante sismica -5% [kN]
X1 X2 X3 Y1 Y2 Y3
Liv. 3 3.79 -0.27 -3.52 33.98 38.04 42.92
Liv. 2 8.07 -0.60 -7.47 69.50 83.63 88.49
Liv. 1 10.61 -0.82 -9.79 88.59 112.51 113.51
Tab. 2.20 – Taglianti alla base +5% [kN]
30
Capitolo 2 – Analisi statica (shear-type)
X1 X2 X3 Y1 Y2 Y3
Liv. 3 -2.03 0.14 1.88 40.63 38.46 35.85
Liv. 2 -4.19 0.31 3.87 83.47 84.54 73.61
Liv. 1 -5.37 0.41 4.96 106.75 113.74 94.12
Tab. 2.21 – Taglianti alla base -5% [kN]
2.5. Commenti
La precedente analisi è stata sviluppata determinando la posizione del baricentro delle masse,
supposto per semplicità di calcolo coincidente con il baricentro geometrico di ogni impalcato,
e il baricentro delle rigidezze, variabile su ogni impalcato. Da notare che questi due punti non
coincidendo tra loro creano effetti torcenti sulla struttura.
In seguito si è imposto l'equilibrio statico ai vari impalcati per ogni grado di libertà
permettendo di determinare, attraverso il contributo delle rigidezze di telaio, l'aliquota di forza
sismica che compete ai vari telai.
Infine la valutazione dei taglianti è stata eseguita sommando sui vari piani il contributo
sismico corrispondente.
Le procedure descritte nel presente capitolo riguardano la più semplice analisi sismica
concepita dal D.M. 14-01-08. Essa è caratterizzata da forti semplificazioni di base:
a) comportamento della struttura shear-type;
b) posizione baricentri approssimata;
c) analisi di tipo statico;
d) periodo di vibrazione della struttura valutato in maniera empirica;
e) modo di vibrazione lineare;
f) massa eccitata 85%.
L'ipotesi a) è utilizzabile generalmente in fase di predimensionamento, per strutture basse e
con rigidezza flessionale delle travi almeno 10 volte quella dei pilastri. Questa condizione in
generale si realizza maggiormente ai piani alti delle strutture. Tale ipotesi semplifica
enormemente la cinematica della struttura e quindi la risoluzione del problema.
L'ipotesi b) consiste nel distribuire uniformemente la massa dell'impalcato e fare affidamento
31
Capitolo 2 – Analisi statica (shear-type)
unicamente alle rigidezze dei pilastri. La prima consegue dal difficile calcolo legato alla
situazione reale, la seconda consegue dal comportamento shear-type attribuito al telaio.
L'ipotesi b) è anch'essa semplificatoria e consiste nell'applicare il valore massimo di progetto
della forzante sismica (che ha natura dinamica) in maniera statica. Tale ipotesi è lecito farla,
soprattutto in fase iniziale di progetto, poiché non è conveniente avviare un'analisi dinamica
vista la scarsa definizione delle caratteristiche della struttura.
L'ipotesi c) consiste nella valutazione del più basso periodo proprio di vibrazione della
struttura. Tale ipotesi è bene usarla per strutture semplici, il più regolari possibili. È una
relazione di prima approssimazione da utilizzare nelle fasi iniziali di progetto.
Le ipotesi d) ed e) sono congiunte. Si usano per strutture semplici e il più regolari possibili.
La struttura vibra con una deformata di tipo lineare, prossima a quella del primo modo di
vibrare, sulla quale viene eccitato il minimo normativo di massa.
Le precedenti considerazioni portano a risultati che si avvicinano a quelli veri ma che tali non
possono essere considerati, poiché sottostimano nettamente i reali taglianti agenti alla base dei
pilastri dei telai.
Inoltre vista la forte approssimazione che riguarda la posizione del baricentro delle masse, la
normativa si cautela attraverso l'eccentricità accidentale.
Paragonando i taglianti alla base delle due precedenti analisi si nota che questi hanno valori
notevolmente differenti tra i vari casi.
32
Capitolo 3 – Analisi statica (condensazione statica)
3.1. Introduzione
In questo capitolo verrà affrontata una tipologia di analisi sismica che assume forzante
sismica applicata in modo statico sulla struttura. Tale analisi è semplificata e ipotizza, su ogni
telaio, un solo grado di libertà alla traslazione per ogni piano e tutti i gradi di libertà
rotazionali. È una semplificazione quella di non considerare la totalità dei cinematismi della
struttura che, per alcune tipologie di strutture, come la presente, effettua errori di modesta
entità.
Questa analisi permette di ripartire attraverso le rigidezze dei telai la forzante sismica su ogni
telaio e livello.
È possibile eseguire un'analisi semplificata che non fa lavorare il telaio secondo il modello
shear-type ma che analizza la struttura nel suo insieme. Questa ipotizza:
- inestensibilità per tutti gli elementi ( EA ∞ );
- infinita rigidezza tagliante per tutti gli elementi ( GA ∞ ).
Ne consegue che ogni telaio potrà solo subire una traslazione orizzontale in corrispondenza
dei traversi e una rotazione in corrispondenza ogni nodo. I gradi di libertà vengono
schematizzati nella seguente fig. 3.1.
Trasportando queste ipotesi direttamente sulla struttura 3D si va verso una risoluzione
complicata del problema vista la grande quantità di incognite. Perciò si preferisce semplificare
il tutto lavorando sulla matrice di rigidezza condensata.
33
Capitolo 3 – Analisi statica (condensazione statica)
Vengono eliminati sui telai tutti i gradi di libertà sui quali forze e coppie non fanno lavoro.
Nel nostro caso tutti i contributi dovuti alle rotazioni vengono inglobati nei contributi dovuti
alle traslazioni.
Algebricamente il problema si scrive come:
dove:
- ktt è la matrice di rigidezza che lega forze a spostamenti del singolo telaio;
- ktφ è la matrice di rigidezza che lega forze a rotazioni del singolo telaio;
- kφt è la matrice di rigidezza che lega momenti a spostamenti del singolo telaio;
- kφφ è la matrice di rigidezza che lega momenti a rotazioni del singolo telaio;
- ut è il vettore delle traslazioni del singolo telaio;
- uφ è il vettore delle rotazioni del singolo telaio;
- Ft è l'aliquota di forzante sismica sul singolo telaio.
Sviluppando il sistema di equazioni si arriva alla scrittura della matrice condensata (kt).
Si ottengono tante matrici delle rigidezze condensate quanti sono i telai, nel nostro problema
si ottiene ktx1, ktx2, ktx3, kty1, kty2, kty3.
34
Capitolo 3 – Analisi statica (condensazione statica)
Come fatto nelle precedenti analisi, trasportando le stesse ipotesi dei telai sulla struttura 3D e
quindi sull'impalcato, ne consegue che esso è infinitamente rigido e potrà solo subire due
traslazioni e una rotazione secondo l'asse verticale.
Si conclude che ogni struttura 3D, sotto queste ipotesi, ha 3 gradi di libertà per impalcato.
La matrice di rigidezza dell'intera struttura si scrive come:
dove:
- ;
- ;
- ;
- ;
- ;
- ;
- ;
- ;
- ;
35
Capitolo 3 – Analisi statica (condensazione statica)
K
204751 -107893 7502 0.00 0.00 0.00 37874 -19363 1389
-107893 153267 -57414 0.00 0.00 0.00 -19363 26266 -9701
7502 -57414 50376 0.00 0.00 0.00 1389 -9701 8432
0.00 0.00 0.00 192144 -103423 10513 32165 -17262 2936
0.00 0.00 0.00 -103423 139466 -53083 -17262 19386 -6986
0.00 0.00 0.00 10513 -53083 43996 2936 -6986 4586
37874 -19363 1389 32165 -17262 2936 6537064 -3509996 232593
-19363 26266 -9701 -17262 19386 -6986 -3509996 5128635 -1955992
1389 -9701 8432 2936 -6986 4586 232593 -1955992 1737054
Tab. 3.3 – Matrice spaziale delle rigidezze [kN/m, kN, kNm]
in cui:
36
Capitolo 3 – Analisi statica (condensazione statica)
Dagli spostamenti degli impalcati si sono ottenuti i contributi di spostamento relativi ad ogni
telaio.
La ripartizione della forzante sismica è stata fatta, per sisma in direzione x e y, come:
In ultimo i taglianti alla base dei pilastri sono stati valutati come la sommatoria delle forzanti
sismiche superiori alla base considerata.
Per quanto precedentemente detto si è ottenuto:
Sisma lungo x
ηX1 ηX2 ηX3 ηY1 ηY2 ηY3
Liv. 3 0.0102 0.0099 0.0097 -0.0003 0.0000 0.0003
Liv. 2 0.0074 0.0072 0.0071 -0.0002 0.0000 0.0002
Liv. 1 0.0039 0.0038 0.0037 -0.0001 0.0000 0.0001
Sisma lungo y
ηX1 ηX2 ηX3 ηY1 ηY2 ηY3
Liv. 3 0.0001 0.0000 -0.0001 0.0122 0.0123 0.0125
Liv. 2 0.0001 0.0000 -0.0001 0.0088 0.0090 0.0091
Liv. 1 0.0001 0.0000 -0.0001 0.0044 0.0045 0.0046
Tab. 3.4 – Spostamenti telai ai vari livelli [m]
37
Capitolo 3 – Analisi statica (condensazione statica)
Sisma lungo x
X1 X2 X3 Y1 Y2 Y3
Liv. 3 39.27 38.03 37.64 -1.05 -0.03 1.07
Liv. 2 41.41 46.28 38.99 -1.11 0.00 1.11
Liv. 1 19.62 35.36 18.01 -0.89 -0.11 1.00
Sisma lungo y
X1 X2 X3 Y1 Y2 Y3
Liv. 3 0.32 -0.02 -0.30 49.81 14.52 50.61
Liv. 2 0.40 -0.03 -0.37 60.38 6.02 60.29
Liv. 1 1.27 -0.11 -1.16 7.26 55.43 10.30
Tab. 3.5 – Forzante sismica [kN]
Sisma lungo x
X1 X2 X3 Y1 Y2 Y3
Liv. 3 39.27 38.03 37.64 -1.05 -0.03 1.07
Liv. 2 80.68 84.31 76.63 -2.16 -0.03 2.19
Liv. 1 100.30 119.67 94.64 -3.05 -0.13 3.19
Sisma lungo y
X1 X2 X3 Y1 Y2 Y3
Liv. 3 0.32 -0.02 -0.30 49.81 14.52 50.61
Liv. 2 0.73 -0.06 -0.67 110.18 20.54 110.90
Liv. 1 2.00 -0.17 -1.83 117.44 75.97 121.20
Tab. 3.6 – Taglianti alla base [kN]
3.4. Commenti
38
Capitolo 3 – Analisi statica (condensazione statica)
rigidezze condensate dei telai, l'aliquota di forza sismica che gli compete.
Infine la valutazione dei taglianti è stata eseguita sommando sui vari piani il contributo
sismico corrispondente.
La procedura descritta nel presente capitolo riguarda una semplice e anche accurata analisi
sismica concepita dal D.M. 14-01-08. Essa è caratterizzata da alcune semplificazioni di base:
a) condensazione statica della matrice di rigidezza;
b) analisi di tipo statico;
c) periodo di vibrazione della struttura valutato in maniera empirica;
d) modo di vibrazione lineare;
e) massa eccitata 85%.
L'ipotesi a) è utilizzabile generalmente in fase di predimensionamento, per strutture che hanno
un notevole numero di gradi di libertà. Tale ipotesi semplifica nettamente la cinematica della
struttura e quindi la risoluzione del problema.
L'ipotesi b) è anch'essa semplificatoria e consiste nell'applicare il valore massimo di progetto
della forzante sismica (che ha natura dinamica) in maniera statica. Tale ipotesi è lecito farla,
soprattutto in fase iniziale di progetto, poiché non è conveniente avviare un'analisi dinamica
vista la scarsa definizione delle caratteristiche della struttura.
L'ipotesi c) consiste nella valutazione del più basso periodo proprio di vibrazione della
struttura. Tale ipotesi è bene usarla per strutture semplici, il più regolari possibili. È una
relazione di prima approssimazione da utilizzare nelle fasi iniziali di progetto.
Le ipotesi d) ed e) sono congiunte. Si usano per strutture semplici e il più regolari possibili.
La struttura vibra con una deformata di tipo lineare, prossima a quella del primo modo di
vibrare, sulla quale viene eccitato il minimo normativo di massa.
Le precedenti considerazioni portano a risultati che si avvicinano a quelli veri ma che tali non
possono essere considerati, poiché sottostimano i reali taglianti agenti alla base dei pilastri dei
telai.
Paragonando i taglianti alla base delle due precedenti analisi si nota che questi hanno valori
differenti tra loro. È semplice concludere che l'analisi per condensazione statica ha un grado
di accuratezza maggiore rispetto all'analisi statica con shear-type vista la minore
approssimazione del problema.
39
Capitolo 4 – Analisi dinamica (pseudo-dinamica)
4.1. Introduzione
In questo capitolo verrà affrontato un tipo di analisi che determina la massima forzante
sismica attraverso l'analisi dinamica e l'assume applicata in modo statico sulla struttura. È una
metodologia più accurata delle precedenti, che però non tiene in conto della alternata
applicazione dei versi della forzante sismica e conseguentemente dei cicli a fatica.
Tutte le grandezze statiche che questa analisi utilizza sono le medesime riportate nel capitolo
precedente.
Inizialmente si sono analizzati i modi di vibrare della struttura. Per far ciò si sono studiate le
vibrazioni libere non smorzate, cioè con forzante nulla applicata e assenza di dissipazione.
L'equazione alla base di questo problema dinamico è:
dove:
- M è la matrice delle masse della struttura spaziale;
- K è la matrice di rigidezza della struttura spaziale;
- u è il vettore spostamento della struttura spaziale.
È stato necessario costruire la matrice M. Per convenienza si è scelto di costruirla non
ponendo il sistema di riferimento di ogni impalcato nel baricentro delle rigidezze C, come
fatto precedentemente. Se valutata rispetto al baricentro delle masse G, la matrice delle masse
risulta diagonale poiché si azzerano i termini di momento statico rispetto G. I termini non
40
Capitolo 4 – Analisi dinamica (pseudo-dinamica)
nulli sono riferiti alle masse degli impalcati per le accelerazioni ü e v̈ e alle inerzie
rotazionali per le accelerazioni ̈ .
Ipotizzando una soluzione per u(t) ne uscirebbe un problema singolare, risolvibile solo
attraverso un problema agli autovalori.
Ponendo nullo il determinante della matrice tra parentesi si ottiene la soluzione non banale
dell'autovettore Φ e dell'autovalore ω.
Il primo è una matrice quadrata di dimensioni pari ai gradi di libertà della struttura. Essa
rappresenta i modi di vibrare della struttura, una forma adimensionalizzata della deformata.
Il secondo è un vettore di dimensioni pari ai gradi di libertà della struttura. Esso rappresenta la
frequenza circolare dei rispettivi modi di vibrare. Da questo è possibile ricavare il periodo di
vibrazione modale come:
2
T i=
i
Nelle seguenti tabelle vengono riportati i termini trattati finora:
M
115.65 0 0 0 0 0 0 0 0
0 109.98 0 0 0 0 0 0 0
0 0 68.72 0 0 0 0 0 0
0 0 0 115.65 0 0 0 0 0
0 0 0 0 109.98 0 0 0 0
0 0 0 0 0 68.72 0 0 0
0 0 0 0 0 0 1946.75 0 0
0 0 0 0 0 0 0 1851.36 0
0 0 0 0 0 0 0 0 1156.73
Tab. 4.1 – Matrice spaziale delle masse [kN·s2/m]
41
Capitolo 4 – Analisi dinamica (pseudo-dinamica)
Φ
0.0001 0.0319 0.0012 0.0004 -0.0570 -0.0017 -0.0019 0.0660 0.0031
0.0003 0.0608 0.0026 0.0002 -0.0384 -0.0023 0.0017 -0.0625 -0.0028
0.0003 0.0831 0.0037 -0.0005 0.0818 0.0038 -0.0008 0.0304 0.0014
0.0307 -0.0001 0.0001 0.0568 0.0004 0.0002 -0.0668 -0.0020 0.0025
0.0610 -0.0003 0.0010 0.0395 0.0003 -0.0023 0.0616 0.0018 -0.0022
0.0836 -0.0004 0.0018 -0.0813 -0.0006 0.0024 -0.0306 -0.0009 0.0011
-0.0002 -0.0004 0.0081 -0.0003 0.0006 -0.0143 0.0002 -0.0005 0.0156
-0.0003 -0.0006 0.0149 -0.0001 0.0003 -0.0087 -0.0008 0.0008 -0.0156
-0.0003 -0.0008 0.0200 0.0005 -0.0009 0.0200 0.0007 -0.0006 0.0079
Tab.4.2 – Autovettori analisi modale
ω
132.87 163.80 324.53 908.51 1047.25 2094.31 2524.35 2678.75 5222.40
Tab. 4.3 – Autovalori analisi modale [rad/s]
Ti
0.545 0.491 0.349 0.208 0.194 0.137 0.125 0.121 0.087
Tab. 4.4 – Periodo di vibrazione modale [s]
Siccome non tutti i modi sono fondamentali ad un'analisi dinamica, attraverso la valutazione
delle masse efficaci nelle due direzioni, è possibile studiare i modi che eccitano
maggiormente la struttura.
42
Capitolo 4 – Analisi dinamica (pseudo-dinamica)
mEX
0.005 258.689 0.458 0.001 26.993 0.037 0.007 8.138 0.019
0.00 % 87.89 % 0.16 % 0.00 % 9.17 % 0.01 % 0.00 % 2.76 % 0.01 %
mEY
256.462 0.006 0.064 28.449 0.002 0.004 9.335 0.008 0.015
87.13 % 0.00 % 0.02 % 9.76 % 0.00 % 0.00 % 3.17 % 0.00 % 0.01 %
2
Tab. 4.5 – Masse efficaci modali [kN·s /m]
In seguito sono state analizzate le forzanti sismiche associate ai modi di vibrare della struttura.
Per far ciò si sono studiate le vibrazioni forzate, cioè con forzante non nulla applicata.
L'equazione alla base di questo problema dinamico è:
dove:
- ü g t è l'accelerazione al suolo del sisma;
- i è il vettore di trascinamento nelle due direzioni.
La soluzione di questa equazione del moto è:
dove:
Sai
- Sdi è lo pseudo-spostamento spettrale modale, determinabile come Sdi = ;
2i
43
Capitolo 4 – Analisi dinamica (pseudo-dinamica)
7,00
6,00
5,00
4,00
Sa [m/s2]
elastico
inelastico
3,00
2,00
1,00
0,00
0,00 0,50 1,00 1,50 2,00 2,50 3,00 3,50 4,00
To [s]
Sisma lungo x
0.069 16.084 0.677 0.034 -5.195 -0.193 -0.082 2.853 0.136
Sisma lungo y
16.014 -0.079 0.254 5.334 0.040 -0.065 -3.055 -0.090 0.124
Tab. 4.6 – Coefficiente di partecipazione modale
44
Capitolo 4 – Analisi dinamica (pseudo-dinamica)
Sa
1.1498 1.2560 1.2560 1.2560 1.2560 1.5610 1.6911 1.7299 2.0961
Sd
0.0087 0.0077 0.0039 0.0014 0.0012 0.0007 0.0007 0.0006 0.0004
Tab. 4.7 – Pseudo-accelerazione e spostamento spettrale modale [m/s2, m]
1 2 3 4 5 6 7 8 9
8.45E-08 0.0039 3.21E-06 1.81E-08 0.0004 2.43E-07 1.02E-07 0.0001 1.67E-07
1.54E-07 0.0075 6.74E-06 1.15E-08 0.0002 3.37E-07 -9.33E-08 -0.0001 -1.56E-07
2.06E-07 0.0102 9.62E-06 -2.49E-08 -0.0005 -5.45E-07 4.33E-08 0.0001 7.55E-08
1.83E-05 -1.71E-05 3.88E-07 2.71E-06 -2.50E-06 -3.40E-08 3.66E-06 -3.63E-06 1.37E-07
3.63E-05 -3.74E-05 2.62E-06 1.88E-06 -1.98E-06 3.32E-07 -3.38E-06 3.35E-06 -1.18E-07
4.97E-05 -5.31E-05 4.84E-06 -3.88E-06 3.75E-06 -3.38E-07 1.68E-06 -1.67E-06 5.76E-08
-1.00E-07 -4.32E-05 2.11E-05 -1.38E-08 -3.71E-06 2.05E-06 -1.36E-08 -1.01E-06 8.54E-07
-1.62E-07 -7.82E-05 3.90E-05 -4.68E-09 -2.09E-06 1.24E-06 4.44E-08 1.47E-06 -8.50E-07
-2.03E-07 -1.05E-04 5.24E-05 2.53E-08 5.47E-06 -2.87E-06 -3.94E-08 -1.03E-06 4.33E-07
Tab. 4.8 – Spostamento modale impalcato per sisma in x [m]
1 2 3 4 5 6 7 8 9
1.97E-05 -1.93E-05 1.20E-06 2.80E-06 -2.74E-06 8.17E-08 3.80E-06 -3.86E-06 1.52E-07
3.58E-05 -3.68E-05 2.53E-06 1.77E-06 -1.84E-06 1.13E-07 -3.48E-06 3.65E-06 -1.41E-07
4.80E-05 -5.03E-05 3.61E-06 -3.85E-06 3.92E-06 -1.83E-07 1.62E-06 -1.78E-06 6.85E-08
0.0043 8.39E-08 1.46E-07 0.0004 1.93E-08 -1.14E-08 0.0001 1.15E-07 1.24E-07
0.0085 1.83E-07 9.82E-07 0.0003 1.52E-08 1.12E-07 -0.0001 -1.06E-07 -1.07E-07
0.0116 2.61E-07 1.82E-06 -0.0006 -2.89E-08 -1.14E-07 0.0001 5.31E-08 5.23E-08
-2.33E-05 2.12E-07 7.92E-06 -2.14E-06 2.86E-08 6.90E-07 -5.09E-07 3.19E-08 7.75E-07
-3.77E-05 3.84E-07 1.46E-05 -7.23E-07 1.61E-08 4.19E-07 1.66E-06 -4.67E-08 -7.71E-07
-4.74E-05 5.14E-07 1.96E-05 3.91E-06 -4.21E-08 -9.66E-07 -1.47E-06 3.25E-08 3.93E-07
Tab. 4.9 – Spostamento modale impalcato per sisma in y [m]
Come fatto nel par. 2.5 è conveniente calcolare i contributi di spostamento modali relativi ad
ogni telaio.
La ripartizione della forzante sismica modale è stata fatta, per sisma in direzione x e y, come:
45
Capitolo 4 – Analisi dinamica (pseudo-dinamica)
In ultimo i taglianti modali alla base dei pilastri sono stati valutati come la sommatoria delle
forzanti sismiche superiori alla base considerata.
I taglianti totali da considerare sono stati ottenuti come composizione di quelli modali
principali, secondo la relazione:
T xj= Txj T xj T yi = Tyi T yi
2 2 52 i 12 42
Modo 2°
ηX1 ηX2 ηX3 ηY1 ηY2 ηY3
Liv. 3 0.0107 0.0102 0.0098 -0.0006 -0.0001 0.0005
Liv. 2 0.0078 0.0075 0.0071 -0.0005 -0.0001 0.0004
Liv. 1 0.0041 0.0039 0.0037 -0.0003 0.0000 0.0002
Modo 5°
ηX1 ηX2 ηX3 ηY1 ηY2 ηY3
Liv. 3 -0.0005 -0.0005 -0.0005 0.0000 0.0000 0.0000
Liv. 2 0.0002 0.0002 0.0002 0.0000 0.0000 0.0000
Liv. 1 0.0004 0.0004 0.0003 0.0000 0.0000 0.0000
Tab. 4.10 - Spostamenti modali telai ai vari livelli per sisma in x [m]
Modo 1°
ηX1 ηX2 ηX3 ηY1 ηY2 ηY3
Liv. 3 0.0003 0.0000 -0.0002 0.0113 0.0116 0.0119
Liv. 2 0.0002 0.0000 -0.0001 0.0082 0.0084 0.0087
Liv. 1 0.0001 0.0000 -0.0001 0.0041 0.0042 0.0044
Modo 4°
ηX1 ηX2 ηX3 ηY1 ηY2 ηY3
Liv. 3 0.0000 0.0000 0.0000 -0.0006 -0.0006 -0.0006
Liv. 2 0.0000 0.0000 0.0000 0.0003 0.0003 0.0003
Liv. 1 0.0000 0.0000 0.0000 0.0004 0.0004 0.0004
Tab. 4.11 – Spostamenti modali telai ai vari livelli per sisma in y [m]
46
Capitolo 4 – Analisi dinamica (pseudo-dinamica)
Modo 2°
X1 X2 X3 Y1 Y2 Y3
Liv. 3 40.28 38.09 36.95 -2.40 -0.13 1.93
Liv. 2 45.22 49.10 40.66 -2.83 -0.05 2.21
Liv. 1 20.67 36.18 17.75 -1.60 -0.48 1.76
Modo 5°
X1 X2 X3 Y1 Y2 Y3
Liv. 3 -12.80 -12.18 -11.70 0.80 0.08 -0.60
Liv. 2 9.82 8.64 9.09 -0.55 -0.06 0.37
Liv. 1 13.94 16.57 12.53 -0.94 -0.13 0.77
Tab. 4.12 – Forzante sismica modale per sisma in x [kN]
Modo 1°
X1 X2 X3 Y1 Y2 Y3
Liv. 3 0.77 0.10 -0.43 45.60 13.10 47.14
Liv. 2 0.99 0.13 -0.59 57.83 6.90 58.86
Liv. 1 1.55 0.02 -1.26 5.53 50.96 8.97
Modo 4°
X1 X2 X3 Y1 Y2 Y3
Liv. 3 -0.43 -0.05 0.24 -14.73 -7.10 -15.59
Liv. 2 0.23 0.04 -0.10 11.74 5.21 12.17
Liv. 1 0.57 0.07 -0.35 15.16 12.63 16.25
Tab. 4.13 – Forzante sismica modale per sisma in y [kN]
Modo 2°
X1 X2 X3 Y1 Y2 Y3
Liv. 3 40.28 38.09 36.95 -2.40 -0.13 1.93
Liv. 2 85.50 87.19 77.61 -5.23 -0.17 4.14
Liv. 1 106.17 123.37 95.36 -6.83 -0.65 5.90
Modo 5°
X1 X2 X3 Y1 Y2 Y3
Liv. 3 -12.80 -12.18 -11.70 0.80 0.08 -0.60
Liv. 2 -2.99 -3.54 -2.61 0.25 0.02 -0.23
Liv. 1 10.95 13.03 9.91 -0.69 -0.11 0.54
Tab. 4.14 – Taglianti modali alla base per sisma in x [kN]
47
Capitolo 4 – Analisi dinamica (pseudo-dinamica)
Modo 1°
X1 X2 X3 Y1 Y2 Y3
Liv. 3 0.77 0.10 -0.43 45.60 13.10 47.14
Liv. 2 1.76 0.23 -1.01 103.42 19.99 106.00
Liv. 1 3.30 0.24 -2.27 108.95 70.95 114.97
Modo 4°
X1 X2 X3 Y1 Y2 Y3
Liv. 3 -0.43 -0.05 0.24 -14.73 -7.10 -15.59
Liv. 2 -0.19 -0.01 0.14 -2.99 -1.89 -3.42
Liv. 1 0.38 0.06 -0.21 12.16 10.74 12.83
Tab. 4.15 – Taglianti modali alla base per sisma in y [kN]
Sisma lungo x
X1 X2 X3 Y1 Y2 Y3
Liv. 3 42.27 39.99 38.76 2.53 0.15 2.03
Liv. 2 85.55 87.26 77.65 5.23 0.17 4.15
Liv. 1 106.73 124.06 95.88 6.86 0.66 5.92
Sisma lungo y
X1 X2 X3 Y1 Y2 Y3
Liv. 3 0.88 0.12 0.49 47.92 14.90 49.65
Liv. 2 1.77 0.23 1.02 103.46 20.08 106.06
Liv. 1 3.33 0.25 2.28 109.63 71.76 115.68
Tab. 4.16 – Taglianti alla base [kN]
4.4. Commenti
48
Capitolo 4 – Analisi dinamica (pseudo-dinamica)
49