Nella Fig. 6.2 sono riportati i risultati delle prove triassiali ... - Padis

144 9 Implementazione numerica di equazioni costitutive in codici agli elementi finiti
L’abbassamento u2 è stato imposto in una serie di passi temporali successivi,
a ciascuno dei quali è assegnata la stessa deformazione verticale, che nel caso di n
passi è pari a ɛ n 2 = 5
n %.
In una prima fase, al fine di valutare l’influenza del numero dei passi sulla
soluzione, sono state svolte diverse simulazioni, riportate in Tab. 9.3, nelle quali
lo stesso abbassamento stato imposto in un numero variabile di passi.
Tabella 9.3. Numero di steps di carico nelle diverse simulazioni.
TEST A B C D
N. PASSI 10 20 50 1000
Nella Fig. 9.2 sono riportate le curve ɛv : ɛs e q : ɛs relative alle quattro
simulazioni.
Figura 9.2. Andamento delle curve: (a) ɛv : ɛs ; (b) q : ɛs nelle quattro prove.
Le curve riproducono il comportamento tipico di un terreno normalmente
consolidato sottoposto ad un percorso di carico TX–CD.
Si noti che tutte le curve relative alle quattro simulazioni sono praticamente
coincidenti, a testimonianza dell’elevata accuratezza dell’algoritmo anche per
ampiezze elevate del passo di carico.
La efficacia della procedura di linearizzazione numerica è stata valutata calcolando
la velocità di convergenza dell’algoritmo di integrazione ad un fissato passo
temporale. Ricordando la definizione di residuo 8.58 la convergenza della procedura
iterativa per la soluzione del problema discreto di equilibrio risulta quadratica
se:
R (k+1)
= c
(9.34)
R (k)
2