Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio...

85
1 Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del Lazio Meridionale Anno accademico 2014-2015

Transcript of Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio...

Page 1: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

1

Appunti di

MECCANICA DELLE STRUTTURE

Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del Lazio Meridionale

Anno accademico 2014-2015

Page 2: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

2

1 Richiami di Meccanica del continuo

1.1 Analisi della deformazione

Si consideri un mezzo continuo Ω che nel tempo cambi configurazione. Detta 0C la

configurazione del corpo Ω al tempo iniziale del moto 0t t= , sia C la configurazione di

Ω al generico istante 0t t> . In Figura 1-1 è riportato schematicamente il cambiamento di

configurazione del corpo Ω .

Figura 1-1: Cambiamento di configurazione del corpo Ω .

Si indicano nel seguito con:

x uy vz w

ξηζ

= = =

x y u (1.1)

il vettore posizione del generico punto materiale di Ω al tempo 0t , il vettore posizione

dello stesso punto materiale al tempo t ed il vettore spostamento del punto materiale, tale che:

= −u y x (1.2)

Nel seguito si utilizza la notazione di Voigt, per cui le matrici a due indici sono organizzate come vettori. In tale spirito, si introduce il vettore di deformazione infinitesima

T

x y z xy yz xzε ε ε γ γ γ=ε , definito dalla relazione:

u = y - x

C , t

P

Ω

o o

O

o

x y C , t

P

Page 3: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

3

=ε Su (1.3)

con

,

,

,

, ,

, ,

, ,

0 00 00 0

00

0

x

y

z

y x

z y

z x

• • •

= • • • • • •

S (1.4)

avendo indicato con la virgola la derivata parziale.

Le componenti x y zε ε ε del vettore di deformazione infinitesima rappresentano le

dilatazioni lineari lungo gli assi di riferimento. Le componenti xy yz xzγ γ γ

rappresentano gli scorrimenti angolari tra le direzioni.

I valori delle deformazioni principali e le direzioni principali si determinano risolvendo il problema di autovalori ed autovettori:

( )

( )

( )

1 1 02 2

1 1 02 21 1 02 2

x x xy y xz z

xy x y y yz z

xz x yz y z z

n n n

n n n

n n n

ε λ γ γ

γ ε λ γ

γ γ ε λ

− + + =

+ − + =

+ + − =

(1.5)

dove T

x y zn n n=n è il versore della direzione principale. Per ammettere soluzione

diversa dalla banale, il sistema di equazioni (1.5) non deve avere rango massimo, ovvero λ deve soddisfare la seguente equazione caratteristica:

3 21 2 3 0A A Aλ λ λ− + + = (1.6)

dove 1A , 2A e 3A sono invarianti della deformazione.

1.2 Analisi della tensione

Si suppone che sul corpo possano agire due tipi di forze: p forze superficiali (o di

contatto) e b forze di massa (o di volume, o azioni a distanza).

Il vettore t è la tensione nel punto P agente sul generico piano π di normale n . Vale il seguente teorema di azione e reazione o di reciprocità:

( ) ( )= − −t n t n (1.7)

Page 4: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

4

La dipendenza della tensione dalla normale è fornita dal teorema del tetraedro di Cauchy:

( )x x xy y xz z

xy x y y yz z

xz x yz y z z

n n nn n nn n n

σ τ ττ σ ττ τ σ

+ + = + + + +

t n (1.8)

dove

T

x y z xy yz xzσ σ σ τ τ τ=σ (1.9)

è il vettore delle tensioni. Introducendo la matrice dei coseni direttori N , la formula (1.8) assume la forma:

0 0 0

0 0 00 0 0

x y z

y x z

z y x

n n nn n n

n n n

= =

t Nσ N (1.10)

Indicando con T

x y zb b b=b il vettore delle forze per unità di volume nel solido,

valgono le seguenti equazioni indefinite di equilibrio:

, , ,

, , ,

, , ,

000

x x xy y xz z x

xy x y y yz z y

xz x yz y z z z

bbb

σ τ τ

τ σ τ

τ τ σ

+ + + =

+ + + =

+ + + =

(1.11)

ovvero, in forma compatta:

T + =S σ b 0 (1.12)

Le direzioni e le tensioni principali si determinano risolvendo il problema di autovalori ed autovettori:

( )( )

( )

0

0

0

x x xy y xz z

xy x y y yz z

xz x yz y z z

n n n

n n n

n n n

σ λ τ τ

τ σ λ τ

τ τ σ λ

− + + =

+ − + =

+ + − =

(1.13)

Per ammettere soluzione diversa dalla banale, il sistema di equazioni (1.13) non deve avere rango massimo, ovvero λ deve soddisfare la seguente equazione caratteristica:

3 21 2 3 0J J Jλ λ λ− + + = (1.14)

dove 1J , 2J e 3J sono invarianti della tensione.

Page 5: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

5

1.3 Principio dei lavori virtuali

Si consideri un corpo Ω con frontiera f u∂Ω = ∂ Ω ∪ ∂ Ω . Siano assegnate le forze di

volume b in Ω , le forze di superficie p su f∂ Ω e gli spostamenti u su u∂ Ω .

Si considera un campo di tensioni in equilibrio con le forze attive e reattive agenti sul corpo; ovvero si assume che siano soddisfatte le equazioni (1.11) e le seguenti equazioni al contorno:

su x x xy y xz z x

xy x y y yz z y f

xz x yz y z z z

n n n pn n n pn n n p

σ τ τ

τ σ τ

τ τ σ

+ + =

+ + = ∂ Ω

+ + =

(1.15)

con x y zp p p le componenti del vettore p .

Si assegni poi un campo di spostamenti u congruente con un campo di deformazioni ε , ovvero si suppone che le deformazioni soddisfino l'equazione di congruenza (1.3); inoltre sia assume ˆ=u u su u∂ Ω .

Sotto tali ipotesi, si dimostra che il lavoro virtuale esterno veL è uguale a quello interno viL

per qualunque insieme di spostamenti virtuali (infinitesimi) compatibili con la continuità del corpo:

ve viL L= (1.16)

con

( ) ( ) ( )

ˆ

ˆ ˆ ˆ

ˆ ˆ ˆ

f u

f u

f u

T T Tve

x x x

y y y

z z z

x y z x y z x y z

L dV dA dA

b p tu v w b dV u v w p dA u v w t dA

b p t

b u b v b w dV p u p v p w dA t u t v t w dA

Ω ∂ Ω ∂ Ω

Ω ∂ Ω ∂ Ω

Ω ∂ Ω ∂ Ω

= + +

= + +

= + + + + + + + +

∫ ∫ ∫

∫ ∫ ∫

∫ ∫ ∫

u b u p u t

(1.17)

( )

x

y

zTvi x y z xy yz xz

xy

yz

xz

x x y y z z xy xy yz yz xz xz

L dV dV

dV

σσσ

ε ε ε γ γ γτττ

σ ε σ ε σ ε τ γ τ γ τ γ

Ω Ω

Ω

= =

= + + + + +

∫ ∫

ε σ (1.18)

Page 6: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

6

Si evidenzia che non è stata supposta alcuna relazione tra il sistema di forze ed il campo di spostamenti. In altri termini, il sistema delle forze e quello degli spostamenti non sono legati fra loro da alcun nesso di tipo causa-effetto.

Vale inoltre quanto segue.

• Il soddisfacimento delle equazioni di equilibrio e delle equazioni di congruenza implicano la validità dell’equazione (1.16), ovvero del principio dei lavori virtuali.

• Se sono soddisfatte le equazioni di congruenza, l’equazione del principio dei lavori virtuali conduce all’equilibrio delle tensioni con le forze esterne.

• Se sono soddisfatte le equazioni di equilibrio, l’equazione del principio dei lavori virtuali conduce alla congruenza degli spostamenti con le deformazioni.

1.4 Legame costitutivo

Le equazioni di legame forniscono una relazione tra la deformazione e la tensione; tale relazione dipende specificamente dal materiale che costituisce il corpo Ω .

In particolare si considera un materiale elastico, ovvero un materiale capace di restituire tutta l’energia spesa per deformarlo. Si assume inoltre che il materiale abbia un comportamento lineare, ovvero che la relazione tensione – deformazione sia di tipo lineare:

=σ Cε (1.19)

dove C è la matrice (6×6 simmetrica) elastica del materiale.

Nel caso di materiali elastici lineari il tensore elastico C è definito da 21 costanti elastiche. Si vogliono ora classificare i materiali in funzione di particolari proprietà, chiamate simmetrie materiali.

1.5 Simmetrie materiali

Si consideri l’intorno ( )I P di un punto P . Sia x il vettore posizione che individua il

generico punto ( )Q I P∈ ed u sia lo spostamento del punto Q .

Nel seguente paragrafo si farà uso della notazione tensoriale, per semplicità degli sviluppi. A tale scopo si indicano con [ ]ε e [ ]σ il tensore delle deformazioni e quello delle tensioni.

La deformazione in Q , per l’equazione di congruenza (1.3) scritta in forma tensoriale,

vale:

[ ] 1 ( )2

T= ∇ + ∇ε u u (1.20)

Page 7: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

7

Si ruoti poi l’intorno ( )I P con centro in P , così da individuare un intorno ( )I P∗ ruotato

di ( )I P tramite un tensore di rotazione Q . Il generico punto ( )Q I P∗ ∗∈ sia individuato

dal vettore ∗x ruotato di x , i.e. ∗ =x Qx . Si evidenzia che si è ruotato solo l’intorno geometrico ( )I P , e non un insieme di punti materiali del corpo con le sue proprietà. Siano

i punti di ( )I P∗ soggetti a spostamenti ∗u ruotati di u , i.e. ∗ =u Qu . La deformazione

associata agli spostamenti ∗u vale:

( ) ( ) ( )k mij i ik k ik k ik ik km jm

j j j m j

u xu Q u Q u Q Q Qx x x x x

∗ ∗∗ ∗ ∗ ∗

∂ ∂∂ ∂ ∂∇ = = = = = ∇

∂ ∂ ∂ ∂ ∂u u (1.21)

da cui si ricava:

[ ] T∗ = ε Q ε Q (1.22)

Per l’equazione costitutiva, alla deformazione [ ]ε è associato il tensore delle tensioni

[ ] [ ]=σ ε , con tensore costitutivo del quarto ordine, mentre alla deformazione ruotata

corrisponde la tensione:

[ ][ ]T∗ ∗ = = σ ε Q ε Q (1.23)

D’altra parte, è possibile anche definire il ruotato del tensore di tensione [ ]σ , come:

[ ] T♦ = σ Q σ Q (1.24)

In generale la tensione associata alla deformazione ruotata ∗ σ è differente dal ruotato

della tensione ♦ σ .

Si dice che Q è una trasformazione di simmetria per il materiale se accade che la tensione corrispondente alla deformazione ruotata eguaglia la ruotata della tensione associata alla deformazione iniziale:

[ ] [ ]* [ ] [ ]T T♦ = ⇔ = ∀ σ σ εQ ε Q Q ε Q (1.25)

L’insieme di tutti i tensori Q che, per un dato materiale, soddisfano la proprietà appena enunciata prende il nome di gruppo di simmetria per il materiale e si indica con G .

Si evidenzia che G soddisfa le seguenti proprietà di gruppo:

1 T−∈ ⇒ = ∈

, ∈ ⇒ ∈Q Q QQ P QP

G GG G

(1.26)

Tali proprietà si verificano come segue:

1. 1 T−∈ ⇒ = ∈Q Q QG G

Page 8: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

8

Dalla relazione (1.25) si ricava: [ ] [ ][ ]TT =Q Q εQ ε Q

e quindi

[ ] [ ] [ ][ ]TT T T= =Q Q Q Q ε Q Q εQ ε Q

allora, posto [ ] [ ]T = ηQ ε Q , ovvero [ ] [ ]T=ε Q η Q , si ha:

[ ] [ ] [ ][ ] [ ]TTT T T= =Q η Q Q Q Q η Qη QQQQ

da cui si deduce che 1T −= ∈Q Q G .

2. , ∈ ⇒ ∈Q P QPG G Applicando la (1.25) si ha:

[ ] [ ][ ]TT=ε P PP ε P

ed anche

[ ] [ ][ ] [ ( ) ]T TT T= Q Q Q QP ε P P ε P

Combinando le ultime due relazioni si ottiene:

( ) ( ) ( ) ( )[ ] [ ( ) ] [ ]T TT T TT= =C ε P Q C Q Q QP QP C QP ε QP QPPεP

da cui si deduce che ∈QP G .

Il più piccolo sottogruppo di tensori ortogonali che sono trasformazioni di simmetria per

un materiale contiene almeno , −I I ed è relativo ai materiali anisotropi.

1.5.1 Materiale monoclino Si esamina ora il caso di materiali che ammettono un piano di simmetria: essi sono detti monoclini. Si assuma come piano di simmetria il piano individuato da due versori 1e e 2e

ortogonali ad un versore 3e . In tal caso, il gruppo di simmetria materiale è:

3

= ,− ,I I PG

dove

3 1 1 2 2 3 3= ⊗ + ⊗ − ⊗P e e e e e e

così che nel riferimento 1 2 3, ,e e e si ha:

3

1 0 00 1 00 0 1

= −

P

Si calcoli allora il tensore di deformazione ∗ε nel caso che 3=Q P :

Page 9: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

9

[ ]11 1112 13 12 13

3 312 1222 23 22 23

13 1323 33 23 33

1 0 0 1 0 00 1 0 0 1 00 0 1 0 0 1

Tε ε ε ε ε εε ε ε ε ε εε ε ε ε ε ε

− = = = − − − − −

ε P ε P (1.27)

Le componenti del tensore delle tensioni ∗σ sono allora:

11 1111 11 1122 22 1133 33 1123 23 1113 13 1112 12

22 2211 11 2222 22 2233 33 2223 23 2213 13 2212 12

33 3311 11 3322 22 3333 33 3323 23 3313 13 3312 12

23 2311 11 23

2 2 22 2 22 2 2

C C C C C CC C C C C CC C C C C CC C

σ ε ε ε ε ε ε

σ ε ε ε ε ε ε

σ ε ε ε ε ε ε

σ ε

= + + − − +

= + + − − +

= + + − − +

= + 22 22 2333 33 2323 23 2313 13 2312 12

13 1311 11 1322 22 1333 33 1323 23 1313 13 1312 12

12 1211 11 1222 22 1233 33 1223 23 1213 13 1212 12

2 2 22 2 22 2 2

C C C CC C C C C CC C C C C C

ε ε ε ε ε

σ ε ε ε ε ε ε

σ ε ε ε ε ε ε

+ − − +

= + + − − +

= + + − − +

(1.28)

D’altra parte, il ruotato tramite 3P del tensore delle tensioni σ è:

[ ]11 12 13

3 312 22 23

13 23 33

Tσ σ σσ σ σσ σ σ

− = = −

− −σ P σ P (1.29)

che per il legame costitutivo, in esplicito fornisce:

11 1111 11 1122 22 1133 33 1123 23 1113 13 1112 12

22 2211 11 2222 22 2233 33 2223 23 2213 13 2212 12

33 3311 11 3322 22 3333 33 3323 23 3313 13 3312 12

23 2311 11 2

2 2 22 2 22 2 2

C C C C C CC C C C C CC C C C C C

C C

σ ε ε ε ε ε ε

σ ε ε ε ε ε ε

σ ε ε ε ε ε ε

σ ε

= + + + + +

= + + + + +

= + + + + +

= − +( )( )

322 22 2333 33 2323 23 2313 13 2312 12

13 1311 11 1322 22 1333 33 1323 23 1313 13 1312 12

12 1211 11 1222 22 1233 33 1223 23 1213 13 1212 12

2 2 2

2 2 2

2 2 2

C C C C

C C C C C C

C C C C C C

ε ε ε ε ε

σ ε ε ε ε ε ε

σ ε ε ε ε ε ε

+ + + +

= − + + + + +

= + + + + +

(30)

Se 3P appartiene al gruppo di simmetrie materiali, deve verificarsi l’eguaglianza ∗ ♦ = σ σ , che in componenti fornisce:

1111 11 1122 22 1133 33 1123 23 1113 13 1112 12

1111 11 1122 22 1133 33 1123 23 1113 13 1112 12

2 2 22 2 2

C C C C C CC C C C C C

ε ε ε ε ε εε ε ε ε ε ε

+ + − − += + + + + +

2211 11 2222 22 2233 33 2223 23 2213 13 2212 12

2211 11 2222 22 2233 33 2223 23 2213 13 2212 12

2 2 22 2 2

C C C C C CC C C C C C

ε ε ε ε ε εε ε ε ε ε ε

+ + − − += + + + + +

3311 11 3322 22 3333 33 3323 23 3313 13 3312 12

3311 11 3322 22 3333 33 3323 23 3313 13 3312 12

2 2 22 2 2

C C C C C CC C C C C C

ε ε ε ε ε εε ε ε ε ε ε

