Teoria della consolidazione 1D - people.unica.it · 2019. 4. 15. · 2 Teoria di Terzaghi...

Corso di Laurea in Ingegneria CivileCorso di Integrato di Geotecnica e Costruzioni di strade – Modulo di Geotecnica

Docente: ing. Giuseppe Tropeano

‐ 1 ‐

TEORIADELLACONSOLIDAZIONEDEITERRENIAGRANAFINE

Schemariassuntivo

1) Premessa

a. Descrizione e definizione della fenomenologia

2) Teoria di Terzaghi (consolidazione monodimensionale)

a. Ipotesi di base

ip.1) Flusso e deformazioni 1D

ip.2) Fluido non comprimibile

ip.3) Mezzo completamente saturo

ip.4) Incomprimibilità dei granuli solidi (i.e. granuli con rigidezza elevata; ρs=costante)

ip.5) Validità della legge di Darcy

ip.6) Comportamento elastico lineare dello scheletro solido

ip.7) Coefficiente di permeabilità costante nel tempo

ip.8) Carico applicato in maniera istantanea e costante (i.e. σz costante nel tempo)

b. Equazioni di base

1. Conservazione delle massa solida 11 0

sn vnt z

2. Conservazione della massa liquida 0

wn vn

t z

3. Continuità della massa (mezzo bifase) w s s

n v v vz z

4. Legge di Darcy f zhv kz

5. Legame costitutivo 1

z zed

d dE

6. Equilibrio dello scheletro solido 0zt

c .Equazione finale 2

2

z edv v

w

k Eu uc ct z

d. Integrazione della equazione del modello



‐ 2 ‐

1. Adimensionalizzazione 2 vczZ T t

H H

e. Soluzione

1. Isocrone

2. Grado di consolidazione

f. Valutazione sperimentale del coefficiente di consolidazione verticale



‐ 3 ‐

1 Premessa

1.1 Descrizioneedefinizionedellafenomenologia

Per consolidazione si intende la progressiva deformazione di un terreno prodotta da fenomeni idrodinamici connessi alla dissipazione delle sovrappressioni interstiziali, ū, indotte da un incremento della tensione totale, Δσz, dovuta a carichi esterni, q. Poiché questo fenomeno dipende dalla resistenza al moto della fase fluida nei terreni a grana fine, ha senso pratico parlare di consolidazione solo in riferimento a terreni con conducibilità idraulica bassa.

2 TeoriadiTerzaghi(consolidazionemonodimensionale)

2.1 Ipotesidibase

Si considera un terreno a grana fine omogeneo, le ipotesi di base del modello matematico della consolidazione 1D sono:

ip.1) Flusso e deformazioni 1D ip.2) Fluido non comprimibile ip.3) Mezzo completamente saturo ip.4) Incomprimibilità dei granuli solidi (i.e. granuli con rigidezza elevata; ρs=costante) ip.5) Validità della legge di Darcy ip.6) Comportamento elastico lineare dello scheletro solido ip.7) Coefficiente di permeabilità costante nel tempo ip.8) Carico applicato in maniera istantanea e costante (i.e. σz costante nel tempo)

2.2 Equazionidibase

2.2.1 Conservazionedellemassasolida

l’eq. di conservazione della massa è data dalla espressione:

massa uscente massa entrante massa iniziale - massa finaleunità di tempo unità di tempo unità di tempo

u e f im m m m (1)



‐ 4 ‐

Figura 1: Schema per la definizione dell’equazione di continuità della massa solida nel caso 1D.

in cui la massa entrante, me, nella faccia dx . dy è pari a:

1 e s sm n v dxdy dt (2)

dove: ρs è la massa volumica (densità) dei granuli del terreno; n è la porosità dell’elemento di volume dV = dx . dy . dz; vs velocità del flusso solido, dt istante temporale;

la massa uscente, mu:

1

1

s su s s

n vm n v dz dx dy dt

z (3)

la massa inizialmente presente nell’elemento di volume, mi:

1 i sm n dxdy dz (4)

la massa presente nell’elemento di volume dopo l’intervallo di tempo dt, mf:

1

1

sf s

nm n dt dx dy dz

t (5)

Sostituendo le eqq. (2), (3) , (4) e (5) nella (1) si ottiene:

11 1

11 1

s ss s s s

ss s

n vn v dz dx dy dt n v dx dy dt

z

nn dt dx dy dz n dx dy dz

t



‐ 5 ‐

Facendo gli opportuni passaggi algebrici si ottiene infine:

1 1

s s sn v nz t

Essendo il terreno omogeneo, la densità del solido non varia con la profondità, inoltre, per la (ip.4), è costante nel tempo, ottenendo infine:

1 1

sn v n

z t (6)

2.2.2 Conservazionedellemassafluida

