Appunti di trasporto solido

appunti di
TRASPORTO SOLIDO

Estratto da appunti
a cura di: prof. Silvio Franzetti
prof. Alberto Guadagnini ing.
Francesco Ballio
Appunti di trasporto solido

1 STABILITÀ E MOVIMENTO DEI GRANELLI

1.1 INTRODUZIONE
Nel calcolo delle correnti a pelo libero la geometria dell'alveo (planimetria, altimetria, sezioni)
è solitamente considerata come un dato di ingresso del problema, nota sulla base di rilievi
topografici. In effetti, nel caso generale in cui l'alveo sia scavato in materiale incoerente, la
geometria non può essere considerata a priori fissa nel tempo, a causa dei fenomeni di erosione
e deposito del materiale costitutivo del fondo e delle sponde, mosso dalla corrente. Sotto
l'azione delle spinte idrodinamiche, infatti, i grani solidi possono essere messi in movimento e
trasportati verso valle; si differenziano inoltre il fenomeno del trasporto di fondo, in cui i
granelli si muovono sul fondo, con moti rotatori e/o di strisciamento, più o meno intermittenti,
dal fenomeno del trasporto in sospensione, in cui il solido percorre lunghi tratti trascinato
all'interno della corrente, e solidale ad essa.
In conseguenza dei processi descritti l'alveo di un corso d'acqua può subire mutamenti
progressivi nel tempo, in relazione alle sue caratteristiche planimetriche, e quindi al suo
percorso, altimetriche (pendenze), alla forma delle sezioni e alla distribuzione granulometrica
del materiale di fondo, cui sono legati i valori di scabrezza. Le modifiche dell'alveo possono
essere conseguenza di eventi limitati nel tempo (piene, alluvioni), ovvero al regime normale del
corso d'acqua: in ogni caso, terminato il transitorio evolutivo, si possono avere situazioni di
equilibrio statico, in cui i fenomeni di trasporto solido sono assenti, ovvero di equilibrio
dinamico tra erosione e deposito del materiale movimentato.
E' inoltre opportuno distinguere i fenomeni distribuiti, che coinvolgono l'intera sezione d'alveo
su tratti estesi longitudinalmente, dai fenomeni concentrati, limitati a zone specifiche della
sezione, ovvero a tratti limitati dell'alveo, e legati alla presenza di singolarità geometriche nel
corso d'acqua (curve, opere idrauliche, ...). Nel seguito si considereranno essenzialmente i
fenomeni di tipo distribuito, accennando solamente a quelli concentrati in un paragrafo finale.

Il fenomeno di movimentazione del fondo è essenzialmente legato alle:

1. caratteristiche della corrente: date la geometria e la portata vengono calcolate mediante le
usuali metodologie dell'idraulica;
2. caratteristiche dei sedimenti, ed in particolare la distribuzione granulometrica e le
caratteristiche di coesività. I fenomeni di trasporto solido per sedimenti coesivi richiedono
una trattazione specifica, e non vengono qui considerati.

1.2 SCHEMI CONCETTUALI DI TRASPORTO - INCIPIENTE MOVIMENTO

Da un punto di vista fenomenologico, il moto di fondo delle particelle solide a causa delle forze
esercitate dalla corrente può essere inquadrato con diversi schemi concettuali. In ogni caso si
deve valutare la stabilità della particella di fronte alle forze attive di trascinamento della
corrente e alle forze resistive dovute al peso immerso della particella medesima ed
eventualmente al conseguente attrito sul fondo. A parità di ogni altro parametro, all'aumentare
della velocità le azioni fluidodinamiche riescono a vincere le azioni stabilizzanti, mettendo in
movimento la particella. Le condizioni di flusso in cui la particella si trova al limite del
movimento vengono indicate come "critiche", o di incipiente movimento.
Appunti di trasporto solido

Strisciamento (attrito)

La particella striscia sul fondo; le forze

idrodinamiche devono superare la reazione
vincolare di attrito radente.

Ribaltamento (senza attrito)

La particella è parzialmente incastrata sul fondo; le

forze idrodinamiche devono farla ruotare rispetto
al punto di appoggio.

"Saltelli"

La particella è parzialmente incastrata sul fondo; le

forze idrodinamiche devono sollevarla rispetto al
fondo.

Condizione di incipiente movimento, schema di attrito radente

Sulla base dello schema di attrito radente si cercano le condizioni, e quindi i legami fra i
parametri del sistema, in condizione di incipiente movimento. In figura sono indicate le forze
che agiscono sulla particella.
Appunti di trasporto solido

L = Lift (portanza)
Umax D = Drag (resistenza)
G = spinta di Archimede
P = peso proprio
u(ζ)
δ R = reazione vincolare normale
A = reazione vincolare tangenziale
ζ ud G L (Attrito)

D d = dimensione del “granello”

A θ
R P

In condizioni di equilibrio statico le reazioni vincolari devono valere:

R = (P - G) cosθ - L
A = (P - G) sinθ + D

Il massimo valore della forza tangenziale è: Amax = R tanβ (β = angolo di attrito).