+ + − − += + + + + +

2311 11 2322 22 2333 33 2323 23 2313 13 2312 12

2311 11 2322 22 2333 33 2323 23 2313 13 2312 12

2 2 22 2 2

C C C C C CC C C C C Cε ε ε ε ε ε

ε ε ε ε ε ε

+ + − − +

=− + + + + +

Page 10: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

10

1311 11 1322 22 1333 33 1323 23 1313 13 1312 12

1311 11 1322 22 1333 33 1323 23 1313 13 1312 12

2 2 22 2 2

C C C C C CC C C C C Cε ε ε ε ε ε

ε ε ε ε ε ε

+ + − − +

=− + + + + +

1211 11 1222 22 1233 33 1223 23 1213 13 1212 12

1211 11 1222 22 1233 33 1223 23 1213 13 1212 12

2 2 22 2 2

C C C C C CC C C C C C

ε ε ε ε ε εε ε ε ε ε ε

+ + − − += + + + + +

che possono verificarsi se e solo se:

1123 1113 2223 2213 3323 3313 1223 1213 0C C C C C C C C= = = = = = = = (1.31)

In definitiva, si osservi che per essere il materiale monoclino, devono essere nulle le 8 costanti elastiche specificate nelle relazioni (1.31), se ne deduce allora che per materiali monoclini le costanti elastiche necessarie alla definizione di C si riducono a 13.

1.5.2 Materiale triclinoclino o ortotropo Si esamini ora il caso in cui il materiale ammetta due piani di simmetria e si supponga che il gruppo di simmetria sia fornito da:

3 2

= ,− , ,I I P PG

dove

2 1 1 2 2 3 3= ⊗ − ⊗ + ⊗P e e e e e e

così che nel riferimento 1 2 3, ,e e e si ha:

2

1 0 00 1 00 0 1

= −

P

Per la seconda proprietà dei gruppi, si ha che se 2 ∈P G e 3 ∈P G anche 2 3 ∈P P G , ma

2 3 1 1 2 2 3 3 1 1 2 2 3 3

1 1 2 2 3 3 1

= ⊗ − ⊗ + ⊗ ⊗ + ⊗ − ⊗

= ⊗ − ⊗ − ⊗ =

P P e e e e e e e e e e e e

e e e e e e P

così che nel riferimento 1 2 3, ,e e e è:

1

1 0 00 1 00 0 1

= − −

P

Se ne deduce che se il materiale ammette due piani di simmetria ne ammette anche un terzo; in questo caso il materiale è detto triclino o ortotropo. Per i materiali triclini deve accadere che:

Page 11: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

11

1123 1113 2223 2213

3323 3313 1223 1213

1112 2212 3312 1323

0 C C C CC C C CC C C C

= = = == = = == = = =

(1.32)

ovvero il numero delle costanti elastiche necessarie alla definizione di C si riduce ulteriormente, passando da 13 a 9.

In notazione Voigt, l’inversa della matrice elastica per il materiale ortotropo è del tipo:

1

1 0 0 0

1 0 0 0

1 0 0 0

10 0 0 0 0

10 0 0 0 0

10 0 0 0 0

yx zx

x y z

xy zy

x y z

yzxz

x y z

xy

yz

xz

E E E

E E E

E E E

G

G

G

ν ν

ν ν

νν−

− −

− − − −

=

C (1.33)

dove xE , yE e zE sono i moduli elastici lungo le direzioni x , y e z ; ijν è il rapporto di

Poisson, definito come il rapporto tra la deformazione lungo la j − esima direzione quando

è applicata ina tensione lungo la i − esima direzione, con , , ,i j x y z= .

A causa della simmetria della matrice elastica C e quindi della sua inversa 1−C , deve accadere:

; ;yx xy zy yzzx xz

y x z x z yE E E E E Eν ν ν νν ν

= = = (1.34)

per cui le costanti indipendenti sono solo 9.

Imponendo ulteriori condizioni di simmetria il numero delle costanti elastiche indipendenti si riduce sempre più.

1.5.3 Materiale trasversalmente isotropo Per i materiali trasversalmente isotropi accade che i gruppo di simmetrie materiali contenga tutti i tensori di rotazione intorno ad un asse. In tal case le costanti elastiche indipendenti sono solo 5 e l’inversa della matrice elastica diventa:

Page 12: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

12

( )

( )

1

1 0 0 0

1 0 0 0

1 0 0 0

2 10 0 0 0 0

2 10 0 0 0 0

10 0 0 0 0

z

z

z

z

z z

z z z

xz

E E E

E E E

E E E

E

E

G

νν

νν

ν ν

ν

ν

− −

− − − −

= + +

C (1.35)

In questo caso il piano ortogonale all’asse di tutte le rotazioni che appartengono a G è il

piano x y− .

1.5.4 Materiale isotropo Il più ampio gruppo di simmetria possibile, cui corrisponde il minimo numero di costanti elastiche, è costituito dall’insieme di tutti i tensori di rotazione. In tal caso il materiale è detto isotropo, ed ha uguale risposta per ogni direzione.

Un materiale è isotropo se si comporta nello stesso modo in tutte le direzioni, ovvero le proprietà meccaniche del materiale non dipendono dalla direzione. In questo caso, la matrice elastica dipende da solo due costanti elastiche, E modulo di Young e ν rapporto di Poisson.

Tornando alla notazione di Voigt, la relazione tensione-deformazione si scrive come:

( )( )

( )

( )

( )

1 0 0 01 0 0 0

1 0 0 010 0 0 1 2 0 021 1 2

10 0 0 0 1 2 02

10 0 0 0 0 1 22

E

ν ν νν ν νν ν ν

νν ν

ν

ν

− − − −= + − −

C (1.36)

e, nella forma inversa,

Page 13: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

13

( )( )

( )

1

1 0 0 01 0 0 0

1 0 0 010 0 0 2 1 0 00 0 0 0 2 1 00 0 0 0 0 2 1

E

ν νν νν ν

νν

ν

− − − − − −

= = + +

+

A C (1.37)

1.6 Problema delle Equilibrio Elastico

Dato un corpo Ω con frontiera f u∂Ω = ∂ Ω ∪ ∂ Ω , il problema dell’equilibrio elastico

consiste in:

• assegnate o le forze di volume b in Ω , o le forze di superficie p su f∂ Ω ,

o gli spostamenti u su u∂ Ω ;

• determinare lo stato elastico definito da o il campo di spostamenti u , o il campo di deformazioni ε , o il campo di tensioni σ ,

che soddisfano le equazioni di campo:

in T

=+ = Ω

=

ε SuS σ b 0σ Cε

(1.38)

e le condizioni al contorno:

su

ˆ su f

u

= ∂ Ω

= ∂ Ω

Nσ pu u

(1.39)

1.7 Energia Potenziale Totale

Il problema dell’equilibrio elastico può essere affrontato e risolto ricorrendo ad un approccio basato sulla definizione dell’energia coinvolta nel sistema.

Per formulare in modo semplice il problema in forma energetica, si considera inizialmente un caso molto semplice, il cui sistema meccanico è costituito da una molla ed una massa.

Page 14: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

14

La molla di lunghezza L è vincolata ad un estremo e sull’estremo opposto è presente una massa m che da luogo ad una forza F m g= , come illustrato in Figura 1-2.

Figura 1-2: Molla vincolata ad un estremo e caricata da una forza F .

Per effetto della forza applica F , nella molla nasce uno sforzo normale N ; d’altra parte, per la presenza dello sforzo normale N si ha un allungamento elastico della molla e, di conseguenza, si ha uno spostamento u dell’estremo su cui è applicata la forza. Le equazioni che governano il problema sono le seguenti:

congruenza

0 equilibriolegame costitutivo

uL

N FN k

ε

ε

=

− ==

(1.40)

dove k denota l’elasticità (per unità di lunghezza) della molla.

Dopo semplici sostituzioni, le equazioni (1.40) forniscono:

F K u= (1.41)

dove /K k L= è la rigidezza (globale) della molla. Risolvendo la (1.41) si ottiene:

FuK

= (1.42)

Energia interna

Tenuto conto della terza delle equazioni (1.40), la densità di energia interna immagazzinata dalla molla vale:

uF

configurazione iniziale

configurazione finale

Page 15: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

15

( ) 2

0 0

12

N k kε ε

ϕ ε δε ε δε ε= = =∫ ∫ (1.43)

ovvero si ricava come area sottesa dalla curva che definisce la relazione tra la sollecitazione interna N e la deformazione ε . L’energia immagazzinata in tutta la molla è:

( ) ( ) 2 2

0 0

1 12 2

L Ldx k dx k Lε ϕ ε ε εΦ = = =∫ ∫ (1.44)

Tenendo conto delle (1.40)1 e (1.40)3, la (1.44) fornisce:

21 12 2

uk u K uL

Φ = = (1.45)

Potenziale delle forze

Si definisce poi energia potenziale delle forze applicate la quantità:

F uΛ = − (1.46)

Energia potenziale totale

L’energia potenziale totale del sistema è ottenuta come la somma dell’energia potenziale interna della molla (1.45) e dell’energia potenziale delle forze (1.46):

21( )2

u K u F uΠ = Φ + Λ = − (1.47)

Si evidenzia che Π è funzione dello spostamento dell’estremo libero u della molla; in particolare, Π è una funzione quadratica dello spostamento.

In Figura 1-3 è schematicamente illustrato il grafico dell’energia potenziale totale.

Figura 1-3: Grafico dell’energia potenziale totale.

( )uΠ

u/F K

Page 16: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

16

Il teorema di stazionarietà dell’energia potenziale totale assicura che: condizione necessaria e sufficiente affinché lo spostamento u che soddisfi l’equazione di equilibrio (1.41), è che l’energia potenziale totale (1.47) attinga un punto di stazionarietà proprio in u . Infatti:

0 K u Fu

∂Π= = −

∂ (1.48)

ovvero si ottiene l’equazione (1.41).

Il teorema di stazionarietà dell’energia potenziale totale vale non solo per il semplice sistema considerato, costituito da una molla caricata da una forza, ma per ogni sistema elastico. Infatti, si può dimostrare che la soluzione di un problema dell’equilibrio elastico si può determinare imponendo la condizione di stazionarietà dell’energia potenziale totale. Nel seguito si introdurrà l’energia potenziale per specifici modelli strutturali.

Si considera ora il problema di un corpo rigido connesso ad un piano rigido tramite due molle parallele di differente rigidezza, come schematicamente illustrato in Figura 1-4.

Figura 1-4: Corpo rigido vincolato ad un piano tramite due molle elastiche.

La deformata del corpo rigido è governata dai due parametri cinematici u e ϕ ; gli

spostamenti dei punti di connessione delle molle con il corpo e lo spostamento del punto di applicazione della forza valgono rispettivamente:

1 1

2 2

F

u u du u du u e

ϕϕ

ϕ

= += −= +

(1.49)

configurazione iniziale

configurazione finale

F

2K1K

1d2d

eu

ϕ

Page 17: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

17

L’energia interna del sistema vale allora:

( ) ( )

2 21 1 2 2

2 21 1 2 2

1 12 21 12 2

K u K u

K u d K u dϕ ϕ

Φ = +

= + + − (1.50)

mentre, l’energia potenziale del carico vale:

( )FF u

F u eϕΛ = −

= − + (1.51)

In definitiva, l’energia potenziale totale è:

( ) ( ) ( )2 2

1 1 2 21 12 2

K u d K u d F u eϕ ϕ ϕ

Π = Φ + Λ

= + + − − + (1.52)

La condizione di stazionarietà dell’energia potenziale totale fornisce:

( ) ( )

( ) ( )

1 1 2 2

1 1 1 2 2 2

0

0

K u d K u d Fu

K u d d K u d d F e

ϕ ϕ

ϕ ϕϕ

∂Π= = + + − −

∂∂Π

= = + − − −∂

(1.53)

ovvero, in forma matriciale:

1 2 1 1 2 22 2

1 1 2 2 1 1 2 2

K K K d K d u FK d K d K d K d F eϕ

+ − = − +

(1.54)

La matrice a primo membro della (1.54) è detta matrice di rigidezza; essa è sempre simmetrica e definita positiva.

Page 18: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

18

2 Modellazione della trave: Problema assiale

2.1 Il modello di trave

Le equazioni che governano il problema assiale della trave sono di seguito riportate.

Spostamenti

Il campo degli spostamenti si scrive come segue:

( ) spostamento lungo w w z z= (2.1)

Deformazione

La deformazione si calcola come:

'wε = (2.2)

avendo indicato con l’apice la derivata rispetto a z .

Tensione

La tensione normale alla sezione retta si determina tramite l’equazione del legame costitutivo:

'E E wσ ε= = (2.3)

Sforzo normale

La presenza della tensione normale σ fa nascere uno sforzo normale:

A A

N dA E dA E Aσ ε ε= = =∫ ∫ (2.4)

essendo A l’area della sezione retta della trave.

Equilibrio

Dall’equilibrio del tratto infinitesimo di trave, rappresentato in Figura 2-1, si ricava l’equazione:

'N f= − (2.5)

Problema dell’equilibrio elastico

In definitiva, le equazioni che governano il problema assiale della trave sono:

' 'w N E A N fε ε= = = − (2.6)

e quindi

''E Aw f= − (2.7)

Page 19: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

19

Facendo riferimento allo schema rappresentato in Figura 2-2, il problema dell’equilibrio elastico è completato dalle seguenti condizioni al contorno:

0 0 0 in 00 inL L

N P zN P z L

− = =− = =

(2.8)

Figura 2-1: Equilibrio del tratto di trave infinitesimo.

2.2 Energia potenziale totale

Si consideri un filo di area dA e lunghezza dz . L’energia interna in questo caso si calcola come:

212

E dAdzφ ε= (2.9)

Poiché il tratto di trave di lunghezza dz può essere visto come l’insieme di tanti fili di area dA lunghezza dz , l’energia interna si può determinare come:

2 21 12 2A

A

E dAdz EA dzφ ε ε= =∫ (2.10)

ovvero

( )21 '2A EA w dzφ = (2.11)

Ne consegue che l’energia interna di una trave di lunghezza L vale:

( )2

0

1 '2

L

EA w dzΦ = ∫ (2.12)

Il potenziale dei carichi, con riferimento alla Figura 2-2, vale:

0 00

L

L Lf wdz P w P w

Λ = − + − ∫ (2.13)

N N N+ ∆

z∆

f

Page 20: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

20

Figura 2-2: Trave soggetta a carichi assiali distribuiti ed a forze sulle sezioni di estremità.

L’energia potenziale totale è:

( )20 0

0 0

1 '2

L L

L LEA w dz f wdz P w P w

Π = Φ + Λ = − + −

∫ ∫ (2.14)

Si evidenzia che Π dipende dalla funzione ( )w z ; l’integrale che definisce Π è un valore

scalare che dipende da una funzione. Tale speciale applicazione viene chiamata in matematica funzionale. Il concetto di funzionale può caratterizzarsi come la generalizzazione del concetto classico di funzione. In altri termini, un funzionale è una funzione a valore scalare nella quale la variabile indipendente è, a sua volta, una funzione.

Imporre la stazionarietà dell’energia, come già fatto nel caso del semplice sistema di una molla caricata che ha condotto alla semplice equazione (1.48), non risulta immediato nel caso in cui si debba imporre la stazionarietà di un funzionale come quello definito dall’equazione (2.14) per la trave.

La variazione prima di un funzionale è definita come:

( )0

d w ud α

δ α αα =

Π = Π + (2.15)

essendo u wα δ= una funzione che rappresenta una possibile variazione della funzione di

spostamento assiale della trave w .

Applicando la (2.15), la variazione prima di Π si determina come:

( ) ( )

( ) ( )( )

2

0 0

0 0 00

1 ' '2

L L

L L L

d EA w u dz f w u dzd

P w u P w uα

δ α α αα

α α=

Π = + − +

+ − +

∫ ∫ (2.16)

Effettuando la derivata rispetto ad α e calcolandone il valore per 0α = , la (2.16) fornisce:

( )0 00 0

' 'L L

L LEAw w dz f wdz P w P wδ δ δ δ δΠ = − − −∫ ∫ (2.17)

La condizione di stazionarietà dell’energia potenziale totale wδ∀ fornisce:

f LP0P

Page 21: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

21

( )0 00 0

0 ' 'L L

L LEAw w dz f wdz P w P w wδ δ δ δ δ δ= Π = − − − ∀∫ ∫ (2.18)

Integrando per parti, l’equazione (2.18) diventa:

( )

( ) ( ) ( )

0 0

0 0 0 0

0 0 00

0 '' '

'

'' ' '

L L

L L

L L

L

L L L

EAw wdz f wdz EAw w

EAw w P w P w

EAw f wdz EAw P w EAw P w w

δ δ δ

δ δ δ

δ δ δ δ

= − − +

− − −

= − + + − − − ∀

∫ ∫

(2.19)

Vista l’arbitrarietà di wδ , la (2.19) fornisce:

( )

0 0

'' 0 per 0,' 0 in 0' 0 inL L

EAw f z LEAw P zEAw P z L

+ = ∈− = =− = =

(2.20)

Le equazioni (2.20) coincidono le (2.6) e (2.7). Ne consegue che se w rappresenta un punto di stazionarietà dell’energia potenziale totale, ovvero se è soddisfatta la (2.18) e quindi le (2.20), allora sono soddisfatte le (2.6) e (2.7). Al contrario se w risolve le equazioni (2.6) e (2.7), allora sono soddisfatte le (2.20) e quindi w rappresenta un punto di stazionarietà dell’energia potenziale totale.

Page 22: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

22

3 Modellazione della trave: Trave inflessa di Eulero Bernoulli

3.1 Il modello di trave

Il modello di trave di Eulero Bernoulli (EB) è basato sull’ipotesi che la generica sezione retta a deformazione avvenuta resti piana ed ortogonale alla deformata della linea d’asse della trave, come schematicamente illustrato in Figura 3-1. Nel seguito viene ignorato il problema assiale, trattato nel paragrafo precedente.

Figura 3-1: Trave inflessa, ipotesi di Eulero Bernoulli.

Spostamenti

Il campo degli spostamenti si può scrivere come segue:

( ) spostamento lungo ( ) spostamento lungo

v v z yw y z zϕ

==

(3.1)

dove ( )zϕ rappresenta la rotazione della sezione retta, che per il modello di Eulero-

Bernoulli deve risultare:

'vϕ = − (3.2)

Deformazione

La deformazione si calcola come:

' 'w yε ϕ= = (3.3)

Ricordando la relazione (3.2), la deformazione (3.3) diventa:

''y vε = − (3.4)

inoltre, introducendo la curvatura come 'χ ϕ= , ovvero anche ''vχ = − , si ha:

v

Page 23: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

23

yε χ= (3.5)

Tensione

La tensione normale alla sezione retta si determina tramite l’equazione del legame costitutivo:

E E yσ ε χ= = (3.6)

essendo E il modulo elastico o modulo di Young.

Caratteristiche della sollecitazione

La presenza della tensione normale σ può far nascere sia sforzo normale che momento flettente:

2

A A

A A

N dA E y dA E S

M y dA E y dA E I

σ χ χ

σ χ χ

= = =

= = =

∫ ∫

∫ ∫ (3.7)

essendo A , S e I rispettivamente l’area, il momento statico e il momento d’inerzia della sezione retta della trave. Si evidenzia che scegliendo, come spesso accade, un sistema di riferimento baricentrico nella sezione retta, il momento statico è nullo e quindi lo sforzo normale è nullo: 0 0S N= ⇒ = .

Risolvendo la seconda delle (3.7) rispetto a χ e sostituendo il risultato nella (3.6), si

ottiene la formula di Navier:

M yI

σ = (3.8)

Equilibrio

Dall’equilibrio del tratto infinitesimo di trave, rappresentato schematicamente in Figura 3-2, si ricavano le equazioni:

'

'M TT q

== −

(3.9)

Figura 3-2: Equilibrio del tratto di trave infinitesimo.

T

T T+ ∆

M M M+ ∆

z∆

q

Page 24: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

24

Problema dell’equilibrio elastico

In definitiva, le equazioni che governano il problema della trave inflessa sono:

' '' '

'' ''

v T qM T

v M E I M q

ϕχ ϕ

χ χ

= − = −= =⇓ ⇓= − = = −

(3.10)

e quindi

''''E Iv q= (3.11)

Il problema dell’equilibrio elastico scritto nella forma (3.10) o (3.11) si completa con opportune condizioni al contorno.

3.2 Energia potenziale totale

Si consideri la trave rappresentata in Figura 3-3. Le equazioni che governano il problema sono le equazioni di campo (3.10), ovvero:

( )

''

per 0,'

'

v

M E I z LM TT q

ϕχ ϕ

χ

= −=

= ∈

== −

(3.12)

completate dalle condizioni al contorno:

0 0

0 0

0in 0

00

in0

L L

L L

T Fz

M mT F

z LM m

− ==

− =− =

=− =

(3.13)

Si consideri ora un filo di area dA lunghezza dz della trave. L’energia interna in questo caso si calcola come:

212

E dAdzφ ε= (3.14)

Poiché il tratto di trave di lunghezza dz può essere visto come l’insieme di tanti fili di area dA e lunghezza dz , tenuto conto della formula (3.5), l’energia interna si determina come:

2 2 2 21 1 12 2 2A

A A

E dAdz E y dAdz EI dzφ ε χ χ= = =∫ ∫ (3.15)

Page 25: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

25

ovvero

( )21 ''2A EI v dzφ = (3.16)

Ne consegue che l’energia interna di una trave inflessa di lunghezza L vale:

( )2

0

1 ''2

L

EI v dzΦ = ∫ (3.17)

Il potenziale dei carichi, con riferimento alla Figura 3-3, vale:

0 0 0 00

L

L L L Lq v dz F v m F v mϕ ϕ

Λ = − + + − − ∫ (3.18)

Figura 3-3: Trave soggetta a carico distribuito ed a carichi concentrati sulle sezioni di estremità.

Ne consegue che l’energia potenziale totale è:

( )20 0 0 0

0 0

1 ''2

L L

L L L LEI v dz q v dz F v m F v mϕ ϕ

Π = Φ + Λ = − + + − −

∫ ∫ (3.19)

Si evidenzia che Π dipende dalla funzione ( )v z ; l’integrale che definisce Π è un valore

scalare che dipende da una funzione, ovvero è un funzionale.

Come già vista nel paragrafo 2.2, la variazione prima di un funzionale è definita come:

( )0

d v ud α

δ α αα =

Π = Π + (3.20)

che nel caso specifico della trave EB fornisce:

q

m0mL

F0

FL

Page 26: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

26

( ) ( )

( ) ( ) ( ) ( )( )

2

0 0 0

0 0 0 0 0 00

1 '' ''2

' ' ' '

L L

L L L L L L

d EI v u dz q v u dzd

d F v u m v u F v u m v ud

α

α

δ α α αα

α α α α αα

=

=

Π = + − +

− + − + − + + +

∫ ∫(3.21)

essendo uα una funzione che rappresenta una possibile variazione della funzione di

inflessione v , come illustrato in Figura 3-4.

Figura 3-4: Variazione della funzione

La condizione di stazionarietà dell’energia potenziale totale diventa allora:

( )

( )

0 0 0 00 0

0 0 0 00 0

'' '' ' ' '

'' '' ' ' ' 0

L L

L L L L

L L

L L L L

EI v u dz qu dz F u m u F u m u

EI v v dz q v dz F v m v F v m v

δ α α

δ δ δ δ δ δ

Π = − − − − +

= − − − − + =

∫ ∫

∫ ∫ (3.22)

essendo v uδ α= .

La formula (3.22) si può riscrivere come:

( )0 0 00 0

'' '' ' ' 'L L

L L L L LEI v v dz q v dz F v m v F v m vδ δ δ δ δ δ= + − − +∫ ∫ (3.23)

che rappresenta la ben nota equazione del principio dei lavori virtuali. Integrando per parti, si ottiene:

( )

0 00

0 0 00

''' ' '' ' '' '

' ' '

L

L L

L

L L L L L

EI v v dz EI v v EI v v

q v dz F v m v F v m v

δ δ δ

δ δ δ δ δ

− + −

= + − − +

∫ (3.24)

ed integrando ancora:

Page 27: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

27

( )

0 0 0 00

0 0 00

'''' '' ' '' ' ''' '''

' ' '

L

L L L L

L

L L L L L

EI v v dz EI v v EI v v EI v v EI v v

q v dz F v m v F v m v

δ δ δ δ δ

δ δ δ δ δ

+ − − +

= + − − +

∫ (3.25)

ovvero, riorganizzando i termini:

( ) ( ) ( )

( ) ( )

0 0 00

0 0 0

'''' '' ' '' '

''' ''' 0

L

L L L

L L L

EI v q v dz EI v m v EI v m v

EI v F v EI v F v

δ δ δ

δ δ

− + + − +

− + + + =

∫ (3.26)

Dovendo la formula (3.26) valere per ogni vδ , deve accadere:

( )0 0

0

'''' 0 per 0,'' 0 in 0'' 0 in''' 0 in 0''' in

L L

L L

EI v q z LEI v m zEI v m z LEI v F zEI v F z L

− = ∈+ = =+ = =+ = =+ =

(3.27)

Le equazioni (3.27) coincidono le (3.12) e (3.13). Ne consegue che se v rappresenta un punto di stazionarietà dell’energia potenziale totale, ovvero se è soddisfatta la (3.23) e quindi le (3.27), allora sono soddisfatte le (3.12) e (3.13). Al contrario se v risolve le equazioni (3.12) e (3.13), allora sono soddisfatte le (3.27) e quindi v rappresenta un punto di stazionarietà dell’energia potenziale totale.

3.3 Metodo variazionale approssimato

Si assume una forma di rappresentazione per la funzione incognita v :

1

( ) ( )n

i ii

v z c zφ=

= ∑ (3.28)

essendo

( ) 1,2,..., funzioni approssimanti1,2,..., coefficienti incogniti

i

i

z i nc i nφ =

= In altre parole, la funzione incognita v si assume combinazione lineare di funzioni note. Le funzioni approssimanti sono scelte, mentre i coefficienti devono essere determinati.

L’energia potenziale totale di formula (3.19) diventa:

Page 28: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

28

2

1 10 0

0 01 1 1 1

1 ''2

( ) '( ) (0) '(0)

L Ln n

i i i ii i

n n n n

L i i L i i i i i ii i i i

EI c dz q c dz

F c L m c L F c m c

φ φ

φ φ φ φ

= =

= = = =

Π = − +

− − +

∑ ∑∫ ∫

∑ ∑ ∑ ∑ (3.29)

La condizione di stazionarietà dell’energia potenziale totale conduce alle n equazioni:

1

2

0

0

...

0n

c

c

c

∂Π=

∂∂Π

=∂

∂Π=

(3.30)

Il sistema di n equazioni (3.30) nelle n incognite 1c , 2c , .., nc risulta risolvibile se le

funzioni approssimanti iφ sono scelte linearmente indipendenti tra loro.

Poiché lo schema di trave analizzato illustrato in Figura 3-3 è labile, il sistema di equazioni (3.30) non risulta risolvibile, a meno che il sistema di forze non sia equilibrato; in tal caso, il sistema ammette soluzione a meno di moti rigidi. Per ottenere un’unica soluzione è necessario eliminare i moti rigidi introducendo opportuni vincoli. In tal caso, è possibile verificare che all’aumentare del numero delle funzioni scelte la soluzione approssimata tende alla soluzione esatta purché lo spazio delle funzioni approssimanti sia completo. Nel caso di polinomi approssimanti, un polinomio di grado p è completo se sono presenti tutti

i termini del polinomio con grado minore di p .

Per risolvere il problema di una trave non labile è necessario scegliere delle funzioni di approssimazione che soddisfano le condizioni di vincolo imposte sulla trave, cosa non sempre semplice da effettuare. Allo scopo di evitare la determinazione di funzioni di approssimazione che implicitamente soddisfino le condizioni di vincolo, nel paragrafo successivo viene presentato un metodo semplice ed efficace.

3.4 Metodo di penalizzazione (penalty)

La tipica condizione di vincolo può essere rilassata sostituendo i vincoli rigidi con vincoli cedevoli elasticamente. In Figura 3-5 è illustrato il caso di un appoggio. I vincoli elastici sono caratterizzati da una rigidezza che governa la relazione tra la reazione del vincolo e l’abbassamento. Incrementando il valore della rigidezza vk lo spostamento v che subirà la

sezione vincolata si riduce; per un valore della rigidezza che tende a infinito, si tende

Page 29: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

29

quindi al vincolo rigido, per il quale lo spostamento della sezione vincolata è nullo. La forza che nasce nel vincolo elastico vale:

vR k v= − (3.31)

Figura 3-5: Comportamento del vincolo elastico.

L’energia interna immagazzinata nel vincolo elastico vale:

21 12 2v vR v k vΦ = − = (3.32)

Nota che nell’espressione (3.32) compare il segno meno, in quanto R e v hanno versi opposti.

L’energia espressa dalla formula (3.32) deve essere introdotta nel calcolo dell’energia potenziale totale del sistema. In definitiva si ottiene:

21ˆ2v vk vΠ = Φ + Λ + Φ = Φ + Λ + (3.33)

La stazionarietà del nuovo funzionale Π definito dalla relazione (3.33) fornisce:

ˆvk v vδ δ δ δΠ = Φ + Λ + (3.34)

3.5 Moltiplicatori di Lagrange

Ricordando l’espressione (3.31), l’equazione di equilibrio (3.34) diventa:

ˆ R vδ δ δ δΠ = Φ + Λ − (3.35)

Tale equazione si può ricavare introducendo il funzionale:

L R vΠ = Φ + Λ − (3.36)

detto di Lagrange, dove R è il moltiplicatore di Lagrange del vincolo 0v = . Infatti, LΠ

dipende da R e da v , per cui la condizione di stazionarietà diventa:

v

R

vR

R

vk

Page 30: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

30

0

0v L

R L

R vv R

δ δ δ δδ δ

Π = Φ + Λ − =Π = − =

(3.37)

La seconda delle (3.37) rappresenta la condizione di vincolo 0v = .

3.6 Trave su suolo elastico

Un caso di notevole interesse ingegneristico è quello della trave appoggiata su un letto di molle distribuite (alla Winkler) indipendenti tra loro. Soggetto a inflessione positiva, il letto di molle reagisce in modo proporzionale al valore dell’abbassamento ( )v z in verso

opposto all’inflessione; in altre parole, la trave oltre ad essere soggetta al carico esterno assegnato q è soggetta anche alla reazione delle molle r K v= − , così che l’equazione

(3.11) diventa:

''''E I v q r= + (3.38)

ovvero

''''E I v K v q+ = (3.39)

Posto:

44KE I

β= (3.40)

l’equazione (3.39) diventa:

4'''' 4 qv vEI

β+ = (3.41)

La soluzione dell’omogena associata all’equazione (3.41) assume la forma:

( ) ( )( )( ) ( )( )

1 2

3 4

cos sin

cos sin

z

z

v e C z C z

e C z C z

β

β

β β

β β

−= +

+ + (3.42)

L’energia interna immagazzinata dal letto di molle deformato congruentemente con la trave dell’inflessione ( )v z , si calcola come:

2

0 0

1 12 2

L L

W r v dz K v dzΦ = − =∫ ∫ (3.43)

Per tenere quindi conto della presenza del letto di molle alla Winkler, nell’espressione dell’energia potenziale della trave (3.19) si deve introdurre anche il termine WΦ fornito dalla (3.43).

Page 31: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

31

4 Stabilità di travi inflesse Si consideri una trave soggetta ad un carico assiale uniformemente distribuito tale da indurre un andamento lineare dello sforzo normale N di compressione lungo l’asse della trave, come riportato in Figura 4-1.

Figura 4-1: schema di una trave soggetta a carico assiale uniformemente distribuito.

Le ipotesi su cui la trattazione presentata si fonda prevedono che la trave sia:

- perfettamente rettilinea, - soggetta a sforzo di compressione N perfettamente assiale, - omogenea, - sottoposta a tensioni entro il limite elastico.

La trattazione è svolta assumendo che la trave si trovi in una configurazione deformata, ma che quest'ultima sia sufficiente vicina ma non coincida con l'indeformata.

4.1 Approccio statico

Facendo riferimento all'elemento di trave di lunghezza infinitesima dz riportato in Figura 4-1, le equazioni di equilibrio sono:

- equilibrio alla traslazione verticale

' 0T T dT T= + → = (4.1)

L

N

dz

M+dM

N+dN

M T

f

T+dT

f

zQ0 QL

Page 32: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

32

- equilibrio alla rotazione

' 'M N dv Tdz M dM M Nv T+ + = + → − = (4.2)

combinando le equazioni (4.1) e (4.2), si ottiene:

'' '' ' ' 0M Nv N v− − = (4.3)

Nelle ipotesi che la curvatura della trave sia sufficientemente piccola e che quindi valga ancora la classica relazione tra curvatura e momento flettente:

''M EI EI vχ= = − (4.4)

l’equazione differenziale che descrive il problema è:

'''' '' ' ' 0EI v N v N v+ + = (4.5)

Dalla condizione di equilibrio in direzione dell’asse della trave, si ha inoltre:

0 'dN f dz N f− = → = (4.6)

Per 0f ≠ , l'equazione (4.5) è un’ equazione differenziale del quarto ordine a coefficienti

non costanti; per tale equazione, in generale, non è possibile determinare una soluzione analitica in forma chiusa.

Qualora si pone 0f = , lo sforzo di compressione N risulta essere costante lungo l’asse

della trave, e l’equazione differenziale (4.5) si particolarizza in:

'''' '' 0EI v N v+ = (4.7)