Analogamente al caso precedente, l’eq. di conservazione della massa fluida è data dalla espressione:

massa uscente massa entrante massa iniziale - massa finaleunità di tempo unità di tempo unità di tempo

u e f im m m m (7)

Figura 2: Schema per la definizione dell’equazione di continuità della massa fluida nel caso 1D.

in cui la massa entrante, me, nella faccia dx . dy è pari a:

e w wm nv dx dy dt (8)

dove: ρw è la massa volumica (densità) del fluido interstiziale (p. es. acqua); n è la porosità dell’elemento di volume dV = dx . dy . dz; vw velocità assoluta del flusso, dt istante temporale;

la massa uscente, mu, è data da:



‐ 6 ‐

w wu w w

n vm n v dz dxdy dt

z (9)

la massa di fluido inizialmente presente nell’elemento di volume, mi:

massa iniziale i wm ndx dy dz (10)

la massa di fluido presente nell’elemento di volume dopo l’intervallo di tempo dt, mf:

massa finale

wf w

nm n dt dx dy dz

t (11)

Sostituendo le eqq. (8), (9), (10) e (11) nella (7) si ottiene:

w ww w w w

ww w

n vn v dz dx dy dt n v dx dy dt

z

nn dt dx dy dz n dx dy dz

t

(12)

Facendo gli opportuni passaggi algebrici si ottiene infine:

w w wn v nz t

(13)

Per le (ip.2) e (ip.3) la densità del fluido è costante sia con la profondità sia col tempo, si ottiene:

wn v nz t

(14)

2.2.3 Continuitàdellamassa(mezzobifase)

L’equazione di continuità del mezzo bifase si ottiene sommando le eqq. (6) e (14):

1 (1 ) 0

s w sw sn v n v n vn v vn n n nz z t t z z t t

quindi alla fine:

w s s

n v v vz z

(15)

Introducendo la velocità di filtrazione, vf, che è la velocità relativa del flusso del fluido che entra nella faccia dx dy rispetto al flusso del terreno nella medesima faccia:



‐ 7 ‐

f w sv n v v

si ottiene infine:

f sv vz z

(16)

2.2.4 LeggediDarcy

La legge di Darcy in condizioni 1D è data da:

f zhv kz

dove kz è la conducibilità idraulica (permeabilità) nella direzione z.

Essendo h pari a:

0 0derivando per

w w w w w

u uu u h uh zz z z

sostituendo:

0

f zw w

u uv kz z

Se il fluido nelle condizioni iniziali è in regime idrostatico si ha:

0 0

w

uz

quindi alla fine:

z

fw

k uvz (17)

2.2.5 Legamecostitutivo

Per la (ip.6) le deformazioni lungo la direzione z, dz, sono espresse da:

1

z zed

d dE

(18)



‐ 8 ‐

dove Eed è il modulo di rigidezza edometrica (essendo in condizioni di deformazioni 1D, cioè di deformazione laterale impedita)

2.2.6 Equilibriodelloscheletrosolido

Per la (ip.8) le tensioni totali, σz, agenti sull’elemento di volume, dopo l’applicazione del carico rimangono costanti nel tempo ovvero:

0zt

(19)

2.3 Equazionefinale

Derivando l’eq. (17) rispetto a z si ottiene:

f z

w

v k uz z z

Per le (ip.2) e (ip.7), w e kz non sono variabili con z, inoltre, per la (ip.7), kz non varia nel tempo, pertanto:

2

2

f z

w

v k uz z

Sostituendo il primo membro con l’equazione di continuità della massa per mezzo bifase (eq.(16));

2

2

s z

w

v k uz z

(20)

Essendo per definizione:

s s s zv s sz z t t z t

(21)

dove con ss si è indicato lo spostamento e con εz la deformazione (il segno negativo è per convenzione: una variazione positiva dello spostamento corrisponde ad un allungamento).

Sostituendo la eq. (21) nelle (20) si ottiene:

2

2

z z

w

k ut z

(22)

Introducendo il legame costitutivo (eq. (18)), la deformazione εz si può riscrivere convenientemente come:



‐ 9 ‐

1

zz

edt E t (23)

Dall’equilibrio dello scheletro solido (eq. (19)) si ottiene inoltre:

0 0

z z zu u ut t t t

(24)

Sostituendo le eqq. (23) e (24) nella (22) si ottiene:

2 2

2 21

z edz

ed w w

k Eku u u uE t z t z

ovvero:

2

2

vu uct z

(25)

dove cv è il coefficiente di consolidazione verticale, definito come:

2

z edvw

k E LcT

2.4 Integrazionedellaequazionedelmodello