La condizione di equilibrio limite (incipiente movimento) viene raggiunta per A = Amax
I termini possono essere esplicitati come:
D = (0.5 ρ ud2) CD S
L = (0.5 ρ ud2) CL S = KL D
P - G = (ρs - ρ) g W
dove: ud velocità della corrente (vedi successivamente)
S Superficie frontale del granello: S = α1 d2
W volume del granello: W = α3 d3
CD, CL coefficienti di Drag, Lift: dipendono dalla forma del granello e dal numero di
Reynolds; in prima approssimazione si assume costante il rapporto L/D = KL
I coefficienti αi sono legati alla forma dei granelli, così come CD e CL, e si considerano costanti
e caratteristici della granulometria dei sedimenti, benché a priori incogniti.
Imponendo che A = Amax (limite di attrito, incipiente movimento) si ottiene:
(P - G) sinθ + D = [(P - G) cosθ - L] tanβ ⇒ D + L tanβ = (P - G) cosθ [tanβ - tanθ]
Sostituendo le espressioni per D, L, P-G:
(1 + KL tanβ) CD α1 d2 (0.5 ρ ud2) = (ρs - ρ) g α3 d3 cosθ [tanβ - tanθ]
2
2α
ud,c 3
= tanβ cosθ (1 - tanθ / tanβ)
g [(ρs - ρ)/ρ] d (1 + KL tanβ) CD α1

La velocità ud rappresenta la velocità efficace agli effetti delle forze fluidodinamiche, e può
essere considerata pari alla velocità della corrente ad una opportuna distanza dal fondo,
proporzionale alla dimensione del granello. Il pedice "c" si riferisce alle condizioni critiche.
Sia:
Appunti di trasporto solido

ud = u(ζ = α4 d)
Esprimendo il profilo di velocità (strato limite) secondo una forma di similitudine del tipo
u(ζ) / u* = f(ζ/d) (paragrafo 1.2) si ottiene evidentemente ud/u* = f(α4) = α5.
Tramite la forma del profilo di velocità si è riusciti a legare le condizioni idrodinamiche alla
velocità u*, rappresentante degli sforzi sul fondo - noti o facilmente calcolabili per una
corrente, anziché alla meno pratica ud. La relazione di incipiente movimento risulta pertanto:

uc*2 2 α3 tanβ
= cosθ (1 - tanθ / tanβ)
g [(ρs - ρ)/ρ] d (1 + KL tanβ) CD α1 α52

La frazione del termine di destra dipende, oltre che dalla granulometria, dal numero di
Reynolds tramite CD, KL, ed eventualmente α5. In sintesi risulta:

uc*2
φc = = f(Re) cosθ (1 - tanθ / tanβ)
g [(ρs - ρ)/ρ] d

Il termine adimensionale:
u *2  ρ − ρ
φ= = τ/ ρ ∆ = s
g∆d g∆d  ρ 
è denominato numero indice di Shields, e può essere interpretato come il rapporto fra le forze
di trascinamento della corrente (proporzionali a τ d2) e le forze stabilizzanti (peso
immerso, proporzionale a g(ρs-ρ)d3). Date le caratteristiche della corrente (e quindi ρ, τ) e
quelle dei sedimenti (∆, ρs, d), per valutare la stabilità delle particelle si deve confrontare il
valore di φ con la soglia critica individuata, in cui deve essere ancora definita la funzione
incognita f( ). Si avrà in sostanza:
per φ < φc granello stabile
per φ > φc granello in movimento
Il valore di φ per la corrente si calcola sulla base degli sforzi sul fondo, ricordando (paragrafo
1.2) che questi ultimi sono legati alla cadente:
τ/ ρ RJ
φ= = [ τ = γRJ ] =
g∆d ∆d
Si noti che il valore di φ è legato alla dimensione d dei granuli: il fenomeno di trasporto è
selettivo, agendo diversamente sui granelli di diversa dimensione.

Tipici valori per la densità dei materiali rispetto alla densità dell'acqua

Materiale ρs/ρ ∆
Quarzo (sabbia) 2.6 ÷ 2.7 1.6 ÷ 1.7
Argilla 2.6 ÷ 2.8 1.6 ÷ 1.8
Basalto 2.7 ÷ 2.9 1.7 ÷ 1.9
Carbone 1.3 ÷ 1.4 0.3 ÷ 0.4
Plastica 1.0 ÷ 1.5 0.0 ÷ 0.5
Appunti di trasporto solido

Condizione di incipiente movimento, schema di rotolamento

L = Lift (portanza)
Umax D = Drag (resistenza)
T = L + D (vettoriale!)
u(ζ) G = spinta di Archimede
δ P = peso proprio
ζ G L
ud d = dimensione del “granello”

D
P θ

Si impone l’equilibrio alla rotazione rispetto al punto C di contatto. Siano bT e bGP i bracci di rotazione
di T e G + P rispetto a tale punto. Risulta, in condizione limite di ribaltamento:
(P - G) bGP = T bT
I termini bx sono espressi medianti opportuni coefficienti legati alla forma e disposizione dei granelli:
bD = β 1 d bL = β 1 d bGP = β3 d
Esplicitate i termini D, L, P, G analogamente a quanto precedentemente descritto, e sostituendo le
espressioni nella relazione di equilibrio alla rotazione, si ottiene:
(ρs - ρ) g α3 d3 bGP = CD α1 d2 (0.5 ρ ud2) (1 + K L2)0.5 bT
da cui:
u d2 2 α3 bGP
=
g [(ρs - ρ)/ρ] d CD α1 bT (1 + K L2)0.5