Tale equazione si differenzia dalla classica equazione della linea elastica di Eulero-Bernoulli per la presenza del termine ''N v responsabile dell’equilibrio nella

configurazione inflessa.

Ponendo:

2 NEI

α = (4.8)

l'equazione (4.7) si riscrive nella forma equivalente:

2'''' '' 0v vα+ = (4.9)

L’equazione (4.9) è un’ equazione differenziale del quarto ordine a coefficienti costanti,

che ammette una soluzione generale del tipo:

( ) sin( ) cos( )v z A z B z Cz Dα α= + + + (4.10)

Imponendo opportune condizioni al contorno è possibile scrivere un sistema omogeneo di 4 equazioni nelle 4 incognite A, B, C e D. Per ottenere una soluzione diversa dalla banale si

Page 33: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

33

impone che il determinate dei coefficienti delle incognite sia nullo. In tal modo si determina il valore più piccolo di α e, per la (4.8), il valore del carico critico.

4.1.1 La trave semplicemente appoggiata Facendo riferimento alla classica trave appoggiata e appoggiata, le equazioni che esprimono le condizioni al contorno forniscono:

2 2

(0) 0 0(0) 0 0( ) 0 sin( ) cos( ) 0( ) 0 sin( ) cos( ) 0

v B DM Bv L A L B L CL D

M L A L B Lα α

α α α α

= → + == → == → + + + == → + =

(4.11)

Tali condizioni forniscono in definitiva:

0, 0, 0, sin( ) 0B C D A Lα= = = = (4.12)

Se anche A fosse nulla, la soluzione sarebbe banale, ovvero l’equilibrio si avrebbe nella

configurazione indeformata.

Poiché si intende determinare la condizione di equilibrio nella configurazione deformata, si

deve porre:

sin( ) 0 nL L nLπα α π α= → = → = (4.13)

Tenendo conto della posizione (4.8), la terza delle relazioni (4.13) fornisce:

2

2crN EILπ

= (4.14)

avendo assunto 1n = , ossia di quel valore di N superato il quale si attiva il fenomeno dell’instabilità dell’elemento.

4.2 Approccio energetico

L'energia potenziale si scrive come somma dell'energia interna Φ e del potenziale dei carichi esterni Λ :

Π = Φ + Λ (4.15)

L'energia interna, dovuta all'inflessione della trave assume la forma (3.17):

( )2

0

1 ''2

L

EI v dzΦ = ∫ (4.16)

mentre l'energia dei carichi, con riferimento alla Figura 4-1, vale:

Page 34: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

34

0

L

L LQ w f wdzΛ = − − ∫ (4.17)

dove w rappresenta lo spostamento assiale della trave nella generica sezione e Lw è il

valore di w per z L= ; si evidenzia inoltre che lo spostamento è stato posto uguale a zero per 0z = , i.e. 0 0w = .

Figura 4-2: spostamento assiale dovuto all’inflessione.

Con riferimento allo schema illustrato in Figura 4-2 , lo spostamento assiale si calcola come:

( ) ( ) ( ) ( )22

0 0 0 0

1 11 cos '2 2

z z z z

w z d d d d v dζ ξ ϕ ζ ϕ ζ ζ= − = − ≈ ≈∫ ∫ ∫ ∫ (4.18)

In definitiva, l'energia potenziale totale assume la forma specifica:

( ) ( ) ( )2

2 2

0 0 0 0

1 1 1'' ' '2 2 2

L L L z

LEI v dz Q v dz f v d dzζ

Π = − −

∫ ∫ ∫ ∫ (4.19)

La stazionarietà dell'energia (4.19) fornisce la condizione di equilibrio della trave, ovvero conduce all’equazione (4.5).

4.3 Metodo variazionale approssimato

L'approccio all'energia potenziale totale può considerarsi un'ottima base per lo sviluppo di tecniche numeriche capaci di determinare soluzioni approssimate del carico critico. Tale considerazione è particolarmente importante quando l'equazione differenziale che governa il problema ha una soluzione complessa oppure non ammette soluzione analitica.

Assumendo una forma specifica per il campo di spostamenti da inflessione del tipo:

( ) ( ) ( )1 1 2 2( ) .... n nv z c z c z c zφ φ φ= + + + (4.20)

dw

Configurazione indeformata

Configurazione deformata

Page 35: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

35

con 1φ , 2φ , ..., nφ n funzioni assegnate e 1c , 2c , ..., nc n parametri incogniti, l'equazione

(4.19) assume la forma:

( )

( )

( )

2

1 1 2 20

21 1 2 2

0

21 1 2 2

0 0

1 '' '' ... ''2

1 '' '' ... ''2

1 '' '' ... ''2

L

n n

L

L n n

L z

n n

EI c c c dz

Q c c c dz

f c c c d dz

φ φ φ

φ φ φ

φ φ φ ζ

Π = + + +

− + + +

− + + +

∫ ∫

(4.21)

Le funzioni assegnate 1φ , 2φ , ..., nφ sono scelte in modo da soddisfare la compatibilità

cinematica dei vincoli introdotti nella trave. In alternativa si può utilizzare un metodo tipo penalty.

Imponendo la stazionarietà dell'energia (4.21), si ottiene il seguente sistema di equazioni algebriche:

( ) 0 , 1,2,...,ij L ij ij ji

K Q G f M c i j nc

∂Π= − − = =

∂ (4.22)

con

0 0 0 0

'' '' ' ' ' 'L L L z

ij i j ij i j ij i jK EI dz G dz M d dzφ φ φ φ φ φ ζ= = =∫ ∫ ∫ ∫ (4.23)

Si ottiene in definitiva un sistema omogeneo di n equazioni nelle n incognite 1c , 2c , ..., nc ;

per ottenere una soluzione diversa dalla banale si impone che il determinate dei coefficienti delle incognite sia nullo:

( )det 0LQ f− − =K G M (4.24)

in tal modo si determina il più piccolo valore di LQ che rappresenta il valore del carico

critico.

4.4 Trave su suolo elastico

Si supponga ora che la trave sia adagiata su un suolo elastico e si schematizzi il terreno come un mezzo continuo elastico, ossia come costituito da un letto di molle dotate tutte della medesima rigidezza.

Si pensi ad esempio alle rotaie ferroviarie, in cui un aumento di temperatura può provocare forze assiali di compressione così elevate da indurre fenomeni di svergolamento se la rigidezza del collegamento con il suolo non è sufficiente. Altro esempio è quello dei lunghi

Page 36: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

36

pali di fondazione portanti di punta. La reazione del terreno sarà, in virtù di tale posizione, proporzionale in ciascuna sezione retta della trave alla deformazione impressa.

Per determinare il carico critico in tale circostanza si adotta un criterio di natura energetica. Si indica con K il modulo di rigidezza della fondazione, ovvero la rigidezza per unità di lunghezza.

4.4.1 Soluzione approssimata tramite approccio energetico Si affronta il problema imponendo la stazionarietà dell'energia potenziale totale. In questo caso è necessario introdurre l’ulteriore termine nella formula (4.19) che definisce l'energia elastica delle molle trasversali alla trave. Assumendo il carico assiale distribuito nullo, cioè

0f = , l'energia potenziale diventa:

( ) ( )2

22

0 0 0

1'' '2 2 2

L L LEI Kv dz v dz N v dzΠ = + −∫ ∫ ∫ (4.25)

L’espressione della deformata della trave è approssimata dalla polinomiale:

2 3 41 2 3 1... n

nv c z c z c z c z−= + + + + (4.26)

Sostituendo l'espressione (4.26) nella formula (4.25), ed imponendo la stazionarietà dell'energia si ottiene:

( ) 0 , 1,2,...,ij ij ij jK K H N G c i j n+ − = = (4.27)

che rappresenta un sistema omogeneo di n equazioni nelle n incognite 1c , 2c , ..., nc ; per

ottenere una soluzione diversa dalla banale si impone che il determinate dei coefficienti delle incognite sia nullo:

( )det 0K N+ − =K H G (4.28)

in tal modo si determina il più piccolo di valore N che rappresenta il valore del carico critico.

Page 37: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

37

5 Modellazione della trave: Trave di Timoshenko

5.1 Il modello di trave

Il modello di trave di Timoshenko è basato sull’ipotesi che la generica sezione retta a deformazione avvenuta resti piana ma non necessariamente ortogonale alla deformata della linea d’asse della trave. In altri termini, rispetto al modello di trave di Eulero Bernoulli, la rotazione della sezione retta ( )zϕ non risulta definita come l’opposto della derivata prima

dell’inflessione, ma è una funzione indipendente dalla ( )v z .

Il modello di trave di Timoshenko (T) viene spesso anche denominato modello con deformazione a taglio del primo ordine.

Spostamenti

Il campo degli spostamenti si può scrivere come segue:

( ) spostamento lungo ( ) spostamento lungo

v v z yw y z zϕ

==

(5.1)

Deformazione

La deformazione si calcola come:

'

'

w yzv w vz y

ε ϕ

γ ϕ

∂= =

∂∂ ∂

= + = +∂ ∂

(5.2)

avendo indicato con l’apice la derivata rispetto a z . Si evidenzia che, al contrario del modello EB, il valore dello scorrimento angolare γ può essere diverso da zero.

Introducendo la curvatura come 'χ ϕ= , si ha:

'

w yzv w vz y

ε χ

γ ϕ

∂= =

∂∂ ∂

= + = +∂ ∂

(5.3)

Tensioni

La tensione normale e tangenziale nella sezione retta della trave si determinano tramite le equazioni del legame costitutivo:

Page 38: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

38

( )'

'E E y E y

G G vσ ε χ ϕτ γ ϕ

= = =

= = + (5.4)

essendo E il modulo elastico o modulo di Young e G il modulo a taglio.

Caratteristiche della sollecitazione

La presenza della tensione normale σ e della tensione tangenziale τ fa insorgere momento flettente e taglio:

( ) ( )

2 ' '

' 'A A

sA A

M y dA E y dA E I

T dA G v dA GA v

σ ϕ ϕ

τ ϕ ϕ

= = =

= = + = +

∫ ∫

∫ ∫ (5.5)

essendo sA ed I rispettivamente l’area a taglio (concetto meglio precisato nel seguito) ed

il momento d’inerzia della sezione retta della trave. L’area a taglio viene spesso definita come il prodotto dell’area geometrica A moltiplicata per il fattore di deformazione a taglio κ , sA Aκ= .

Risolvendo la prima delle (5.5) rispetto a 'ϕ e sostituendo il risultato nella prima delle

(5.4), si ottiene la formula di Navier:

M yI

σ = (5.6)

Equilibrio

Dall’equilibrio del tratto infinitesimo di trave, rappresentato in Figura 3-2, si ricavano le equazioni:

'

'M TT q

== −

(5.7)

come nel caso della trave EB.

Problema dell’equilibrio elastico

In definitiva, le equazioni che governano il problema della trave inflessa sono:

' '' '

sv T G A T qM E I M T

γ ϕ γχ ϕ χ

= + = = −= = =

(5.8)

e quindi

( )

( )'' '

'' 's

s

E I G A v

G A v q

ϕ ϕ

ϕ

= +

+ = − (5.9)

ovvero, derivando la prima delle (5.9) e tenendo conto della seconda delle (5.9):

Page 39: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

39

( )''''' 's

E I qE I G A v

ϕϕ ϕ

= −

= + (5.10)

Facendo riferimento allo schema di trave illustrato in Figura 3-3, il problema dell’equilibrio elastico della trave inflessa è completato dalle condizioni al contorno:

0 0

0 0

0in 0

00

in0

L L

L L

T Fz

M mT F

z LM m

− ==

− =

− ==

− =

(5.11)

5.2 Determinazione del fattore di correzione a taglio

Il fattore di deformazione a taglio per il modello di trave T nasce dalla seguente esigenza.

Le tensioni tangenziali nella generica sezione retta della trave si possono determinare tramite l’equazione del legame costitutivo, la seconda delle (5.4), oppure tramite l’equazione di equilibrio:

0y zτ σ∂ ∂

+ =∂ ∂

(5.12)

Assumendo inizialmente per semplicità che la sezione retta della trave sia rettangolare di dimensioni b h× , l’equazione di equilibrio (5.12) si può integrare, fornendo:

/2

y

hd

zστ η

∂= −

∂∫ (5.13)

Ricordando la prima delle (5.4), si ottiene:

2 2

2

/2/2

1'' '' ''2 2 4

yy

hh

hE d E E yητ η ϕ η ϕ ϕ−

= − = − = − −

∫ (5.14)

In definitiva, si giunge al seguente risultato:

( )

22

' legame costitutivo

1'' equilibrio2 4

G v

hE y

τ ϕ

τ ϕ

= +

= −

(5.15)

In Figura 5-1 è illustrato il profilo delle tensioni tangenziali ottenuti utilizzando le due formule riportate per il calcolo (5.15).

Page 40: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

40

Figura 5-1: Profilo delle tensioni tangenziali nella sezione retta della trave.

Si evidenzia che le tensioni da legame costitutivo sono costanti nella sezione e, per la simmetria delle tensioni tangenziali, fanno emergere tensioni tangenziali non nulle sull’estradosso ed intradosso della trave. I valori delle tensioni tangenziali che emergono sono funzioni del taglio nella sezione e quindi si determinano come parte della soluzione del problema della trave. In realtà, il valore delle tensioni all’estradosso ed intradosso della trave devono essere prescritte, ovvero rappresentano un dato del problema e non possono quindi risultare come parte della soluzione del problema. Spesso, le tensioni tangenziali all’estradosso ed intradosso della trave sono assunte nulle. Quindi, le tensioni che emergono sono inaccettabili.

Al contrario, le tensioni che si determinano tramite le equazioni di equilibrio hanno un andamento parabolico e sono nulle all’estradosso ed intradosso della trave. Inoltre, è possibile dimostrare che i valori che si determinano tramite la seconda delle equazioni (5.15) approssima molto bene la soluzione esatta del taglio. Si evidenzia la seconda delle equazioni (5.15) rappresenta un caso particolare della ben nota formula di Jourawski.

E’ allora possibile derivare due differenti espressioni della tensione tangenziale in funzione della caratteristica della sollecitazione tagliante. Infatti, assumendo la tensione tangenziale costante in tutta la sezione della trave, si ottiene:

Tbh

τ = (5.16)

D’altra parte si calcoli la risultante delle tensioni tangenziali, ovvero la caratteristica di sollecitazione di taglio, integrando le tensioni tangenziali ottenute tramite la seconda delle equazioni (5.15):

Travetensioni tangenziali da legame costitutivo

tensioni tangenziali da equilibrio

Emergono tensioni tangenziali inaccettabili

Distribuzione ottimale

Page 41: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

41

22

2/2 /2 2

/2 /2

/22 3

/2

3 33

1 ''2 4

1 ''2 4

1 ''2 4 3

1 1'' ''2 4 12 12

A A

h b

h b

h

h

hT dA E y dA

hE y dA

h yE b y

h hE b E bh

τ ϕ

ϕ

ϕ

ϕ ϕ

− −

= = −

= −

= −

= − =

∫ ∫

∫ ∫ (5.17)

da cui si ricava:

3

12'' TEbh

ϕ = (5.18)

Sostituendo la formula (5.18) nella (5.14), si ottiene:

( )2 23

3 42

T h ybh

τ = − (5.19)

Si determina quindi l’energia associata alla deformazione provocata dalla presenza delle tensioni tangenziali nella generica sezione retta della trave. L’espressione dell’energia è:

12 A

dAψ τ γ= ∫ (5.20)

Tenuto conto che si è pervenuti a due differenti espressioni della tensione tangenziale e, di conseguenza, dello scorrimento angolare, è possibile determinare due differenti formule per l’energia di deformazione:

2 2

22

1 12 2 2lc

A A

T TdA dAG G A GA

ψ τκ κ κ

= = =∫ ∫ (5.21)

( )2

2 2 23

1 1 3 42 2 2eq

A A

TdA h y dAG G bh

ψ τ = = − ∫ ∫ (5.22)

dove i pedici lc e eq denotano l’energia calcolata utilizzando l’espressione delle tensioni tangenziali da legame costitutivo e da equilibrio. Calcolando l’integrale al secondo membro dell’equazione (5.22), si ottiene:

Page 42: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

42

( )2

4 4 2 22 6

25

2 6

2 25 5

2 6 2 6

2 25

2 6

1 9 16 82 4

1 9 16 812 4 5*16 3*41 9 1 2 1 9 15 3 101

2 4 5 3 2 4 151 9 8 3

2 4 15 5

eqA

T h y h y dAG b h

T bhG b h

T Tbh bhG b h G b h

T TbhG b h G bh

ψ = + −

= + −

+ − = + − =

= =

(5.23)

Imponendo l’eguaglianza dell’energia calcolata con l’espressione delle tensioni tangenziali da legame costitutivo e da equilibrio, si ha:

2 23 5

2 5 6lc eqT TGA G bh

ψ ψ κκ

= ⇒ = ⇒ = (5.24)

Nel caso più generale di trave con sezione generica, per la determinazione del fattore di correzione a taglio è necessario determinare il profilo delle tensioni tangenziali tramite l’equazione di equilibrio.

Con riferimento alla Figura 5-2, si considera l’equilibrio della parte di trave definita tra le sezioni z e z z+ ∆ e dal piano passante per i punti A e B della sezione retta della trave. Su tale parte di trave agiscono tensioni nomali sulla sezione retta e tensioni tangenziali sia sulla sezione retta che, per la simmetria delle tensioni tangenziali, anche sulla faccia che divide la trave in senso longitudinale.

Figura 5-2: Tratto di trave soggetto a flessione e taglio.

Per l’equilibrio lungo l’asse della trave, cioè in direzione z , si può scrivere:

( ) ( )*

0A

z z z dA b zσ σ τ+ ∆ − + ∆ = ∫ (5.25)

∆z

M(z) M (z+ ∆z)A*

τ ττ A B

Page 43: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

43

dove *A è l’area della parte di sezione al di sopra della corda AB e b è la dimensione della corda. Dividendo tutti i termini della formula (5.25) per z∆ e facendo il limite per

z∆ che tende a zero, si ottiene:

*

' 0A

dA bσ τ+ =∫ (5.26)

Ricordando la formula di Navier (5.6), l’equazione (5.26) diventa:

* *

'0A A

M Ty dA b y dA bI I

τ τ= + = +∫ ∫ (5.27)

essendo, per la prima delle (5.7), il taglio la derivata del momento flettente. Risolvendo l’equazione (5.27), la tensione tangenziale vale:

*T SI b

τ = − (5.28)

nota come formula di Jourawksi per il taglio; nell’equazione (5.28), *S rappresenta il

momento statico della parte di sezione retta *A . Si evidenzia che sia la corda b che il

momento statico *S sono funzioni di y .

L’espressione dell’energia diventa allora:

22 *

212eq

A

T S dAG I b

ψ

=

∫ (5.29)

Imponendo l’eguaglianza dell’energia calcolata con l’espressione delle tensioni tangenziali da legame costitutivo e da equilibrio:

22 2 * 2

22 *2 2lc eqA

A

T T S IdAGA G I b SA dA

b

ψ ψ κκ

= ⇒ = ⇒ =

∫∫

(5.30)

Sezione circolare

Nel caso della sezione circolare, si utilizza per il calcolo del fattore di correzione a taglio un sistema di coordinate polari. Con riferimento alla Figura 5-3, per la corda individuata dall’anomalia θ , si ha:

cossin

x Ry R

θθ

==

(5.31)

Page 44: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

44

Figura 5-3: Sezione circolare, calcolo del fattore di correzione a taglio.

La lunghezza della corda risulta allora:

2 2 cosb x R θ= = (5.32)

L’area ed il momento d’inerzia della sezione circolare si calcolano rispettivamente come:

( )/2 2 2

/22 cos

R

RA

A dA bdy R d Rπ

πθ θ π

− −= = = =∫ ∫ ∫ (5.33)

( ) ( )/2 2 22 4

/2

12 sin cos4A

I y dA R R d Rπ

πθ θ θ π

−= = =∫ ∫ (5.34)

Analogamente, il momento statico dell’area *A vale:

( ) ( )*

2 3*

/2

22 sin cos cos3A

S y dA R R d Rϕ

πθ θ θ ϕ

−= = = −∫ ∫ (5.35)

Si calcola anche la quantità:

( )( )

( )

23

2* /2 2

/2

/2 66 6

/2

2 cos3 2 cos

2 cos

2 5cos9 72

A

RS dA R db R

R d R

π

π

π

π

ϕϕ ϕ

ϕ

ϕ ϕ π

=

= =

∫ ∫

(5.36)

Sostituendo le espressioni (5.33), (5.34) e (5.36) nella formula (5.30), si ottiene:

*Aθ

ϕ

Page 45: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

45

24

2 6

194

10 109 16

R

R R

πκ

π π

= =

(5.37)

5.3 Soluzione analitica

Il sistema costituito dalle due equazioni differenziali (5.10) nelle due funzioni incognite ( )v z e ( )zϕ si risolve come segue. Dalla prima delle (5.10) si ottiene:

1

21 2

3 21 2 3