L’eq. (25) è una equazione alle derivate parziali di tipo parabolico (possiede quindi una famiglia di curve caratteristiche). Può essere facilmente integrata analiticamente in condizioni monodimensionali ad esempio per il dominio di integrazione mostrato nella Fig. 3, dove H prende il nome di “percorso di filtrazione” e rappresenta la massima distanza che una particella di acqua deve percorrere per raggiungere una delle superfici drenanti (nel caso in Fig. 3a, H è pari a metà dello spessore dello strato di terreno che consolida essendoci 2 superfici drenanti, nell’altro caso ‐singola superficie drenante ‐ H coincide con l’intero spessore del banco).

Il dominio di integrazione è in questo caso è: t > 0; 0 z 2H.

Le condizioni al contorno sui bordi del dominio (Γ0, Γ1, Γ2) nel caso in cui la falda sia inizialmente idrostatica e ci sia una distribuzione costante delle ū(z, 0) pari all’incremento delle tensioni totali σz (distribuzione “rettangolare” di ū(z, 0)):

0: ū(z, 0) = σz = q Condizione iniziale 1: ū(0, t) = 0 Condizione di drenaggio libero

(interfaccia superiore) 2: ū(2H, t) = 0 Condizione di drenaggio libero

(interfaccia inferiore)



‐ 10 ‐

(a) (b)

Figura 3: dominio di integrazione dell’equazione di consolidazione 1D

2.4.1 Adimensionalizzazione

Per poter risolvere in maniera analitica l’eq. (25) è conveniente effettuare una trasformazione delle variabili indipendenti:

2

(fattore spessore)

(fattore tempo)

v

zZHcT tH

l’eq. (25) si riscrive come problema di Cauchy in “grande” (trattazione “in grande”):

2

2

,0 0 2; 00, 0 0; 02, 0 2; 0

z

u uT Z

u Z Z Tu T Z Tu T Z T

(26)

2.4.2 Soluzioneanalitica(APPROFONDIMENTO)

Una generica soluzione per l’eq. (26) è del tipo:

, u Z T Z T (27)

ovvero:



‐ 11 ‐

22

2 2

quindi:

Z Zu T TZ Z Z Z

Tu TT T

2 2

2 2

1 1separando le variabili T TZ Z

Essendo il primo membro funzione solo di T e il secondo solo di Z, come soluzione di primo tentativo si otterrà una costante :

2

22

2

1 ( )

1 ( )

aT

bZ

(28)

L’equazione (28)a è a variabili separabili si può integrare facilmente:

22 2 2ln ln

TdT T T T e

dove ln (costanti di integrazione).

L’equazione (28)b è quella dell’oscillazione armonico (non smorzato), la soluzione è una combinazione lineare di armoniche semplici del tipo:

1 2cos sin Z C Z C Z

Quindi una soluzione “particolare” è il prodotto tra le due:

2sin cos Tu A Z B Z e dove: 1A C e 2B C (29)

Più in generale, una famiglia di soluzioni è data dalla combinazione lineare di funzioni del tipo (29) ovvero:

21

sin cos

mTm m m mm

u A Z B Z e

Le costanti Am, Bm, m si calcolano a partire dalle condizioni al contorno

20, 0 0; 0 0 0 mTm mu T Z T B e B



‐ 12 ‐

2

2, 0 2; 0 sin 2 0

sin 2 02

m Tm m

m m

u T Z T A em

1

,0 0 2; 0 sin2

z m zm

u Z Z T A m Z

moltiplicando ambo i membri per sin 2 j Z con j indice generico, e integrando nell’intervallo 0 Z 2, si ottiene:

2 2

0 01

sin sin sin2 2 2z mmj Z dZ A j Z m Z dZ

Indicando con I1 e I2 gli integrali:

2

1 0

0 2 (indice pari)sin 4 2 1 (indice dispari)2

j m

I j Z dZj m

j

2

2 0

0sin sin

12 2

j m

I j Z m Z dZj m

si ottengono i coefficienti Am:

1

21

0 2 (indice pari)4 2 1 (indice dispari)

zm

m

j mIA

j mI j

La soluzione finale è:

2 22 1

4

1

4 1 2 1sin2 1 2

mT

zm

mu Z em

(30)

dove (2m‐1) è un intero dispari.

2.4.3 Soluzionegrafica

La soluzione analitica del problema di Cauchy in “grande” (26), nel caso in cui ū(z,0) = Δσz è del tipo:



‐ 13 ‐

2

2(2 1)4

1

4 1 2 1sin2 1 2

m T

zm

mu Z em

(31)

in cui compare una sommatoria (convergente) di termini. Al fine di rendere agevole l’uso della eq. (31), si preferisce fare ricorso alla soluzione grafica attraverso l’abaco riportato in Fig. 4, in cui la soluzione è espressa attraverso la variabile adimensionale U definita come “grado di consolidazione”:

0

1 1

v

v v

u uUu