La relazione è sostanzialmente analoga a quella ottenuta dallo schema di attrito radente; la dipendenza
da θ entra nel braccio bGP (bGP = 0 per θ = 0), mentre bT è, in linea di principio, indipendente da θ.

1.3 LEGGI DI INCIPIENTE MOVIMENTO / 1 (abaco di Shields)

Nel precedente paragrafo si è visto come, sulla base di opportuni schemi fenomenologici, la
stabilità di un granello di fondo può essere valutata sulla base del confronto dello sforzo
tangenziale τ (opportunamente normalizzato) esercitato dalla corrente sul fondo rispetto ad
una soglia critica di cui rimane da determinare la dipendenza dal numero di Reynolds (ed
eventualmente dalla granulometria dei sedimenti).
Ci si ponga sotto le condizioni:
1. i=0 (θ = 0)
2. granulometria uniforme (d = cost)
I valori per la soglia critica sono stati studiati sperimentalmente da Shields (1936) in funzione
del numero di Reynolds Re* = u*d/ν, ottenendo l'andamento rappresentato in figura (l'apice
"0" apposto a φc sta ad indicare la condizione i = 0). In altre parole, i risultati della figura
rappresentano la definizione sperimentale della funzione f( ), la cui dipendenza dal numero di
Reynolds non era stata specificata.
Appunti di trasporto solido

φ
oc

granelli in movimento

0.1

granelli stabili
0.01
1.E-1 1.E+0 1.E+1 1.E+2 1.E+3 1.E+4
Re*

Abaco di Shields

In relazione ai risultati dell'abaco di Shields si possono fare alcune osservazioni.

1) La soglia φc di incipiente movimento risulta dipendere dal numero di Reynolds soprattutto ai
valori inferiori di Re*. Agli alti numeri di Reynolds il fenomeno risulta di fatto indipendente
da Re* (la viscosità non ha più alcun effetto).
Peraltro, nella maggior parte dei casi reali, si ha Re* > 100, e conseguentemente è lecito
assumere φ*c = cost. Si vedano a tal proposito i due esempi sotto riportati.

Canale Fiume di montagna

d = 0.001 m d = 0.01 m
Ks = 50 Ks = 25
i = 0.1 % i=1%

1:2 1:2
5m 3m

Q = 30 m3/s Q = 145 m3/s

⇓ ⇓
ho = 2.4 m U = 2.0 m/s ho = 2.5 m U = 5.2 m/s
R = 1.4 m u* = 0.12 m/s R = 1.8 m u* = 0.59 m/s
C = 53 m1/2/s C = 28 m1/2/s
⇒ Re* = 120 ⇒ Re* = 5900

2) Nel campo ai maggiori numeri di Reynolds, la maggior parte dei punti sperimentali è
compresa in φc = 0.04÷0.06, ma si arriva anche a φc = 0.1. Scarti del 100% !
3) Molteplici cause giustificano la notevole dispersione dei punti:
Appunti di trasporto solido

• Sono stati di fatto trascurati i parametri caratteristici della forma e della disposizione dei
sedimenti (cfr. modello fenomenologico); tali parametri sono in effetti non determinabili,
e del tutto casuali per il singolo granello, sebbene si possa pensare di identificarne valori
"efficaci" rappresentativi dell'insieme dei sedimenti.
• Imperfetta similitudine dei profili di velocità, implicitamente contenuti nel confronto fra φ
e φc.
• Non univocità del criterio per identificare le condizioni di incipiente movimento, da
valutarsi sulla base del moto di una quantità finita di granelli.
Nel seguito sono discusse alcune estensioni ai risultati di Shields, per condizioni al di fuori dei
limiti sperimentali relative ai dati sopra discussi.

Granulometria non uniforme

Se i sedimenti non hanno dimensione costante si ritiene di potere confrontare il parametro φ di
mobilità con il valore di soglia all’incipiente movimento selettivamente per ogni frazione
granulometrica:
φ(d) = RJ / ∆d ↔ φc(Re*)
Data una condizione di flusso (e quindi fissati i valori di R e J) la
mobilità dei sedimenti è inversamente correlata a d; esiste pertanto
un valore minimo per d al di sotto del quale i granelli vengono
trasportati dalla corrente.
A causa del fenomeno descritto si determina una situazione di
erosione selettiva dei sedimenti a minor diametro, e di conseguenza
una variazione della distribuzione granulometrica nello strato ⇓
superficiale, caratterizzato da un aumento della dimensione media
dei granelli; si è soliti definire tale fenomeno con il termine
corazzamento, ad indicare una aumentata resistenza superficiale del
fondo all'erosione. Si noti peraltro che, qualora lo strato superiore
di sedimenti venga rimosso, la resistenza all'erosione cala
bruscamente, ritornando esposti alle azioni idrodinamiche gli strati
inferiori caratterizzati dalla granulometria originale.
Oltre al fenomeno del corazzamento si deve segnalare che, volendo applicare il criterio di
incipiente movimento selettivamente per ogni frazione granulometrica d, la soglia di mobilità
deve essere in effetti corretta sulla base del rapporto d/d50: la presenza di particelle di elevate
dimensioni tende a “proteggere” quelle più piccole, e viceversa una particella grossa isolata fra
particelle di dimensioni minori risulta più esposta che nel caso di granulometria uniforme. Per
la definizione quantitativa del problema si rimanda ai testi specializzati.