'''

''

'2

16 2

EI qq z c

EIq z c z cEIq z c z c z cEI

ϕ

ϕ

ϕ

ϕ

= −

= − +

= − + +

= − + + +

(5.38)

sostituendo poi l’espressione della ϕ e della ''ϕ nella seconda delle (5.10) si ha:

3 21 1 2 3

3 21 2 1 3

' ''

16 2

16 2

s

s

s

EIvGAEI q qz c z c z c z c

GA EI EI

q EI q EIz c z c z c cEI GA EI GA

ϕ ϕ

χ

= −

= − + − − + + +

= − − + + −

(5.39)

che integrata fornisce:

4 3 21 2 1 3 4

1 124 6 2 s s

q EI q EIv z c z c z c c z cEI GA EI GA

= − − + + − +

(5.40)

Le costanti di integrazione si determinano imponendo opportune condizioni al contorno.

5.4 Energia potenziale totale

Per il calcolo dell’energia interna, nel caso del modello di trave di Timoshenko, si deve tener conto non solo dell’energia che dipende dall’effetto flessionale, ma anche di quella dovuta alla presenza di tensioni tangenziali che lavorano per gli associati scorrimenti angolari.

Ne consegue che l’energia interna della trave di lunghezza L vale:

Page 46: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

46

( ) ( )2 2

0 0

1 1' '2 2

L L

sEI dz GA v dzϕ ϕΦ = + +∫ ∫ (5.41)

Il potenziale dei carichi è ancora espresso dalla formula (3.18), per cui l’energia potenziale totale del sistema rappresentato in Figura 3-3, vale:

( ) ( ) ( )2 20 0 0 0

0 0 0

1 1, ' '2 2

L L L

s L L L Lv EI dz GA v dz q v dz F v m F v mϕ ϕ ϕ ϕ ϕ

Π = + + − + + − −

∫ ∫ ∫ (5.42)

La condizione di stazionarietà dell’energia potenziale totale fornisce allora:

( ) ( )

( )

0 0 0

0 0 00

' ' ' ' '

' 0

L L L

s s

L

L L L L L

EI dz GA v dz GA v v dz

q v dz F v m F v m

δ ϕ δϕ ϕ δϕ ϕ δ

δ δ δϕ δ δϕ

Π = + + + +

− − + − − =

∫ ∫ ∫

∫ (5.43)

ovvero si ottengono le due equazioni di stazionarietà:

( ) ( )0 00 0

' ' ' 0L L

s L LEI dz GA v dz m mϕδ ϕ δϕ ϕ δϕ δϕ δϕΠ = + + − − =∫ ∫ (5.44)

( ) ( )0 00 0

' ' 0L L

v s L LGA v v dz q v dz F v F vδ ϕ δ δ δ δΠ = + − − − =∫ ∫ (5.45)

Integrando per parti l’espressione (5.44), si ottiene:

( ) ( ) ( )0 0 00

0 '' ' ' 'L

s L L LEI GA v dz EI m EI mϕ ϕ δϕ ϕ δϕ ϕ δϕ= − + + + − − − ∫ (5.46)

che, dovendo essere valida per ogni possibile δϕ , si localizza in:

( ) ( )

0 0

'' ' 0 per 0,' in 0' in

s

L L

EI GA v z LEI m zEI m z L

ϕ ϕϕϕ

− + = ∈= == =

(5.47)

Analogamente, integrando per parti l’espressione (5.45), si ottiene:

( )

( ) ( )0

0 0 0 0

0 '' '

' '

L

s

s L L L L s

GA v q v dz

GA v F v GA v F v

ϕ δ

ϕ δ ϕ δ

= − + +

+ + − − + −

∫ (5.48)

che, dovendo essere valida per ogni possibile vδ , si localizza in:

( ) ( )( )( )

0 0 0

'' ' 0 per 0,' in 0' in

s

s

s L L L

GA v q z LGA v F zGA v F z L

ϕϕϕ

+ + = ∈+ = =+ = =

(5.49)

Page 47: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

47

Le equazioni (5.47) e (5.49) equivalgono alle già determinate (5.8)-(5.10) e (5.11).

5.5 Metodo variazionale approssimato

Si assumono le seguenti forme di rappresentazione per le funzioni incognita v e ϕ :

1

1

( ) ( )

( ) ( )

n

i ii

m

i ii

v z c z

z g z

φ

ϕ ψ

=

=

=

=

∑ (5.50)

essendo

( )

( )( ) 1,2,..., funzioni approssimanti per

1,2,..., coefficienti incogniti per i

i

z i n v zc i n v zφ =

=

( )

( )( ) 1,2,..., funzioni approssimanti per

1,2,..., coefficienti incogniti per i

i

z i m zg i m zψ ϕ

ϕ==

L’energia potenziale totale di formula (5.42) diventa:

2 2

1 1 1 10 0 0

0 01 1 1 1

1 1' '2 2

( ) ( ) (0) (0)

L L Lm n m n

i i s i i i i i ii i i i

n m n m

L i i L i i i i i ii i i i

EI g dz GA c g dz q c dz

F c L m g L F c m g

ψ φ ψ φ

φ ψ φ ψ

= = = =

= = = =

Π = + + − +

+ − −

∑ ∑ ∑ ∑∫ ∫ ∫

∑ ∑ ∑ ∑ (5.51)

La condizione di stazionarietà dell’energia potenziale totale conduce alle n m+ equazioni:

1

2

1

2

0

0

...

0

0

0

...

0

n

m

c

c

c

g

g

g

∂Π=

∂∂Π

=∂

∂Π=

∂∂Π

=∂∂Π

=∂

∂Π=

(5.52)

Page 48: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

48

Il sistema di n m+ equazioni (5.52) nelle n m+ incognite 1c , 2c , .., nc , 1g , 2g , .., mg

risulta risolvibile se le funzioni approssimanti iφ sono scelte linearmente indipendenti tra

loro e le funzioni approssimanti iψ sono scelte linearmente indipendenti tra loro.

Page 49: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

49

6 Problemi bidimensionali: lastra e piastra

Gli elementi strutturali lastra e piastra sono oggetti tridimensionali in cui una dimensione è significativamente più piccola delle altre due. Si indica con Ω il corpo tridimensionale di

spessore costante h e piano medio M , così che ,2 2h h Ω = − ×

M .

Nel seguito si indica con:

• lastra, la struttura intesa bidimensionale sulla quale il carico agisce nel piano stesso della struttura;

• piastra, la struttura intesa bidimensionale sulla quale il carico agisce ortogonalmente al piano stesso della struttura.

6.1 Lastra

6.1.1 Cinematica Il campo di spostamenti è descritto da vettore:

( )( )

,,

u x yv x y

=

u (6.1)

Le deformazioni non nulle valgono:

(6.2)

che organizzate in un vettore forniscono:

0

0x

y

xy

x

y

y x

εεγ

∂ ∂

∂ = = = ∂ ∂ ∂ ∂ ∂

ε Su S (6.3)

6.1.2 Legame costitutivo Tenuto conto delle (1.19) e (1.37), il legame costitutivo 3D isotropo assume la forma:

, , , ,x x y y xy y xu v u vε ε γ= = = +

Page 50: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

50

( ) ( ) ( )1 1 1

1 1 1

x x y z y y x z z z x y

xy xy yz yz xz xz

E E E

G G G

ε σ ν σ ν σ ε σ ν σ ν σ ε σ ν σ ν σ

γ τ γ τ γ τ

= − − = − − = − −

= = = (6.4)

Si esamina inizialmente il caso di stato piano di tensione, definito dalle seguenti condizioni:

0z yz xzσ τ τ= = = (6.5)

In tal caso le equazioni di legame costitutivo (6.4) si semplificano in:

( ) ( ) ( )1 1 1

1 0 0

x x y y y x z x y

xy xy yz xz

E E E

G

ε σ ν σ ε σ ν σ ε ν σ ν σ

γ τ γ γ

= − = − = − −

= = = (6.6)

Invertendo le (6.6), si ottiene:

( )

( )

2

2

1

1

x x y

y y x

xy xy

E

E

G

σ ε ν εν

σ ε ν εν

τ γ

= +−

= +−

=

(6.7)

ovvero in forma matriciale:

( )

2 2

2 2

1 01 1

1 01 1

10 02 1

E

νν ν

νν ν

ν

− − = = − − +

σ Dε D (6.8)

Lo stato piano di deformazione è definito dalle seguenti condizioni:

0z yz xzε γ γ= = = (6.9)

In tal caso le equazioni di legame costitutivo (6.4) si semplificano in:

Page 51: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

51

( )

( )

( )

1

1

10

1

10

10

x x y z

y y x z

z x y

xy xy

yz

xz

E

E

E

G

G

G

ε σ ν σ ν σ

ε σ ν σ ν σ

σ ν σ ν σ

γ τ

τ

τ

= − −

= − −

= − −

=

=

=

(6.10)

Invertendo le(6.10), si ottiene:

( ) ( ) ( )( )

( ) ( ) ( )( )

11 1 2

11 1 2

x x y

y y x

xy xy

E

E

G

σ ν ε ν εν ν

σ ν ε ν εν ν

τ γ

= − ++ −

= − ++ −

=

(6.11)

ovvero in forma matriciale:

1 01 2 1 2

1 01 1 2 1 2

10 02

E

ν νν ν

ν νν ν ν

− − −

− = = + − −

σ Dε D (6.12)

Inoltre si ottiene anche:

( )z x yσ ν σ σ= + (6.13)

6.1.3 Equazioni di equilibrio Si introducono le caratteristiche membranali della sollecitazione come:

/2

/2/2

/2/2

/2

h

x x xhh

y y yhh

xy xy xyh

N dz h

N dz h

N dz h

σ σ

σ σ

τ τ

= =

= =

= =

∫∫∫

(6.14)

essendo xσ , yσ e xyτ funzioni solo di x e y .

Le equazioni di equilibrio si ricavano con riferimento alla Figura 6-1, sono:

Page 52: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

52

( ) ( )( ) ( )

0

0x x x xy xy xy x

xy xy xy y y y y

N N N y N N N x b x y

N N N y N N N x b x y

+ ∆ − ∆ + + ∆ − ∆ + ∆ ∆ =

+ ∆ − ∆ + + ∆ − ∆ + ∆ ∆ = (6.15)

dove xb e yb sono le componenti di un vettore b che rappresenta un carico distribuito

agente sul piano medio della lastra.

Figura 6-1: Equilibrio dell’area elementare per la lastra.

Dividendo tutti i membri delle equazioni (6.15) per x y∆ ∆ , nel limite per 0x∆ → e

0y∆ → si ottiene:

, ,

, ,

00

x x xy y x

xy x y y y

N N bN N b

+ + =

+ + = (6.16)

6.1.4 Energia potenziale totale e formulazione variazionale approssimata L’energia potenziale totale per una lastra vale:

( ) ( ) ( )Π = Φ + Λu u u (6.17)

dove l’energia interna ( )Φ u è:

x

Nx Nx+∆Nx

Ny

Ny+∆Ny

Nxy

Nxy

Nxy+∆Nxy

Nxy+∆Nxy∆y

∆x

y

Page 53: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

53

( )

( )

1 12 212

T T

T

dV h dA

h dA

ΩΦ = =

=

∫ ∫

u ε σ ε Dε

S u DS u

M

M

(6.18)

e il potenziale dei carichi ( )Λ u è:

( ) T TdA ds∂

Λ = − −∫ ∫u u b u pM M

(6.19)

dove p è il vettore che rappresenta un carico distribuito agente sul contorno del piano

medio della lastra.

La formulazione variazionale approssimata del problema della lastra si ottiene considerando una forma di rappresentazione del campo di spostamenti, come:

( ) ( )

( ) ( )

1

1

, ,

, ,

x

y

n

x i x iin

y i y ii

u x y c x y

v x y c x y

φ

φ

=

=

=

=

∑ (6.20)

essendo

( , ) 1,2,..., funzioni approssimanti di

1,2,..., coefficienti incogniti per l'approssimazione di xx i

xx i

x y i n u

c i n u

φ =

=

( , ) 1,2,..., funzioni approssimanti di

1,2,..., coefficienti incogniti per l'approssimazione di yy i

yy i

x y i n v

c i n v

φ =

=

Si evidenzia che il numero xn e yn delle funzioni approssimanti i campi u e v e le

funzioni approssimanti ( , )x i x yφ e ( , )y i x yφ possono essere fra loro diverse.

Introducendo le approssimazioni scelte per u e v (6.20) nell’espressione dell’energia potenziale totale (6.17), tenuto conto delle espressioni (6.18) e (6.19), ed imponendo la condizione di stazionarietà di Π rispetto a tutti i coefficienti incogniti x ic e y ic si ottiene

un sistema di x yn n+ equazioni in x yn n+ incognite. Scelte bene le funzioni approssimanti

il sistema è risolvibile e fornisce un’unica soluzione approssimata del problema.

6.1.5 Applicazione A titolo di esempio si può studiare il problema della lastra rappresentata in Figura 6-2, determinando una soluzione approssimata e confrontandola con le soluzioni fornita dai modelli trave Eulero-Bernoulli e Timoshenko.

Page 54: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

54

Figura 6-2: Lastra soggetta a carico verticale uniforme.

Page 55: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

55

6.2 Piastra di Kirchhoff-Love

Vengono utilizzate le seguenti notazioni:

• Spostamenti

uvw

=

u (6.21)

• Deformazioni

,

,

,

, ,

, ,

, ,

0 0

0 0

0 0

0

0

0

x x

y y

z z

xy y x

yz z y

xz z x

x

u yvw z

u vy xv w

u wz y

z x

εεεγγγ

∂ ∂

∂ ∂ ∂ ∂ = = = = ∂ ∂+

∂ ∂ + ∂ ∂+

∂ ∂ ∂ ∂ ∂ ∂

ε Su S (6.22)

• Tensioni

x

y

z

xy

yz

xz

σσστττ

=

σ (6.23)

• Legame costitutivo o Ortotropo

Page 56: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

56

1

1

1

1

1

1

yx zxx x y z

x y z

xy zyy x y z

x y z

yzxzz x y z

x y z

xy xyxy

yz yzyz

xz xzxz

E E E

E E E

E E E

G

G

G

ν νε σ σ σ

ν νε σ σ σ

ννε σ σ σ

γ τ

γ τ

γ τ

= − −

= − + −

= − − +

=

=

=

(6.24)

ovvero

1 0 0 0

1 0 0 0

1 0 0 0

10 0 0 0 0

10 0 0 0 0

10 0 0 0 0

yx zx

x y z

xy zy

x y z

yzxz

x y z

xy

yz

xz

E E E

E E E

E E E

G

G

G

ν ν

ν ν

νν

− −

− −

− − = =

ε Aσ A (6.25)

con

yx xy zy yzzx xz

y x z x z yE E E E E Eν ν ν νν ν

= = = (6.26)

per cui per definire completamente il legame elastico lineare ortotropo è necessario determinare il valore di 9 costanti elastiche;

o Isotropo

Page 57: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

57

( )

( )

( )

( )

( )

( )

2 11

2 11

1 2 1

xy xyx x y z

y y x z yz yz

z z x y xz xz

EE

E E

E E

νγ τε σ ν σ σ

νε σ ν σ σ γ τ

νε σ ν σ σ γ τ

+= = − +

+ = − + =

+ = − + =

(6.27)

ovvero

( )( )

( )

1 0 0 01 0 0 0

1 0 0 010 0 0 2 1 0 00 0 0 0 2 1 00 0 0 0 0 2 1

E

ν νν νν ν

νν

ν

− − − − − −

= = + + +

ε Aσ A (6.28)

6.2.1 Ipotesi di base Il modello di piastra di Kirchhoff-Love è basato sulle seguenti ipotesi:

• la dilatazione lineare in direzione dello spessore della piastra è trascurabile: la fibra ortogonale al piano medio della piastra risulta inestensibile;

• gli scorrimenti angolari tra il piano della piastra e la fibra ortogonale a tale piano sono trascurabili: la fibra ortogonale al piano medio della piastra risulta rettilinea ed ortogonale alla superficie media della piastra a deformazione avvenuta;

• la tensione normale in direzione dello spessore della piastra è trascurabile.

Si introduce un sistema di riferimento, come rappresentato in Figura 6-3.

Figura 6-3: Piastra di spessore costante.

x

y z

qh: spessore della piastra

Page 58: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

58

In formule, le ipotesi appena introdotte si traducono in:

0

00

z

yz xz

z

εγ γ

σ

== =

=

(6.29)

La (6.29)1 risulta:

,0 0z zwε = ⇒ = (6.30)

che integrata nello spessore fornisce:

( , )w w x y= (6.31)

Le ipotesi sugli scorrimenti angolari (6.29)2 risultano:

, ,

, ,

0 00 0

yz z y

xz z x

v wu w

γ

γ

= ⇒ + =

= ⇒ + = (6.32)

che integrati nelle spessore forniscono:

,

,

= −

= −x

y

u z wv z w

(6.33)

In definitiva, sulla base delle ipotesi cinematiche (6.29)1 e (6.29)2 si ottiene la seguente forma di rappresentazione degli spostamenti:

,

,

( , )( , )

( , )

x

y

u z w x yv z w x yw w x y

= −

= −

=

(6.34)

Sulla base del campo di spostamenti (6.34) e delle equazioni di congruenza (6.22), le deformazioni non nulle valgono:

,

,,

,, ,

00

x x x xx

y y y yy

xy xy xy y x

wz

w

ε χ χε χ χγ χ χ

•− = = = = • − • •

ε S S (6.35)

dove il vettore , ,

T

x yw w− − è l’opposto del gradiente dello spostamento trasversale w ,

che rappresenta la rotazione (infinitesima) della fibra ortogonale al piano medio della piastra. Inoltre le quantità xχ , yχ e xyχ sono le curvature della piastra:

,

,

, , ,2

x xx

y yy

xy xy yx xy

www w w

χ

χ

χ

= −

= −

= − − = −

(6.36)

Page 59: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

59

6.2.2 Equazioni di equilibrio tramite PLV Le tensioni che compiono lavoro per le deformazioni non nulle della piastra di Kirchhoff-Love sono solo le xσ , yσ e xyτ , per cui il lavoro virtuale interno vale:

T TviL dV dA

Ω= =∫ ∫ε σ χ M

M (6.37)

avendo introdotto il vettore dei momenti M come:

/2

/2/2

/2/2

/2

h

xhx

h

y yhhxy

xyh

z dzMM z dzM

z dz

σ

σ

τ

= =

∫∫∫

M (6.38)

In Figura 6-4 sono illustrati i momenti agenti nel piano medio della piastra; si evidenzia che:

• xM è il momento flettente agente nella direzione y ;

• yM è il momento flettente agente nella direzione x− ;

• xyM è il momento torcente.

Figura 6-4: Momenti agenti nel piano medio della piastra.

In forma esplicita, l’espressione (6.37) è:

( )

( ), , , ,

vi x x y y xy xy

x xx xy yx y yy xy xy

L M M M dA

M w M w M w M w dA

χ χ χ= + +

= − + + +

∫∫

M

M

(6.39)

Integrando per parti, si ottiene:

xyM

xyM

xM

yM

Page 60: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

60

( )

( ) ( ), , , , , , , ,

, , , ,

vi x x x xy x y y y y xy y x

x x xy y x y y xy x y

L M w M w M w M w dA

M w M w n M w M w n ds∂

= + + +

− + + +

∫∫

M

M

(6.40)

integrando ancora per parti si ha:

( )( ) ( )( ) ( )

, , ,

, , , ,

, , , ,

2vi x xx xy xy y yy

x x xy y x xy x y y y

x x xy y x y y xy x y

L M M M wdA

M M n M M n wds

M w M w n M w M w n ds∂

= − + +

+ + + +

− + + +

∫∫∫

M

M

M

(6.41)

Si introducono i tagli trasversali, definiti come:

, , , ,= + = +x x x xy y y xy x y yQ M M Q M M (6.42)

per cui la (6.41) diventa:

( )

( ) ( )

, , ,

, , , ,

2

= − + +

+ +

− + + +

∫∫∫

vi x xx xy xy y yy

x x y y

x x xy y x y y xy x y

L M M M wdA

Q n Q n wds

M w M w n M w M w n ds

M

M

M

(6.43)

In relazione ai termini di bordo dell’espressione (6.43), si consideri sulla curva ∂M un riferimento locale ,n t , definito dalla normale n e dalla tangente t alla curva, come

rappresentato in Figura 6-5. Si evidenzia che il momento nM è diretto secondo t .

Figura 6-5: Tipico tratto della frontiera del piano medio della piastra.

I termini di bordo si possono riscrivere allora come:

( ) ( )[ ] ( ) ( )

( )

, , , ,

, , , ,

, ,

x x y y x x xy y x y y xy x y

n n t t n n nt t n t t nt n t

n n n nt t

Q n Q n wds M w M w n M w M w n ds

Q n Q n wds M w M w n M w M w n ds

Q wdA M w M w ds

∂ ∂

∂ ∂

∂ ∂

+ − + + + = + − + + +

= − +

∫ ∫∫ ∫∫ ∫

M M

M M

M M

(6.44)

Page 61: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

61

essendo 1nn = e 0tn = .

Lungo la curva ∂M , il termine ,nt tM w può ancora essere integrato per parti, per cui si

ottiene:

( )

( ) [ ], ,

, ,

∂ ∂

∂∂∂ ∂

− +

= + − −

∫ ∫∫ ∫

n n n nt t

n nt t n n nt

Q wdA M w M w ds

Q M wdA M w ds M wM M

MM M

(6.45)

avendo indicato con ∂∂M i punti che rappresentano l’inizio e la fine della curva ∂M . Il termine ,+n nt tQ M è detta taglio di Kirchhoff.

Il lavoro virtuale interno, in definitiva si scrive come:

( )( ) [ ]

, , ,

, , , ,

2vi x xx xy xy y yy

n n nt t nt t n n nt

L M M M wdA

M M M wds M w ds M w∂∂∂ ∂

= − + +

+ + + − −

∫∫ ∫

M

MM M

(6.46)

Per consistenza con l’espressione derivata per il lavoro virtuale interno, il lavoro virtuale esterno deve essere compiuto da enti che lavorano in M per w , su ∂M per w e per ,nw e

su ∂∂M per w . Ne consegue che il lavoro virtuale esterno vale:

[ ],ve nL q wds f wds m w ds F w∂∂∂ ∂

= + − −∫ ∫ ∫ MM M M (6.47)

dove q è il carico per unità di superficie, trasversale al piano medio della piastra, che

compie lavoro per lo spostamento w in M ; f è il carico per unità di lunghezza,

trasversale al piano medio della piastra, che compie lavoro per lo spostamento w su ∂M ; m è una coppia per unità di lunghezza che compie lavoro per la rotazione ,nw− su ∂M ;

F è una forza concentrata agente su punti di discontinuità della frontiera di ∂M , che compie lavoro per lo spostamento w .

Eguagliando il lavoro virtuale esterno a quello interno si ottiene:

( )( ) [ ]

[ ]

, , ,

, , , ,

,

2x xx xy xy y yy

n n nt t nt t n n nt

n

M M M wdA

M M M wds M w ds M w

q wdA f wds m w ds F w w

∂∂∂ ∂

∂∂∂ ∂

− + +

+ + + − −

= + − − ∀

∫∫ ∫

∫ ∫ ∫

M

MM M

MM M M

(6.48)

da cui si ricavano le equazioni di equilibrio di campo e sulla frontiera:

, , ,

,

2 0 in su su in

+ + + =+ = ∂= ∂= ∂∂

x xx xy xy y yy

n nt t

n

nt

M M M qQ M fM mM F

MMMM

(6.49)

Page 62: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

62

Le equazioni di equilibrio di campo si possono determinare anche seguendo un ragionamento più meccanico, facendo cioè riferimento al tipico elemento di piastra.

Facendo riferimento alla Figura 6-6, l’equazione di equilibrio nella direzione z fornisce:

( ) ( ) 0x x x y y yQ Q Q y Q Q Q x q x y+ ∆ − ∆ + + ∆ − ∆ + ∆ ∆ = (6.50)

ovvero nel limite per 0x∆ → e 0y∆ → si ottiene:

, , 0x x y yQ Q q+ + = (6.51)

Figura 6-6: Tipico elemento di piastra soggetto a carico distribuito ed alle sollecitazioni di taglio trasversale.

Figura 6-7: Tipico elemento di piastra soggetto a carico distribuito ed alle sollecitazioni flessionali, torsionali e di taglio trasversale.

y yQ Q+ ∆

x xQ Q+ ∆

xQ

yQ

qx

QxQx+∆Qx

Qy+∆Qy

Qy

∆y∆x

y

x

Mx

Mx+∆Mx

My

My+∆My

Mxy

Mxy

Mxy+∆Mxy

Mxy+∆Mxy∆y∆x

y

Qx

Qx+∆Qx

Qy+∆Qy

Qy

Page 63: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

63

Facendo riferimento alla Figura 6-7, l’equazione di equilibrio alla rotazione intorno all’asse y ed intorno all’asse x , rispettivamente forniscono:

( ) ( ) ( )

( ) ( ) ( )

2

2

02

02

x x x xy xy xy x x

y y y xy xy xy y y

x yM M M y M M M x Q Q x y q

x yM M M x M M M y Q Q x y q

∆ ∆+ ∆ − ∆ + + ∆ − ∆ − + ∆ ∆ ∆ − =

∆ ∆+ ∆ − ∆ + + ∆ − ∆ − + ∆ ∆ ∆ − =

(6.52)

Semplificando, dividendo tutti i membri per x y∆ ∆ , nel limite per 0x∆ → e 0y∆ → si

ottiene:

, ,

, ,

x x x xy y

y xy x y y

Q M MQ M M

= +

= + (6.53)

Sostituendo le espressioni trovate per xQ e yQ nell’equazione (6.51), si perviene alla:

, , ,2 0x xx xy xy y yyM M M q+ + + = (6.54)

ovvero alla prima delle (6.49).

Figura 6-8: Calcolo del momento flettente e torcente per la generica direzione.

Le equazioni di equilibrio per l’elemento di piastra intorno agli assi definiti da t e da n si scrivono come:

x

Mx Mn

My

Mxy

MxyMnt

y

nt

Page 64: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

64

( ) ( )( ) ( )

n y xy xy x

nt y xy xy x

M M x M x M y M y

M M x M x M y M y

+ ∆ + ∆ − + ∆ + ∆ − =

− + ∆ + ∆ − + ∆ + ∆ − =

t i j i j 0

n i j i j 0

(6.55)

ovvero, moltiplicando scalarmente la prima per t e la seconda per n :

00

n y x xy y xy x x y

nt y x xy y xy x x y

M M x t M x t M y t M y tM M x n M x n M y n M y n

+ ∆ − ∆ + ∆ − ∆ =

− + ∆ − ∆ + ∆ − ∆ =

(6.56)

dividendo le due equazioni per e ricordando che:

x y x yy x x yn n t t∆ ∆ ∆ ∆

= = = − =

(6.57)

si ottiene:

( ) ( )

2 2

2 2

2n x y y x xy x y

nt y x x y xy x y

M M t M t M t t

M M M n n M n n

= + −

= − + − (6.58)

6.2.3 Legame costitutivo Ricordando l’ipotesi sul valore nullo della tensione normale nello spessore della piastra,

0zσ = , le prime due e la quarta delle equazioni costitutive per il problema 3D, riportate

nelle (6.27), diventano:

( )

1

1

2 1

x x y

y y x

xy xy

E

E

E

ε σ νσ

ε σ νσ

νγ τ

= −

= −

+=

(6.59)

che invertite forniscono:

( )

2

2

1

1

2 1

x x y

y y x

xy xy

E

E

E

σ ε ν εν

σ ε ν εν

τ γν

= + −

= + −

=+

(6.60)

ovvero in forma matriciale:

Page 65: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

65

( )

2 2

2 2

1 01 1

1 01 1

10 02 1

E

νν ν

νν ν

ν

− − = = − − +

σ Dε D (6.61)

Nel caso si consideri un legame costitutivo ortotropo, si ha:

( )

01 1

01 1

10 02 1

yx xx

xy yx xy yx

xy y y

xy yx xy yx

EE

E E

νν ν ν ν

νν ν ν ν

ν

− −

= = − −

+

σ Dε D (6.62)

Tenuto conto dell’equazione (6.35), si ha:

z=σ Dχ (6.63)

e, tramite la (6.38), si determina il vettore dei momenti come:

3/2 /2 2

/2 /2 12h h

h h

hz dz z dz− −

= = = =∫ ∫M σ Dχ Dχ D D (6.64)

che in esplicito, per il caso isotropo, fornisce:

( ) ( )( ) ( )

( )

, ,

, ,

,1 1

2

x x y xx yy

y y x yy xx

xy xy xy

M D D w w

M D D w w

M D Dw

χ ν χ ν

χ ν χ ν

ν χ ν

= + = − +

= + = − +

−= = − −

(6.65)

avendo posto:

( )3

21 12E hDν

=−

(6.66)

In Figura 6-9 sono schematicamente rappresentati gli andamenti delle tensioni normali e tangenziali nello spessore della piastra.

Page 66: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

66

Figura 6-9: Andamento delle tensioni normali e tangenziali nello spessore della piastra.

Sostituendo le espressioni (6.65) nell’equazione di equilibrio (6.49)1, si ottiene l’equazione di campo della piastra di Eulero-Bernoulli:

( ) ( ) ( ), , , , ,, ,0 2 1 in xx yy xyxy yy xxxx yy

D w w Dw D w w qν ν ν= + + − + + − M (6.67)

ovvero:

, , ,2 in xxxx xyxy yyyyq qw w w wD D

+ + = ⇒ ∆∆ = M (6.68)

6.2.4 Energia potenziale totale L’energia potenziale totale per una piastra vale:

(6.69)

dove l’energia interna è:

( )

2

,

, , , ,

,

1 1 12 2 2

1 22

2

T T T

xx

xx yy xy yy

xy

dV z dV dA

ww w w w dA

w

Ω ΩΦ = = =

− = − − − − −

∫ ∫ ∫

u ε σ χ Dχ χ Dχ

D

M

M

(6.70)

e il potenziale dei carichi è:

( ) [ ],nq wdA f wds m w ds F w∂∂∂ ∂

Λ = − + − − ∫ ∫ ∫uMM M M

(6.71)

( ) ( ) ( )Π = Φ + Λu u u

( )Φ u

( )Λ u

Page 67: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

67

La formulazione variazionale approssimata del problema della piastra si ottiene considerando una forma di rappresentazione del campo degli spostamenti, come:

( ) ( )1

, ,n

i ii

w x y c x yφ=

= ∑ (6.72)

essendo

( , ) 1,2,..., funzioni approssimanti di

1,2,..., coefficienti incogniti per l'approssimazione di i

i

x y i n wc i n wφ =

=

Introducendo l’approssimazione scelta per w (6.72) nell’espressione dell’energia potenziale totale (6.69), tenuto conto delle espressioni (6.70) e (6.71), ed imponendo la condizione di stazionarietà di rispetto ai coefficienti incogniti ic si ottiene un sistema di

n equazioni in n incognite. Scelte bene le funzioni approssimanti il sistema è risolvibile e fornisce un’unica soluzione approssimata del problema.

6.2.5 Applicazione A titolo di esempio si può studiare il problema della piastra rappresentata in Figura 6-10, determinando una soluzione approssimata.

Figura 6-10: Piastra appoggiata sui quattro lati.

Π

Page 68: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

68

6.3 Piastra di Reissner-Mindlin

6.3.1 Ipotesi di base Il modello di piastra di Reissner-Mindlin (FSDT, first-order shear deformation theory) è basato sulle seguenti ipotesi:

• la dilatazione lineare in direzione dello spessore della piastra è trascurabile: la fibra ortogonale al piano medio della piastra risulta inestensibile;

• gli scorrimenti angolari tra il piano della piastra e la fibra ortogonale a tale piano sono trascurabili: la fibra ortogonale al piano medio della piastra risulta rettilinea ma non necessariamente ortogonale alla superficie media della piastra a deformazione avvenuta;

• la tensione normale in direzione dello spessore della piastra è trascurabile.

Si introduce un sistema di riferimento, come rappresentato in Figura 6-3.

In formule, le ipotesi appena introdotte si traducono in:

( ) ( )0

, ,0

z

yz y xz x

z

g x y g x yεγ γ

σ

=

= =

=

(6.73)

Sulla base delle ipotesi cinematiche (6.73)1 e (6.73)2 si ottiene la seguente forma di rappresentazione degli spostamenti:

( , )( , )

( , )

x

y

u z x yv z x yw w x y

ϕϕ

==

=

(6.74)

Sulla base del campo di spostamenti (6.74), le deformazioni non nulle valgono:

,

,

, ,

0 00 00

x x x x

y y y x y

xy xy xy y y x

wz

ε χ χε χ χ ϕγ χ χ ϕ

• = = = = • • •

ε S S (6.75)

,

,

1 00 1

xz xx

yz yy

ϕγ

ϕ

• = = = •

γ T T (6.76)

dove le quantità xχ , yχ e xyχ sono le curvature della piastra:

,

,,

,, ,

x xx x

y yy y

x y y x

ww

ϕϕ

ϕϕ

ϕ ϕ

+ = = + +

χ γ (6.77)

Page 69: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

69

6.3.2 Equazioni di equilibrio tramite PLV Le tensioni che compiono lavoro per le deformazioni non nulle della piastra di Reissner-

Mindlin sono la , ,T

x y xyσ σ τ=σ e ,T

xz yzτ τ=τ per cui il lavoro virtuale interno vale:

( ) ( )T T T TviL dV dA

Ω= + = +∫ ∫ε σ γ τ χ M γ Q

M (6.78)

avendo introdotto il vettore dei momenti M ed il vettore dei tagli Q come:

/2

/2/2/2 /2

/2/2/2 /2

/2

h

x hhx

xzh hxy y hh

xyzhxy h

xyh

z dzM dzQM z dz

Q dzMz dz

στ

στ

τ

= = = =

∫∫

∫∫

M Q (6.79)

In forma esplicita, l’espressione (6.78) fornisce:

( )( )

( ) ( ), , , ,

, ,

vi x x y y xy xy x xz y yz

x x x xy y x y y y xy x y

x x x y y y

L M M M Q Q dA

M M M M dA

Q w Q w dA

χ χ χ γ γ

ϕ ϕ ϕ ϕ

ϕ ϕ

= + + + +

= + + +

+ + + +

∫∫∫

M

M

M

(6.80)

Integrando per parti, si ottiene:

( )( ) ( )

( ) ( )( )

, , , ,

, ,

vi x x x xy x y y y y xy y x

x x y y x x y y

x x xy y x y y xy x y

x x y y

L M M M M dA

Q Q wdA Q Q dA

M M n M M n ds

Q n Q n wds

ϕ ϕ ϕ ϕ

ϕ ϕ

ϕ ϕ ϕ ϕ∂

= − + + +

− + + +

+ + + +

+ +

∫∫ ∫∫∫

M

M M

M

M

(6.81)

ovvero

( ) ( )( )

( ) ( )( )

, , , ,

, ,

vi x x xy y x x xy x y y y y

x x y y

x x xy y x y y xy x y

x x y y

L M M Q M M Q dA

Q Q wdA

M M n M M n ds

Q n Q n wds

ϕ ϕ

ϕ ϕ ϕ ϕ∂

= − + − + + −

− + +

+ + + +

+ +

∫∫∫∫

M

M

M

M

(6.82)

In relazione ai termini di bordo dell’espressione (6.82), si consideri sulla curva ∂M un riferimento locale ,n t , definito dalla normale e dalla tangente; tali termini si possono riscrivere come:

Page 70: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

70

( ) ( )( )( ) ( )

( )

x x xy y x y y xy x y

x x y y

n n nt t n t t nt t t

n n

n n nt t n n

M M n M M n ds

Q n Q n wds

M M n M M n ds

Q n wds

M M ds Q n wds

ϕ ϕ ϕ ϕ

ϕ ϕ ϕ ϕ

ϕ ϕ

∂ ∂

+ + +

+ +

= + + +

+

= + +

∫∫∫∫∫ ∫

M

M

M

M

M M

(6.83)

essendo 1nn = e 0tn = .

Il lavoro virtuale interno, in definitiva si scrive come:

( ) ( )( )( ) ( )

, , , ,

, ,

vi x x xy y x x xy x y y y y

x x y y

n n nt t n n t t

L M M Q M M Q dA

Q Q wdA

M M ds Q n Q n wds

ϕ ϕ

ϕ ϕ∂ ∂

= − + − + + −

− + +

+ + + +

∫∫∫ ∫

M

M

M M

(6.84)

Per consistenza con l’espressione derivata per il lavoro virtuale interno, il lavoro virtuale esterno deve essere compiuto da enti che lavorano in M per w , per xϕ e per yϕ e su

∂M per w , per nϕ e per tϕ . Ne consegue che il lavoro virtuale esterno vale:

( ) ( )ve x x y y n n nt tL q wdA c c dA f wds m m dsϕ ϕ ϕ ϕ∂ ∂

= + + + + +∫ ∫ ∫ ∫M M M M (6.85)

dove q è il carico per unità di superficie, trasversale al piano medio della piastra, che

compie lavoro per lo spostamento w in M ; xc e yc sono coppie distribuite per unità di

superficie, che compiono lavoro per le rotazioni spostamento xϕ e yϕ in M ; f è il

carico per unità di lunghezza, trasversale al piano medio della piastra, che compie lavoro per lo spostamento w su ∂M ; nm e ntm sono coppie per unità di lunghezza che

compiono lavoro per le rotazioni nϕ e tϕ su ∂M .

Eguagliando il lavoro virtuale esterno a quello interno si ottiene:

( ) ( )( )

( ) ( )

( )

, , , ,

, ,

0 , , in ; , , in

x x xy y x x x xy x y y y y y

x x y y

n n n nt nt t

n n x y n t

M M Q c M M Q c dA

Q Q q wdA

M m M m ds

Q n f wds w w

ϕ ϕ

ϕ ϕ

ϕ ϕ ϕ ϕ∂

+ − + + + − +

+ + +

− − + −

− − = ∀ ∀ ∂

∫∫∫∫

M

M

M

MM M

(6.86)

da cui si ricavano le equazioni di equilibrio di campo e sulla frontiera:

Page 71: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

71

, ,

, ,

, ,

0 in 0 in

0 in 0 su 0 su 0 su

x x xy y x x

xy x y y y y

x x y y

n n

nt nt

n n

M M Q cM M Q cQ Q qM mM mQ n f

+ − + =+ − + =

+ + =− = ∂− = ∂− = ∂

MMM

MMM

(6.87)

6.3.3 Legame costitutivo Ricordando l’ipotesi sul valore nullo della tensione normale nello spessore della piastra,

0zσ = , le equazioni costitutive diventano:

( )

( )

2 2

2 2

1 01 1

1 01 1

10 02 1

1 00 12 1

E

E

νν ν

νν ν

ν

ν

− − = = − − +

= = +

σ Dε D

τ Gγ G

(6.88)

Tenuto conto delle equazioni (6.75) e (6.88), si ha:

z= =σ Dχ τ Gγ (6.89)

per cui i vettori dei momenti e dei tagli trasversali si determinano come:

3/2 /2 2

/2 /2

/2 /2

/2 /2

12h h

h h

h h

h h

hz dz z dz

dz dz hκ

− −

− −

= = = =

= = = =

∫ ∫

∫ ∫

M σ Dχ Dχ D D

Q τ G γ G γ G G (6.90)

essendo κ il fattore di correzione a taglio, di cui si è già discusso per la trave. Per la piastra si ha 5 / 6κ = .

In esplicito, per il caso isotropo, fornisce:

( ) ( )( ) ( )

( )( )( )

, ,

, ,

, ,

,

,

1 12 2

x x y x x y y

y y x y y x x

xy xy x y y x

x xz x x

y yz y y

M D D

M D D

M D D

Q H H w

Q H H w

χ ν χ ϕ ν ϕ

χ ν χ ϕ ν ϕ

ν νχ ϕ ϕ

γ ϕ

γ ϕ

= + = +

= + = +

− −= = +

= = +

= = +

(6.91)

avendo posto:

Page 72: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

72

( ) ( )3

2

56 2 11 12

E h ED H hνν

= =+−

(6.92)

Sostituendo le espressioni (6.92) nelle equazioni di equilibrio (6.87)1, (6.87)2, e (6.87)3 si ottengono le equazioni di campo della piastra di Mindlin:

( ) ( )

( ) ( )( )

, , , , , ,

, , , , , ,

, , , ,

2 0 in 2

2 0 in 2

0 in

x xx x yy y xy y yx x yy x x x

x yx y xx y yy x xy y xx y y y

xx yy x x y y

D H w c

D H w c

H w w q

ϕ ϕ ϕ ν ϕ νϕ ϕ

ϕ ϕ ϕ ν ϕ νϕ ϕ

ϕ ϕ

+ + + − − + + =

+ + + − − + + =

+ + + + =

M

M

M

(6.93)

6.3.4 Energia potenziale totale L’energia potenziale totale per una piastra vale:

(6.94)

dove l’energia interna è:

( ) ( ) ( )

, ,, ,

, ,, ,

, , , ,

1 12 2

12

T T T T

TTx x x x

x x x xy y y y

y y y yx y y x x y y x

dV dA

w wdA

w w

ϕ ϕϕ ϕ

ϕ ϕϕ ϕ

ϕ ϕ ϕ ϕ

ΩΦ = + = +

+ + = + + + + +

∫ ∫

u ε σ γ τ χ Dχ γ Gγ

D G

M

M

(6.95)

e il potenziale dei carichi è:

( ) ( ) ( )x x y y n n nt tq wdA c c dA f wds m m dsϕ ϕ ϕ ϕ∂ ∂

Λ = − + + + + + ∫ ∫ ∫ ∫uM M M M

(6.96)

La formulazione variazionale approssimata del problema della piastra si ottiene considerando una forma di rappresentazione del campo di spostamenti w , del campo delle rotazioni xϕ e del campo delle rotazioni yϕ come:

( ) ( )

( ) ( )

( ) ( )

1

1

1

, ,

, ,

, ,

w

x

y

nw w

i ii

nx x

x i iin

y yy i i

i

w x y c x y

x y c x y

x y c x y

φ

ϕ φ

ϕ φ

=

=

=

=

=

=

(6.97)

essendo

( ) ( ) ( )Π = Φ + Λu u u

( )Φ u

( )Λ u

Page 73: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

73

( ), 1,2,..., funzioni approssimanti di 1,2,..., coefficienti incogniti per l'approssimazione di

wi w

wi w

x y i n wc i n wφ =

=

( ), 1,2,..., funzioni approssimanti di 1,2,..., coefficienti incogniti per l'approssimazione di

xi x x

xi x x

x y i nc i nφ ϕ

ϕ==

( ), 1,2,..., funzioni approssimanti di

1,2,..., coefficienti incogniti per l'approssimazione di

yi y y

yi y y

x y i nc i nφ ϕ

ϕ==

Introducendo l’approssimazione scelta per w , per xϕ e per yϕ nell’espressione (6.97)

dell’energia potenziale totale (6.94), tenuto conto delle espressioni (6.95) e (6.96), ed

imponendo la condizione di stazionarietà di rispetto ai coefficienti incogniti wic , x

ic e y

ic si ottiene un sistema di w x yn n n+ + equazioni in w x yn n n+ + incognite. Scelte bene

le funzioni approssimanti il sistema è risolvibile e fornisce un’unica soluzione approssimata del problema.

6.3.5 Applicazione A titolo di esempio si può studiare il problema della piastra rappresentata in Figura 6-10, determinando una soluzione approssimata.

Π

Page 74: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

74

7 Soluzione di problemi strutturali non lineari

7.1 Introduzione

Nella meccanica delle strutture, spesso si assume che gli spostamenti siano piccoli rispetto alla dimensione caratteristica della struttura e che le deformazioni siano infinitesime, così che la misura di deformazione è ottenuta come la parte simmetrica del gradiente degli spostamenti e che le equazioni di equilibrio sono scritte direttamente sulla configurazione indeformata del corpo. D’altra parte, si assume anche molto frequentemente che la relazione tra le tensioni e le deformazioni sia di tipo elastico lineare. Sotto le ipotesi sopra riportate il problema dell’equilibrio elastico è lineare. Ne consegue che implementando il metodo degli elementi finiti per un problema dell’equilibrio elastico lineare, si giunge ad un sistema algebrico di tipo lineare.

Talvolta però non è possibile con buona approssimazione assumere valide ipotesi alla base della linearità del problema strutturale. Può infatti accadere che la struttura soggetta alle dovute condizioni di carico si deformi in modo tanto significativo che non è appropriato considerare come misura di deformazione quella infinitesima, ovvero la parte simmetrica del gradiente degli spostamenti, e di conseguenza non è lecito scrivere le equazioni di equilibrio sulla configurazione indeformata della struttura. E’ necessario allora tenere in conto di misure non lineari della deformazione ed imporre l’equilibrio sulla configurazione deformata, che non è nota a priori e costituisce un’incognita del problema strutturale. In tal caso il problema dell’equilibrio strutturale diventa non lineare e tale non linearità è generalmente denotata come non linearità geometrica.

I materiali che costituiscono le costruzioni presentano comportamento di tipo elastico lineare in un intervalli di tensioni o deformazioni limitato. Superato un certo valore della deformazione il materiale si comporta in modo tale che la relazione costitutiva è non lineare. Questo è il caso della plasticità, del danneggiamento e del comportamento unilaterale. Anche in tal caso il problema dell’equilibrio strutturale diventa non lineare tale non linearità è generalmente denotata come non linearità materiale.

Nel seguito si descrivono alcune procedure classiche utilizzate per risolvere problemi algebrici non lineari. La discussione è limitata al caso di non linearità materiale.

Page 75: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

75

7.2 Materiale con comportamento elastico non lineare

Un materiale è detto elastico quando assegnata che sia la deformazione, la tensione è univocamente determinata ed, inoltre, l’energia interna non dipende dal particolare percorso deformativo seguito. Per una struttura costitutiva da materiale elastico, lo stato deformativo e tensionale non dipende dalla particolare storia di carico seguita ma solo dal valore attuale delle forze agenti su essa.

In Figura 7-1 è riportata schematicamente la risposta meccanica in termini di uno spostamento s di una struttura soggetta ad un sistema di forze (F). La relazione F-s risulta non lineare. Per una struttura costituita da materiale elastico, in fase di carico e di scarico si percorre comunque la stessa curva Λ(s).

Figura 7-1: Risposta strutturale schematica in termini di forza-spostamento.

Il problema che si intende affrontare è allora di determinare una strategia di calcolo, ovvero una procedura numerica che permetta di determinare lo spostamento che si ha una volta che sia assegnato che sia un valore della forza. In realtà il problema si pone nella seguente forma: dato un sistema strutturale, determinare il vettore degli spostamenti (nodali) s di una struttura che si hanno quando agisce un sistema di carico definito dal vettore F. In altra forma, si tratta di risolvere il problema struttura non lineare:

( ) =Λ s F (7.1)

L’equazione (7.1) si riscrive nella forma residuale:

( ) ( )= − =R s Λ s F 0 (7.2)

F

s

Λ(s)

Page 76: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

76

essendo R(s) il residuo da annullare.

Le procedure numeriche descritte nel seguito sono algoritmi capaci di trasformare il problema della ricerca della soluzione di un sistema di equazioni algebriche non lineari in una successione di problemi lineari. La procedura converge alla soluzione nel caso in cui la successione di soluzioni converge alla soluzione del problema iniziale non lineare.

7.3 Metodo di Newton-Raphson

Assumendo che la funzione residuo R(s) sia sufficientemente regolare, si effettua lo sviluppo in serie di R(s) a partire dal vettore spostamento so:

( )

( )

( ) ( )

( )

o

o

o o

o o

=

=

∂= + −

∂= − + − =

s s

s s

RR s R s s ss

ΛΛ s F s s 0s

(7.3)

Si evidenzia che la derivate del vettore Λ(s) rispetto al vettore s è una matrice detta matrice tangente e denotata con Kt:

, o

o o

t iij

jij

Ks= =

∂Λ∂= = ∂ ∂

s

s s s s

Λs

(7.4)

Assumendo so=0, è possibile determinare lo spostamento s1 risolvendo l’equazione (7.3):

( ) ( ) ( )1 11 ,0 ,0( ) ( )t t− −= − − = −s K Λ 0 F K R 0 (7.5)

Figura 7-2: schema iterativo secondo il metodo di Newton-Raphson

F

s

Λ(s) F*

s* 0 s1

R(s1

Kt0

R(s2

Kt1

Page 77: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

77

Nota che sia la soluzione s1 è possibile determinare il residuo

1 1( ) ( )= −R s Λ s F (7.6)

che, in generale, risulta diverso da zero, come illustrato in Figura 7-2. A questo punto risolvendo l’equazione (7.3), con so=s1, si ottiene:

( ) ( ) ( )1 12 1 ,1 1 1 ,1 1( ) ( )t t− −= − − = −s s K Λ s F s K R s (7.7)

Noto s2, in modo analogo, si determina s3.

In definitiva il metodo di Newton-Raphson ovvero della tangente, consiste nella determinazione della successione di soluzioni:

( ) ( ) ( )1 11 , ,( ) ( )k k t k k k t k k− −+ = − − = −s s K Λ s F s K R s (7.8)

La successione converge quando il residuo R(sk) tende a zero. Naturalmente, numericamente non è possibile giungere al valore nullo per cui si assumerà soddisfacente la soluzione sk per la quale accade che il residuo R(sk) è sufficientemente piccolo. In particolare si definisce errore k-esimo la quantità:

( )k

ke =R s

F (7.9)

Si dimostra che tale procedura numerica quando converge, converge alla soluzione e converge con velocità quadratica. Di conseguenza l’algoritmo è caratterizzato da una buona velocità di convergenza. Si evidenzia che qualora la curva che rappresenta la risposta strutturale illustrata in Figura 7-2 tende ad una retta orizzontale, la matrice tangente tende a diventare singolare e, di conseguenza, il sistema di equazioni (7.8) tende a ad essere mal posto.

7.4 Metodo di Newton-Raphson modificato

Il metodo di Newton-Raphson modificato si ottiene a partire dal metodo di Newton-Raphson assumendo la matrice tangente costante nel processo iterativo, e pari a quella corrispondente al punto s=0. In altri termini, l’algoritmo consiste nella costruzione della successione di soluzioni ottenute risolvento l’equazione:

Page 78: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

78

( ) ( ) ( )1 11 ,0 ,0( ) ( )k k t k k t k− −+ = − − = −s s K Λ s F s K R s (7.10)

In Figura 7-3 è rappresentato lo schema iterativo del metodo di Newton-Raphson modificato. Si nota che anche quando la curva che rappresenta la risposta strutturale tende ad una retta orizzontale, l’algoritmo non entra in crisi, in quanto la matrice da invertire è sempre quella iniziale. In generale, tale procedura risulta molto stabile nelle iterazioni e raramente entra in crisi. Un vantaggio di tale metodo rispetto al metodo di Newton-Raphson è che la matrice che governa l’equazione (7.10) non varia durante le iterazioni e deve essere calcolata solo all’inizio della procedura. D’altra parte si nota che la convergenza di tale metodo risulta in generale più lenta rispetto a quella che si ha per il metodo della tangente esposto nel paragrafo precedente.

Figura 7-3: schema iterativo secondo il metodo di Newton-Raphson modificato

7.5 Metodo della secante

Il metodo della secante, detto anche di Picard, è basato sull’ipotesi, molto spesso soddisfatta, che il vettore delle forze interne Λ(s) sia esprimibile come il prodotto di una matrice K(s), che dipende dal valore degli spostamenti, moltiplicato per il vettore degli spostamenti, in formula:

( ) ( )=Λ s K s s (7.11)

F

s

Λ(s) F*

s* 0 s1

R(s1

Kt0

s2 s4

Page 79: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

79

così che l’equazione di equilibrio (7.1) diventa:

( ) =K s s F (7.12)

La matrice tangente si determina come:

[ ] [ ], ( ) ( )( )

o

o o o

t o

= = =

∂ ∂∂= = = +

∂ ∂ ∂s

s s s s s s

K s s K sΛK K s ss s s

(7.13)

La matrice secante è ottenuta dalla relazione (7.13) trascurando l’ultimo termine:

, ( )os o=sK K s (7.14)

L’algoritmo della matrice secante è basato sulla soluzione dell’equazione (7.12), assumendo noto il vettore degli spostamenti per la valutazione della matrice secante. Sostituendo la matrice secante al posto della tangente nella formula (7.8) si ha:

( ) ( )( ) ( )

( )

11 ,

1, ,

1,

( )k k s k k

k s k s k k

s k

−+

= − −

= − −

=

s s K Λ s F

s K K s F

K F

(7.15)

In Figura 7-4 è rappresentato lo schema iterativo del metodo della secante. In generale, tale procedura risulta abbastanza stabile nelle iterazioni. La velocità di convergenza è minore di quella ottenibile con il metodo di Newton-Raphason ma è maggiore di quella con il metodo di Newton-Raphason modificato. Come per il metodo di Newton-Raphason la matrice che governa l’equazione (7.15) varia durante le iterazioni e deve essere calcolata ed invertita ad ogni iterazione.

Page 80: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

80

Figura 7-4: schema iterativo secondo il metodo della secante

7.6 Sezione in cemento armato

A scopo esemplificativo, ma anche per illustrare un problema di notevole interesse applicativo, si sviluppano nel seguito i tre algoritmi introdotti nel paragrafo precedente per risolvere il problema della verifica di una sezione in cemento armato. I metodi sopra esposti, naturalmente, si possono applicare al problema della verifica agli stati limite per una sezione generica. Per semplicità nel seguito viene trattato il caso della verifica alle tensioni ammissibili di una sezione rettangolare in cemento armato.

Si considera la sezione riportata in Figura 7-5, soggetta alle sollecitazioni di sforzo normale e momento flettente. Si assume che il materiale conglomerato cementizio sia non reagente a trazione, mentre presenti un comportamento elastico lineare a compressione, come illustrato in Figura 7-6. In particolare, si assume:

0 se 0

se 0Eσ εσ ε ε

= >= ≤

(7.16)

La cinematica della sezione è caratterizzata dalla deformazione assiale ε0 e dalla curvatura χ, così che la deformazione nel generico punto della sezione è definita dalla classica relazione:

0 yε = ε + χ (7.17)

F

s

Λ(s) F*

s* 0 s1 s2 s4

K(s1) s1

Page 81: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

81

Figura 7-5: sezione in c.a.

Figura 7-6: modello di materiale non reagente a trazione per conglomerato

cementizio

Nella sezione si individua l’asse neutro, come:

00: 0n n ny y y εε χ

χ= + ⇒ = − (7.18)

x1

x2

Af2

Af1

yf2

yf1

h

b

ε

σ

Page 82: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

82

Si distinguono due possibili casi:

caso (a) χ ≥ 0

caso (b) χ < 0

Si definiscono gli assi individuati da y = y1, y = y2, tali che:

1 2

1 2

caso (a) / 2 min / 2,max- / 2, caso (b) min / 2,max- / 2, / 2

n

n

y h y h h yy h h y y h

= − == =

(7.19)

così che la parte di sezione compresa tra y1 ed y2 risulta soggetta ad una tensione di compressione che varia linearmente;

Ne consegue che, in termini di risultanti, si ha:

( ) ( ) ( )

( ) ( ) ( )

2

1

2

1

0 2 2 1 0 1 2 0 2 0

2 0 2 2 1 1 0 1 2 2 0 2 0

y

f f f f f fy

y

f f f f f f f fy

N Eb x dx E A y E A y A S

M x Eb x dx y E A y y E A y S J

ε χ ε χ ε χ ε χ

ε χ ε χ ε χ ε χ

= + + + + + = +

= + + + + + = +

∫(7.20)

dove

( ) ( )( ) ( )

( ) ( )

2 1 1 2

2 22 1 1 1 2 2

3 3 2 22 1 1 1 2 2

2

3

f f f

f f f f f

f f f f f

A Eb y y E A A

EbS y y E y A y A

EbJ y y E y A y A

= − + +

= − + +

= − + +

(7.21)

Si evidenzia che l’area, il momento statico ed il momento di inerzia sono funzioni di y1 ed y2, ovvero di yn, e quindi, tenendo conto delle (7.19) e della (7.18), sono funzioni della deformazione assiale ε0 e della curvatura χ.

In definitiva si ottiene un sistema di equazioni, che scritto in forma matriciale diventa:

0( ) ( )( ) ( )

n n

n n

A y S yNS y J yM

εχ

=

(7.22)

Page 83: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

83

7.6.1 Algoritmo di Newton-Raphson L’algoritmo consiste nei seguenti passi:

1. All’inizio della procedura si pone y1=-h/2 e y2=h/2 e si determinano le proprietà geometriche tramite le equazioni (7.21).

2. Si calcola la matrice delle proprietà geometriche:

0 0

0 0

A SS J

=

G (7.23)

3. Alla (k+1)-iterazione si risolve il problema (7.8) ottenendo una nuova coppia di valori per la deformazione assiale e per la curvatura:

( ) ( )( ) ( )

( ) ( )

110 0 0

10

0 0 0

0 0

1 10 0

0

, ,, ,

0

k k k k k k k

k k k k k k kn n

k kk k k k k

k

k

NA S A S Ndd d dMd dy d dyS J S J M

A S NR

S J M

R RR

NM

ε ε ε χε χχ χ

ε χ ε χ εε χ ε χ χ

ε ε χ χε χ

−+

+

+ +

= − + −

= −

∂ ∂

= + − + −∂ ∂

=

( ) ( )( ) ( ) ( )

( ) ( )( ) ( ) ( )

( )

0 0 0 10 0

0 00

0 0 0 1

0 0

10 0

0

, ,, ,

, ,, ,

10

k k

k k

k k k k kk k

k k k k k

A SNS JM

A SS J

NN A S A SMM S J S J

ε χ ε χ εε ε

ε χ ε χ χε

ε χ ε χ εχ χ

ε χ ε χ χχ

ε εε

+

+

+

∂− + − ∂

∂+ − ∂ ∂

= − + − + ∂ ( )

( ) ( )

( )

0 10 0

01 1

010 0

0

01

10

k k

k k k kk k k k

k k k k

k k k k kk k

k k k k k

A S A SS J S J

NN A S A SMM S J S J

εε ε

χ

εχ χ χ χ

χχ

εε ε

χε

+

+ +

+

∂+ − + − ∂ ∂

= − + − + ∂ ( )

( ) ( )

10 0

01 101

k k

k k k kk k k k

k k k k

A S A SS J S J

ε ε

εχ χ χ χ

χχ

+

+ +

∂+ − + − ∂

(7.24)

4. Si determina la nuova posizione dell’asse neutro dall’equazione (7.18). 5. Si definisce la nuova sezione reagente tramite le (7.19). 6. Noti che siano alla (k+1)-esima iterazione y1

k+1 e y2k+1, si determinano le proprietà

geometriche tramite le equazioni (7.21).

Page 84: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

84

7. In base alla soluzione disponibile si valutano le sollecitazioni ad esse associate come:

1 1 1 1

01 1 1 1

k k k k

k k k k

N A SM S J

εχ

+ + + +

+ + + +

=

(7.25)

8. Si determina il residuo:

1 1

1 1

k kN

k kM

R N NR M M

+ +

+ +

= −

= − (7.26)

9. Si valuta l’errore:

( ) ( )2 21 1

2 2

k kN Mk

R Re

N M

+ ++=

+ (7.27)

10. Se ek < tolleranza esci dalla procedura, altrimenti si torna al punto 2.

7.6.2 Algoritmo di Newton-Raphson modificato L’algoritmo consiste nei seguenti passi:

1. All’inizio della procedura si pone y1=-h/2 e y2=h/2 e si determinano le proprietà geometriche tramite le equazioni (7.21).

2. Si calcola la matrice delle proprietà geometriche:

0 0

0 0

A SS J

=

G (7.28)

3. Alla (k+1)-iterazione si risolve il problema (7.10) ottenendo una nuova coppia di valori per la deformazione assiale e per la curvatura:

1

10 01

k k k

k k k

NNMM

ε εχ χ

+−

+

= − −

G (7.29)

4. Si determina la nuova posizione dell’asse neutro dall’equazione (7.18). 5. Si definisce la nuova sezione reagente tramite le (7.19). 6. Noti che siano alla (k+1)-esima iterazione y1

k+1 e y2k+1, si determinano le proprietà

geometriche tramite le equazioni (7.21). 7. In base alla soluzione disponibile si valutano le sollecitazioni ad esse associate

come:

1 1 1 1

01 1 1 1

k k k k

k k k k

N A SM S J

εχ

+ + + +

+ + + +

=

(7.30)

8. Si determina il residuo;

1 1

1 1

k kN

k kM

R N NR M M

+ +

+ +

= −

= − (7.31)

Page 85: Appunti di MECCANICA DELLE STRUTTURE Prof. Elio Sacco ... · MECCANICA DELLE STRUTTURE . Prof. Elio Sacco Dipartimento di Ingegneria Civile e Meccanica Università di Cassino e del

85

9. si valuta l’errore:

( ) ( )2 21 1

2 2

k kN Mk

R Re

N M

+ ++=

+ (7.32)

10. Se ek < tolleranza esci dalla procedura, altrimenti si torna al punto 2.

7.6.3 Algoritmo della secante L’algoritmo consiste nei seguenti passi:

1. All’inizio della procedura si pone y1=-h/2 e y2=h/2 e si determinano le proprietà geometriche tramite le equazioni (7.21).

2. Alla (k+1)-iterazione si risolve il problema (7.22) ottenendo una nuova coppia di valori per la deformazione assiale e per la curvatura:

11

01

k k k

k k k

NA SMS J

εχ

−+

+

=

(7.33)

3. Si determina la nuova posizione dell’asse neutro dall’equazione (7.18). 4. Si definisce la nuova sezione reagente tramite le (7.19). 5. Noti che siano alla (k+1)-esima iterazione y1

k+1 e y2k+1, si determinano le proprietà

geometriche tramite le equazioni (7.21). 6. In base alla soluzione disponibile si valutano le sollecitazioni ad esse associate

come:

1 1 1 1

01 1 1 1

k k k k

k k k k

N A SM S J

εχ

+ + + +

+ + + +

=

(7.34)

7. Si determina il residuo;

1 1

1 1

k kN

k kM

R N NR M M

+ +

+ +

= −

= − (7.35)

8. si valuta l’errore:

( ) ( )2 21 1

2 2

k kN Mk

R Re

N M

+ ++=

+ (7.36)

9. Se ek < tolleranza esci dalla procedura, altrimenti si torna al punto 2.