che rappresenta l’incremento delle tensioni efficaci, rapportato all’incremento delle tensioni totali (costante), in funzione della posizione relativa, Z, e al variare del fattore tempo, T (isocrone).

L’evoluzione nel tempo del processo di consolidazione dell’intero banco di terreno di spessore 2H

è espresso attraverso la variabile adimensionale “grado di consolidazione medio”, U , definita come:

2 2

0 02

0

, ,1 1

2

H H

Hv

v

u z t dz u z t dzU

Hdz

Figura 4: Abaco di consolidazione.



‐ 14 ‐

Figura 5: definizione del grado di consolidazione medio del banco consolidante.

ovvero con riferimento alla Fig. 5:

area sottesa dall'isocrona (t) area campita1area rettangolo 2 area totale

v

UH

Il valore di U in funzione di T, è riportato nell’abaco di Fig. 6 e rappresenta l’evoluzione nel tempo del processo di consolidazione del banco di terreno.

Tenuto conto della linearità del legame costitutivo:

2 2 2

0 0 0

1 1 1, 2 , H H H

z z zed ed ed

w t dz u z t dz H u z t dzE E E

Figura 6. Variazione del grado di consolidazione medio con il fattore tempo T.



‐ 15 ‐

2 2

,0 0

1 1 1 2H H

f z fin z zed ed ed

w dz dz HE E E

quindi alla fine:

22

00

1 2 ( , ) ( , )( ) 1 (c.v.d.)1 22

H

Hz

ed

f zz

ed

H u z t dz u z t dzEw t Uw HH

E

il Grado di Consolidazione Medio corrisponde al cedimento al generico istante di tempo t rapportato al cedimento finale.

La soluzione (31) e gli abachi delle Figg. 4 e 6 valgono anche nel caso in cui la superficie drenante sia solamente una, come mostrato nella Fig. 7, a parità di ipotesi iniziali (isocrona iniziale rettangolare ovvero nel caso in cui l’incremento delle tensioni totali, σz, ‐ quindi quello delle pressioni interstiziali, ū0, al tempo t = 0 ‐ è costante con la profondità).

Figura 7: Corrispondenza tra le soluzioni con una o due superfici drenanti a parità di H.



‐ 16 ‐

Figura 8: Costruzione grafica per la determinazione del coefficiente di consolidazione.

2.4.4 Valutazionesperimentaledelcoefficientediconsolidazioneverticale

Il coefficiente di consolidazione, cv, è determinato sperimentalmente sulla base delle curve cedimento‐tempo (in scala semilogarilmica) ottenute attraverso prova edometrica, interpretate alla luce della teoria di Terzaghi (Fig. 8).

Un metodo tipicamente utilizzato fa riferimento alla procedura grafica di Casagrande (Fig. 8):

1) si stima il cedimento iniziale, w0, valutando la differenza, Δw tra il valore del cedimento misurato al tempo t, w1(t), e quello misurato al tempo 4t, w2(4 t), nelle prime letture della prova;

2) il cedimento w0 si assume pari a w1 – Δw, e individua l’inizio della consolidazione (U = 0%);

3) si individua il punto di flesso della curva e si traccia la retta tangente a tale punto; 4) si valuta l’asintoto per t (per la teoria della consolidazione, tale asintoto dovrebbe

essere orizzontale: in realtà si osserva sperimentalmente un andamento lineare del cedimento dovuto a fenomeni viscosi e/o di creep che prendono il nome di “consolidazione secondaria”);

5) si individua il valore del cedimento, w3, corrispondente alla fine della consolidazione (U = 100%) come l’intersezione tra la tangente al punto di flesso e l’asintoto (obliquo): il cedimento di consolidazione wc è dato da w3‐w0;

6) si legge da grafico il valore dell’istante di tempo t50, in corrispondenza del cedimento, w4 = w0+wc/2, che identifica il valore di w quando U = 50%;

(U=0%)

t 4t

w1w2

ww

w0

wc

w3

t50

w4 (U=50%)

(U=100%)

log t

w



‐ 17 ‐

7) dalla curva U in funzione del fattore tempo di Fig. 6, fissato 50%U si ottiene T50 da cui:

2 2

50 5050 2

50 50

0.197 v vc t T H HT cH t t

dove H è il “percorso di filtrazione” o che nel caso di presenza due superfici drenanti (come nella prova edometrica) si può assumere paria a metà dell’altezza iniziale del provino, h0 (in genere di 20 mm).

il valore di cv varia da 10‐4 10‐3 cm2/s per le argille normalconsolidate e aumenta al crescere di OCR.

Teoria della consolidazione 1D - people.unica.it · 2019. 4. 15. · 2 Teoria di Terzaghi...

Documents

Transcript of Teoria della consolidazione 1D - people.unica.it · 2019. 4. 15. · 2 Teoria di Terzaghi...