Correzione per pendenza diversa da 0

Valori di pendenza del fondo diversi da zero riducono la stabilità dei granelli, riducendosi
l'effetto stabilizzante della gravità. Nel caso limite in cui θ = β (inclinazione del fondo pari
all'angolo di natural declivio dei sedimenti) risulta φc = 0 (i granelli si muovono per φ = 0
ovverosia per τ = 0).
Il correttivo da applicarsi a φo c per tenere conto della pendenza è
Appunti di trasporto solido


φ = φ o cos 1 − tan θ 
c c
θ 
 tan β 
L'andamento di tale fattore correttivo è diagrammato in figura.

i = 0.01 i = 0.1 i = 1.
↓ ↓ ↓
1.00

φc/φ c
o
0.75

0.50

0.25

0.00
0.10 1.0 10. 100.

θ [gradi]

Correttivo per la pendenza del fondo (β = 45°)

Fattori correttivi differenti dalla (1.4.1) diversi sono proposti in letteratura, sulla base di schemi
concettuali diversi (si veda il paragrafo 1.7); gli andamenti risultano comunque qualitativamente
e quantitativamente paragonabili a quello di figura 3.

h/d ridotti (basso grado di sommergenza)

I risultati dell'abaco di Shields valgono se le dimensioni dei granelli sono piccole rispetto a
quelle della corrente, ovverosia per valori non piccoli del grado di sommergenza h/d. Per
h/d < 6 la soglia di incipiente equilibrio deve essere corretta, risultando questa inversamente
correlata con h/d. Tale effetto, non evidenziato nel modello concettuale proposto nel paragrafo
precedente, è dovuto alla deformazione dei profili di velocità: il legame ud ↔u* risulta
diverso per h/d ridotti rispetto alle leggi "universali" valide per i maggiori rapporti di
sommergenza, e quindi risulta differente il legame fra gli sforzi sul fondo (ovvero le perdite di
carico) e la mobilità dei granelli; si riscontra, di fatto, che questa si riduce progressivamente al
diminuire di h/d. Si noti che bassi valori di h/d sono tipici dei torrenti montani, caratterizzati da
granulometrie non fini e pendenze elevate, che inducono limitate altezze d'acqua.
Si riportano di seguito due correttivi proposti in letteratura per tenere conto degli effetti sopra
discussi; i relativi andamenti sono confrontati in figura.
−0.266
φc = 0.0851  h 
o
Suszka (1991)
 d50 
o  d
Armanini e Scotton (1994) φc = 0.06 1+ 0.67  (1.4.3)
 h
Appunti di trasporto solido

0.14 Shields
Armanini e Scotton
φoc Suszka

0.08
0.06
0.04
0.02
0.00
0.1 1.0 10. 100. 1000.
h/d50

Correttivi per h/d ridotti

Entrambi i correttivi determinano un aumento di φo c al diminuire di h/d, in conseguenza della

Appunti di trasporto solido
diminuita mobilità dei granelli al diminuire del grado di sommergenza. E' peraltro evidente che
la formula di Suszka non è tarata per valori elevati di h/d, indicando una diminuzione indefinita
della soglia critica all'aumentare del grado di sommergenza, laddove tale soglia presenta un
asintoto (orizzontale) definito dai risultati dell'abaco di Shields.

1.4 LEGGI DI INCIPIENTE MOVIMENTO / 2

Nei paragrafi pr e c e de nt i le condizioni di incipiente movimento dei sedimenti di un alveo
sono state valutate sulla base del confronto fra lo sforzo esercitato dalla corrente sul fondo
(τ = γRJ) e del valore di soglia per un granello per essere mosso, risultando:
τc = f(ρ, ρs, g, d, ν, θ)
Le condizioni critiche possono evidentemente essere definite sulla base di grandezze diverse
dagli sforzi; in generale si può esprime una qualsiasi grandezza (U, u*, τ, i, Q, ...) in funzione
delle altre, sempre in stato di incipiente movimento, per ottenere relazioni del tipo:
Uc = ... τc = ... ic = ... qc = ...

In figura sono poste a confronto alcune delle formule e s i s t e n t i . Per permettere il

confronto, tutte le formule sono state trasformate in termini di φc; quando necessario il
legame sforzi - velocità è stato esplicitato mediante la legge di resistenza. Tutte le formule
sono riferite a sezione "larga" (in cui sia R ≅ h), in condizioni di moto uniforme, per elevati
numeri di Reynolds.
Osservazioni
1. Per i minori valori di i, ovvero elevati valori di h/d, la maggior parte delle formule portano
a valori approssimativamente costanti per φc, distribuiti all'interno della fascia di valori in
cui si disperdono i punti dell'abaco di Shields (φc = 0.03÷0.06; Re* > 70). Fanno eccezione
alcune formule che appaiono evidentemente inadeguate a stimare la soglia di incipiente
movimento per i minori di valori di pendenza.
2. Per i > 0.2, h/d < 6, ovverosia al di fuori del campo di validità dell'abaco di Shields, si
presentano andamenti contrastanti: alcune formule indicano andamenti crescenti della soglia
di mobilità, in corrispondenza della riduzione del grado di sommergenza, salvo invertire tale
andamento per le pendenze più elevate qualora sia presente un correttivo specifico per la
pendenza; per altre formule prevale senz'altro quest'ultimo effetto, cosicché la soglia di
incipiente movimento è indicata essenzialmente decrescente al crescere di i e quindi di h/d.

11
Appunti di trasporto solido

(5) (5') (2) 0.14

(2') (3) (1)
(1') (11) (12) 0.12
(14) (14') (21)
(22) campo validità di Shields
0.10 φc

0.08

0.06

0.04

0.02

0.00
1.00000 0.10000 0.01000 0.00100 0.00010 0.00001
i

0.14 (5) (5') (2)

(2') (3) (1)
(1') (11) (12)
0.12 (14) (14') (21)
φc (22) campo validità di Shields
0.10

0.08

0.06

0.04

0.02

0.00
1.E-1 1.E+0 1.E+1 1.E+2 1.E+3 1.E+4
h/d50

Figura 8: confronto fra le formule di incipiente movimento.

12
Appunti di trasporto solido

Riferimento Formula Note

Shields φc = 0.05 h/d > 6 (5)

 tan θ
Shields φc = 0.05 cos θ1 −  h/d > 6 (5’)
 tan β 

Armanini e  d
Scotton φoc = 0.06 1+ 0.67  (2)
 h

Armanini e  d  tan θ
Scotton φoc = 0.06 1+ 0.67  cos θ1−  (2’)
 h   tan β 

Graf e Suszka φ = 0.042 102.2 i i = 0.02÷0.2 (3)

c

 −0.266 h/d = 0.5÷50 (1)

Suszka φoc = 0.0851  h 
 d50 

−0 .266
Suszka φco = 0.0851  h  cos θ tan β − ρ s tan θ h/d = 0.5÷50 (1’)
 d 50   ρs − ρ 

Bettes qc 0.134 1.221 i = 0.002÷0.1 (21)

= ln (h/d = 1÷100)
gd 3 i i
qc
Bathurst et al. = 0.15 i −1.12 i = 0.002÷0.1 (22)
3
gd (h/d = 1÷100)

13
Appunti di trasporto solido

1.5 TRASPORTO DI FONDO

Fino ad ora si è affrontato il problema di determinare la soglia di movimento incipiente dei
sedimenti in funzione delle caratteristiche della sovrastante corrente fluida. Il passo successivo
è necessariamente quello di determinare, nel caso tale soglia sia superata, la quantità di solido
movimentato. Il fenomeno del trasporto solido (di fondo) può essere schematizzato secondo
diversi approcci, che non devono essere interpretati come escludentesi a vicenda, bensì
complementari nel descrivere una fenomenologia complessa.

"Strato limite" di solido in moto

Trascinato dagli sforzi tangenziali, uno spessore più o meno

profondo di particelle solide viene messo in movimento; Q
all'aumentare della profondità gli strati di solido sono
caratterizzati da velocità progressivamente decrescenti.

Moto superficiale

Solo uno strato superficiale di solido è considerato in

movimento; le particelle rotolano / strisciano sulle sottostanti, Q
ferme.

"Saltelli"

Il moto delle particelle non è schematizzato come fenomeno Q

continuo, bensì come successione di "salti" finiti delle
particelle superficiali, seguiti da tempi di riposo.

A partire dai diversi schemi concettuali si determinano varie formule che legano la quantità di
solido trasportata ai parametri del flusso. Non ci si addentra in alcuna di tali analisi; basti dire
che tutte conducono a formule del tipo:

14
Appunti di trasporto solido

f(U/Uc) ; f(U - Uc)

qs = f(i - ic)
f(τ/τc) ; f(τ - τc)
..... analoghe
in cui la portata solida qs (portata volumetrica di solido per unità di larghezza, vuoti esclusi) è
legata all'eccedenza fra un parametro caratteristico della corrente e la corrispondente soglia di
incipiente movimento (eccedenza espressa come differenza o come rapporto). Ovviamente
deve risultare qs = 0 in condizioni di incipiente movimento. La portata volumetrica effettiva,
comprensiva dei vuoti, si ottiene dividendo qs per (1-n), essendo n la porosità del materiale.
La portata solida qs viene generalmente adimensionalizzata secondo il numero indice Φ:
q
Φ= 3
s

d g∆
Nella successiva tabella sono riportate alcune formule per il trasporto di fondo. Le formule
devono essere considerate valide solo al di sopra della soglia critica.

i1.5 ( q − q )
2.5
Schocklitsch (1962) q = (A1)
s
ρs / ρ c

Φ= 8 ( φ* −φc ) 1.5
o
Bathurst et al. (1987)
(Meyer-Peter e Müller modificato) ρu*2 (B1)
φ* =
γ∆d cos θ (tan β − tan θ)

Suszka (1991)  φc 
Φ= 10.4 φ3/ 21−
o 5/2
(C1)
 φ

Meyer-Peter e Müller (1948) Φ = 13.3 ( φ − φ )c 1.5 (D1)

qs = d τ 2  τ 1/6  τ 
5/3

ρ 3  ρ  τ 
Pezzoli (1979)    −1 (E1)
c

τ = (ρs − ρ)g d φ

La formula più diffusa è probabilmente la (D1); la (B1) ne rappresenta una variante

specificamente tarata per pendenze non ridotte; anche la (A1) tiene esplicitamente in conto la
dipendenza dalla pendenza del fondo.
Secondo quanto discusso, una formula di trasporto deve sempre essere associata ad una
valutazione delle condizioni di incipiente movimento; negli esperimenti di laboratorio tale
soglia viene misurata; qualora si vogliano utilizzare i risultati in applicazioni pratiche,
raramente si è in grado di misurare la soglia critica, che deve essere pertanto valutata mediante
una delle formule discusse in precedenza. E' allora chiaro che alla inevitabile dispersione dei
risultati relativi alle formule di trasporto si assomma la dispersione dei risultati per l'incipiente
movimento. Per questa ragione, benché la scelta della formula di incipiente movimento da
abbinare ad una formula per il trasporto sia in linea di principio del tutto libera, è opportuno,
quando possibile, utilizzare formule fra loro omogenee dal punto di vista concettuale e della

15
Appunti di trasporto solido

fonte (ovvero dei dati su cui sono state tarate), al fine di minimizzare la dispersione del
risultato. Sono di seguito elencati alcuni possibili abbinamenti per le formule precedentemente
proposte.
qc
= 0.15 i −1.12
Schocklitsch → (A2)
gd 3

Bathurst et al. → φoc = 0.047 (B2)

−0.266
Suszka → φoc = 0.0851 h  (C2)
d 
 d
φoc = 0.06 1+ 0.67 
 h
Meyer-Peter e Müller, Pezzoli →  tan θ
(D2) , (E2)
φ c = φoc cos θ1 − 
 tan β

Le formule di trasporto sono in generale determinate per pendenze trascurabili o nulle, e

rapporti di sommergenza elevati. Se non specificamente indicati (formula (B1)) correttivi per
pendenze elevate e ridotti valori di h/d sono da apportarsi non tanto alla formula di trasporto,
quanto al termine di soglia.
In figura 9 le formule sono poste a confronto per alcuni valori significativi dei parametri, per il
caso di sezioni "larghe" e condizioni di moto uniforme per la corrente; le perdite di carico sono
state valutate tramite la formula di Chézy sulla base dei coefficienti di Strickler indicati in
figura. Le portate solide qs e liquide q sono indicate per unità di larghezza; per le diverse
portate è anche indicata la corrispondente altezza di moto uniforme (asse di destra). Sebbene i
risultati presentati non abbiano carattere di generalità, essendo riferiti a valori specifici per i
parametri, sono comunque significative alcune osservazioni.
1. La dispersione dei risultati è estremamente elevata. Al di là dell'ovvio trend crescente di qs
all'aumentare di q le formule danno valori e gradienti alquanto dispersi, e variabili
relativamente le une alle altre in funzione dei parametri del sistema.
2. Nel confrontare i risultati è opportuno considerare le differenze delle soglie di incipiente
movimento, quando visibili. Parte della dispersione delle curve è infatti da imputarsi ad una
diversa valutazione delle condizioni critiche.
3. Si noti in particolare che per il caso "a" (canale, d = 0.1 mm, i = 0.1%) si ha trasporto anche
per i valori minori della portata (per q = 1 m2/s risulta ho = 0.8m, cui corrisponde una
velocità media U = 1.25 m/s). Solo per valori estremamente ridotti della pendenza (e quindi
elevate altezze d'acqua e velocità ridotte) si riesce ad ottenere stabilità dei sedimenti (caso
"b").
4. Analoga osservazione può essere svolta per i casi "c" e "d", in cui i sedimenti di dimensione
inferiore (d = 1 cm) sono mossi anche per valori relativamente ridotti della portata, mentre
solo sedimenti di dimensione dell'ordine della decina di cm risultano relativamente stabili. Si
notino i diversi ordini di grandezza per qs per il canale e per il fiume di montagna.

16
Appunti di trasporto solido

Canale: confronto delle formule di trasporto per due valori di pendenza

Trasporto solido - d = 0.001 m i = 0.001 Ks = 50

3.0E-3 2.50
(D1)+(D2)
(B1)+(B2)
2.5E-3 (A1)+(A2)
(C1)+(C2) 2.00
(E1)+(E2)
2.0E-3 ho

qs 1.50
[m2/s] ho
1.5E-3
[m]

1.00
1.0E-3

0.50
0.5E-3

0. 0.00
0.00 1.00 2.00 3.00 4.00 5.00 6.00 7.00

q [m2/s]

Trasporto solido - d = 0.001 m i = 0.0001 Ks = 50

6.E-5 2.50
(D1)+(D2)
(B1)+(B2)
5.E-5 (A1)+(A2)
(C1)+(C2) 2.00
(E1)+(E2)
4.E-5 ho

qs 1.50
[m2/s] ho
3.E-5
[m]

1.00
2.E-5

0.50
1.E-5

0. 0.00
0.00 0.50 1.00 1.50 2.00 2.50

q [m2/s]

17
Appunti di trasporto solido

Fiume di montagna:
confronto delle formule di trasporto per due valori di dimensione dei grani

Trasporto solido - d = 0.01 m i = 0.010 Ks = 35

9.E-2 2.50
(D1)+(D2)
8.E-2 (B1)+(B2)
(A1)+(A2)
(C1)+(C2) 2.00
7.E-2
(E1)+(E2)
6.E-2 ho

qs 1.50
2 5.E-2
[m /s] ho
[m]
4.E-2
1.00
3.E-2

2.E-2
0.50

1.E-2

0.E+0 0.00
0.00 2.00 4.00 6.00 8.00 10.00 12.00 14.00 16.00

q [m2/s]

Trasporto solido - d = 0.1 m i = 0.010 Ks = 35

3.5E-2 2.50
(D1)+(D2)
(B1)+(B2)
3.0E-2 (A1)+(A2)
(C1)+(C2) 2.00

2.5E-2 (E1)+(E2)
ho

qs 1.50
2.0E-2
[m2/s] ho
[m]
1.5.E-
1.00

1.0E-2

0.50
0.5E-2

0. 0.00
0.00 2.00 4.00 6.00 8.00 10.00 12.00 14.00 16.00

q [m2/s]

18
Appunti di trasporto solido

In figura (fonte: Bathurst et al., 1987) è riportato il confronto fra un ampio campione di
dati sperimentali di laboratorio ed il modello di trasporto (A1)+(A2). Tali formule si sono
dimostrate le migliori per il confronto con i dati sperimentali considerati, caratterizzati da
valori elevati della dimensione dei sedimenti e delle pendenze. Il confronto fra il valore
misurato per qs e il corrispondente valore calcolato mostra un accordo di massima fra formule
e valori sperimentali, ma le discrepanze quantitative possono essere anche di notevole entità,
specie per i valori minori delle portate solide (fino al 1000%)3. Ben peggiore è in tal senso il
risultato di figura 11, in cui il modello di trasporto è confrontato con dati di campo (eccesso di
portata liquida rispetto alla soglia critica in ascissa, portata solida in ordinata). Lo scarto fra i
valori misurati e la previsione del modello arriva anche a due ordini di grandezza per fiumi di
grandi dimensioni; nella valutazione del trasporto in fiumi di dimensioni ridotte l'errore
connesso al modello risulta ancora più elevato.
Si noti che le notevoli discrepanze tra modello e dati sperimentali sono da associarsi alle
condizioni estreme scelte (bassi h/d, i elevate), caratteristiche dei fiumi e torrenti di montagna,
laddove i modelli di trasporto sono normalmente concepiti e tarati per situazioni tipiche dei
fiumi e dei canali di pianura, con sedimenti di dimensione ridotta e pendenze non eccessive.

Sulla base dei risultati mostrati si potrebbe concludere che, almeno per situazioni relativamente
estreme, non si possa ottenere, con un accettabile grado di confidenza, neppure una stima
dell'ordine di grandezza del trasporto solido. In realtà non si deve rinunciare alla possibilità di
stimare fenomeni quali quelli in discussione, purché si tenga ben presente che:
1. i risultati devono essere valutati con estrema cautela, e con la coscienza che, nel migliore dei
casi, si può arrivare alla stima degli ordini di grandezza connessi al fenomeno, e non certo a
valori precisi;
2. i modelli di calcolo devono necessariamente essere calibrati ed adattati ai valori di campo
relativi al tratto fluviale e al bacino di interesse, o al limite di situazioni confrontabili nelle
vicinanze. Sulla base dei dati storici relativi all'evoluzione planimetrica e altimetrica
dell'alveo, delle caratteristiche dei sedimenti, dei volumi solidi trasportati si deve scegliere il
modello di calcolo che si riveli più adatto alla situazione specifica ed eventualmente apporre
correzioni ragionate ai parametri relativi alla corrente e ai sedimenti da inserire nel modello
medesimo. I principali elementi di scelta e/o taratura sono:
a) la granulometria dei sedimenti, le caratteristiche delle sezioni trasversali, il coefficiente di
resistenza;
b) le formule per la valutazione della soglia di incipiente movimento;
c) le formule per la valutazione dell'entità del trasporto.

3
Si sarebbe potuto ridurre la dispersione dei valori qualora si fosse utilizzato il valore misurato per la soglia di
incipiente movimento, anziché quello calcolato mediante la (A2).

19
Appunti di trasporto solido

Confronto della formula di Schocklitsch con i dati sperimentali

qs misurato
[m2/s]
Figura 10:
dati di laboratorio

qs calcolato [m2/s]

qs
[m2/s]

Figura 11:
dati di campo

q - qc [m2/s]

20
Appunti di trasporto solido

1.9 FORME DI FONDO

Quale diretta conseguenza del trasporto solido si possono verificare, sul fondo degli alvei, le cosiddette
"forme di fondo". I sedimenti, muovendosi, tendono infatti a conformare il letto dell'alveo secondo
ondulazioni di varia misura, forma e disposizione (ripples, dune, antidune). Possibili significativi esempi
sono le ondulazioni della sabbia su di una spiaggia, le dune del deserto (create dal vento).
Le forme di fondo si muovono. La loro natura risulta funzione dei parametri della corrente, secondo
relazioni più o meno bene individuate. I loro principali effetti sono:
1. le forme di fondo creano macroscabrezze che influiscono sulla resistenza dell'alveo; tali effetti
devono essere tenuti in conto con opportune modifiche delle formule del paragrafo 1.2 che legano la
scabrezza alla sola granulometria dei sedimenti nell'alveo.
2. Nelle formule di trasporto Φ = f(φ, φc) la φ va corretta per tenere conto degli effetti cui al punto
(1): la quota parte di perdite di carico attribuibili alle macroscabrezze deve essere depurata dalla
valutazione degli sforzi attivi per il trasporto solido.
3. Le ondulazioni di fondo, per profondità d'acqua ridotte, inducono ondulazioni sulla superficie libera
della corrente; l'interazione fra i due fenomeni può portare ad un effetto di mutua esaltazione degli
stessi.
Si noti che, in generale, il fenomeno delle forme di fondo è poco significativo per i torrenti; da non
trascurarsi per i fiumi.

Jansen pag. 320

Figura 12: Esempio di forme di fondo - esperimento di laboratorio in canaletta idraulica

21
Appunti di trasporto solido

1.10 TRASPORTO IN SOSPENSIONE

A causa delle componenti di agitazione turbolenta presenti nelle correnti è possibile che le
particelle solide vengano sollevate dal fondo e trasportate, anche per tratti di notevole estensione, in
sospensione all'interno della corrente. Un esempio caratteristico è quello della polvere trasportata dal
vento. Il fenomeno, possibile solo per sedimenti di piccola dimensione, si differenzia da quelli fino ad
ora discussi per il fatto che le particelle, una volta in sospensione, assumono la velocità di traslazione
della corrente, mentre le particelle in moto sul fondo si muovono con una velocità (media) inferiore a
quella del fluido sovrastante .
Numerosi studi si sono occupati del trasporto in sospensione; la quantità di solido trasportato è
legato alla capacità della corrente di sollevare le particelle e al tempo di permanenza delle particelle in
sospensione, prima di cadere nuovamente sul fondo (e quindi alla concentrazione di solido
all'interno della corrente), mentre, come detto, è nota la velocità di traslazione del solido.
Il fenomeno risulta di notevole importanza per alcuni corsi d'acqua (Nilo, Gange, ...) ma è di
poco interesse in relazione alla situazione italiana.

2 FENOMENI LOCALIZZATI
Oltre ai fenomeni erosivi distribuiti, caratterizzanti l'intera sezione d'alveo e determinati
dall'azione di trascinamento della corrente nel suo insieme, sono di importanza nei problemi
idraulici i fenomeni localizzati a zone ristrette dell'alveo, dovuti ad accelerazioni locali create a
loro volta da singolarità geometriche (curve, opere idrauliche, ...).
L'argomento non viene approfondito in questa sede; per una definizione quantitativa degli
effetti erosivi è necessaria una trattazione specifica per ogni situazione di campo di moto. Si
descrivono brevemente di seguito tre casi tipici, rimandando ai testi specializzati per i dettagli
quantitativi.

Erosione alla base di pile e spalle.

La presenza delle pile e spalle dei ponti all'interno della corrente è causa dell'insorgere di moti
vorticosi, con conseguenti aumenti di velocità in prossimità del fondo. Quale effetto di tali
accelerazioni locali si instaurano processi erosivi alla base delle strutture, che arrivano
facilmente a profondità pari a due volte le dimensioni trasversali delle strutture medesime;
vengono così scoperte le fondazioni, mettendo a rischio la stabilità dell'opera.

Erosione alla base di una briglia.

La lama d'acqua in caduta da una briglia crea un campo di moto fortemente turbolento alla
base dell'opera, caratterizzato da notevole capacità erosiva. Anche in questo caso la buca di
erosione che si può creare qualora il fondo non sia opportunamente protetto può
compromettere le fondazioni e quindi la stabilità della struttura.

Erosione ai piedi di uno scivolo.

La corrente, veloce lungo lo scivolo, ritorna lenta attraverso un risalto. Tale fenomeno,
fortemente dissipativo, è accompagnato da intensi scambi energetici (ovverosia moti vorticosi)
che esercitano sul fondo azioni erosive anche di notevole entità.

22
Appunti di trasporto solido

Erosione alla base

delle pile dei ponti

Erosione alla base

delle spalle dei ponti

23
Appunti di trasporto solido

Erosione al piede
di una briglia

Erosione a valle
di uno scivolo

24