UNIVERSIDAD DE LOS ANDES

Grupo de Investigación en Geotecnia

Facultad de Ingenierı́a
Departamento de Ingenierı́a Civil y Ambiental

Comparación del ángulo de fricción crı́tico con el

ángulo de reposo. Análisis de la influencia de
algunos factores en la determinación del
ángulo de reposo

Diana Paola Solaque Guzmán

Maestrı́a en Ingenierı́a Civil

Bogotá D.C., Enero 2008

Contenido

1. Estado del Conocimiento 4

1.1. Resistencia de Materiales . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
1.1.1. Fricción entre cuerpos sólidos . . . . . . . . . . . . . . . . . . . . . . 4
1.1.2. Fricción de masas granulares - sin cohesión . . . . . . . . . . . . . . 7
1.2. Ángulo de reposo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9

2. Estado Crı́tico 35
2.1. Comportamiento del suelo en ensayos de corte . . . . . . . . . . . . . . . . . 36
2.2. Estado pico, último y residual del suelo . . . . . . . . . . . . . . . . . . . . 38
2.3. Estado crı́tico del suelo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
2.4. Superficies de fluencia . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42
2.5. Resistencia del suelo en el estado crı́tico medido en ensayos triaxiales . . . . 43

3. Propiedades de los suelos granulares 48

3.1. Tamaño y Distribución de Partı́culas . . . . . . . . . . . . . . . . . . . . . . 49
3.2. Forma de la Partı́cula . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 52
3.3. Análisis Fractal . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57

4. Sensibilidad del modelo hipoplástico con ϕc 63

4.1. Modelo constitutivo hipoplástico . . . . . . . . . . . . . . . . . . . . . . . . 63
4.2. Resultados de ensayos triaxiales en arena de Karlsruhe . . . . . . . . . . . . 68
4.3. Resultados de ensayos triaxiales en arena del Guamo . . . . . . . . . . . . . 69
4.4. Simulaciones con ElementTest de los ensayos triaxiales en arena de Karls-
ruhe y del Guamo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71
4.5. Simulación de ensayos triaxiales con parámetros de la arena del Guamo
variando ϕc - Densidad media (eo = 0,9) . . . . . . . . . . . . . . . . . . . . 73

I
 4.6. Simulación de ensayos triaxiales con parámetros de la arena del Guamo
variando ϕc - Estado suelto (eo = 1,1) . . . . . . . . . . . . . . . . . . . . . 78
4.7. Simulación de ensayos triaxiales con parámetros de la arena de Karlsruhe -
Estado suelto (eo = 0,9) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 82
4.8. Simulación de ensayos triaxiales con parámetros de la arena de Karlsruhe
variando ϕc - Estado denso (eo = 0,6) . . . . . . . . . . . . . . . . . . . . . 86
4.9. Resumen de las simulaciones de sensibilidad del modelo con ϕc . . . . . . . 90
4.10. Simulaciones ensayo triaxial variando eo . . . . . . . . . . . . . . . . . . . . 90
4.11. Ensayos Triaxiales en arenas . . . . . . . . . . . . . . . . . . . . . . . . . . . 91

5. Ensayos de Laboratorio 93
5.1. Ángulo de Reposo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 93
5.1.1. Métodos seleccionados para determinar el ángulo de reposo . . . . . 95
[Link]. Método: embudo . . . . . . . . . . . . . . . . . . . . . . . . 95
[Link]. Método: Cilindro “Santamarina” . . . . . . . . . . . . . . . 97
[Link]. Método: Cilindro “Train” . . . . . . . . . . . . . . . . . . . 98
[Link]. Método: Cilindro abierto . . . . . . . . . . . . . . . . . . . 99
5.1.2. Influencia del método de ensayo y la cantidad de material (masa) en
el ángulo de reposo . . . . . . . . . . . . . . . . . . . . . . . . . . . . 100
[Link]. Método: embudo . . . . . . . . . . . . . . . . . . . . . . . . 100
[Link]. Método: cilindro “Santamarina” . . . . . . . . . . . . . . . 101
[Link]. Método: cilindro “Train” . . . . . . . . . . . . . . . . . . . 102
[Link]. Método: cilindro abierto . . . . . . . . . . . . . . . . . . . . 103
[Link]. Cuadro Resumen - influencia del método y la cantidad del
material en el ángulo de reposo . . . . . . . . . . . . . . . . 104
[Link]. Relación lineal entre el ángulo de reposo y la redondez de
la partı́cula (Santamarina & Cho [26]) . . . . . . . . . . . . 105
5.1.3. Influencia de la velocidad del desplazamiento del embudo en el ángulo
de reposo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 106
5.1.4. Influencia del diámetro de salida del embudo en el ángulo de reposo 107
5.1.5. Influencia del tamaño y distribución de las partı́culas en el ángulo
de reposo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 108
5.1.6. Variación del ángulo de reposo con la forma de las partı́culas . . . . 110
5.1.7. Influencia de la saturación del material en el ángulo de reposo . . . . 112
5.2. Ángulo de Fricción Crı́tico . . . . . . . . . . . . . . . . . . . . . . . . . . . . 112
5.3. Celda Hele Shaw . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 116
5.4. Equipo para medir el ángulo de reposo . . . . . . . . . . . . . . . . . . . . . 118

II
6. Conclusiones 120

III
Indice de figuras

1.1. Fricción y Criterio de deslizamiento . . . . . . . . . . . . . . . . . . . . . . 5

1.2. Métodos descritos por Train [28] para determinar el ángulo de reposo . . . . 10
1.3. Condición para que las partı́culas soportadas sobre 3 o 4 esferas rueden en
la superficie del talud [28] . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
1.4. Métodos para determinar el ángulo de reposo - Kalman et al. [16] . . . . . . 13
1.5. Influencia de algunos factores en el ángulo de reposo [16] . . . . . . . . . . . 14
1.6. Ángulo de reposo dinámico [16] . . . . . . . . . . . . . . . . . . . . . . . . . 15
1.7. Métodos para deteminar el ángulo de fricción (Atkinson [3]) . . . . . . . . . 16
1.8. Montaje para determinar el ángulo Carr de Reposo. (A) contador de tiempo;
(B) mecanismo de vibración; (C) medidor de amplitud; (D) Reostato; (E)
mecanismo golpeador - ASTM [8] . . . . . . . . . . . . . . . . . . . . . . . . 17
1.9. Forma de la Pila - ASTM [8] . . . . . . . . . . . . . . . . . . . . . . . . . . 17
1.10. Determinación de ϕc del ángulo de reposo [11] . . . . . . . . . . . . . . . . . 18
1.11. Variación de L, Zh y W en la celda Hele Shaw [21] . . . . . . . . . . . . . . 20
1.12. Superficies de pilas construidas con esferas de vidrio [21] . . . . . . . . . . . 21
1.13. Aumento de la pendiente entre la parte lineal y la cola de la pila [?] . . . . 21
1.14. Influencia de Zh y W en el ángulo de reposo (θ = ϕrep ) [21] . . . . . . . . . 22
1.15. Método simplificado para determinar el ángulo de reposo [25] . . . . . . . . 24
1.16. Ángulo de reposo [7] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
1.17. Resultados del Ángulo de reposo - Dodds, 2003 [7] . . . . . . . . . . . . . . 27
1.18. Análisis de fuerzas - Talud infinitamente largo [30] . . . . . . . . . . . . . . 28
1.19. Análisis por cuña [30] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
1.20. Rugosidad de la superficie base [30] . . . . . . . . . . . . . . . . . . . . . . . 30
1.21. Influencia de la rugosidad de la base en el ángulo de reposo [30] . . . . . . . 31
1.22. Análisis de la rugosidad de la base [30] . . . . . . . . . . . . . . . . . . . . . 32
1.23. Variación de los ángulos β, δ y ϕ con el porcentaje de finos en la masa de
suelo [30] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34

IV
2.1. Curva tı́pica esfuerzo-deformación de metales. . . . . . . . . . . . . . . . . . 36
2.2. Comportamiento tı́pico del suelo en ensayo de corte drenado (Atkinson, 1993
[3]) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
2.3. Compresión y dilatancia durante el corte (Atkinson, 1993 [3]). . . . . . . . . 38
2.4. Estado pico, último y residual (Atkinson, 1993 [3]) . . . . . . . . . . . . . . 39
2.5. Estado crı́tico del suelo (Atkinson, 1993 [3]). . . . . . . . . . . . . . . . . . . 39
0 0
2.6. Lı́nea del estado crı́tico en el espacio q − p − ν (Atkinson & Bransby [4]). . 41
2.7. Superficies de fluencia (Britto & Gunn [19]) . . . . . . . . . . . . . . . . . . 43
2.8. Ensayo triaxial no drenados en muestras normalmente consolidadas . . . . . 44
2.9. Ensayo triaxial drenados en muestras normalmente consolidadas (Atkinson
& Bransby [4]) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
2.10. Lı́nea del estado crı́tico obtenida de ensayos triaxiales (Atkinson [3]) . . . . 45
2.11. Cı́rculo de Mohr . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46

3.1. Diámetro de una partı́cula [24] . . . . . . . . . . . . . . . . . . . . . . . . . 49

3.2. Tamaño de diferentes granos en un suelo [18]. . . . . . . . . . . . . . . . . . 50
3.3. Curva granulométrica - arena del Guamo . . . . . . . . . . . . . . . . . . . 51
3.4. Caracterización de la forma de la partı́cula dependiendo de la escala [18] . . 52
3.5. Forma de las partı́culas [18]. . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
3.6. Redondez de una partı́cula [24] . . . . . . . . . . . . . . . . . . . . . . . . . 53
3.7. Forma de las partı́culas - arena del Guamo . . . . . . . . . . . . . . . . . . 55
3.8. Angulo de reposo en función de la redondez de las partı́culas [26] . . . . . . 56
3.9. Dimensión fractal . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
3.10. Análisis fractal, (a) perfil elipsoidal (b) perfil rugoso [29] . . . . . . . . . . . 58
3.11. Relación entre Nn y rn Vallejo, 1995 [29] . . . . . . . . . . . . . . . . . . . . 59
3.12. Dimensión fractal partı́cula N o 1 arena del Guamo . . . . . . . . . . . . . . 60
3.13. Dimensión fractal partı́cula N o 2 arena del Guamo . . . . . . . . . . . . . . 60
3.14. Dimensión fractal partı́cula N o 3 arena del Guamo . . . . . . . . . . . . . . 60
3.15. Relación entre Nn y rn arena del Guamo . . . . . . . . . . . . . . . . . . . . 61

4.1. Representación geométrica de las invariantes tan ψ y cos 3ϑ en el espacio de

los esfuerzos principales (Tomado de [31]) . . . . . . . . . . . . . . . . . . . 66
4.2. Rango de las posibles relaciones de vacı́os como función del nivel de esfuerzos
[31] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67
4.3. Ensayos triaxiales CD con arena de Karlsruhe, presión de cámara variable
y eo = 0,53 (Herle & Gudehus [11]). . . . . . . . . . . . . . . . . . . . . . . 69
4.4. Ensayos triaxiales CD con arena del Guamo, presión de cámara variable y
eo = 0,56 (Arias [2]). . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 70

V
4.5. Ensayo triaxial CD con arena de Karlsruhe, presión de cámara variable y
eo = 0,53. (a) Simulación Herle & Gudehus. (b) Simulación con ElementTest. 71
4.6. Ensayo triaxial CD con arena de Karlsruhe, presión de cámara variable y
eo = 0,53. (a) Simulación Herle & Gudehus. (b) Simulación con ElementTest. 71
4.7. Ensayo triaxial CD con arena del Guamo, presión de cámara variable y
eo = 0,56. (a) Simulación realizada por Arias. (b) Simulación con ElementTest. 72
4.8. Ensayo triaxial CD con arena del Guamo, presión de cámara variable y
eo = 0,56. (a) Simulación realizada por Arias. (b) Simulación con ElementTest. 72
4.9. Simulación ensayo triaxial CU con arena del Guamo variando ϕc - densidad
media. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 74
4.10. Simulación ensayo triaxial CD con arena del Guamo variando ϕc - densidad
media. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 76
4.11. Simulación de ensayos triaxiales CU y CD con arena del Guamo, densidad
media . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77
4.12. Simulación ensayo triaxial CU con arena del Guamo variando ϕc - estado
suelto . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 78
4.13. Simulación ensayo triaxial CD con arena del Guamo variando ϕc - estado
suelto . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 80
4.14. Simulación de ensayos triaxiales CU y CD con arena del Guamo, estado suelto 81
4.15. Simulación ensayo triaxial CU con arena de Karlsruhe variando ϕc - estado
suelto . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 82
4.16. Simulación ensayo triaxial CD con arena de Karlsruhe variando ϕc - estado
suelto . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 84
4.17. Simulación de ensayos triaxiales CU y CD con arena de Karlsruhe, estado
suelto . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 85
4.18. Simulación ensayo triaxial CU con arena de Karlsruhe variando ϕc - estado
denso . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 86
4.19. Simulación ensayo triaxial CD con arena de Karlsruhe variando ϕc - estado
denso . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 88
4.20. Simulación ensayos triaxiales CU y CD con arena de Karlsruhe - estado denso 89
4.21. Simulación ensayo triaxial CU y CD con arena del Karlsruhe variando eo . 91
4.22. Ensayos triaxiales no drenados en arenas [14] . . . . . . . . . . . . . . . . . 92
4.23. Ensayos triaxiales no drenados en arenas [14] . . . . . . . . . . . . . . . . . 92

5.1. Formas de las pilas dependiendo de la rugosidad de base . . . . . . . . . . . 94

5.2. Método Embudo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 95
5.3. Flujo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 97
5.4. Método: cilindro - Santamarina . . . . . . . . . . . . . . . . . . . . . . . . . 98

VI
5.5. Método: cilindro “Train” . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 99
5.6. Método: cilindro abierto . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 100
5.7. Influencia de la masa en el ángulo de reposo, método: embudo - Material
pasa tamiz #16, retiene tamiz #30 . . . . . . . . . . . . . . . . . . . . . . . 101
5.8. Influencia de la masa en el ángulo de reposo, método: cilindro “Santamarina”
(material seco) - Material pasa tamiz #16, retiene tamiz #30 . . . . . . . . 102
5.9. Influencia de la masa en el ángulo de reposo, método: Train - Material pasa
tamiz #16, retiene tamiz #30 . . . . . . . . . . . . . . . . . . . . . . . . . . 103
5.10. Influencia del método y la masa en el ángulo de reposo . . . . . . . . . . . . 104
5.11. Influencia de diámetro de salida del embudo en ϕrep . . . . . . . . . . . . . 108
5.12. Variación de ϕrep con el tamaño y distribución de las partı́culas . . . . . . . 109
5.13. Imágenes tomadas con estereoscopio con 10 aumentos. Escala: 1/10 mm . . 111
5.14. Material seco y sumergido . . . . . . . . . . . . . . . . . . . . . . . . . . . . 112
5.15. Ensayo triaxial arena del Guamo . . . . . . . . . . . . . . . . . . . . . . . . 113
5.16. Ensayo triaxial CU arena del Guamo, plano q − 1 . . . . . . . . . . . . . . 114
5.17. Ensayo triaxial CU arena del Guamo, plano ∆u − 1 . . . . . . . . . . . . . 115
0
5.18. Ensayo triaxial CU arena del Guamo, plano q − p . . . . . . . . . . . . . . 115
5.19. Celda Hele Shaw - medidas en cm. . . . . . . . . . . . . . . . . . . . . . . . 117
5.20. Ángulo de reposo con material saturado - medidas en cm. . . . . . . . . . . 117
5.21. Ángulo de reposo - adaptado Norma Din 18126. 1. base, 2. guı́a para man-
tener centrado el embudo, 3. embudo, 4. sistema de poleas, 5. manivela. . . 118
5.22. Ángulo de reposo - adaptado Norma Din 18126. . . . . . . . . . . . . . . . . 119

VII
Indice de tablas

1.1. ángulo de fricción del material variando el porcentaje de finos [30] . . . . . 33

3.1. Forma de partı́cula, tamaño y superficie especı́fica para diferentes materiales

[24]. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
3.2. Métodos para medir tamaños de partı́culas [24]. . . . . . . . . . . . . . . . . 50
3.3. Clasificación de la forma de una partı́cula [24]. . . . . . . . . . . . . . . . . 52
3.4. Métodos para caracterizar la forma de la partı́cula [18] . . . . . . . . . . . . 54
3.5. Clasificación de partı́culas - Arena del Guamo . . . . . . . . . . . . . . . . . 55
3.6. Dimensión fractal calculada para las partı́culas de la Figura 3.10 [29] . . . . 59
3.7. Dimensión fractal - arena del Guamo . . . . . . . . . . . . . . . . . . . . . . 61

4.1. Sensibilidad del modelo hipoplástico con la variación de ϕc . Ensayo triaxial

CU con arena del Guamo, densidad media (eo = 0,9) . . . . . . . . . . . . . 75
4.2. Sensibilidad del modelo hipoplástico con la variación de ϕc . Ensayo triaxial
CD con arena del Guamo, densidad media (eo = 0,9) . . . . . . . . . . . . . 77
4.3. Sensibilidad del modelo hipoplástico con la variación de ϕc . Ensayo triaxial
CU con arena del Guamo, estado suelto (eo = 1,1) . . . . . . . . . . . . . . 79
4.4. Sensibilidad del modelo hipoplástico con la variación de ϕc . Ensayo triaxial
CD con arena del Guamo, estado suelto (eo = 1,1) . . . . . . . . . . . . . . 81
4.5. Sensibilidad del modelo hipoplástico con la variación de ϕc . Ensayo triaxial
CU con arena de Karlsruhe, estado suelto (eo = 0,9) . . . . . . . . . . . . . 83
4.6. Sensibilidad del modelo hipoplástico con la variación de ϕc . Ensayo triaxial
CD con arena de Karlsruhe, estado suelto (eo = 0,9) . . . . . . . . . . . . . 85
4.7. Sensibilidad del modelo hipoplástico con la variación de ϕc . Ensayo triaxial
CU con arena de Karlsruhe, estado denso (eo = 0,6) . . . . . . . . . . . . . 87
4.8. Sensibilidad del modelo hipoplástico con la variación de ϕc . Ensayo triaxial
CD con arena de Karlsruhe, estado suelto (eo = 0,9) . . . . . . . . . . . . . 89

VIII
 MIC 2008-I-34

4.9. Resumen de sensibilidad del modelo hipoplástico con la variación de ϕc . . 90

5.1. Flujo de las partı́culas dependiendo del embudo y altura del material dentro
del embudo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 96
5.2. Influencia de la masa en el ángulo de reposo, método: embudo - Material
pasa tamiz #16, retiene tamiz #30 . . . . . . . . . . . . . . . . . . . . . . . 101
5.3. Influencia de la masa en el ángulo de reposo, método: cilindro “Santamarina”
(material seco) - Material pasa tamiz #16, retiene tamiz #30 . . . . . . . . 102
5.4. Influencia de la masa en el ángulo de reposo, método: Train - Material pasa
tamiz #16, retiene tamiz #30 . . . . . . . . . . . . . . . . . . . . . . . . . . 103
5.5. Influencia del método y la masa en el ángulo de reposo . . . . . . . . . . . . 104
5.6. Ángulo de reposo de la arena del Guamo determinado a partir de la relación
lineal entre ϕrep y la redondez de la partı́cula . . . . . . . . . . . . . . . . . 106
5.7. Variación de la velocidad con que se levanta el embudo - Material pasa T16
retiene T30 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 107
5.8. Influencia del diámetro de salida del embudo en ϕrep - Material pasa tamiz
#16, retiene tamiz #30 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 107
5.9. Descripción de los materiales 1 y 2 . . . . . . . . . . . . . . . . . . . . . . . 108
5.10. Variación de ϕrep con el tamaño y distribución de las partı́culas . . . . . . . 109
5.11. Influencia de la forma de las partı́culas . . . . . . . . . . . . . . . . . . . . 111
5.12. Variación de ϕrep con la saturación del material . . . . . . . . . . . . . . . . 112
5.13. Espesor mı́nimo de la celda Hele Shaw (20φpart ) teniendo en cuenta el diáme-
tro máximo de la partı́cula . . . . . . . . . . . . . . . . . . . . . . . . . . . 116

1
Introducción

El modelo hipoplástico es una ley constitutiva incremental no lineal desarrollada para ma-
teriales granulares, considerando que los materiales granulares son inelásticos (deformacio-
nes permanentes) y no lineales (rigidez variable). Este modelo describe el comportamiento
esfuerzo-deformación del suelo mediante una formulación matemática - ecuación tensorial.

El modelo hipoplástico es formulado en términos de variables de estado y parámetros

del suelo. Las variables de estado son caracterı́sticas del material que cambian durante el
proceso de deformación o carga. Los parámetros del suelo son propiedades del material que
no cambian durante este proceso.

Los parámetros del modelo hipoplástico son 8: ángulo de fricción crı́tico (ϕc ), dureza del
esqueleto granular (hs), sensibilidad del esqueleto granular a la presión (n), exponente de
picnotropı́a (α), exponente de barotropı́a (β) y las relaciones de vacı́os máxima posible
(ei0 ), mı́nima posible (ed0 ) y crı́tica (ec0 ) para presiones nulas de confinamiento.

El ángulo de fricción crı́tico (ϕc ) es un parámetro del suelo porque es un valor que una vez
alcanzado no varı́a durante el proceso de deformación o de carga. Ésto no sucederı́a si se
trabajara con el ángulo de fricción del material (ϕ), porque éste se considera una variable
de estado (depende del esfuerzo aplicado, densidad, entre otros). ϕc es una relación de los
esfuerzos principales en el estado crı́tico y es importante para el modelo hipoplástico ya
que define el estado crı́tico o último propio de cada material. El estado crı́tico o último
de un material está definido como el estado en el cual el material continúa deformándose
indefinidamente sin presentarse cambios en el esfuerzo efectivo aplicado y en el volumen.

Para determinar que tanto influye ϕc en el comportamiento de un material se realizaron

varias simulaciones de ensayos triaxiales con Elementtest (ley constitutiva: Niemunis y

2
 MIC 2008-I-34

código del programa: Cudmani)y con los parámetros de la arena del Guamo y de Karlsruhe,
variando únicamente ϕc . Con base en los resultados preliminares se podrı́a decir que el
ángulo de fricción crı́tico interviene significativamente en el comportamiento de un material.

Tradicionalmente se considera el ángulo de fricción crı́tico igual que el ángulo de reposo

de un material (ϕc = ϕreposo ). Herle y Gudehus, 1999 [11] validan experimentalmente esta
consideración. ϕc se obtiene mediante ensayos triaxiales monotónicos o de ensayos de corte
simple. ϕreposo se obtiene a partir de diferentes métodos (Train, 1958 [28]). Uno de los
métodos aplicables a materiales granulares secos es la construcción de una pila cónica por
medio de un embudo. ϕc es remplazado por el ϕreposo dado lo dispendioso y complejo
de los ensayos necesarios para determinar ϕc comparados con los métodos para obtener
ϕreposo y la similitud de los resultados obtenidos. Esta simplificación se hace debido a que
el objetivo de estos modelos es obtener los parámetros del suelo de la forma más sencilla
posible. Pero algunos autores como Kalman et al. 1993 [16] y Chik y Vallejo 2005 [30] en
sus artı́culos presentan que factores como rugosidad de la base, tamaño de la pila, tamaño
de partı́culas y su distribución, densidad de la pila, efectos dinámicos, entre otros, influyen
significativamente en el valor de ϕrep .

Dos preguntas resultan entonces:

1. ¿Realmente el estado crı́tico de un suelo es único independientemente del estado

inicial?.

2. ¿A pesar de todos los factores que influyen en la determinación del ϕreposo , es correcto
seguir empleando esta simplificación (ϕc = ϕreposo )?. ¿Cómo se debe realizar el ensayo
para utilizar esta simplificación?

Resolver estas dos preguntas son el objetivo de este trabajo.

3
Capı́tulo 1

Estado del Conocimiento

Todos los materiales poseen resistencia debido a sus propiedades internas. Por esta razón
las pendientes de taludes naturales y artificiales se mantienen en equilibrio debido a la
resistencia al corte que posee el suelo o la roca de que estén formados. Parte de la resistencia
de los materiales se debe a la fricción del material. A continuación se presenta un resumen
de la información recopilada desde 1958 sobre el ángulo de fricción, ángulo de reposo y
ángulo de fricción crı́tico de un material.

1.1. Resistencia de Materiales

1.1.1. Fricción entre cuerpos sólidos

La fricción (ϕ) es una propiedad del material de la cual depende, entre otras, la fuerza de
fricción (Pf ) que se desarrolla entre dos cuerpos. A continuación se presenta un ejemplo
basado en la estática de cuerpos rı́gidos para aclarar el concepto de ϕ y Pf y determinar
la interacción entre ellos [27].

En la Fig. 1.1 (a) un cuerpo descansa sobre una superficie horizontal. Sobre el cuerpo actúa
una fuerza vertical Pn que incluye el peso del cuerpo; la reacción a Pn es la fuerza Pr . Estas
fuerzas verticales proporcionan una fuerza de fricción Pf , la cual puede expresarse como:

Pf = Pn tanϕ (1.1)

donde, ϕ es el ángulo de fricción (propiedad de los materiales en contacto). Para la mayorı́a

de los materiales esta propiedad es aproximadamente constante e independiente de las

4
CAPÍTULO 1. ESTADO DEL CONOCIMIENTO MIC 2008-I-34

Figura 1.1. Fricción y Criterio de deslizamiento

5
CAPÍTULO 1. ESTADO DEL CONOCIMIENTO MIC 2008-I-34

fuerzas que actúan. La fuerza de fricción (Pf ) no actúa a menos que sea requerida para
resistir una fuerza aplicada.

En la Fig 1.1 (b) una pequeña fuerza horizontal Ps0 es aplicada al cuerpo. La resultante de
Pn y Ps0 es la fuerza P 0 , la cual forma un ángulo α con respecto a la normal al plano sobre
el cual descansa el cuerpo. El ángulo α se denomina ángulo de oblicuidad de la fuerza P 0
y depende de las fuerzas que actúan y no de las propiedades de los materiales en contacto.
Al aplicar la fuerza horizontal Ps0 surge una fuerza resistente Pf0 . Esta fuerza resistente
depende de las propiedades de los materiales y de las fuerzas aplicadas. En este caso, la
fuerza Pf0 es menor que Pf y el ángulo α es menor que ϕ, por lo tanto no se presentará des-
lizamiento o desplazamiento entre el cuerpo y el plano, es decir, la fuerza aplicada es menor
a la fuerza resistente.

En la Fig 1.1 (c) la fuerza horizontal aplicada Ps0 es igual a la fuerza de fricción Pf , entonces,
toda la fricción de la fuerza Pf actúa para resistir a Ps0 . Debido a que la fuerza actuante es
igual a la fuerza resistente (α = ϕ), se inicia el deslizamiento hacia la derecha. En la Fig
1.1 (d) se puede observar esta misma situación pero el desplazamiento ocurre en dirección
contraria ya que en este caso la dirección de fuerza aplicada es hacia la izquierda. La fuerza
de fricción (fuerza resistente) siempre actúa en dirección contraria al movimiento.

Analizando los casos anteriormente mencionados se pueden establecer dos principios fun-
damentales:

La fuerza de fricción depende directamente de la fuerza aplicada y del ángulo de

fricción. Si alguno de los dos es igual a cero no habrá fuerza de fricción.

Si el ángulo de oblicuidad de la fuerza aplicada es menor que el ángulo de fricción, sólo

actuará una parte de la fricción y por lo tanto no hay posibilidad de deslizamiento.
El ángulo de fricción es el limitante de la oblicuidad. Cuando la oblicuidad es igual
a ϕ empezará a presentarse deslizamiento.

Estos principios fundamentales pueden extenderse a otros casos de fuerzas aplicadas como:

En la Fig 1.1 (e) un cuerpo sobre el cual actúa sólo la gravedad, descansa sobre un plano
inclinado, la pendiente de este plano se aumentó gradualmente hasta llegar la pendiente
crı́tica icr en la cual el deslizamiento ocurre. Esta condición se presenta cuando icr = ϕ ,
esta relación se mantendrá solo si el peso y la reacción son las únicas fuerzas que actúan.

Sobre un cuerpo con peso W actúa una fuerza P , esta fuerza forma un ángulo β con la
normal al plano sobre el cual reposa el cuerpo (Figura 1.1f). Este plano tiene una inclinación

6
CAPÍTULO 1. ESTADO DEL CONOCIMIENTO MIC 2008-I-34

i. La fuerza resultante de P y W forma un ángulo α con la normal al plano inclinado, en

este caso se obtiene:

W sen(i) − P sen(β)
tan α = (1.2)
W cos(i) + P cos(β)

Para esta condición el deslizamiento ocurrirá si el ángulo α es igual a ϕ. Si la componente

del peso en el eje paralelo al plano es mayor que la componente de P sobre este mismo plano,
el deslizamiento es hacia abajo. Si la componente de P es mayor que W el deslizamiento
es hacia arriba.

Para los cuerpos sólidos que están en contacto y no presentan adhesión entre ellos, la fricción
es sinónimo de resistencia al corte. En la mayorı́a de los suelos la fricción representa una
parte importante de la resistencia al corte.

1.1.2. Fricción de masas granulares - sin cohesión

En arenas y otros materiales sin cohesión, la resistencia al deslizamiento sobre cualquier

plano a través del material es similar a la fricción entre cuerpos sólidos (la resistencia de-
pende de la presión directa que existe sobre el plano y del ángulo de fricción). Sin embargo,
la resistencia friccionante de las arenas es más compleja; ésta resistencia está compuesta
por:

resistencia friccional entre los granos, ésta es una combinación de fricción deslizante
y fricción rodante.

entrelace de las parı́culas.

El entrelace de las partı́culas es un fenómeno que contribuye a la resistencia en arenas den-

sas. Éste fenómeno no ocurre en arenas sueltas. La pérdida gradual de resistencia después
de que se ha pasado el punto máximo puede ser atribuido a una disminución gradual del
entrelace, lo cual sucede debido a que la muestra disminuye en densidad.

La densidad de la arena influye apreciablemente en el ángulo de fricción del material (ϕ)

y el entrelace de las partı́culas. ϕ también varı́a con la magnitud del esfuerzo aplicado
durante el ensayo.

Los principios fundamentales sobre los cuales está basada la teorı́a de Mohr no son invalida-
das por la ocurrencia del entrelace; las envolventes de Mohr solamente tendrán ordenadas

7
CAPÍTULO 1. ESTADO DEL CONOCIMIENTO MIC 2008-I-34

más grandes y pendientes más pronunciadas para suelos densos que para suelos sueltos
[27].

De acuerdo con Taylor, 1962 [27] el procedimiento más antiguo y sencillo para obtener el
ángulo de fricción de un suelo granular seco, es el ángulo de reposo de una pequeña pila
de material.

El ángulo de reposo es el ángulo de fricción del material bajo un esfuerzo prácticamente

igual a cero. Una pila cualquiera de material está en equilibrio cuando los granos menos
estables en su superficie estén en equilibrio, por lo tanto, el ángulo de reposo es determinado
por los granos menos estables.

El ángulo de fricción del material se refiere a condiciones internas. Este ángulo es el

promedio de todos los ángulos fricción del suelo analizado, por esta razón tiende a ser un
poco mayor que el ángulo de reposo.

A continuación se citan algunas consideraciones mencionadas por Taylor [27] que se deben
tener en cuenta para no obtener el ángulo de fricción por medio del ángulo de reposo:

La superficie de la pila de arena estará posiblemente en un estado más suelto que

en su interior, por lo tanto, el ángulo de fricción en la zona de la superficie donde el
ángulo de reposo es determinado, tenderá a ser menor que aquel de toda la masa.

El ángulo de reposo no incluye ningún esfuerzo debido al entrelace.

Las arenas casi nunca son completamente limpias o secas, generalmente tienen cierta
cohesión. Esta pequeña cohesión puede ser de magnitud inapreciable para la resis-
tencia al cortante, pero en la construcción de una pila puede hacer que se obtenga
un mayor ángulo de reposo.

Por estas y otras razones el ángulo de reposo debe ser tomado solo como una aproximación
del ángulo de fricción de un material. En suelos sin cohesión por lo general es menor el
ángulo de reposo que el ángulo de fricción.

De acuerdo con Taylor [27] el ángulo de fricción no es un parámetro del suelo ya que como
se lista a continuación hay muchos factores que influyen en el valor de este ángulo:

El ángulo de fricción de un suelo granular determinado mediante ensayos convencio-

nales depende de la relación de vacı́os inicial (picnotropia).

8
CAPÍTULO 1. ESTADO DEL CONOCIMIENTO MIC 2008-I-34

Algunas propiedades de los granos como la forma, influyen en el valor del ángulo
de fricción. Partı́culas angulares tienen un mayor ángulo de fricción que partı́culas
redondeadas.

El contenido mineralógico de los granos influye en el valor del ángulo de fricción.

El tamaño y distribución de las partı́culas de una masa de suelo influyen en el ángu-

lo de reposo, arenas uniformes tienen ángulos de fricción menores que arenas con
variedad de tamaños de granos.

La vieja creencia de que el agua actuando como lubricante, hace que el ángulo de fricción sea
menor cuando el suelo está sumergido, fue rechazada al comparar los resultados de ensayos
realizados para determinar ϕ en arena saturada y completamente seca. Estos ensayos se
realizaron garantizando que la velocidad del ensayo no generara exceso de presión de poros
en las muestras. Sin embargo, en masas de suelo sumergido las presiones intergranulares son
reducidas por la flotación. Ésto hace que la resistencia del suelo sea menor para la condición
sumergida, ya que la resistencia depende directamente de las presiones intergranulares.

A continuación se hace una recopilación de la información existente del ángulo de reposo

de forma cronológica, describiendo el problema tratado en cada artı́culo y las preguntas
abiertas:

1.2. Ángulo de reposo

• Train, 1958 [28] realizó trabajos experimentales para determinar el ángulo de reposo
de polvos. Estos materiales, al igual que los suelos granulares, presentan resistencia al
movimiento entre partı́culas cuando son sometidos a fuerzas externas, resistencia debida a
la fricción entre partı́culas. En este artı́culo enfocado a fármacos, al igual que lo descrito
anteriormente, el ángulo de reposos de masas de polvo sueltas es un método para determinar
la fricción entre partı́culas.

Train presenta cuatro métodos para determinar el ángulo de reposo (ϕrep ) que se describen
a continuación. En su artı́culo Train simboliza el ángulo de reposo con α. De acuerdo con
la nomenclatura que se utiliza a lo largo de este documento α = ϕrep .

Método I: El embudo empleado en la conformación de la pila de polvo se mantiene

a una altura constante H y el material se vierte cuidadosamente dentro del embudo

9
CAPÍTULO 1. ESTADO DEL CONOCIMIENTO MIC 2008-I-34

para que la pila quede perfectamente cónica. La pila se construye hasta que toque el
orificio de salida del embudo (Figura 1.2a). Luego se mide el diámetro de la base de
la pila (2R) y se determina el ángulo de reposo de la siguiente forma:

H
tan ϕrep = tan α =
R
Método II: La pila se construye levantando lentamente el embudo desde el centro de
una base circular de cualquier diámetro (Figura 1.2b). Cuando se alcanza la altura
máxima y se cubre completamente la superficie de la base, se mide la altura de la
pila y se calcula el ϕrep de la misma forma que en el método I.

Método III: Una caja rectangular completamente llena del material se levanta lenta-
mente desde una esquina hasta que el material empieza a deslizarse. El ángulo que
forma la superficie de la caja con la horizontal es el ángulo de reposo (Figura 1.2c).

Método IV: Un cilindro se llena aproximadamente hasta la mitad con el material y

se sella. Luego se gira hasta que el material dentro del cilindro alcance una forma
plana (Figura 1.2d). El máximo ángulo que este plano forma con la horizontal es el
ángulo de reposo.

(a) Método I (b) Método II

(c) Método III (d) Método IV

Figura 1.2. Métodos descritos por Train [28] para determinar el ángulo de reposo

Los resultados obtenidos por cada uno de los métodos fueron comparados y se encontró que
el método usado influye en los resultados. El ángulo de reposo obtenido por los métodos
I y II es inferior a aquel obtenido con los métodos III y IV. También se encontró que el
ángulo de reposo disminuye a medida que el tamaño de la pila aumenta.

10
CAPÍTULO 1. ESTADO DEL CONOCIMIENTO MIC 2008-I-34

Train también observó que en el proceso de formación de las pilas, las esferas ruedan
unas sobre otras en la superficie de la pila, constituyéndose esto en el principal fenómeno
observado. El autor analizó este fenómeno de dos formas: 1. con una esfera ubicada en el
espacio que forman 3 esferas (Figura 1.3a y b), y 2. la esfera ubicada ahora en el espacio
que forman 4 esferas (Figura 1.3c y d). De este análisis Train determinó el ángulo mı́nimo
que debe existir entre el plano que forman las esferas de la base y la horizontal, para que
la esfera ubicada en el espacio pueda rodar.

Para que la esfera n ruede de acuerdo con la configuración que se muestra en la Figura
1.3a es necesario que el plano que forman las esferas o y p tenga una inclinación mı́nima de
19.3o , mientras que para que la esfera q ruede sobre la esfera s (Figura 1.3b) la inclinación
mı́nima es de 35.1o .

(a) (b)

(c) (d)

Figura 1.3. Condición para que las partı́culas soportadas sobre 3 o 4 esferas rueden en la
superficie del talud [28]

Para que una esfera ruede sobre 4 partı́culas distribuidas como se muestra en la Figura
1.3c es necesario una inclinación mı́nima de 35,3o . Y para la configuración de la figura 1.3d,
donde la esfera w debe rodar por encima de la esfera z, la inclinación es de 45,0o .

Train concluye que para que unas partı́culas rueden sobre otras, el menor ángulo de la pen-
diente es 19.3o y el mayor ángulo posible es 45o , sin que se presente movimiento diferente
a que los granos rueden.

11
CAPÍTULO 1. ESTADO DEL CONOCIMIENTO MIC 2008-I-34

• En 1966 Pilpel N. [23] presenta una revisión de trabajos enfocados más a intereses
de la ingenierı́a quı́mica, concernientes con el flujo de polvos finos y sólidos granulares en
producción. En este trabajo Pilpel define el ángulo de reposo como el ángulo de inclinación
de una pila de polvo cuando ésta empieza a deslizarse sobre sı́ misma. El método que
plantéa para medir el ángulo de reposo es utilizando un cilindro abierto en el fondo, de
7.62 cm de alto y 5.08 cm de diámetro aproximadamente. Sobre una base horizontal se
llena el cilindro con el material y se levanta lentamente formando una pila cónica de radio
r e hipotenusa h. El ángulo de reposo ϕrep es determinado por:

r
cos ϕrep = (1.3)
h

Pilpel denomina este ángulo como el ángulo de reposo estático del material. Este mismo
ángulo, pero medido bajo condiciones dinámicas es llamado “ángulo de reposo cinético”.
Éste es obtenido por rotación del polvo en el disco rotatorio. Usualmente el ángulo cinéti-
co de reposo es menor que el ángulo estático de reposo entre 1 y 5o , dependiendo de la
naturaleza del polvo (distribución de tamaños de las partı́culas). Pilpel también presenta
ecuaciones que relacionan el ángulo de reposo con el tamaño de las partı́culas, estas ecua-
ciones fueron desarrolladas para materiales quı́micos como magnesio.

• En 1993 Kalman et al. [16] publicaron su trabajo de investigación sobre la influencia

de diferentes factores en el ángulo de reposo. Los factores analizados fueron:

Rugosidad de la base

Cantidad de material

Tamaño de las partı́culas y distribución de tamaños

Efecto dinámico

Otros factores que también influyen en el ángulo de reposo (rugosidad de las partı́culas,
velocidad y altura de caı́da del material) no fueron analizados en este trabajo.

Los Autores consideran el ángulo de reposo como un parámetro caracterı́stico de polvos y

materiales granulares y presentan diferentes métodos para determinarlos.

1. Método de Inyección (Figura 1.4a). Existen dos formas para determinar el ángulo
de reposo por este método: a) Altura constante y b) diámetro de la base constante.

12
CAPÍTULO 1. ESTADO DEL CONOCIMIENTO MIC 2008-I-34

Estos dos métodos son los mismos métodos I y II planteados por Train [28] explicados
anteriormente.

2. Método de descarga (Figura 1.4b)

3. Método de inclinación (Figura 1.4c). El ángulo de reposo se determina con un cilindro

o una caja como lo describe Train [28] (métodos III y IV)

(a) Método de Inyección (b) Método de descarga

(c) Método de inclinación

Figura 1.4. Métodos para determinar el ángulo de reposo - Kalman et al. [16]

El método de inyección fue empleado para determinar la influencia de los factores mencio-
nados anteriormente. En los ensayos realizados observaron tres mecanismos fundamentales
que gobiernan el comportamiento en el proceso de formación de las pilas: a) las partı́culas
ruedan en la superficie de la pila, b) las partı́culas se deslizan en la base de la pila, y c)
las partı́culas se deslizan y ruedan al interior de la pila para llenar los vacı́os. En la cons-
trucción de las pilas, las partı́culas se depositaron lenta y cuidadosamente para reducir el
efecto del impacto de las partı́culas en la pila.

La influencia de la rugosidad de la base se determinó construyendo pilas sobre dos tipos

de superficie, una superficie con un coeficiente de fricción bajo (vidrio) y la otra con un

13
CAPÍTULO 1. ESTADO DEL CONOCIMIENTO MIC 2008-I-34

coeficiente de fricción alto (vidrio cubierto con pegamento). En estos ensayos se observó que
en las pilas construidas sobre vidrio las partı́culas se deslizan en la base y en las construidas
sobre la base de fricción alta, las partı́culas ruedan en la superficie de la pila. De éste
análisis puede concluirse que el ángulo de reposo varı́a aproximadamente 5o dependiendo
del coeficiente de fricción de la base. En la Figura 1.5a se presentan los resultados de este
análisis.

Para analizar la influencia de la cantidad de material (Figura 1.5b) construyeron pilas

variando la masa de material. De este análisis encontraron que para algunos métodos el
ϕrep disminuye cuando la cantidad de material aumenta, hasta que un valor del ángulo
constante es alcanzado.

(a) Rugosidad de la base (b) Cantidad de material

(c) tamaño y distribución de partı́culas

Figura 1.5. Influencia de algunos factores en el ángulo de reposo [16]

La influencia del tamaño de las partı́culas en el ángulo de reposo es mayor para las pilas
construidas sobre superficies de baja rugosidad (aproximadamente 3o ). Independientemente
de la rugosidad de la base el ϕrep es mayor para partı́culas de menor tamaño como se observa

14
CAPÍTULO 1. ESTADO DEL CONOCIMIENTO MIC 2008-I-34

en la Figura 1.5c.

Aún adicionando con cuidado las partı́culas para evitar efectos dinámicos, se observó oca-
sionalmente colapso (deslizamiento de las partı́culas en la superficie de la pila). Ésto se
presentó en los dos tipos de superficies empleados en la base para construir las pilas (vidrio
y pegamento), pero más en la base de vidrio. Kalman et al. [16] al igual que Pilpel [23]
realizaron ensayos para determinar el ángulo de reposo dinámico. Este ángulo fue obtenido
sometiendo la pila construida a vibraciones durante 60 s. El ángulo de reposo dinámico
(Fig 1.6) es el ángulo que forma la superficie de la pila sometida a vibración con el plano
horizontal después de que ésta adquirió una forma estable.

Figura 1.6. Ángulo de reposo dinámico [16]

Durante la construcción de las pilas con mezclas de tamaños se observó que las partı́culas
de mayor tamaño ruedan más fácilmente en la superficie de la pila. Este fenómeno debe
ser estudiado más a fondo.

• Atkinson, 1993 [3] en su libro establece una investigación experimental simple del
estado crı́tico. Los parámetros en el estado crı́tico son intrı́nsecos de las partı́culas de
suelo y dependen sólo de la naturaleza del suelo (mineralogı́a, forma, textura, etc.). Las
propiedades intrı́nsecas están relacionadas con los parámetros de clasificación, por lo tanto,
los parámetros de resistencia del suelo en el estado crı́tico se pueden obtener de ensayos de
clasificación.

El ángulo de fricción crı́tico (ϕc ) depende de la naturaleza del suelo. Para suelos finos ϕc
aumenta con el decrecimiento de la plasticidad y para suelos granulares ϕc depende en
general de la forma y rugosidad de los granos y de la gradación del material (bien gradado
o pobremente gradado).

15
CAPÍTULO 1. ESTADO DEL CONOCIMIENTO MIC 2008-I-34

Atkinson plantea que el ángulo de reposo (ic ) de un talud fallado es igual al ángulo de
fricción crı́tico (ϕc ) y presenta 3 formas para determinarlo:

1. Cilindro horizontal (Figura 1.7a)

2. Pila de material

3. Pila sumergida (Figura 1.7b)

(a) (b)

Figura 1.7. Métodos para deteminar el ángulo de fricción (Atkinson [3])

El método del cono sumergido consiste en verter la arena en un contenedor lleno de agua,
luego se drena el contenedor y se mide el ángulo de la pila.

• La Sociedad America para ensayos y materiales (ASTM) en 1999 publicó la Norma D

6393-99 donde establece los procedimientos y aparatos para medir algunas propiedades
del suelo. Estas propiedades las llama ı́ndices de Carr y el ángulo Carr de reposo es uno
de estos ı́ndices.

En la Figura 1.8 se observa el montaje para determinar el ángulo Carr de reposo. Este
montaje está diseñado para clasificar el material (tamizarlo si es necesario) antes de hacer-
lo pasar por el embudo para construir la pila, y para vibrar la pila después de construida
para determinar el ángulo Carr de caı́da. El procedimiento de este ensayo prácticamente es
una combinación de los métodos I y II planteados por Train, 1958 [28] descritos anterior-
mente (altura del embudo y diámetro de la base constantes). La pila se construye sobre la
plataforma circular ubicada centrada debajo del embudo; la tasa de vertimiento debe ser
tal que tome de 15 a 30 seg formar la pila cónica.

16
CAPÍTULO 1. ESTADO DEL CONOCIMIENTO MIC 2008-I-34

Figura 1.8. Montaje para determinar el ángulo Carr de Reposo. (A) contador de tiempo;
(B) mecanismo de vibración; (C) medidor de amplitud; (D) Reostato; (E) mecanismo
golpeador - ASTM [8]

El ángulo Carr de reposo se calcula por medio de la ecuación 1.4. Si la pila no queda de
forma completamente cónica, se debe remover y repetir el ensayo.

H
Angulo Carr de Reposo = tan−1 (1.4)
R

donde H es la altura de la pila y R es el radio de la plataforma circular. Se debe indicar la

forma de la pila: cóncava hacia arriba (Fig 1.9A), cóncava hacia abajo (Fig 1.9B) y recta
(Fig 1.9C). Si la pila queda de forma irregular se debe hacer tres veces el ensayo y sacar
un promedio del ángulo.

Figura 1.9. Forma de la Pila - ASTM [8]

En el año 2000 la ASTM publicó la Norma ASTM C1444-00 para determinar el ángulo
de reposo de polvos que fluyen libremente. Este método de ensayo se retiró de la nor-

17
CAPÍTULO 1. ESTADO DEL CONOCIMIENTO MIC 2008-I-34

matividad, sin sustitución en diciembre de 2005 debido a la limitación de su uso por la

industria.

• Herle & Gudehus, 1999 [11] realizan una estimación simple del ángulo de fricción
crı́tico (ϕc ) por medio del ángulo de reposo (ϕrep ) de un material granular seco (Figura
1.10). ϕrep lo determinan por medio de la construcción de una pila y recomiendan hacer
una pequeña excavación en la pata del talud de la pila, con el fin de borrar alguna influencia
generada durante el proceso de construcción.

Figura 1.10. Determinación de ϕc del ángulo de reposo [11]

En el ensayo realizado para determinar ϕrep un estado cuasi-crı́tico es alcanzado debido al

flujo constante de los granos. A pesar de que la forma de la pila no es perfectamente cónica y
la medida del ángulo puede ser levemente influenciada por el método de ensayo, los valores
de ϕc determinados de esta forma se aproximan bastante bien a los resultados de ensayos
de corte. Ésto lo concluyen después de comparar el ángulo de reposo (valor promedio de 6
mediciones) y el ángulo de fricción crı́tico determinado por ensayos de corte.

A continuación se mencionan algunos análisis que Herle & Gudehus presentan en el artı́culo:

Al rodar las partı́culas en el proceso de formación de la pila, los granos grandes

tienden a superar esfuerzos más altos. Esto puede explicar un incremento de ϕc al
aumentar el tamaños de los granos.

Un incremento pequeño del ϕc observado para granos con diámetro d50 < 0,2 mm
puede deberse al agua capilar resultado de la humedad del aire.

A mayor angularidad de los granos, mayor ángulo de reposo.

Cuando decrece el tamaño de los granos, la influencia de varios fenómenos fisico-quimicos

(capilaridad, fuerzas de van der Waals) en la interacción entre los granos se hace más noto-
ria. El ángulo de reposo se vuelve poco realista si el porcentaje de granos con diámetro

18
CAPÍTULO 1. ESTADO DEL CONOCIMIENTO MIC 2008-I-34

< 0,1 mm es grande. En este caso los ensayos de corte son necesarios para determinar ϕc ;
los ensayos triaxiales con volumen constante en muestras sueltas son las más apropiados.
ϕc puede ser obtenido por la extrapolación de ángulos de fricción pico ϕp para varias den-
sidades, o por la extrapolación de la relación observada entre ϕp y el ángulo de dilatancia
ν de ensayos de corte directo con ν = 0.

• Investigaciones experimentales de la forma de las pilas construidas en una celda Hele

Shaw o silos bi-dimensional son presentados por Grasselli & Herrmann, 1999 [21].
Ellos definen 3 regiones diferentes en la superficie de una pila: Una parte redondeada en la
cima debido al impacto de las partı́culas cayendo; una región lineal central donde el ángulo
de reposo es comúnmente definido; y otra parte redondeada en la base (curvatura opuesta
a la de la cima) conocida como la cola de la pila. De estas 3 partes la más estudiada es la
región lineal - el ángulo de reposo - debido a la intención de relacionarlo directamente con las
propiedades microscópicas de las partı́culas (tamaño, fricción y coeficiente de restitución).

El problema tratado por Grasselli & Herrmann [21] es la influencia de la energı́a inicial y
el flujo de las partı́culas en la forma de una pila granular. Ésta influencia la determinaron
por medio de estudios experimentales en la celda Hele Shaw, usando granos de vidrio y
arena (diámetro promedio 250µm). También mencionan otros trabajos numéricos y teóricos
recientes, sobre la dependencia de la forma de una pila:

Karolyi et al. utilizaron un autómata celular para calcular el perfil de la superficie

mientras se va llenando un silo con un flujo débil de partı́culas, como una función del
coeficiente de restitución de las partı́culas.

Boutreux et al. hicieron un estudio teórico de la influencia del flujo de las partı́culas
en el ángulo de reposo, en una pila construida en un silo 2D.

Los experimentos fueron realizados en una celda Hele Shaw vertical de 20 x 30 cm y espesor
fijo de 5 mm. Este espesor (mayor a 20 veces el diámetro de los granos) permite que la
movilidad local del grano sea esencialmente tridimensional. Con la relación del espesor
mayor a 20 veces el diámetro de las partı́culas, las paredes rı́gidas no juegan un papel
importante en el ángulo de reposo. El sistema en general es bidimensional porque no hay
grado de libertad en la dirección perpendicular a la pared.

Los ensayos fueron realizados por los autores variando la longitud de la celda (L), la
energı́a inicial de las partı́culas (Zh ) y el flujo de las partı́culas (W ), Figura 1.11. El flujo
W se controla variando el diámetro del embudo. La energı́a Zh se varió de 0 a 10 cm. Las

19
CAPÍTULO 1. ESTADO DEL CONOCIMIENTO MIC 2008-I-34

superficies experimentales de las pilas son recuperadas por análisis de imagen y mediante
fotografı́as se obtuvo el ángulo de reposo en función de Zh y W .

Figura 1.11. Variación de L, Zh y W en la celda Hele Shaw [21]

Los autores observaron en estos experimentos que el incremento de Zh o W genera una

superficie alargada y un incremento del tamaño de la cola. También observaron que la forma
de la parte redondeada en la cima de la pila es independiente de Zh y W . Lo contrario
ocurre en la región lineal de la pila, el ángulo de reposo depende de Zh , W y L. En la
Figura 1.12 se presentan las imágenes de estos ensayos experimentales.

Otro fenómeno que observaron fue el incremento local de la pendiente cerca de la cola de
la pila (Figura 1.13). Este efecto aparece para grandes L y pequeñas Zh . Los autores dan
dos posibles explicaciones a este incremento de la pendiente: 1. aumento de la densidad de
la pila y 2. deformación plástica de la pila del material debido a su propio peso.

20
CAPÍTULO 1. ESTADO DEL CONOCIMIENTO MIC 2008-I-34

Figura 1.12. Superficies de pilas construidas con esferas de vidrio [21]

Figura 1.13. Aumento de la pendiente entre la parte lineal y la cola de la pila [?]

En los ensayos con esferas de vidrio, los autores obtuvieron que el ángulo de reposo decrece

21
CAPÍTULO 1. ESTADO DEL CONOCIMIENTO MIC 2008-I-34

con el incremento de la energı́a inicial de las partı́culas Zh (Figura 1.14a). En pilas muy
grandes, el ángulo de reposo podrı́a no depender de la energı́a inicial de las partı́culas, ya
que las partı́culas pueden rebotar en la superficie de la pila hasta que liberen su energı́a
inicial y alcancen un estado de energı́a fijo. El ángulo de reposo también decrece al incre-
mentar el flujo de las partı́culas W (Figura 1.14b). Para altos valores de W la forma del
talud cambia continuamente desde la cima hasta la base, en este caso la región lineal no
puede ser claramente definida y el ángulo local de la superficie de la pila decrece a medida
que se viaja desde la cima hasta la base.

(a) (b)

Figura 1.14. Influencia de Zh y W en el ángulo de reposo (θ = ϕrep ) [21]

Los Autores realizaron los mismos ensayos experimentales mencionados anteriormente pero
ahora con arena de tamaño promedio igual a 250µm. El efecto de la energı́a inicial (Zh )
y el flujo (Z) observado en las pilas de vidrio, es el mismo efecto en las pilas de arena.
Sin embargo, el decrecimiento del valor del ángulo de reposo es bajo. Por ejemplo para
variación del flujo entre W = 0,32g/s y W = 2,2g/s sólo disminuyó 1 grado el ángulo de
reposo, y cerca de 2 a 3 grados para una variación del flujo entre Zh = 2cm y Zh = 8cm. El
pequeño decrecimiento observado para arenas se debe al bajo coeficiente de restitución de
los granos en comparación con las esferas de vidrio. Para arenas, el flujo de las partı́culas
en la superficie es débil y el fenómeno de erosión tiende a dar superficies alargadas.

Los autores concluyen que el ángulo de reposo definido en la región lineal de la pila decrece,
mientras que el tamaño de la cola se incrementa cuando aumenta la energı́a y el flujo de las
partı́culas. Si estos parámetros experimentales no son muy grandes, la superficie de la cola

22
CAPÍTULO 1. ESTADO DEL CONOCIMIENTO MIC 2008-I-34

de la pila puede ser descrita por un termino logarı́tmico: γ ln(L/x), donde γ representa el
tamaño de la cola.

• Herrmann et al., 2000 [22] estudiaron experimentalmente la forma de la cima de una

pila granular construida por vertimiento del material dentro de una celda vertical Hele
Shaw y analizaron el efecto de la energı́a de impacto de los granos.

El estudio de materiales granulares ha captado gran atención en los últimos años. Las
pilas de arena son uno de los sistemas granulares más estudiados experimentalmente y
teóricamente. A pesar de su aparente simplicidad, las pila de arena presentan fenómenos
interesantes como avalanchas, textura interna, segregación y estratificación. Diversos mo-
delos teóricos han sido empleados para describir los perfiles estáticos y dinámicos de una
pila de arena.

Las pilas granulares pueden ser generadas por vertimiento lento de granos desde una altura
dada sobre una base rı́gida (pilas 3D) o dentro de una celda Hele-Shaw (pilas 2D). Las
pilas tı́picamente muestran un talud bien definido, correspondiente al ángulo de reposo.
El ángulo de reposo depende de las caracterı́sticas del material como densidad, humedad,
historia de empaquetamiento y condiciones de borde.

Se analizó la forma de la cima de una pila granular, enfocándose en el efecto de la energı́a de

impacto de las partı́culas. Con este fin se elaboraron una serie de experimentos con esferas
de vidrio vertidas a una rata de velocidad muy lenta dentro de la celda Hele-Shaw y se
estudiaron los cambios del perfil cuando los granos eran vertidos desde diferentes alturas.
Los autores proponen un modelo continuo simple para describir las propiedades del estado
de los granos fluyendo en la cima de la pila, proporcionando una expresión explı́cita para
el perfil de la pila. El modelo trata explı́citamente la disipación de la energı́a cinética de los
granos fluyendo durante la colisión inelástica y proporciona un buen ajuste con los datos
experimentales.

Este modelo no fue estudiado en este trabajo.

• Santamarina & Cho, 2001 [25] establecen un procedimiento simple para determinar
los parámetros del estado crı́tico en suelos arenosos. En este trabajo de investigación sólo
se mencionará el procedimiento para determinar el ángulo de fricción crı́tico, objetivo de
esta tesis.

El ensayo triaxial convencional es usado para la determinación de los parámetros en el

estado crı́tico. Sin embargo, éste requiere de tiempo y de otros ensayos. Los autores de-
sarrollan un procedimiento de ensayos simplificados para determinar la lı́nea del estado

23
CAPÍTULO 1. ESTADO DEL CONOCIMIENTO MIC 2008-I-34

crı́tico en suelos arenosos (procedimiento confiable, económico y rápido [25]). Los resul-
tados obtenidos por medio de este procedimiento fueron comparados con parámetros del
estado crı́tico determinados en ensayos triaxiales convencionales. La comparación muestra
buena concordancia en los resultados.

El comportamiento del suelo en el estado crı́tico (SC) ha sido investigado generalmente en

especı́menes dilatantes y ensayos drenados controlando la velocidad de deformación. Ésto
porque la resistencia en el estado crı́tico puede ser alcanzada a un nivel de deformación
relativamente bajo (Been et al. 1991; Lee 1995)[25].

Dentro del procedimiento que Santamarina y Cho, 2001 describen en el artı́culo establecen
el siguiente método para determinar el ángulo de reposo:

“Determine el ángulo de fricción en el estado crı́tico (ϕcs ) vertiendo el suelo en un cilindro

graduado (1000mL) lleno de agua (un contenedor rectangular transparente es preferible).
Inclı́nelo a más de 60o y devuélvalo nuevamente a la posición vertical. Mida el ángulo de
reposo en la mitad de la región de la pendiente (Fig. 1.15; Conforth 1973; Bolton 1986;
Schofield 1999).

Figura 1.15. Método simplificado para determinar el ángulo de reposo [25]

La medida del ángulo de fricción bajo agua es recomendada por el efecto de las fuerzas
de capilaridad entre partı́culas. Sin embargo, arena limosa o suelos con una fracción de
finos, pueden ser licuados durante la medida del ángulo de fricción, reduciendo a un valor
menor que ϕcs . En este caso, se recomienda secar en el horno el suelo y después seguir el
procedimiento simplificado sin agua”.

• Dodds J., 2003 [7] en su tesis de maestrı́a realiza un análisis de diferentes métodos para
determinar el ángulo de reposo y compara los resultados obtenidos de ϕrep con el ángulo
de fricción crı́tico (ϕc ).

24
CAPÍTULO 1. ESTADO DEL CONOCIMIENTO MIC 2008-I-34

El ángulo de fricción en el estado crı́tico es relacionado con el ángulo de fricción entre

partı́culas y la habilidad de las partı́culas a rotar (Cho, 2001), ϕc comparte esta carac-
terı́stica con es el ángulo de reposo (ϕrep ).

El ángulo de fricción crı́tico y el ángulo de reposo se incrementan si disminuye la redondez

e incrementa la rugosidad de las partı́culas. Este incremento de los dos ángulos se presenta
por el impedimento en la rotación de las partı́culas. Las partı́culas que son más angulares
se traban “interlock” y es necesario tener una mayor energı́a de activación para iniciar el
movimiento.
La dilatancia y el arreglo de fábrica no afectan estos dos ángulos. Esto se debe a que ambos
ángulos se presentan bajo una relación de vacı́os constante. Sin embargo, a diferencia del
ángulo de fricción en el estado crı́tico, un suelo puede experimentar diferencias en los
ángulos de reposo, de acuerdo con el método utilizado. La pregunta que surge es: cuál es
la mejor forma de determinar ϕrep ?.

La forma del ángulo de reposo puede variar (convexa, cóncava o plana) de acuerdo con el
método empleado. El ángulo convexo es creado por la caı́da de arena (pluviación) sobre una
superficie plana (Figura 1.9a); mientras el ángulo cóncavo es el ángulo interior creado por
remover un tapón del centro de la base , permitiendo que el material salga (Figura 1.9b).
El ángulo de reposo plano (Figura 1.9c) es creado por rotación de una celda Hele-Shaw
más de 45o y luego se regresa lentamente a la posición vertical Figura 1.16.

Figura 1.16. Ángulo de reposo [7]

Dos posibles fuentes de error en la medición del ángulo de reposo plano son:

1. la tendencia de la arena a segregarse en la celda Hele-Shaw de acuerdo con el tamaño

de las partı́culas

2. la tendencia del ángulo medido a variar hasta que un ángulo estable sea alcanzado
después de 2 o 3 rotaciones (Herrmann, 1998)

25
CAPÍTULO 1. ESTADO DEL CONOCIMIENTO MIC 2008-I-34

El ángulo de reposo plano de 34 arenas fue determinado. Estos resultados fueron compa-
rados con ensayos convexos y cóncavos (realizados por Narsilio 2003) y ángulos de fricción
crı́ticos determinados por medio de ensayos triaxiales (reportados por Cho 2001). El ángulo
de fricción plano de arenas con gran contenido de finos (cerca del 8 %) no fueron determi-
nados porque los granos tienden a adherirse creando un ángulo inconsistente.

El procedimiento recomendado por Santamaria y Cho, 2001 [25] fue utilizado por Dodds
para determinar el ángulo de reposo plano. Una celda Hele-shaw es sustituida por un
cilindro graduado. Previo a los ensayos las arenas son secadas en hornos por 24 horas y
el enfriamiento se hace en empaque sellados. La arena es vertida en el cilindro hasta una
altura justo encima de la primera de dos lı́neas (Figura 1.16). EL cilindro fue rotado pasado
los 60 grados, luego lentamente retornado a la horizontal, esto haciendolo con cuidado para
evitar que golpes o sacudidas puedan generar una avalancha. Un transportador es usado
para determinar el ángulo entre la lı́nea horizontal y la sección central del talud de arena,
el cuál es el ángulo de reposo. El procedimiento es repetido cambiando la dirección de
rotación (derecha - izquierda - derecha) hasta que un ángulo fijo es determinado.

La medida del ángulo de reposo cóncavo y convexo (Narsilio 2003) se obtiene vertiendo
arena por medio de un embudo, para formar una pila en una superficie plana. Un laser se
pasa a través de la superficie de la pila para medir el ángulo convexo. Luego un pequeño
tapón es removido del centro de la pila, la arena fluye hasta que un ángulo interno estable
es formado. Éste ángulo cóncavo es entonces medido con el laser. El ensayo es repetido 5
veces para cada muestra.

Los resultados del ángulo de reposo determinados de 3 métodos diferentes son presentados
en la Figura 1.17. Se presentan únicamente los valores promedio de las mediciones, la
desviación del promedio es de 3 grados. En la figura los ángulos más altos corresponden
a los medidos con la forma cóncava, los intermedios con la forma plana y los bajos con la
forma convexa y con los medidos en triaxial. La diferencia entre los ángulos medidos en la
forma cóncava y plana es cerca de 8o , plana y convexa es cerca de 1.5o y convexa y triaxial
cerca de 5o . Los ensayos de ángulo de reposo fueron realizados todos sobre las mismas
muestras, mientras que los resultados del ensayo triaxial fueron tomados de ensayos con
arena del mismo nombre, pero esta arena puede no tener las mismas propiedades de las
otras muestras.

La diferencia en el ángulo de reposo entre la superficie del talud plana y convexa se debe
probablemente a la variación en el número y dirección de contacto entre partı́culas. Las
partı́culas en una superficie plana tienen un alto número de coordinación y contactos en
más direcciones de partı́culas que en superficie convexa, lo cuál da una mayor estabilidad.

26
CAPÍTULO 1. ESTADO DEL CONOCIMIENTO MIC 2008-I-34

Las partı́culas que descansan en una superficie cóncava forman un tipo de arco. Como el
ángulo se incrementa y las partı́culas comienzan a moverse hacia el centro, esto hace posible
que formen un anillo de compresión y ocurra arqueamiento incrementando el ángulo de
reposo.

Figura 1.17. Resultados del Ángulo de reposo - Dodds, 2003 [7]

• Chik & Vallejo, 2005 [30] establecen que el ángulo de fricción de un material granular
suelto puede ser medido por el ángulo de reposo (al igual que Taylor,1958 [27]). El ángulo de
reposo es una propiedad única del material y se determina por medio de una pila construida
con granos secos que fluyen libremente. Este ángulo corresponde al ángulo del talud de la
pila con respecto a la base horizontal. Muchos fenómenos naturales como las avalanchas,
estratificación y segregación están relacionados con el ángulo de reposo de los materiales
involucrados. Los autores realizaron una serie de ensayos siguiendo la metodologı́a de la
Norma ASTM C1444-00 y aclaran que en esta norma no se menciona la influencia de la
fricción en la base. Es importante anotar que esta Norma fue retirada en diciembre de 2005
debido a su limitado uso para la industria.

La pregunta de investigación en este artı́culo es la influencia de la rugosidad de la base en el

ángulo de reposo. Para responder esta pregunta los autores trabajaron con dos superficies
diferentes en los ensayos: base rugosa (piedra porosa) y base suave (vidrio).

En este artı́culo los Autores hacen un análisis de falla de talud infinito, con el fin de
demostrar que el ángulo de fricción del suelo bajo ciertas condiciones se puede determinar
con el ángulo de inclinación de una pila:

27
CAPÍTULO 1. ESTADO DEL CONOCIMIENTO MIC 2008-I-34

Cuando se construye un cono de suelo granular sobre una superficie plana dejando caer
libremente los granos dentro del embudo, empiezan a generarse y a fluir delgadas capas de
material que se mueven libremente sobre la superficie inclinada del cono. Esta condición
de falla satisface la condición de falla de un talud infinito. Suponiendo que el ángulo de
fricción de la arena es ϕ y la pendiente del cono es β las fuerzas actuantes dentro del
talud infinito pueden ser analizadas usando ecuaciones de equilibrio estático. Mediante el
análisis de estabilidad de un talud infinito para condiciones drenadas, se demuestra que en
la condición de equilibrio limite β = ϕ.

Figura 1.18. Análisis de fuerzas - Talud infinitamente largo [30]

Los Autores también realizan un análisis de estabilidad de un talud por deslizamiento en

cuña, para correlacionar el ángulo de fricción del material, el ángulo de reposo y el ángulo
de fricción de la interfase material-base:
Si el material de la superficie sobre la cual se construye la pila no es friccionante, entonces
la base del pila puede extenderse en un área muy grande. Esta condición de falla indica
que la fricción de la base puede influenciar en el ángulo de reposo. Suponiendo que una
pila de suelo granular es construida formando un cono perfecto, parte del talud (hacia la
base) está actuando como una cuña triangular (abc) Figura 1.19a. En esta condición la
suposición hecha es que la falla del sistema se debe al deslizamiento de la cuña, causado
por la presión lateral activa de tierras ejercida por el resto de la pila en la cara vertical de
la cuña. Para la cuña actuando como muro de contención del relleno de material granular,
la resistencia a que falle (deslizamiento) depende de un balance de fuerzas entre la presión
lateral de tierras y la fricción de la base resistente al deslizamiento. La magnitud de la
resistencia de la fricción la base se deriva de las fuerzas normales de acuerdo con la teorı́a
de fricción de Coulomb.
Cuando el ángulo de reposo se determina por medio del ensayo con el embudo se obtiene
un talud con inclinación β, un ángulo de fricción interna del material ϕ y el ángulo de

28
CAPÍTULO 1. ESTADO DEL CONOCIMIENTO MIC 2008-I-34

fricción δ entre en material granular y la base. Entonces, por medio de una solución teórica
se puede predecir el valor de β conociendo ϕ y δ. De acuerdo con la teorı́a de presión activa
de tierras de Rankine, la presión que ejerce el relleno sobre la cuña es en dirección paralela
a la pendiente del cono. En la figura 1.19b se observan las fuerzas actuando sobre la cuña.

(a) Desplazamiento lateral (b) Análisis de fuerzas

Figura 1.19. Análisis por cuña [30]

Después de hacer el análisis de deslizamiento por cuña a partir de la presión activa de

tierras de Rankine (Pa ) y de la teorı́a de Coulomb, los autores deducen la ecuación 1.5, que
correlaciona el ángulo de fricción de la interfase δ con el ángulo de reposo β y el ángulo de
fricción del material ϕ.

tan β cos2 ϕ
tan δ = r  2 (1.5)
2+2 1 − cos ϕ
cos β − cos2 ϕ

Los ensayos de laboratorio reportados en este artı́culo fueron realizados con materiales
granulares de diferentes tamaños y 2 tipos de superficie en la base. Los materiales granulares
son: arena de cuarzo gruesa (diámetro promedio 1.16 mm) y arena de cuarzo fina (diámetro
promedio 0.10 mm). Y las superficies: piedra porosa y vidrio, Figura 1.20.

29
CAPÍTULO 1. ESTADO DEL CONOCIMIENTO MIC 2008-I-34

(a) superficie lisa (vidrio) (b) Superficie rugosa (piedra porosa)

Figura 1.20. Rugosidad de la superficie base [30]

El ángulo de reposo por este método se puede determinar de dos formas:

1. Medir directamente del talud el ángulo de reposo (5 mediciones).
2. Medir la longitud de la base en 4 secciones diferentes y la altura vertical del cono, para
calcularlo por medio de la siguiente ecuación:

 
−1 2H
ángulo de reposo = tan (1.6)
DA − d

donde H es la altura del cono (altura fija = 3.81 cm); DA es el promedio de cuatro
mediciones del diámetro de la base; y d es el diámetro interno de salida del embudo (0.51
cm).

De acuerdo con los resultados que presentan Chik & Vallejo en la Figura 1.21b, se puede
analizar que la rugosidad de la base afecta más el ángulo de reposo para partı́culas de
mayor tamaño. Cuando el porcentaje de finos es mayor o igual al 70 % de la masa de la
mezcla (predomina fracción fina), la rugosidad de la base no influye.

30
CAPÍTULO 1. ESTADO DEL CONOCIMIENTO MIC 2008-I-34

(a) Montaje (b) Ensayos

Figura 1.21. Influencia de la rugosidad de la base en el ángulo de reposo [30]

A continuación se describe el análisis que los Autores hacen para dar explicación al fenómeno
que ocurre en cada uno de los ensayos:

En la base de vidrio las partı́culas más pequeñas tienen mayor ángulo de reposo. Ésto
se debe a que hay mayor área de contacto Partı́culas-Base. Este contacto hace que
se genere mayor resistencia friccional en la base y por lo tanto un aumento gradual
de ángulo de reposo, (Figura 1.22a).

31
CAPÍTULO 1. ESTADO DEL CONOCIMIENTO MIC 2008-I-34

(a) Superficie lisa (b) Superficie rugosa

Figura 1.22. Análisis de la rugosidad de la base [30]

En la superficie rugosa las partı́culas gruesas tienen mayor ángulo de reposo debido
al Interlocking entre partı́culas-base. Es decir, Si la mezcla no tiene finos hay mayor
trabazón entre las partı́culas y la base rugosa, por lo tanto las partı́culas llenan los
valles y crestas haciendo que el ángulo de reposo sea mayor.
Si la mezcla tiene partı́culas finas, éstas llenan las depresiones de la base rugosa.
Esta acumulación permite que las partı́culas se muevan más fácilmente en la interfase
permitiendo que la fricción de la base decrezca, por lo tanto se disminuye el ángulo
de reposo. Entonces, a mayor diámetro de partı́culas en la base rugosa se obtiene
mayor ángulo de reposo.

Los autores también realizaron ensayos de resistencia (corte directo) para determinar el
ángulo de fricción del material variando el porcentaje de finos en la masa de suelo. En la
Tabla 1.1 se presentan los resultados de éstos ensayos de corte.

32
CAPÍTULO 1. ESTADO DEL CONOCIMIENTO MIC 2008-I-34

Tabla 1.1. ángulo de fricción del material variando el porcentaje de finos [30]

De estos resultados se concluye que el ángulo de fricción del material es afectado por el
tamaño de las partı́culas en la muestra. A menor diámetro de las partı́culas, menor ángulo
de fricción del material.

Para determinar ángulo de fricción de la interfase entre el material y la superficie base (δ),
los Autores pusieron en contacto el material granular con cada una de las superficies base.
La base se fue inclinando gradualmente y el ángulo al cuál se inició el deslizamiento del
material sobre la base, corresponde al ángulo δ. De estos ensayos los autores determinan
que en la superficie de vidrio a mayor porcentaje de finos en la masa de suelo es mayor el
ángulo de la interfase δ (Figura 1.23a). En este caso la resistencia friccional se debe al área
de contacto de las partı́culas y la base. En la superficie rugosa a mayor contenido de finos,
menor ángulo δ (Figura 1.23b). Ésto se debe a que los finos desarrollan una estructura
suelta encima de la interfase (efecto de las paredes - Pillai 1977; Vallejo 2001). En los
dos tipos de superficie cuando predominan los finos en la mezcla de material, la fricción δ
(base-partı́cula) tiende a ser igual (21 grados).

Al comparar los resultados obtenidos del ángulo δ por medio de la ecuación 1.5 y deter-
minado por los ensayos (inclinación de la base) en la base rugosa, los Autores concluyen
que cuando los granos grandes de arena predominan en la mezcla, el modo de la falla no es
bien modelado por la ecuación 1.5. Ésto se debe a que cuando los granos de arena gruesa
predominan en la mezcla, la falla no tiene lugar en la interfase entre los granos de arena y la
superficie rugosa, sino que la falla se desarrolla en una fina capa localizada en la superficie
libre de la pila de arena (contacto arena-arena).

33
CAPÍTULO 1. ESTADO DEL CONOCIMIENTO MIC 2008-I-34

(a) Superficie lisa

(b) Superfice rugosa

Figura 1.23. Variación de los ángulos β, δ y ϕ con el porcentaje de finos en la masa de

suelo [30]

34
Capı́tulo 2

Estado Crı́tico

Las teorı́as del comportamiento del suelo fueron desarrolladas a partir de los ensayos y
teorı́as del comportamiento de los metales. Para entender y predecir las deformaciones del
suelo por medio de la teorı́a del estado crı́tico (teorı́a desarrollada a partir de la teorı́a
de plasticidad de mecánica de suelos), es necesario entender primero el comportamiento
elástico y plástico de los materiales.

La Figura 2.1 (Britton & Gunn, 1987 [19]) muestra la curva esfuerzo-deformación obtenida
de un ensayo a tensión de una barra de metal. Inicialmente, la relación entre el esfuerzo
y la deformación es lineal (tramo OA). En este tramo el comportamiento del material es
elástico: al descargar la barra desde cualquier punto de este tramo, la trayectoria lineal es
la misma pero en dirección contraria y la deformación se recupera totalmente.
Si la barra es cargada con un esfuerzo un poco mayor al del punto A y se descarga, la
trayectoria también es reversible, pero una parte de la relación esfuerzo-deformación es no-
lineal. El punto donde inicia el comportamiento no-lineal es denominado punto de fluencia
del material (punto B).
Cuando la barra es cargada desde el punto O hasta el punto C y luego se descarga, se sigue
la trayectoria CD. En este proceso se desarrollan dos tipos de deformación, una plástica o
permanente (OD) y una elástica o recuperable (DE).
Al recargar la barra de metal desde el punto D se sigue la trayectoria lineal DC; el punto
C es el nuevo punto de fluencia del material. Luego se sigue una trayectoria no-lineal hasta
que se alcanza el esfuerzo máximo (punto F), en el cual la barra falla. El esfuerzo del punto
F es la resistencia última a tensión del material.

35
CAPÍTULO 2. ESTADO CRÍTICO MIC 2008-I-34

Figura 2.1. Curva tı́pica esfuerzo-deformación de metales.

2.1. Comportamiento del suelo en ensayos de corte

La resistencia de un material es el máximo esfuerzo que puede soportar sin que falle. La
falla puede ser súbita (pérdida completa de la resistencia) o puede presentarse mediante
deformación plástica. Para muchos suelos la falla se presenta con grandes deformaciones
plásticas sin perder completamente la resistencia. Una forma de determinar la resistencia
del suelo es por medio de ensayos de corte (Figura 2.2).

Los suelos normalmente consolidados (NC), arcillas levemente sobreconsolidadas y arenas

sueltas (suelos W definidos por Atkinson [3]), se comprimen cuando el esfuerzo cortante o
desviador aplicado es incrementado (Figura 2.2c). Las arcillas altamente sobreconsolidadas
(SC) y arenas densa (suelos D definidos por Atkinson [3]), dilatan después de experimentar
una pequeña compresión (Figura 2.2c). Las curvas tı́picas esfuerzo-deformación para estos
tipos de suelos, presentadas por Atkinson [3], se muestran en la Figura 2.2b. En esta gráfica
se observa que los dos tipos de suelo al final alcanzan el mismo estado, en el cual el esfuerzo
cortante es constante y no hay más deformación volumétrica. Los suelos sobreconsolidados
o arenas densas, alcanzan un esfuerzo pico antes de llegar a este estado.

La relación entre la deformación volumétrica y la deformación por cortante es definida

como el ángulo de dilatancia ψ (Ecuación 2.1). El signo negativo se introduce en esta
ecuación debido a que la dilatancia (deformación volumétrica negativa) es asociada con
ángulos positivos (Figura 2.2a y c).

dv
tan ψ = − (2.1)
dγ

36
CAPÍTULO 2. ESTADO CRÍTICO MIC 2008-I-34

(a) (b)

(c) (d)

Figura 2.2. Comportamiento tı́pico del suelo en ensayo de corte drenado (Atkinson, 1993
[3])

La variación de la relación de vacı́os con el incremento del esfuerzo cortante se observa

en la Figura 2.2(d). Los dos tipos de suelo (NC o arena suelta y SC o arena densa) en el
estado inicial, bajo el mismo esfuerzo efectivo normal, tienen diferente relación de vacı́os
(eoW > eoD ). Al final del ensayo ambas muestras alcanzan el estado último; estado en el
cual tienen la misma relación de vacı́os (ef W = ef D ).

Los cambios de volumen en el suelo se deben principalmente a la reorganización de las

partı́culas. Por lo tanto, suelos suelto o NC se comprimen, mientras que suelos densos o SC
dilatan. Los granos de los suelos sueltos o normalmente consolidados están espaciados entre
ellos y en la fase de corte ellos pueden moverse dentro de los espacios vacı́os, disminuyendo
el volumen (Figura 2.3). Los suelos densos o SC tienen pocos espacios entre los granos, por
lo tanto al experimentar carga y no tener vacı́os para acomodarse, aumentan el volumen
(dilatan).

37
CAPÍTULO 2. ESTADO CRÍTICO MIC 2008-I-34

Figura 2.3. Compresión y dilatancia durante el corte (Atkinson, 1993 [3]).

2.2. Estado pico, último y residual del suelo

Como se observa en la Figura 2.2b, los suelos SC o densos (D) alcanzan un esfuerzo
cortante pico y después llegan al estado último. El estado pico normalmente es alcanzado
a una deformación del 1 %, mientras que el estado último es alcanzado a una deformación
mayor al 10 % (en algunos suelos el estado último se alcanza a una deformación del 50 %
[3]). El estado pico coincide con el punto de mayor ángulo de dilatancia. Suelos NC o
sueltos (W ) son compresibles, no presentan estado pico y el esfuerzo máximo se presenta
en el estado último.

En la Figura 2.4 se muestra el comportamiento de una arena y una arcilla plástica con
grandes desplazamientos. En el estado último el movimiento de los granos es básicamente
rotacional en ambos tipos de suelos (arena y arcilla). A mayores desplazamientos la defor-
mación comienza a localizarse en diferentes zonas y el esfuerzo cortante aplicado en el suelo
arcilloso decrece. El menor esfuerzo cortante se alcanza después de largos desplazamientos
y esto es llamado estado residual. Los granos en las arcillas se mueven con orientación
paralela en la zona de falla como se muestra en la Figura 2.4. En arenas y otros suelos con
granos redondeados el estado residual es igual que en el estado último. En cierto tipo de
arcillas el esfuerzo en el estado residual puede disminuir al 50 % del esfuerzo en el estado
último (τultimo > τresidual ).

38
CAPÍTULO 2. ESTADO CRÍTICO MIC 2008-I-34

Figura 2.4. Estado pico, último y residual (Atkinson, 1993 [3])

2.3. Estado crı́tico del suelo

El estado crı́tico o último es alcanzado a una deformación mayor al 10 %. En la Figura 2.5

0 0
Atkinson [3] presenta la relación entre el esfuerzo cortante (τ ), el esfuerzo normal (σ ) y la
0
relación de vacı́os (e) del suelo en el estado crı́tico. En el plano e − log σ (Figura 2.5c) se
observa que la lı́nea del estado crı́tico (CSL) es paralela a la lı́nea de normal consolidación
(NCL).

Figura 2.5. Estado crı́tico del suelo (Atkinson, 1993 [3]).

39
CAPÍTULO 2. ESTADO CRÍTICO MIC 2008-I-34

De acuerdo con el criterio de falla de Coulomb, la lı́nea del estado crı́tico está dada por la
ecuación 2.2. En el estado crı́tico la cohesión del suelo es cero (CSL pasa por el origen) y
el ángulo de fricción es el ángulo de fricción crı́tico (ϕc ).

0 0 0
τf = σf tan ϕc (2.2)

0
En el estado crı́tico no hay resistencia cuando el esfuerzo normal efectivo es cero (σ = 0):

0 0 0 0
τf = σf tan ϕc si σf = 0, entonces:
0
τf = 0

En suelos cementados la cohesión es diferente de cero, pero para llegar al estado crı́tico la
deformación es tan grande que esta unión se rompe. Por lo tanto, se sigue cumpliendo que
en el estado crı́tico la cohesión es cero.
0
En el plano e−log σ la lı́nea del estado crı́tico está dada por la ecuación 2.3. Las lı́neas CSL
y NCL son paralelas, por lo tanto, la lı́nea del estado crı́tico se puede definir en términos
de la pendiente de la lı́nea de normal consolidación (Cc ). El parámetro eΓ define la posición
de la lı́nea del estado crı́tico de igual forma que eo define la posición de la lı́nea de normal
consolidación.

0
ef = eΓ − Cc log σf (2.3)

En el estado crı́tico, el suelo continúa deformándose sin cambio de esfuerzo cortante, esfuer-
zo normal, ni relación de vacı́os (Figura 2.2). Atkinson [3] establece que la caracterı́stica
esencial del estado crı́tico es: “durante el corte todos los suelos alcanzan un estado crı́tico
que es independiente del estado inicial”. Es decir, un mismo tipo de suelo con diferente
estado inicial (denso o SC, y suelto o NC) alcanza en el estado crı́tico el mismo esfuerzo
cortante o desviador y llegan a la misma relación de vacı́os.

Algunos autores justifican la existencia de un único estado crı́tico, diciendo que si no se

llegara a un estado constante y la deformación del suelo continuara, la muestra de suelo
estarı́a dilatando o comprimiendo indefinidamente, y ésto no es posible. Otra justificación
es que durante el corte desde el estado inicial hasta el estado crı́tico, se presenta una gran
deformación y el suelo es remoldeado por la deformación cortante. Entonces, el suelo olvida
su estado inicial y el nuevo suelo alcanza un estado único, independiente del estado inicial.

Atkinson sugiere que en los suelos granulares gruesos durante el proceso de carga y corte

40
CAPÍTULO 2. ESTADO CRÍTICO MIC 2008-I-34

se fracturan algunas partı́culas. Por ésto es necesario llevarlos a grandes deformaciones

(> 1000 kPa) para identificar todo el rango del comportamiento.

El concepto fundamental en el estado crı́tico es que una única superficie de falla existe en
0
el espacio q 0 − p − ν(ν = 1 + e) (Figura 2.6). Esta superficie define la falla de un suelo
independientemente de la historia de carga o la trayectoria de esfuerzos seguida. La falla y
el estado crı́tico son sinónimos, la lı́nea de falla = lı́nea del estado crı́tico (CSL). La lı́nea
0
CSL se puede graficar en un plano tridimensional q − p − ν, pero se usan las proyecciones
0 0 0
en el espacio q − p y ν − p por simplicidad. En el espacio p − q la lı́nea del estado crı́tico
es una lı́nea recta de pendiente M , a partir de la cuál se puede determinar el ángulo de
fricción crı́tico.

0 0
Figura 2.6. Lı́nea del estado crı́tico en el espacio q − p − ν (Atkinson & Bransby [4]).

Mooney et al., 1998 [17] es su artı́culo: Un único estado crı́tico para arenas?, describen
0
que los resultados obtenidos en su investigación indican que la relación entre q y p en el
estado crı́tico es única. Sin embargo, no es única la relación de vacı́os en el estado crı́tico
0
para un esfuerzo efectivo medio (p ) dado. Ésto contradice lo dicho por otros autores que
especifican un único estado crı́tico para suelos granulares y suelos fino (relación de vacı́os
única en el estado crı́tico).

41
CAPÍTULO 2. ESTADO CRÍTICO MIC 2008-I-34

2.4. Superficies de fluencia

La superficie de fluencia en el espacio de esfuerzos es la que separa el estado de esfuerzos

que produce deformaciones elásticas, del estado de esfuerzos que produce deformaciones
plásticas. El tamaño inicial o eje mayor de la superficie de fluencia es determinado por el
0
esfuerzo de preconsolidación (p ). La evidencia experimental indica que una superficie de
fluencia elı́ptica es una aproximación razonable para el comportamiento del suelo. El suelo
se comporta elásticamente para esfuerzos que están dentro de la superficie de fluencia y
presenta deformaciones elastoplásticas para esfuerzos por fuera de la superficie de fluencia.

Para los metales Tresca y Von Mises definieron dos criterios por falla elástica. El criterio de
Tresca establece que la fluidización empieza cuando el esfuerzo cortante máximo alcanza
un valor k. Ésto ocurre cuando el esfuerzo principal satisface la siguiente ecuación:

Max(|σa − σb |, |σb − σc |, |σc − σa |) = 2k

Esta ecuación puede ser representada en el espacio de esfuerzos principales como una
superficie prismática con una sección hexagonal centrada en el eje hidrostático (σa = σb =
σc ) (Figura 2.7a). Cuando el estado de esfuerzos es representado como un punto dentro
de esta superficie, el comportamiento del material es elástico; cuando es descrito por un
punto en la superficie, el material fluye. Un estado de esfuerzos fuera de la superficie es
imposible de alcanzar.

El criterio de Von Mises establece que la fluencia plástica inicia cuando se satisface la
siguiente ecuación:
(σa − σb )2 + (σb − σc )2 + (σc − σa )2 = 2σY2

En el espacio de esfuerzos principales esta superficie es equivalente a una superficie cilı́ndri-

ca (Figura 2.7)b. La superficie de Von Mises coincide con la superficie de Tresca en los
bordes.

Las superficies de fluencia de Mohr-Coulomb y Drucker-Prager fueron sugeridas para suelos.

El ingeniero Frances Coulomb introdujo en su análisis de empujes activos en una pared
de retención, la condición de falla para el suelo (usualmente llamada el criterio de falla de
Mohr-Coulomb):
τ = c + σ tan ϕ

Hoy la ingenierı́a geotécnica prefiere escribir esta ecuación en términos de esfuerzos efecti-
vo. Este criterio de fluencia es equivalente a una pirámide hexagonal irregular en esfuerzos
efectivos principales (Figura 2.7c). Drucker y Prager desarrollaron basándose en la super-

42
CAPÍTULO 2. ESTADO CRÍTICO MIC 2008-I-34

ficie de Mohr-Coulomb una superficie de fluencia cónica como la de la Figura 2.7d.

(a) Tresca (b) Von Mises

(c) Mohr-Coulomb (d) Drucker - Prager

Figura 2.7. Superficies de fluencia (Britto & Gunn [19])

2.5. Resistencia del suelo en el estado crı́tico medido en en-

sayos triaxiales

El ensayo más común y usado para determinar el comportamiento del suelo en el estado
crı́tico, es el ensayo triaxial. Este ensayo se realiza con una muestra cilı́ndrica sujeta a es-
fuerzo axial y radial. La presión de poros y el volumen de la muestra pueden ser controlados
y medidos independientemente y la deformación de la muestra puede ser determinada. En
el ensayo triaxial el suelo alcanza el estado crı́tico o último cuando continúa deformándose
con esfuerzo y volumen constante.

43
CAPÍTULO 2. ESTADO CRÍTICO MIC 2008-I-34

Ensayo triaxial no drenado: en muestras de suelo NC o arenas sueltas, durante la fase

0
de corte se generan excesos de presión de poros positiva. Entonces, en el espacio q − p la
trayectoria de esfuerzos efectivos se desplaza hacia la izquierda de la trayectoria de esfuerzos
totales (Figura 2.8a). En este tipo de ensayo durante la fase de corte la relación de vacı́os
0
no cambia y el esfuerzo efectivo (p ) decrece hasta llegar a la lı́nea del estado crı́tico y
luego sigue por la lı́nea del estado crı́tico. En muestras de suelo SC o arenas densas, se
generan excesos de presión de poros negativa; la trayectoria de esfuerzos efectivos en el
0
espacio q − p se desplaza a la derecha de la trayectoria de esfuerzos totales.

(a) (b)

Figura 2.8. Ensayo triaxial no drenados en muestras normalmente consolidadas

La respuesta de un suelo depende de la trayectoria de esfuerzos efectivos. Si en el espacio

0
q−p la pendiente de esta trayectoria es menor que la pendiente de la lı́nea del estado crı́tico
(CSL), el suelo no presentará falla cortante, porque la trayectoria de esfuerzos totales no
intersectarı́a la lı́nea CSL [5].

Ensayo triaxial drenado: en un ensayo triaxial drenado no se generan excesos de presión de

0
poros en la muestra, por lo tanto, la trayectoria en el espacio q−p es una lı́nea con pendiente
1:3 (Figura 2.9a), hasta llegar a la lı́nea del estado crı́tico y continúa su trayectoria por
esta lı́nea. En este tipo de ensayo hay cambio en la relación de vacı́os de la muestra como
se ve en la Figura 2.9.

44
CAPÍTULO 2. ESTADO CRÍTICO MIC 2008-I-34

(a) (b)

Figura 2.9. Ensayo triaxial drenados en muestras normalmente consolidadas (Atkinson

& Bransby [4])

0
La lı́nea del estado crı́tico en el espacio q − p está dada por la ecuación 2.4 y en el espacio
0
ν − ln p por la ecuación 2.5. En la Figura 2.10 se presentan estos dos planos y la lı́nea CSL.

0 0
q f = M pf (2.4)
0
υf = Γ − λ ln pf (2.5)

(a) (b)

Figura 2.10. Lı́nea del estado crı́tico obtenida de ensayos triaxiales (Atkinson [3])

45
CAPÍTULO 2. ESTADO CRÍTICO MIC 2008-I-34

La pendiente M es equivalente al ángulo de fricción crı́tico:

Figura 2.11. Cı́rculo de Mohr

0 0
(σ1 − σ3 )/2
sen ϕc = 0 0
(σ1 + σ3 )/2

0 0 0 0
σ1 − σ3 = σ1 sen ϕc + σ3 sen ϕc

 
0 0 1 + sen ϕc
σ1 = σ3
1 − sen ϕc

0 0 0
q = σ1 − σ3

 
0 0 0 1 + sen ϕc 0
σ1 − σ3 = σ3 − σ3
1 − sen ϕc

 
0 0 2 sen ϕc
q = σ3
1 − sen ϕc
0 0
0 σ1 + 2σ3
p =
3
!
0 1 + sen ϕc 0
σ3 + 2σ3
0
1 − sen ϕc
p =
3
0
0 σ3 (3 − sen ϕc )
p =
3 (1 − sen ϕc )
0
q
M=
p0

46
CAPÍTULO 2. ESTADO CRÍTICO MIC 2008-I-34

!
0 2 sen ϕc
σ3
1 − sen ϕc
M= 0
σ3 (3 − sen ϕc )
3 (1 − sen ϕc )

6 sen ϕc
M=
3 − sen ϕc

3M
sin ϕc = (2.6)
6+M

47
Capı́tulo 3

Propiedades de los suelos

granulares

Las caracterı́sticas de una partı́cula reflejan su origen y proceso de formación, y afectan el

comportamiento del suelo. Por ejemplo, la geometrı́a de un grano depende de su composi-
ción quı́mica y del proceso de formación.

Las caracterı́sticas fı́sicas de los suelos granulares son determinadas por el tamaño de la
partı́cula, forma, textura de la superficie y distribución de tamaños. Estas caracterı́sticas
de los granos influyen en las propiedades esfuerzo-deformación y resistencia del suelo. Por
ejemplo, partı́culas planas y elongadas pueden desarrollar orientaciones y ésto puede incidir
en las propiedades anisotrópicas del suelo [18]. Por otro lado, la composición mineral de
los granos determina la dureza, el clivaje y la resistencia a la degradación fı́sica y quı́mica.

El comportamiento mecánico de materiales granulares es gobernado principalmente por su

estructura y por los esfuerzos efectivos aplicados. La estructura depende del arreglo de las
partı́culas, densidad y anisotropı́a.

El tamaño de las partı́culas, forma y distribución, junto con el arreglo y contacto de los
granos, comprenden la fábrica del suelo.

48
CAPÍTULO 3. PROPIEDADES DE LOS SUELOS GRANULARES MIC 2008-I-34

3.1. Tamaño y Distribución de Partı́culas

El tamaño de una partı́cula puede ser dimensionado por el diámetro de un circulo de igual
área que la partı́cula (Fig 3.1a) ó por medio de dos diámetros frecuentemente usados:
diámetro de Feret (Fig. 3.1b) y diámetro de Martı́n (Fig 3.1c).

Figura 3.1. Diámetro de una partı́cula [24]

El tamaño de los granos es usado para clasificar los suelos en: gravas, arenas, limos o arcillas.
Existen diferentes sistemas de clasificación que definen los tamaños de las partı́culas que
delimitan estos suelos. El tamaño que separa las gravas de las arenas es 4.74 mm (Tamiz
#4 - 4 aberturas/pulgada), y las arenas de los materiales finos es 0.075 mm (Tamiz #200).

En la Tabla 3.1 se presentan los rangos promedios de tamaños de partı́culas para varios
materiales.

Tabla 3.1. Forma de partı́cula, tamaño y superficie especı́fica para diferentes materiales
[24].

Material Forma de la Longitud o diámetro promedio Superficie

partı́cula de la partı́cula [m] especı́fica [m2 /g]
Montmorillonita placa 1x10−9 - 500x10−9 400 - 800
Atapulguita tubo 4 x10−9 - 5 x10−9 ≈ 300
Illita placa 0.1x10−6 - 2x10−6 80 - 100
Kaolinita placa 1 x10−6 - 3 x10−6 10 - 20
Limo 2 x10−6 - 60x10−6 0.04 - 1.1
Arena 60x10−6 - 2 x10−3 (1-40)10−3
Grava > 2x10−3 < 1x10−3

El tamaño de las partı́culas puede ser determinado mediante una base estadı́stica (gra-
nulometrı́a ó hidrómetro), o por medio de una base individual (microscopio). La Tabla
3.2 presenta diferentes métodos para medir de acuerdo con el correspondiente rango de
tamaños de partı́culas.

49
CAPÍTULO 3. PROPIEDADES DE LOS SUELOS GRANULARES MIC 2008-I-34

Tabla 3.2. Métodos para medir tamaños de partı́culas [24].

Método Diámetro Partı́cula [µm]

Tamizado 5 - 1000
Sedimentación 0.1 - 100
Contador Coulter 0.5 - 100
Microscopio de luz >1
Microscopio de escanéo de electrones (SEM) >0.1
Microscopio de transmisión electrones (TEM) >0.001
Dispersión de luz 0.1 - 1000
Rayos X <0.1

Una masa de suelo puede estar compuesta por partı́culas de diferentes tamaños como se
muestra en la Figura 3.2. El análisis granulométrico es una forma de determinar los tamaños
de las partı́culas que conforman una masa de suelo y su distribución. Este análisis se
representa por medio de una curva (curva granulométrica). Pero un análisis granulométrico
no da información sobre la forma de las partı́culas ni la textura.

Figura 3.2. Tamaño de diferentes granos en un suelo [18].

El análisis granulométrico de la arena del Guamo se presenta en la Figura 3.3. En esta

curva se observa que el material tal como llega del rı́o, sin ninguna clasificación, presenta

50
CAPÍTULO 3. PROPIEDADES DE LOS SUELOS GRANULARES MIC 2008-I-34

un porcentaje (aproximadamente el 20 %) de grava; el material restante es arena (fina ≈

20 %, media ≈ 40 % y gruesa ≈ 20 %).

Arena Grava
Limo
Fina Media Gruesa Fina Gruesa
100

80
Porcentaje que pasa [%]

60

40

20

0
0.01 0.1 1 10 100
Diametro de partícula [mm]

Figura 3.3. Curva granulométrica - arena del Guamo

La pendiente de la curva granulométrica es caracterizada por los coeficientes de uniformidad

Cu y de curvatura Cc :
d60
Cu =
d10
d230
Cc =
d10 d60

Los coeficientes de uniformidad y curvatura de la arena del Guamo son:

1,93
Cu = = 6,7
0,29

0,722
Cc = = 0,9
(0,29)(1,93)

De acuerdo con el Sistema Unificado de Clasificación de Suelos (USCS) la arena del Guamo

51
CAPÍTULO 3. PROPIEDADES DE LOS SUELOS GRANULARES MIC 2008-I-34

se clasifica como una arena bien grada (SW).

3.2. Forma de la Partı́cula

La forma de la partı́cula es una caracterı́stica inherente del suelo que influye en su compor-
tamiento mecánico. El análisis de esta caracterı́stica depende de la escala, como se observa
en la Figura 3.4. A mayores escalas la morfologı́a de la partı́cula puede ser descrita como
esférica, redondeada, tubular, aplanada, elı́ptica, alargada, etc. A escalas más pequeñas se
puede analizar la textura de la partı́cula, como la rugosidad de la superficie y redondez de
las esquinas.

Figura 3.4. Caracterización de la forma de la partı́cula dependiendo de la escala [18]

La redondez y esfericidad son usadas para describir la forma de una partı́cula como se
observa en la Tabla 3.3 y Figura 3.5.

Tabla 3.3. Clasificación de la forma de una partı́cula [24].

Clasificación Redondez
Muy angular 0.12 - 0.17
Angular 0.17 - 0.25
Sub angular 0.25 - 0.35
Sub redondeado 0.35 - 0.49
Redondeado 0.49 - 0.70
Bien redondeado 0.70 - 1.00

52
CAPÍTULO 3. PROPIEDADES DE LOS SUELOS GRANULARES MIC 2008-I-34

Figura 3.5. Forma de las partı́culas [18].

La esfericidad de una partı́cula está definida como la relación del diámetro de una esfera
de igual volumen que la partı́cula, con el diámetro de la esfera circunscrita. La redondez
es una medida de los ángulos o filos de los bordes de la partı́cula, y está definida como la
relación entre el radio promedio de curvatura de los bordes de la partı́cula ri , con el radio
de la esfera máxima que puede ser inscrita rmax (Figura 3.6). En la Tabla 3.4 se presentan
diferentes definiciones de esfericidad y redondez. La textura de la superficie se refiere
a los rasgos en escala micro como marcas y cicatrices que aparecen en la superficie de la
partı́cula.

Figura 3.6. Redondez de una partı́cula [24]

53
CAPÍTULO 3. PROPIEDADES DE LOS SUELOS GRANULARES MIC 2008-I-34

Tabla 3.4. Métodos para caracterizar la forma de la partı́cula [18]

La forma de las partı́culas de la arena del Guamo se determinó mediante el análisis de

imágenes obtenidas con estereoscopio (aumento de 10 veces el tamaño de la partı́cula).
En la Figura 3.9 se observan algunas de estas imágenes. Empleando el método descrito
en la Tabla 3.4 - Redondez 1 y estas imágenes, se determinó la redondez de 9 partı́culas
caracterı́sticas y representativas de la arena del Guamo. Posteriormente con este resultado
se clasificaron las partı́culas de acuerdo con la Tabla 3.3. Un resumen de los valores de

54
CAPÍTULO 3. PROPIEDADES DE LOS SUELOS GRANULARES MIC 2008-I-34

redondez y clasificación de cada una de las partı́culas analizadas se presenta en la Tabla

3.5. De estos resultados se puede analizar que la forma de las partı́culas de la arena del
Guamo varı́a desde sub redondeado a muy redondeado.

(a) (b)

(c) (d)

Figura 3.7. Forma de las partı́culas - arena del Guamo

Tabla 3.5. Clasificación de partı́culas - Arena del Guamo

P
Partı́cula # ri /N [mm] rmax [mm] Redondez Clasificación
1 0.103 0.188 0.55 Redondeado
2 0.158 0.206 0.77 Bien Redondeado
3 0.077 0.150 0.52 Redondeado
4 0.184 0.250 0.74 Bien Redondeado
5 0.133 0.049 0.53 Redondeado
7 0.021 0.050 0.42 Sub Redondeado
8 0.016 0.039 0.42 Sub Redondeado
9 0.016 0.038 0.43 Sub Redondeado
10 0.021 0.038 0.56 Redondeado

55
CAPÍTULO 3. PROPIEDADES DE LOS SUELOS GRANULARES MIC 2008-I-34

La forma de las partı́culas influyen en el comportamiento del suelo, a continuación se

presentan varios ejemplos de esta influencia:

Suelos que contienen partı́culas angulares tienden a ser menos densos que suelos con
partı́culas redondeadas.

Materiales débiles y angulares pueden presentar mayor fracturamiento durante la

compresión o compactación.

La relación de vacı́os e está en función de la distribución de tamaños de las partı́culas,

e decrece cuando el rango de tamaños de partı́culas incrementa (Cu incrementa). La
relación de vacı́os también se afectada por la forma de la partı́cula, e incrementa
entre más angular sea la partı́cula.

El ángulo de fricción aumenta cuando se incrementa la angularidad de las partı́culas,

posiblemente esto se debe al incremento en el número de coordinación.

Santamarina y Cho [26] establecieron una relación entre el ángulo de reposo ϕrep y las
redondez de las partı́culas R (Figura 3.8), en donde se puede analizar que a mayor an-
gularidad de las partı́culas (menor redondez) el ángulo de reposo es mayor. El ϕrep se
obtiene por el método de Santamarina [25] y R de la forma que se describe en la Tabla 3.4
- Redondez 1.

Figura 3.8. Angulo de reposo en función de la redondez de las partı́culas [26]

Santamarina y Cho proponen esta relación lineal entre el ángulo de reposo y la redondez
de la partı́cula:

ϕreposo = 42 − 17R (3.1)

56
CAPÍTULO 3. PROPIEDADES DE LOS SUELOS GRANULARES MIC 2008-I-34

3.3. Análisis Fractal

El concepto de dimensión fractal (D) introducido por Mandelbrot (1977) puede explicarse
como una medición de formas (lı́nea, superficie o cuerpo) por medio de elementos cada
vez más pequeños. Por ejemplo para determinar la dimensión D de una lı́nea se empieza
a recorrer ésta con un segmento de cualquier longitud rn (longitud fija) y se establece el
número de segmentos necesarios para recorrerla completamente Nn ; luego se repite esto
con lı́neas de menor dimensión. De acuerdo con Mandelbrot (1977) y Turcotte (1992), si
se desarrolla una relación lineal entre los valores de Nn y rn cuando se grafican en escala
log-log, la partı́cula representa una dimensión fractal D que corresponde al valor absoluto
de la pendiente de esta relación lineal.

Geométricamente un segmento tiene dimensión uno, un cı́rculo o cuadrado tiene dimensión

dos, y una esfera o cubo tiene dimensión tres. Por este motivo la dimensión fractal para una
lı́nea de cualquier forma varı́a entre 1 y 2, y para una superficie entre 2 y 3. En la Figura
3.9 se observa el análisis de una lı́nea, un cuadrado y un cubo por medio de secciones cada
vez más pequeñas para determinar su dimensión fractal.

(a) Dimensión fractal = 1 (1D)

(b) Dimensión fractal = 2 (2D)

(c) Dimensión fractal = 3 (3D)

Figura 3.9. Dimensión fractal

57
CAPÍTULO 3. PROPIEDADES DE LOS SUELOS GRANULARES MIC 2008-I-34

La dimensión fractal puede emplearse para medir la forma de la superficie (rugosidad) de

una partı́cula. Este análisis (1D) se realiza por el contorno de la partı́cula (perı́metro) tal
como lo presenta Vallejo, 1995 [29]:
La Figura 3.10 muestra el perfil de dos partı́culas que tienen área semejante pero diferentes
perfiles. La Figura 3.10(a) muestra en dos dimensiones el perfil de una partı́cula elipsoidal
lisa y la Figura 3.10(b) muestra el perfil de una partı́cula elipsoidal rugosa. La longitud del
perı́metro Ln de estas partı́culas se mide usando una longitud fija rn . Tres segmentos de
longitudes diferentes han sido empleados para medir ambas partı́culas y los resultados se
presentan en la Tabla 3.6. De acuerdo con estos resultados puede concluirse que a menor
longitud del segmento rn mayor es la longitud del perı́metro de la partı́cula Ln , esto se
debe a que la forma del perfil de la partı́cula se detalla mejor entre más pequeña la longitud
del segmento.

Figura 3.10. Análisis fractal, (a) perfil elipsoidal (b) perfil rugoso [29]

58
CAPÍTULO 3. PROPIEDADES DE LOS SUELOS GRANULARES MIC 2008-I-34

Tabla 3.6. Dimensión fractal calculada para las partı́culas de la Figura 3.10 [29]

Perfil Longitud Número de Longitud del

Segmento rn segmentos Nn perfil
1 27.0 27.0
A 2 13.1 26.2
4 6.3 25.2
1 31.4 31.4
B 2 14.5 29.0
4 6.8 27.2

Los valores de Nn y rn para los perfiles de las dos partı́culas son dibujados en la Figura
3.11. En esta figura se observa que las tendencias son dos lı́neas rectas con pendientes
negativas; el valor absoluto de esta pendiente representa la dimensión fractal D. Para el
perfil elipsoidal liso A, el valor de D es 1.0498. Para el perfil rugoso B el valor de D es de
1.1036. La dimensión fractal se incrementa cuando la rugosidad de la partı́cula incrementa.
Partı́culas con un perfil más complejo que el perfil B podrı́an tener dimensión fractal cerca
de 2 (Turcotte, 1992).

Figura 3.11. Relación entre Nn y rn Vallejo, 1995 [29]

La dimensión fractal puede ser usada para medir la rugosidad de una partı́cula granular,
ya que el recorrido con segmentos se hace por el perı́metro de la partı́cula. La dimensión
D no indica angularidad o redondez de la partı́cula, por ejemplo una esfera y un cuadrado
tienen ambos dimensión fractal D=1 pero diferente redondez.

59
CAPÍTULO 3. PROPIEDADES DE LOS SUELOS GRANULARES MIC 2008-I-34

Análisis fractal - arena del Guamo: tres partı́culas de arena del Guamo se analizaron
con dimensión fractal (1D). En las Figuras 3.12, 3.13 y 3.14 se presentan las 3 partı́culas
con los segmentos de diferentes dimensiones para el análisis.

(a) rn =200 (b) rn =100 (c) rn =50

Figura 3.12. Dimensión fractal partı́cula N o 1 arena del Guamo

(a) rn =50 (b) rn =25 (c) rn =10 (d) rn =5

Figura 3.13. Dimensión fractal partı́cula N o 2 arena del Guamo

(a) rn =200 (b) rn =100 (c) rn =50 (d) rn =25

Figura 3.14. Dimensión fractal partı́cula N o 3 arena del Guamo

La Figura 3.15 presenta la relación lineal en escala log-log entre la longitud del segmento
y el número de segmentos necesarios para recorrer el perı́metro de cada una de las tres

60
CAPÍTULO 3. PROPIEDADES DE LOS SUELOS GRANULARES MIC 2008-I-34

partı́culas analizadas. La Tabla 3.7 presentan los resultados de este análisis y la dimensión
fractal de cada grano de arena. Dimensión Fractal - Arena del Guamo

100 Partícula 1 D = 1.019

Partícula 2 D = 1.047
Partícula 3 D = 1.049

Número de segmentos Nn [-]

10

10 100
Longitud del segmento rn [mm]

Figura 3.15. Relación entre Nn y rn arena del Guamo

Tabla 3.7. Dimensión fractal - arena del Guamo

Partı́cula Longitud Número de Dimensión

No Segmento rn segmentos Nn Fractal D
200 8.9
1 100 18.3 1.019
50 36.6
50 5.4
2 25 11.3 1.047
10 29.3
5 60.6
200 6.1
3 100 12.8 1.049
50 26.4
25 54.1

Las imágenes de las partı́culas de la arena del Guamo se obtuvieron por medio de un
estereoscopio con escala de 10 aumentos. Por lo tanto, la dimensión D calculada de cada
partı́cula es una aproximación de la rugosidad (la escala no es suficiente para detallar bien
la superficie de la partı́cula). Para un análisis más exacto de la rugosidad de las partı́culas,
el aumento del estereoscopio debe ser mayor. Con base en los resultados obtenidos y de
acuerdo con Vallejo, 1995 [29], estas partı́culas presentan un perfil no muy complejo puesto

61
CAPÍTULO 3. PROPIEDADES DE LOS SUELOS GRANULARES MIC 2008-I-34

que D esta muy cerca de 1.

62
Capı́tulo 4

Sensibilidad del modelo

hipoplástico con ϕc

Como se dijo anteriormente uno de los objetivos de esta investigación es determinar si el

ángulo de fricción crı́tico (ϕc ) es un parámetro del suelo o una variable de estado. Esta
pregunta surge porque algunos modelos constitutivos como el modelo hipoplástico trabajan
con ϕc como un parámetro del suelo. Por este motivo el primer análisis realizado fue
determinar la sensibilidad del modelo hipoplástico con la variación de ϕc . Para este análisis
se realizaron varias simulaciones de ensayos triaxiales CU (Consolidado - No drenado) y
CD (consolidado - drenado), con los parámetros de la arena del Guamo y Karlsruhe en
estado de densidad suelto, medio y denso. Estas simulaciones se realizaron con el programa
de análisis hipoplástico “ElementTest” (Código del programa: Roberto Cudmani, UMAT:
Andrzej Niemunis).

4.1. Modelo constitutivo hipoplástico

Las ecuaciones hipoplásticas han evolucionado con el tiempo. La primera versión fue publi-
cada por Kolymbas, 1977 [13], y usaba como una única variable de estado, el tensor actual
de esfuerzos de Cauchy T. Esta ecuación demostró ser muy compleja y no muy exitosa en
sus resultados.

Los primeros modelos no tuvieron en cuenta la influencia del nivel de esfuerzos en el

comportamiento volumétrico (barotropı́a) y el de la densidad en los cambios de volumen,
tales como la dilatancia o contractancia (picnotropı́a). Por lo tanto, este modelo no podı́a

63
CAPÍTULO 4. SENSIBILIDAD DEL MODELO HIPOPLÁSTICO CON ϕC MIC 2008-I-34

describir las diferencias en los ángulos de fricción y la rigidez entre materiales sueltos y
densos [14].

Uno de los primeros modelos desarrollados incorporó cuatro términos tensoriales, acom-
pañados de cuatro parámetros de material C1 , C2 , C3 y C4 [32]:

tr(TD) T2 T∗2
T̊ = C1 (trT)D + C2 + C3 √ + C4 √ (4.1)
tr(D) tr(T) trD2 tr(T) trD2

donde T∗ representa el esfuerzo desviador:

1
T∗ = T − (trT)1 (4.2)
3

Una forma de presentación de la ecuación más simplificada pero más general se introdujo,
separando los términos lineales y no lineales de la ecuación. Esta condición de linealidad
está relacionada con el término de velocidad de deformación D:

T̊ = L : D + NkDk (4.3)

donde L es un tensor de cuarto orden (que agrupa los términos lineales y es definido
positivo) asociado a D; N es un tensor de segundo orden asociado con los términos no
√
lineales, y kDk representa la norma Euclidiana del tensor de elongación (kDk = trD2 ).

Esta forma general de la ecuación logra superar las dificultades observadas en los modelos
iniciales como las formas de corazón o de “8” de las envolventes de respuesta que se obtenı́an
al utilizar términos no lineales diferentes a kDk. El término no lineal NkDk modifica la
rigidez, lo cual implica que la ecuación hipoplástica no requiere de definir explı́citamente
un término para la rigidez en carga y otro para la rigidez en descarga.

Se requirieron plantear extensiones adicionales al modelo original para incluir los efectos
de barotropı́a, picnotropı́a y para corregir los frecuentes errores observados durante la
implementación del modelo como el excesivo ratcheting o el gran aumento de la presión de
poros bajo condiciones cı́clicas.

El modelo propuesto por Von Wolffersdorff produce una superficie de fluencia idéntica
a la propuesta originalmente por Matsuoka y Nakai [20]. La derivación de este modelo
comenzó a partir de la ecuación propuesta por Gudehus y Bauer, en la cual se consideraron
dos variables de estado: el esfuerzo actual y la relación de vacı́os (T̊ = h(T, e, D)). Esta
relación constitutiva requiere de cuatro constantes del material que permiten obtener una
descripción realista del comportamiento de los materiales granulares en un amplio rango

64
CAPÍTULO 4. SENSIBILIDAD DEL MODELO HIPOPLÁSTICO CON ϕC MIC 2008-I-34

de esfuerzos y densidades [31].

La velocidad de la segunda variable de estado e se expresa como:

ė = (1 + e)trD (4.4)

La función h puede expresarse como:

h := fb (trT, e)fe (trT, e)(L(T̂, D) + fd (trT, e)N(T̂)kDk) (4.5)

donde T̂ es un tensor adimensional de esfuerzos definido por T̂ = T/trT; L es un término

lineal en D expresado como L = L(T̂)[D], siendo L un tensor de cuarto orden.

La función h presenta dos caracterı́sticas importantes: la primera, la respuesta incremental

del material depende de la relación de vacı́os e; la segunda, la rigidez incremental kFk/kDk
desaparece cuando T = 0 [31]. La Ec. [4.5] separa la barotropı́a y picnotropı́a (que están
descritas por los factores fb , fd y fe ) de la parte tensorial de la función de respuesta.

Un avance importante en el desarrollo de los modelos hipoplásticos fue propiciado por la

inclusión de algunos de los conceptos de la mecánica de suelos en el estado crı́tico. Bauer
modificó la ecuación constitutiva hipoplástica para el estado crı́tico de tal forma que la
superficie de fluencia circunscribiera la pirámide de Mohr–Coulomb [31].

1 h 2 ∗
i
T̊ = fb (trTc , ec ) F D + a2 tr(T̂c D)T̂c + aF (T̂c + T̂c )kDk (4.6)
tr(Tˆ2 )
c

donde T̂c representa la relación de esfuerzos crı́tica con:

√
3(3 − sin ϕc )
a= √ (4.7)
2 2 sin ϕc

El factor 1/tr(Tˆ2c ) considera el hecho que la compresibilidad bajo relaciones de esfuerzos

grandes es evidentemente más alta que bajo esfuerzos isotrópicos [20].

La Ec.[4.6] está caracterizada por el factor de rigidez fd y por un único parámetro de

material, el ángulo de fricción crı́tico ϕc .

F es una función de esfuerzos definida como parte de la función de fluencia que satisface
∗
la condición lı́mite de Matsuoka-Nakai. F es función de las invariantes de T̂ ó T̂ . F es

65
CAPÍTULO 4. SENSIBILIDAD DEL MODELO HIPOPLÁSTICO CON ϕC MIC 2008-I-34

igual a 1 para una compresión triaxial ó para condiciones de esfuerzos isotrópicos.

s
1 2 − tan2 ψ 1
F = tan2 ψ + √ − √ tan ψ (4.8)
8 2 + 2 tan ψ · cos 3ϑ 2 2

con las siguientes invariantes (el significado de cada ángulo se ilustra en la Figura [4.1]):
√
tan ψ = 3kT̂∗c k (4.9)

√ tr(Tˆ∗3 c )
cos 3ϑ = − 6 ∗2 (4.10)
[tr(T̂c )]3/2

Figura 4.1. Representación geométrica de las invariantes tan ψ y cos 3ϑ en el espacio de

los esfuerzos principales (Tomado de [31])

No obstante, la Ec.[4.6] no tenı́a en cuenta el efecto de la picnotropı́a. Por lo tanto, Gudehus

66
CAPÍTULO 4. SENSIBILIDAD DEL MODELO HIPOPLÁSTICO CON ϕC MIC 2008-I-34

y Bauer propusieron una extensión a la ecuación de Wu:

1 n 2 ∗
o
T̊ = fb fe F D + a2 tr(T̂D)T̂ + fd aF (T̂ + T̂ )kDk (4.11)
tr(Tˆ2 )

La ley de compresión para las relaciones de vacı́os caracterı́sticas formulada por Bauer es:

trT n
   
ei ec ed
= = = exp − − (4.12)
ei0 ec0 ed0 hs

donde n es la sensibilidad a la presión del esqueleto granular, hs es una constante conocida

como la dureza granular y ei (relación de vacı́os para la mı́nima densidad), ec (relación de
vacı́os crı́tica) y ed (relación de vacı́os con la máxima densidad) dependen de la presión
media, tal como se muestra en la Figura 4.2.

Figura 4.2. Rango de las posibles relaciones de vacı́os como función del nivel de esfuerzos
[31]

Las funciones de picnotropı́a son:

 e β
c
fe = (4.13)
e

 α
e − ed
fd = (4.14)
ec − ed
donde α y β son los exponentes de picnotropı́a.

Cuando ec > e > ed , el factor de picnotropı́a fd controla la máxima relación de esfuerzos

T̂p en el pico ası́ como la dilatancia asociada a esta relación ; es decir, fd tiene en cuenta la
densidad crı́tica de la arena. El factor de picnotropı́a fe controla la influencia de la relación
de vacı́os e sobre la rigidez incremental [31].

67
CAPÍTULO 4. SENSIBILIDAD DEL MODELO HIPOPLÁSTICO CON ϕC MIC 2008-I-34

La función de barotropı́a está definida como:

β  1−n   −1
√ ei0 − ed0 α
  
hs 1 + ei ei0 trT 2
fb = − 3 + a − 3a (4.15)
n ei ec0 hs ec0 − ed0

donde ec0 , ed0 , ei0 son las relaciones de vacı́os crı́tica, máxima y mı́nima respectivamente
pra un nivel de esfuerzos igual a cero (0). hs es la presión de referencia del esqueleto
granular denominada como la dureza granular y n es el exponente de la ley de compresión.
Estas dos constantes se obtienen de ensayos oedométricos.

El factor escalar fb · fe /tr(Tˆ2 ) en la Ec. 4.11 escala la magnitud de T̊ y consecuentemente

el tamaño de la envolvente de respuesta [12].

De esta forma, el nuevo modelo requiere de cuatro parámetros que son independientes de
T y e: hs , ϕc (ángulo de reposo del material suelto), ec0 y ed0 . También requiere de cuatro
constantes adicionales n, α, β and ei0 .

4.2. Resultados de ensayos triaxiales en arena de Karlsruhe

Herle y Gudehus, 1999 [11] en su artı́culo presentan los resultados de un ensayo triaxial
drenado con arena de Karlsruhe. El ensayo fue realizado con muestras densas (eo = 0,53)
y presión de cámara variable (Figura 4.3a y c). Los Autores presentan también las simu-
laciones de estos ensayos con los parámetros de la arena de Karlsruhe (Figura 4.3b y d).
Los parámetros de esta arena son [6]:

Ángulo de fricción crı́tico - ϕc = 30o

Relación de vacı́os máxima posible - ei = 1.00
Relación de vacı́os crı́tica - ec = 0.84
Relación de vacı́os mı́nima posible - ed = 0.53
Dureza granular del esqueleto - hs = 5800 MPa
Exponente, sensibilidad del esqueleto ante presión - n = 0.28
Exponente, picnotropı́a - α = 0.13
Exponente, barotropı́a - β = 1.05

68
CAPÍTULO 4. SENSIBILIDAD DEL MODELO HIPOPLÁSTICO CON ϕC MIC 2008-I-34

(a) Resultados de ensayos de laboratorio (b) Resultados de simulaciones

(c) Resultados de ensayos de laboratorio (d) Resultados de simulaciones

Figura 4.3. Ensayos triaxiales CD con arena de Karlsruhe, presión de cámara variable y
eo = 0,53 (Herle & Gudehus [11]).

4.3. Resultados de ensayos triaxiales en arena del Guamo

Pablo Arias, 2006 [2] en su tesis de maestrı́a presenta los resultados de un ensayo triaxial
drenado con arena del Guamo. El ensayos fue realizado con muestras densas (eo = 0,56) y
presión de cámara variable (Figura 4.4a y c). También presenta las simulaciones de estos
ensayos con los parámetros de la arena del Guamo (Figura 4.4b y d). Los parámetros de
la arena del Guamo determinados por Arias [2] son:

Ángulo de fricción crı́tico - ϕc = 30o

Relación de vacı́os máxima posible - ei = 1.15
Relación de vacı́os crı́tica - ec = 1.00

69
CAPÍTULO 4. SENSIBILIDAD DEL MODELO HIPOPLÁSTICO CON ϕC MIC 2008-I-34

Relación de vacı́os mı́nima posible - ed = 0.52

Dureza granular del esqueleto - hs = 4000 MPa
Exponente, sensibilidad del esqueleto ante presión - n = 0.27
Exponente, picnotropı́a - α = 0.17
Exponente, barotropı́a - β = 1.00

(a) Resultados de ensayos de laboratorio (b) Resultados de simulaciones

(c) Resultados de ensayos de laboratorio (d) Resultados de simulaciones

Figura 4.4. Ensayos triaxiales CD con arena del Guamo, presión de cámara variable y
eo = 0,56 (Arias [2]).

70
CAPÍTULO 4. SENSIBILIDAD DEL MODELO HIPOPLÁSTICO CON ϕC MIC 2008-I-34

4.4. Simulaciones con ElementTest de los ensayos triaxiales

en arena de Karlsruhe y del Guamo

Con el fin de validar los resultados del programa ElementTest para realizar el análisis
de sensibilidad del modelo constitutivo hipoplástico, se reprodujeron las simulaciones de
los ensayos triaxiales en arena de Karlsruhe realizados por Herle y Gudehus, 1999 [11] y
en arena del Guamo realizados por Arias, 2006 [2]. En las Figuras 4.5, 4.6, 4.7 y 4.8 se
presentan estas simulaciones.
Ensayos simulados Elementtest (Cudmani)
Triaxial CD con arena de Karlsruhe - eo = 0.53

Esfuerzo desviador q [kPa]

3000
1000

800
2000
600
500
400
1000
300
200
100
50
0
0 0.04 0.08 0.12

Deformación axial ε1 [-]

(a) (b)

Figura 4.5. Ensayo triaxial CD con arena de Karlsruhe, presión de cámara variable y
eo = 0,53. (a) Simulación Herle & Gudehus. (b) Simulación con ElementTest.

Ensayos simulados Elementtest (Cudmani)

Triaxial CD con arena de Karlsruhe - eo = 0.53
Deformación volumétrica εV [-]

0.06
1000
800
600
500
400
0.04 300
200
100
50

0.02

0 0.04 0.08 0.12

Deformación axial ε1 [-]

(a) (b)

Figura 4.6. Ensayo triaxial CD con arena de Karlsruhe, presión de cámara variable y
eo = 0,53. (a) Simulación Herle & Gudehus. (b) Simulación con ElementTest.

71
CAPÍTULO 4. SENSIBILIDAD DEL MODELO HIPOPLÁSTICO CON ϕC MIC 2008-I-34

Ensayos simulados ElementTest

Triaxial CD con arena del Guamo

Esfuerzo desviador q [kPa]

1200
eo = 0.56

σ3 =300 kPa

800

σ3 =150 kPa
400
σ3 =75 kPa

0
0 0.04 0.08 0.12 0.16 0.2

Deformación axial ε1 [-]

(a) (b)

Figura 4.7. Ensayo triaxial CD con arena del Guamo, presión de cámara variable y
eo = 0,56. (a) Simulación realizada por Arias. (b) Simulación con ElementTest.

0.08 Ensayos simulados ElementTest

Deformación volumétrica εV [-]

Triaxial CD con arena del Guamo

eo = 0.56

0.06

0.04

0.02 σ3 = 75 kPa
σ3 = 150 kPa
σ3 = 300 kPa
0

0 0.04 0.08 0.12 0.16 0.2

Deformación axial ε1 [-]

(a) (b)

Figura 4.8. Ensayo triaxial CD con arena del Guamo, presión de cámara variable y
eo = 0,56. (a) Simulación realizada por Arias. (b) Simulación con ElementTest.

Al comparar los resultados de las simulaciones realizadas con ElementTest con las simu-
laciones realizadas por Herle & Gudehus y Arias, se validan los resultados del programa
ElementTest para realizar los análisis de sensibilidad del modelo con ϕrep .

72
CAPÍTULO 4. SENSIBILIDAD DEL MODELO HIPOPLÁSTICO CON ϕC MIC 2008-I-34

4.5. Simulación de ensayos triaxiales con parámetros de la

arena del Guamo variando ϕc - Densidad media (eo =
0,9)

Ensayo triaxial CU con parámetros de la arena del Guamo - Densidad

media

Las simulaciones se realizaron para una muestra con una relación de vacı́os inicial igual
a 0.9. La muestra se comprimió isotrópicamente con una presión de cámara (σ3 ) igual a
20 kPa y luego se realizó el corte no drenado hasta una deformación axial del 15 %. Estas
simulaciones se realizaron variando ϕc con incremento de 2o , desde 20o hasta 40o . Los
resultados de estas simulaciones se observan en la Figura 4.9 y Tabla 4.1. La variación
porcentual y numérica de ϕc y q en la Tabla 4.1 se obtuvo comparando estos valores
obtenidos de cada simulación, con los resultados de la simulación realizada con el valor del
parámetro ϕc real de la arena del Guamo (ϕc = 30o ).

En la Figura 4.9a se observa que para las condiciones de este ensayo (triaxial CU; 1 = 15 %)
la trayectoria de esfuerzos en la fase de corte llega hasta un p0 mı́nimo y luego sigue por la
lı́nea del estado crı́tico. Para verificar si realmente esta tendencia es por la lı́nea del estado
crı́tico, se determinó la pendiente M de cada trayectoria y a partir de ésta se calculó ϕc
con la ecuación 4.16.

3M
sin ϕc = (4.16)
6+M

73
CAPÍTULO 4. SENSIBILIDAD DEL MODELO HIPOPLÁSTICO CON ϕC MIC 2008-I-34

120 120

Ensayo triaxial CU - Arena del Guamo

Simulado - Elementtest Cudmani
eo = 0.9 σ3 = 20 kPa

ϕc = 20o
Esfuerzo desviador q [kPa]

ϕc = 20o

Esfuerzo desviador q [kPa]

ϕc = 22o
ϕc = 22o
ϕc = 24o
80 ϕc = 24o 80
ϕc = 26o
ϕc = 26o
ϕc = 28o
ϕc = 28o
ϕc = 30o
ϕc = 30o
ϕc = 32o
ϕc = 32o
ϕc = 34o ϕc = 34o

ϕc = 36o ϕc = 36o

ϕc = 38o ϕc = 38o
40 40
ϕc = 40o ϕc = 40o

Ensayo triaxial CU - Arena del Guamo

Simulado - Elementtest Cudmani
eo = 0.9 σ3 = 20 kPa

0 0
0 40 80 120 0 0.04 0.08 0.12 0.16
Esfuerzo promedio efectivo p' [kPa] Deformación axial ε1

(a) (b)

20

Sensibilidad del modelo constitutivo hipoplástico

Simulado - ElementTest (Cudmani)
Arena del Guamo - triaxial CU
eo = 0.9 σ3 = 20kPa
Esfuerzo desviador q [kPa]

16

12

4
20 30 40
Ángulo de fricción crítico ϕc [ o ]

(c)

Figura 4.9. Simulación ensayo triaxial CU con arena del Guamo variando ϕc - densidad
media.

El ϕc calculado se comparó con el valor de ϕc parámetro del suelo en cada simulación y

las variaciones observadas fueron despreciables.

La sensibilidad del modelo hipoplástico con la variación del parámetro ϕc en un ensayo

triaxial CU con arena del Guamo (densidad media) se observa en la Figura 4.9c. El valor
del esfuerzo desviador q en esta gráfica corresponde al menor valor del esfuerzo promedio

74
CAPÍTULO 4. SENSIBILIDAD DEL MODELO HIPOPLÁSTICO CON ϕC MIC 2008-I-34

p0 en la fase de corte.

Tabla 4.1. Sensibilidad del modelo hipoplástico con la variación de ϕc . Ensayo triaxial
CU con arena del Guamo, densidad media (eo = 0,9)

ϕc q Variación de ϕc Variación de q
[o ] [kPa] [o ] [ %] [kPa] [ %]
20 6.73 10 33 3.7 36
22 7.45 8 27 3.0 29
24 8.23 6 20 2.2 21
26 9.00 4 13 1.5 14
28 9.75 2 7 0.7 7
30 10.46 0 0 0.0 0
32 11.32 2 7 0.9 8
34 11.97 4 13 1.5 14
36 12.79 6 20 2.3 22
38 13.58 8 27 3.1 30
40 14.36 10 33 3.9 37

Ensayo triaxial CD con parámetros de la arena del Guamo - Densidad

media

Las simulaciones se realizaron para una muestra con una relación de vacı́os inicial igual
a 0.9. La muestra se comprimió isotrópicamente con una presión de cámara (σ3 ) igual a
20 kPa y luego se realizó el corte drenado hasta una deformación axial del 15 %. Estas
simulaciones se realizaron variando ϕc con incremento de 2o , desde 20o hasta 40o . Los
resultados de estas simulaciones se observan en la Figura 4.10 y tabla 4.2. La variación
porcentual y numérica de ϕc y q en la Tabla 4.2 se obtuvo comparando estos valores
obtenidos de cada simulación con los resultados de la simulación realizada con el valor del
parámetro ϕc real de la arena del Guamo (ϕc = 30o ).

75
CAPÍTULO 4. SENSIBILIDAD DEL MODELO HIPOPLÁSTICO CON ϕC MIC 2008-I-34

120
120

Ensayo triaxial CD - Arena del Guamo

Ensayo triaxial CD - Arena del Guamo
Simulado - Elementtest Cudmani
Simulado - Elementtest Cudmani
eo = 0.9 σ3 = 20 kPa
eo = 0.9 σ3 = 20 kPa

Esfuerzo desviador q [kPa]

ϕc = 20o
Esfuerzo desviador q [kPa]

ϕc = 22o ϕc = 40o
80
80 ϕc = 24o
ϕc = 38o
ϕc = 26o
ϕc = 28o
ϕc = 36o
ϕc = 30o ϕc = 34o
ϕc = 32o ϕc = 32o
ϕc = 34o
ϕc = 30o
40 ϕc = 28o
40 ϕc = 36o ϕc = 26o
ϕc = 38o ϕc = 24o
ϕc = 22o
ϕc = 40o
ϕc = 20o

0
0 0 0.04 0.08 0.12 0.16 0.2
0 40 80 120
Esfuerzo promedio efectivo p' [kPa] Deformación axial ε1 [-]

(a) (b)

100

Sensibilidad del modelo constitutivo hipoplástico

Simulado - ElementTest (Cudmani)
Arena del Guamo - triaxial CD
eo = 0.9 σ3 = 20kPa
Esfuerzo desviador q [kPa]

80

60

40

20
20 30 40

Ángulo de fricción crítico ϕc [ o ]

(c)

Figura 4.10. Simulación ensayo triaxial CD con arena del Guamo variando ϕc - densidad
media.

La sensibilidad del modelo hipoplástico con la variación del parámetro ϕc en un ensayo

triaxial CD con arena del Guamo (densidad media) se observa en la Figura 4.10c. El valor
del esfuerzo desviador q en esta gráfica corresponde al máximo q en la fase de corte.

76
CAPÍTULO 4. SENSIBILIDAD DEL MODELO HIPOPLÁSTICO CON ϕC MIC 2008-I-34

Tabla 4.2. Sensibilidad del modelo hipoplástico con la variación de ϕc . Ensayo triaxial
CD con arena del Guamo, densidad media (eo = 0,9)

ϕc q Variación de ϕc Variación de q
[o ] [kPa] [o ] [ %] [kPa] [ %]
20 22.23 10 33 22.1 50
22 25.75 8 27 18.6 42
24 29.66 6 20 14.7 33
26 34.02 4 13 10.4 23
28 38.89 2 7 5.5 12
30 44.37 0 0 0.0 0
32 50.57 2 7 6.2 14
34 57.64 4 13 13.3 30
36 65.74 6 20 21.4 48
38 75.11 8 27 30.7 69
40 86.04 10 33 41.7 94

Comparación de ensayos triaxiales CD y CU con parámetros de la arena

del Guamo - Densidad media

120 120

Ensayo triaxial CD - Arena del Guamo Ensayo triaxial CD - Arena del Guamo
Simulado - Elementtest Cudmani Simulado - Elementtest Cudmani
eo = 0.9 σ3 = 20 kPa ϕc = 30o eo = 0.9 σ3 = 20 kPa ϕc = 30o
Esfuerzo desviador q [kPa]
Esfuerzo desviador q [kPa]

CD CD
CU CU
80 CSL 80

40 40

0 0
0 40 80 120 0 0.04 0.08 0.12 0.16 0.2

Deformación axial ε1 [-] Deformación axial ε1 [-]

(a) (b)

Figura 4.11. Simulación de ensayos triaxiales CU y CD con arena del Guamo, densidad
media

La variación de la resistencia del suelo por las condiciones de drenaje en la fase de corte
en un ensayo triaxial con arena del Guamo (densidad media) se observa en la Figura 4.11.

77
CAPÍTULO 4. SENSIBILIDAD DEL MODELO HIPOPLÁSTICO CON ϕC MIC 2008-I-34

Estas gráficas se realizaron con los mismos parámetros del suelo y condiciones del ensayo,
variando únicamente la condición de drenaje en la fase de corte.

4.6. Simulación de ensayos triaxiales con parámetros de la

arena del Guamo variando ϕc - Estado suelto (eo = 1,1)

Ensayo triaxial CU con parámetros de la arena del Guamo - Estado suelto

12
25

Ensayo triaxial CU - arena del Guamo

Ensayo triaxial CU - arena del Guamo Simulado - Elementtest Cudmani
Simulado - Elementtest Cudmani eo = 1.1 σ3 = 20 kPa
20 eo = 1.1 σ3 = 20 kPa

Esfuerzo desviador q [kPa]

Esfuerzo desviador q [kPa]

ϕc = 20
ϕc = 20o
8 ϕc = 30
ϕc = 30o
ϕc = 40
15 ϕc = 40o

10
4

0
0 0 0.04 0.08 0.12 0.16 0.2
0 5 10 15 20 25
Esfuerzo promedio efectivo p' [kPa] Deformación axial ε1 [-]
(a) (b)

Sensibilidad del modelo constitutivo hipoplástico

Simulado - ElementTest (Cudmani)
Arena del Guamo - triaxial CU
Esfuerzo desviador q [kPa]

12
eo = 1.1 σ3 = 20 kPa

4
20 30 40
Ángulo de fricción crítico ϕc
(c)

Figura 4.12. Simulación ensayo triaxial CU con arena del Guamo variando ϕc - estado
suelto

78
CAPÍTULO 4. SENSIBILIDAD DEL MODELO HIPOPLÁSTICO CON ϕC MIC 2008-I-34

Las simulaciones se realizaron para una muestra con una relación de vacı́os inicial igual
a 1.1. La muestra se comprimió isotrópicamente con una presión de cámara (σ3 ) igual a
20 kPa y luego se realizó el corte no drenado hasta una deformación axial del 15 %. Estas
simulaciones se realizaron variando ϕc con incremento de 2o , desde 20o hasta 40o . Los
resultados de estas simulaciones se observan en la Figura 4.12 y Tabla 4.3. La variación
porcentual y numérica de ϕc y q en la Tabla 4.3 se obtuvo comparando estos valores
obtenidos en cada simulación, con los resultados de la simulación realizada con el valor del
parámetro ϕc real de la arena del Guamo (ϕc = 30o ).

La sensibilidad del modelo hipoplástico con la variación del parámetro ϕc en un ensayo

triaxial CU con arena del Guamo (estado suelto) se observa en la Figura 4.12c. El valor
del esfuerzo desviador q en esta gráfica corresponde al máximo q en la fase de corte.

Tabla 4.3. Sensibilidad del modelo hipoplástico con la variación de ϕc . Ensayo triaxial
CU con arena del Guamo, estado suelto (eo = 1,1)

ϕc q Variación de ϕc Variación de q
[o ] [kPa] [o ] [ %] [kPa] [ %]
20 5.14 10 33 2.8 36
22 5.79 8 27 2.3 29
24 6.26 6 20 1.7 22
26 6.83 4 13 1.2 14
28 7.41 2 7 0.6 7
30 7.99 0 0 0.0 0
32 8.57 2 7 0.6 7
34 9.15 4 13 1.2 15
36 9.73 6 20 1.7 22
38 10.31 8 27 2.3 29
40 10.89 10 33 2.9 36

Ensayo triaxial CD con parámetros de la arena del Guamo - Estado suelto

Las simulaciones se realizaron para una muestra con una relación de vacı́os inicial igual
a 1.1. La muestra se comprimió isotrópicamente con una presión de cámara (σ3 ) igual a
20 kPa y luego se realizó el corte drenado hasta una deformación axial del 15 %. Estas
simulaciones se realizaron variando ϕc con incremento de 2o , desde 20o hasta 40o . Los
resultados de estas simulaciones se observan en la Figura 4.13 y Tabla 4.4. La variación
porcentual y numérica de ϕc y q en la Tabla 4.4 se obtuvo comparando estos valores
obtenidos en cada simulación, con los resultados de la simulación realizada con el valor del
parámetro ϕc real de la arena del Guamo (ϕc = 30o ).

79
CAPÍTULO 4. SENSIBILIDAD DEL MODELO HIPOPLÁSTICO CON ϕC MIC 2008-I-34

120
120

Ensayo triaxial CD - arena del Guamo

Ensayo triaxial CD - arena del Guamo Simulado - Elementtest Cudmani
Simulado - Elementtest Cudmani
eo = 1.1 σ3 = 20 kPa
eo = 1.1 σ3 = 20 kPa

Esfuerzo desviador q [kPa]

Esfuerzo desviador q [kPa]

ϕc = 20o ϕc = 40o
ϕc = 22o 80
80
ϕc = 24o ϕc = 38o
ϕc = 26o ϕc = 36o
ϕc = 28o ϕc = 34o
ϕc = 30o ϕc = 32o
ϕc = 32o ϕc = 30o
ϕc = 34o 40 ϕc = 28o
40 ϕc = 26o
ϕc = 36o ϕc = 24o
ϕc = 38o ϕc = 22o
ϕc = 40o ϕc = 20o

0
0 0 0.04 0.08 0.12 0.16 0.2
0 40 80 120
Esfuerzo promedio efectivo p' [kPa] Deformación axial ε1 [-]
(a) (b)

80

Sensibilidad del modelo constitutivo hipoplástico

Simulado - ElementTest (Cudmani)
Arena del Guamo - triaxial CD
eo = 1.1 σ3 = 20 kPa
Esfuerzo desviador q [kPa]

60

40

20

20 30 40

Ángulo de fricción crítico ϕc [ o ]

(c)

Figura 4.13. Simulación ensayo triaxial CD con arena del Guamo variando ϕc - estado
suelto

La sensibilidad del modelo hipoplástico con la variación del parámetro ϕc en un ensayo

triaxial CD con arena del Guamo (estado suelto) se observa en la Figura 4.13c. El valor
del esfuerzo desviador q en esta gráfica corresponde al máximo q en la fase de corte.

80
CAPÍTULO 4. SENSIBILIDAD DEL MODELO HIPOPLÁSTICO CON ϕC MIC 2008-I-34

Tabla 4.4. Sensibilidad del modelo hipoplástico con la variación de ϕc . Ensayo triaxial
CD con arena del Guamo, estado suelto (eo = 1,1)

ϕc q Variación de ϕc Variación de q
[o ] [kPa] [o ] [ %] [kPa] [ %]
20 20.00 10 33 17.7 47
22 22.98 8 27 14.7 39
24 26.22 6 20 11.5 30
26 29.73 4 13 8.0 21
28 33.55 2 7 4.2 11
30 37.72 0 0 0.0 0
32 42.26 2 7 4.5 12
34 47.23 4 13 9.5 25
36 52.68 6 20 15.0 40
38 58.66 8 27 20.9 56
40 65.24 10 33 27.5 73

Comparación de ensayos triaxiales CD y CU con parámetros de la arena

del Guamo - Estado suelto

50
50

Ensayo triaxial - arena del Guamo

Simulado - Elementtest Cudmani
eo = 1.1 σ3 = 20 kPa ϕc = 30o 40
40
Esfuerzo desviador q [kPa]
Esfuerzo desviador q [kPa]

CU
CD
CSL
30 Ensayo triaxial - arena del Guamo
30
Simulado - Elementtest Cudmani
eo = 1.1 σ3 = 20 kPa ϕc = 30o

CU
20
20 CD

10
10

0
0 0 0.04 0.08 0.12 0.16 0.2
0 10 20 30 40 50
Esfuerzo promedio efectivo p' [kPa] Deformación axial ε1 [-]

(a) (b)

Figura 4.14. Simulación de ensayos triaxiales CU y CD con arena del Guamo, estado
suelto

La variación de la resistencia del suelo por las condiciones de drenaje en la fase de corte
en un ensayo triaxial con arena del Guamo (estado suelto) se observa en la figura 4.14.

81
CAPÍTULO 4. SENSIBILIDAD DEL MODELO HIPOPLÁSTICO CON ϕC MIC 2008-I-34

Estas gráficas se realizaron con los mismos parámetros del suelo y condiciones del ensayo,
variando únicamente la condición de drenaje en la fase de corte.

4.7. Simulación de ensayos triaxiales con parámetros de la

arena de Karlsruhe - Estado suelto (eo = 0,9)
Ensayo triaxial CU con parámetros de la arena de Karlsruhe - Estado
suelto
600 400

Ensayo triaxial CU - arena de Karlsruhe

Ensayo triaxial CU - arena de Karlsruhe
Simulado - Elementtest Cudmani ϕc = 20o
Simulado - Elementtest Cudmani
eo = 0.9 σ3 = 500 kPa eo = 0.9 σ3 = 500 kPa ϕc = 25o
ϕc = 28o

Esfuerzo desviador q [kPa]

Esfuerzo desviador q [kPa]

ϕc = 20o 300
ϕc = 30o
ϕc = 25o
400 ϕc = 32o
ϕc = 28o
ϕc = 30o
ϕc = 35o
ϕc = 32o ϕc = 40o
200
ϕc = 35o
ϕc = 40o

200

100

0 0
0 200 400 600 0 0.04 0.08 0.12 0.16

Esfuerzo promedio efectivo p' [kPa] Deformación axial ε1 [-]

(a) (b)

280
Sensibilidad del modelo constitutivo hipoplástico
Simulado - ElementTest (Cudmani)
Arena de Karlsruhe - triaxial CU
Esfuerzo desviador q [kPa]

eo = 0.9 σ3 = 500 kPa

240

200

160

120

20 30 40

Ángulo de fricción crítico ϕc [ o ]

(c)

Figura 4.15. Simulación ensayo triaxial CU con arena de Karlsruhe variando ϕc - estado
suelto

82
CAPÍTULO 4. SENSIBILIDAD DEL MODELO HIPOPLÁSTICO CON ϕC MIC 2008-I-34

Las simulaciones se realizaron para una muestra con una relación de vacı́os inicial igual
a 0.9. La muestra se comprimió isotrópicamente con una presión de cámara (σ3 ) igual a
500 kPa y luego se realizó el corte no drenado hasta una deformación axial del 12 %. Estas
simulaciones se realizaron variando ϕc con incremento de 2o , 5o y 10o . Los resultados de
estas simulaciones se observan en la Figura 4.15 y Tabla 4.5. La variación porcentual y
numérica de ϕc y q en la Tabla 4.5 se obtuvo comparando estos valores obtenidos en cada
simulación, con los resultados de la simulación realizada con el valor del parámetro ϕc real
de la arena de Karlsruhe (ϕc = 30o ).

Tabla 4.5. Sensibilidad del modelo hipoplástico con la variación de ϕc . Ensayo triaxial
CU con arena de Karlsruhe, estado suelto (eo = 0,9)

ϕc q Variación de ϕc Variación de q
[o ] [kPa] [o ] [ %] [kPa] [ %]
20 132.12 10 33 73.1 36
25 168.31 5 17 36.9 18
28 190.39 2 7 14.8 7
30 205.23 0 0 0.0 0
32 220.12 2 7 14.9 7
35 242.51 5 17 37.3 18
40 279.71 10 33 74.5 36

La sensibilidad del modelo hipoplástico con la variación del parámetro ϕc en un ensayo

triaxial CU con arena de Karlsruhe (estado suelto) se observa en la Figura 4.15c. El valor
del esfuerzo desviador q en esta gráfica corresponde al máximo q en la fase de corte.

Ensayo triaxial CD con parámetros de la arena de Karlsruhe variando ϕc

- Estado suelto

Las simulaciones se realizaron para una muestra con una relación de vacı́os inicial igual
a 0.9. La muestra se comprimió isotrópicamente con una presión de cámara (σ3 ) igual a
500 kPa y luego se realizó el corte drenado hasta una deformación axial del 12 %. Estas
simulaciones se realizaron variando ϕc con incremento de 2o , 5o y 10o . Los resultados de
estas simulaciones se observan en la Figura 4.16 y Tabla 4.6. La variación porcentual y
numérica de ϕc y q en la Tabla 4.6 se obtuvo comparando estos valores obtenidos en cada
simulación, con los resultados de la simulación realizada con el valor del parámetro ϕc real
de la arena de Karlsruhe (ϕc = 30o ).

83
CAPÍTULO 4. SENSIBILIDAD DEL MODELO HIPOPLÁSTICO CON ϕC MIC 2008-I-34

ϕc = 40o
1600 1600

Esfuerzo desviador q [kPa]

Esfuerzo desviador q [kPa]

ϕc = 35o
1200 1200

ϕc = 20o ϕc = 32o
ϕc = 25o ϕc = 30o
ϕc = 28o ϕc = 28o
800
800 ϕc = 30o
ϕc = 25o
ϕc = 32o
ϕc = 35o
ϕc = 20o
ϕc = 40o
400
400
Ensayo triaxial CD - arena de Karlsruhe
Ensayo triaxial CD - arena de Karlsruhe Simulado - Elementtest Cudmani
Simulado - Elementtest Cudmani eo = 0.9 σ3 = 500 kPa
eo = 0.9 σ3 = 500 kPa
0
0 0 0.04 0.08 0.12 0.16
0 400 800 1200 1600
Esfuerzo promedio efectivo p' [kPa] Deformación axial ε1 [-]

(a) (b)
2000

Sensibilidad del modelo constitutivo hipoplástico

Simulado - ElementTest (Cudmani)
Arena de Karlsruhe - triaxial CD
Esfuerzo desviador q [kPa]

1600 eo = 0.9 σ3 = 500 kPa

1200

800

400
20 30 40

Ángulo de fricción crítico ϕc [ o ]

(c)

Figura 4.16. Simulación ensayo triaxial CD con arena de Karlsruhe variando ϕc - estado
suelto

La sensibilidad del modelo hipoplástico con la variación del parámetro ϕc en un ensayo

triaxial CD con arena de Karlsruhe (estado suelto) se observa en la Figura 4.16c. El valor
del esfuerzo desviador q en esta gráfica corresponde al máximo q en la fase de corte.

84
CAPÍTULO 4. SENSIBILIDAD DEL MODELO HIPOPLÁSTICO CON ϕC MIC 2008-I-34

Tabla 4.6. Sensibilidad del modelo hipoplástico con la variación de ϕc . Ensayo triaxial
CD con arena de Karlsruhe, estado suelto (eo = 0,9)

ϕc q Variación de ϕc Variación de q
[o ] [kPa] [o ] [ %] [kPa] [ %]
20 504.85 10 33 451.6 47
25 706.70 5 17 249.8 26
28 849.95 2 7 106.5 11
30 955.48 0 0 0.0 0
32 1073.02 2 7 116.5 12
35 1268.92 5 17 312.4 33
40 1660.38 10 33 703.9 74

Comparación de ensayos triaxiales CD y CU con parámetros de la arena

de Karlsruhe - Estado suelto

La variación de la resistencia del suelo por las condiciones de drenaje en la fase de corte
en un ensayo triaxial con arena de Karlsruhe (estado suelto) se observa en la figura 4.17.
Estas gráficas se realizaron con los mismos parámetros del suelo y condiciones del ensayo,
variando únicamente la condición de drenaje en la fase de corte.

1200
1200

Ensayo triaxial - arena de Karlsruhe

Simulado - Elementtest Cudmani
eo = 0.9 σ3 = 500 kPa ϕc = 30o
Esfuerzo desviador q [kPa]
Esfuerzo desviador q [kPa]

ϕc = 30o
ϕc = 40o
800 800
CSL
Ensayo triaxial - arena de Karlsruhe
Simulado - Elementtest Cudmani
eo = 0.9 σ3 = 500 kPa ϕc = 30o

CD
CU

400 400

0 0
0 400 800 1200 0 0.04 0.08 0.12 0.16

Esfuerzo promedio efectivo p' [kPa] Deformación axial ε1 [-]

(a) (b)

Figura 4.17. Simulación de ensayos triaxiales CU y CD con arena de Karlsruhe, estado

suelto

85
CAPÍTULO 4. SENSIBILIDAD DEL MODELO HIPOPLÁSTICO CON ϕC MIC 2008-I-34

4.8. Simulación de ensayos triaxiales con parámetros de la

arena de Karlsruhe variando ϕc - Estado denso (eo = 0,6)

Ensayo triaxial CU con parámetros de la arena de Karlsruhe - Estado

denso

2500 2000

Ensayo triaxial CU - arena de Karlsruhe

Simulado - Elementtest Cudmani
eo = 0.6 σ3 = 500 kPa
2000 1600
ϕc = 20o

Esfuerzo desviador q [kPa]

Esfuerzo desviador q [kPa]

ϕc = 25o
ϕc = 28o
ϕc = 20o 1200
1500
ϕc = 30o
ϕc = 25o
ϕc = 28o ϕc = 32o
ϕc = 30o ϕc = 35o
1000 ϕc = 32o 800 ϕc = 40o
ϕc = 35o
ϕc = 40o

500 400
Ensayo triaxial CU - arena de Karlsruhe
Simulado - Elementtest Cudmani
eo = 0.6 σ3 = 500 kPa

0 0
0 500 1000 1500 2000 2500 0 0.01 0.02 0.03 0.04 0.05

Esfuerzo promedio efectivo p' [kPa] Deformación axial ε1 [-]

(a) (b)
500

Sensibilidad del modelo constitutivo hipoplástico

Simulado - ElementTest (Cudmani)
Arena de Karlsruhe - triaxial CU
Esfuerzo desviador q [kPa]

400 eo = 0.6 σ3 = 500 kPa

300

200

100
20 30 40

Ángulo de fricción crítico ϕc [ o ]

(c)

Figura 4.18. Simulación ensayo triaxial CU con arena de Karlsruhe variando ϕc - estado
denso

86
CAPÍTULO 4. SENSIBILIDAD DEL MODELO HIPOPLÁSTICO CON ϕC MIC 2008-I-34

Las simulaciones se realizaron para una muestra con una relación de vacı́os inicial igual
a 0.6. La muestra se comprimió isotrópicamente con una presión de cámara (σ3 ) igual a
500 kPa y luego se realizó el corte no drenado hasta una deformación axial del 12 %. Estas
simulaciones se realizaron variando ϕc con incremento de 2o , 5o y 10o . Los resultados de
estas simulaciones se observan en la Figura 4.18 y Tabla 4.7. La variación porcentual y
numérica de ϕc y q en la Tabla 4.7 se obtuvo comparando estos valores obtenidos en cada
simulación con los resultados de la simulación realizada con el valor del parámetro ϕc real
de la arena de Karlsruhe (ϕc = 30o ).

Tabla 4.7. Sensibilidad del modelo hipoplástico con la variación de ϕc . Ensayo triaxial
CU con arena de Karlsruhe, estado denso (eo = 0,6)

ϕc q Variación de ϕc Variación de q
[o ] [kPa] [o ] [ %] [kPa] [ %]
20 199.44 10 33 115.2 37
25 248.25 5 17 66.4 21
28 238.78 2 7 30.8 10
30 314.60 0 0 0.0 0
32 329.76 2 7 15.2 5
35 354.86 5 17 40.3 13
40 418.47 10 33 103.9 33

La sensibilidad del modelo hipoplástico con la variación del parámetro ϕc en un ensayo

triaxial CU con arena de Karlsruhe (estado denso) se observa en la Figura 4.18c. El valor
del esfuerzo desviador q en esta gráfica corresponde al menor valor de p0 en la fase de corte.

Ensayo triaxial CD con parámetros de la arena de Karlsruhe - Estado

denso

Las simulaciones se realizaron para una muestra con una relación de vacı́os inicial igual
a 0.6. La muestra se comprimió isotrópicamente con una presión de cámara (σ3 ) igual a
500 kPa y luego se realizó el corte drenado hasta una deformación axial del 12 %. Estas
simulaciones se realizaron variando ϕc con incremento de 2o , 5o y 10o . Los resultados de
estas simulaciones se observan en la Figura 4.19 y Tabla 4.8. La variación porcentual y
numérica de ϕc y q en la Tabla 4.8 se obtuvo comparando estos valores obtenidos en cada
simulación, con los resultados de la simulación realizada con el valor del parámetro ϕc real
de la arena de Karlsruhe (ϕc = 30o ).

87
CAPÍTULO 4. SENSIBILIDAD DEL MODELO HIPOPLÁSTICO CON ϕC MIC 2008-I-34

5000
5000

Ensayo triaxial CD - arena de Karlsruhe

Ensayo triaxial CD - arena de Karlsruhe Simulado - Elementtest Cudmani
Simulado - Elementtest Cudmani eo = 0.6 σ3 = 500 kPa
eo = 0.6 σ3 = 500 kPa 4000
4000

Esfuerzo desviador q [kPa]

Esfuerzo desviador q [kPa]

ϕc = 20o ϕc = 40o
ϕc = 25o
3000 3000
ϕc = 28o
ϕc = 30o
ϕc = 32o
ϕc = 35o
2000
ϕc = 35o 2000
ϕc = 40o
ϕc = 32o
ϕc = 30o
ϕc = 28o
1000 1000
ϕc = 25o
ϕc = 20o

0 0
0 1000 2000 3000 4000 5000 0 0.04 0.08 0.12 0.16

Esfuerzo promedio efectivo p' [kPa] Deformación axial ε1 [-]

(a) (b)
5000

Sensibilidad del modelo constitutivo hipoplástico

4000
Simulado - ElementTest (Cudmani)
Arena de Karlsruhe - triaxial CD
Esfuerzo desviador q [kPa]

eo = 0.6 σ3 = 500 kPa

3000

2000

1000

0
20 30 40

Ángulo de fricción crítico ϕc [ o ]

(c)

Figura 4.19. Simulación ensayo triaxial CD con arena de Karlsruhe variando ϕc - estado
denso

La sensibilidad del modelo hipoplástico con la variación del parámetro ϕc en un ensayo

triaxial CD con arena de Karlsruhe (estado denso) se observa en la Figura 4.19c. El valor
del esfuerzo desviador q en esta gráfica corresponde al mayor valor de q en la fase de corte.

88
CAPÍTULO 4. SENSIBILIDAD DEL MODELO HIPOPLÁSTICO CON ϕC MIC 2008-I-34

Tabla 4.8. Sensibilidad del modelo hipoplástico con la variación de ϕc . Ensayo triaxial
CD con arena de Karlsruhe, estado suelto (eo = 0,9)

ϕc q Variación de ϕc Variación de q
[o ] [kPa] [o ] [ %] [kPa] [ %]
20 676.82 10 33 834.2 55
25 1013.20 5 17 497.8 33
28 1286.22 2 7 224.8 15
30 1511.02 0 0 0.0 0
32 1781.27 2 7 270.3 18
35 2302.66 5 17 791.6 52
40 3671.24 10 33 2160.2 143

Comparación de ensayos triaxiales CD y CU con parámetros de la arena

de Karlsruhe - Estado denso

La variación de la resistencia del suelo por las condiciones de drenaje en la fase de corte
en un ensayo triaxial con arena de Karlsruhe (estado denso)se observa en la figura 4.20.
Estas gráficas se realizaron con los mismos parámetros del suelo y condiciones del ensayo,
variando únicamente la condición de drenaje en la fase de corte.

8000

6000 Ensayo triaxial - arena de Karlsruhe

Simulado - Elementtest Cudmani
eo = 0.6 σ3 = 500 kPa ϕc = 30o
Esfuerzo desviador q [kPa]

Esfuerzo desviador q [kPa]

6000
CD
CU

4000

4000
CU
CD
CSL

2000 Ensayo triaxial - arena de Karlsruhe

Simulado - Elementtest Cudmani 2000

eo = 0.6 σ3 = 500 kPa ϕc = 30o

0 0
0 2000 4000 6000 0 0.04 0.08 0.12 0.16

Esfuerzo promedio efectivo p' [kPa] Deformación axial ε1 [-]

(a) (b)

Figura 4.20. Simulación ensayos triaxiales CU y CD con arena de Karlsruhe - estado

denso

89
CAPÍTULO 4. SENSIBILIDAD DEL MODELO HIPOPLÁSTICO CON ϕC MIC 2008-I-34

4.9. Resumen de las simulaciones de sensibilidad del modelo

con ϕc

En la Tabla 4.9 se resume la variación porcentual del esfuerzo desviador q, obtenida con
la variación del parámetro ϕc en 2o . De esta Tabla resumen puede analizarse que tanto
en ensayos CU y CD en estado de densidad suelto, medio y denso, este parámetro influye
en la resistencia del suelo. Una variación del 10 % en la resistencia del suelo (capacidad
portante) conlleva a un aumento en 5 % del tamaño de la cimentación y del 10 % en los
costo.

Tabla 4.9. Resumen de sensibilidad del modelo hipoplástico con la variación de ϕc

Variación de q [ %]
Arena Densidad relativa Ensayo ϕc − 2o ϕ c + 2o
Guamo Media Triaxial CU 7 8
Guamo Media Triaxial CD 12 14
Guamo Suelto Triaxial CU 7 7
Guamo Suelto Triaxial CD 11 12
Karlsruhe Denso Triaxial CU 10 5
Karlsruhe Denso Triaxial CD 15 18
Karlsruhe Suelto Triaxial CU 7 7
Karlsruhe Suelto Triaxial CD 7 12

4.10. Simulaciones ensayo triaxial variando eo

La influencia de la relación de vacı́os inicial (eo ) en el comportamiento del suelo puede ana-
lizarse en la Figura 4.21. En esta Figura se observan las simulaciones de ensayos triaxiales
CU y CD con parámetros de la arena de Karlsruhe, σ3 = 500 kPa, 1 = 12 % y eo variable.

En las simulaciones de los ensayos triaxiales no drenados se observa que cuando eo > ec
0
la trayectoria en el plano q − p (Figura 4.21b) alcanza el estado crı́tico y continúa por
0
la lı́nea del estado crı́tico pero en dirección al origen (q y p decrecen) . Cuando eo < ec
(triaxial CU) la trayectoria en el plano q−1 no presenta falla (estado pico), por el contrario
continua aumentando el esfuerzo desviador (q). Éste ultimo comportamiento descrito es
caracterı́stico de arenas densas en ensayos no drenado como lo presenta Kolymbas [14], lo
curioso es que el modelo hipoplástico presenta este comportamiento aún en muestras con
densidad media.

90
CAPÍTULO 4. SENSIBILIDAD DEL MODELO HIPOPLÁSTICO CON ϕC MIC 2008-I-34

800
2000
Ensayo Triaxial CU arena de Karlsruhe
Ensayo triaxial CD arena de Karlsruhe Simulado - Elementtest Cudmani
Simulado - Elementtest Cudmani σ3 = 500 kPa ε1 = 12%
σ3 = 500 kPa ε1 = 12%
1600 eo = 0.9

Esfuerzo desviador q [kPa]

Esfuerzo desviador q [kPa]

600 eo = 0.8
eo = 0.7
eo = 0.6
1200 CSL

eo= 0.9
400
eo= 0.8
eo= 0.7
800 eo= 0.6
CSL

200
400

0 0
0 400 800 1200 1600 2000 0 200 400 600 800

Esfuerzo promedio efectivo p' [kPa] Esfuerzo promedio efectivo p' [kPa]

(a) (b)

1600 2000

eo = 0.9 Ensayo Triaxial CU

eo = 0.8 Arena de Karlsruhe
eo = 0.7 Simulado - Elementtest Cudmani
1600 σ3 = 500 kPa ε1 = 12%
eo = 0.6
Esfuerzo desviador q [kPa]
Esfuerzo desviador q [kPa]

1200

1200

800

800
Ensayo Triaxial CD arena de Karlsruhe
Simulado - Elementtest Cudmani
σ3 = 500 kPa ε1 = 12%
400
eo = 0.9 400
eo = 0.8
eo = 0.7
eo = 0.6

0 0
0 0.04 0.08 0.12 0 0.04 0.08 0.12

Deformación axial ε1 [-] Deformación axial ε1 [-]

(c) (d)

Figura 4.21. Simulación ensayo triaxial CU y CD con arena del Karlsruhe variando eo

4.11. Ensayos Triaxiales en arenas

En ensayos triaxiales drenados en arenas las muestras sueltas siguen una trayectoria esfuerzo-
deformación curva cóncava hacia abajo sin un valor máximo pronunciado y en muestras
densas esta trayectoria presenta un máximo a partir del cual la resistencia se reduce [15].
En las Figuras 4.22 y 4.23 Kolymbas [14] describe el comportamiento de muestras de arena

91
CAPÍTULO 4. SENSIBILIDAD DEL MODELO HIPOPLÁSTICO CON ϕC MIC 2008-I-34

en ensayos triaxiales no drenados considerando la misma presión de cámara pero la densi-

dad relativa de las muestras variable. La trayectoria (A) corresponde a una muestra suelta
(decrece el esfuerzo desviador q “softening”) Fig 4.23. La trayectoria (B) corresponde al
comportamiento de una muestra con densidad media (inicialmente presentará ablanda-
miento “softening” y luego endurecimiento ”hardening”) Fig 4.23. Una muestra densa (C)
no presentará ablandamiento (Fig 4.23). La relación de vacı́os de la muestra e y esfuerzo
de confinamiento inicial σ3 definen el estado de la muestra: densa o suelta (Figura 4.22a).
En la Figura 4.22b se puede analizar que en una muestra suelta (A) el exceso de presión
de poros que se genera durante la aplicación del esfuerzo desviador es mayor que para una
muestra densa (C) bajo las mismas condiciones del ensayo. Esto explica las trayectorias de
la Figura 4.23.

(a) (b)

Figura 4.22. Ensayos triaxiales no drenados en arenas [14]

Figura 4.23. Ensayos triaxiales no drenados en arenas [14]

92
Capı́tulo 5

Ensayos de Laboratorio

5.1. Ángulo de Reposo

Después de analizar la información recopilada sobre el ángulo de reposo en el Capı́tulo 1,

se seleccionaron algunos de los métodos plantados por los autores para determinar ϕc . Ésto
con el fin de analizar la influencia de los siguientes factores en el resultado del ángulo:

Método del ensayo

Variación de la cantidad de material (masa) en cada uno de los métodos

Velocidad de salida del material por el embudo

Tamaño del embudo

Tamaño y distribución de las partı́culas (método: embudo)

Forma de las partı́culas (método: embudo)

Humedad (método: cilindro “Santamarina”)

La rugosidad de la superficie sobre la cual se realiza el ensayo, es un factor que influye

en el ángulo de reposo. Ésto ya ha sido estudiado por algunos autores como Kalman et
al., 1993 [16] y Chik & Vallejo, 2005 [30]. La primera referencia concluye: “el ángulo de
reposo varı́a aproximadamente 5o dependiendo del coeficiente de fricción de la base”; y la
segunda: “la rugosidad de la base afecta más el ángulo de reposo para partı́culas gruesas.
Cuando el porcentaje de finos es mayor o igual al 70 % de la masa de la mezcla (predomina

93
CAPÍTULO 5. ENSAYOS DE LABORATORIO MIC 2008-I-34

fracción fina), la rugosidad de la base no influye”. Los resultados de estas investigaciones se

presentan en el Capı́tulo 1. Con base en éstos resultados se descartó esta variable (rugosidad
de la base), para analizar los otros factores mencionados anteriormente.

En este trabajo se realizaron los ensayos sobre una superficie de vidrio. Ésto para no generar
ningún tipo de fuerza debida a la fricción entre el materia granular y la base durante el
proceso de formación de la pila.

Como se dijo anteriormente, en esta investigación no se realizó el análisis de la rugosidad de

la base, pero se construyeron pilas (método: embudo) con arena del Guamo sobre 2 tipos de
superficies: una superficie rugosa (papel lija No 180) y una superficie lisa (vidrio). De ésto
se encontró que la forma de los taludes de las pilas son diferentes. Para la superficie rugosa
se obtiene una pila convexa como lo describe Dodds, 2003 [7] en su tesis (Figura 5.1a).
Para la superficie de vidrio la forma de la pila es como la describen Grasselli & Herrmann,
1997 [22] en sus artı́culos: “el talud está compuesto por una parte redondeada en la parte
superior de la pila, una parte lineal en el centro donde generalmente se determina el ángulo
de reposo y una parte curva en la base del talud (descrita por una función logarı́tmica) ”
(Figura 5.1b).

(a) base rugosa: lija (b) base lisa: vidrio

Figura 5.1. Formas de las pilas dependiendo de la rugosidad de base

Durante el proceso de formación de estas pilas (base: vidrio y papel lija) se observó que
cuando la superficie es rugosa (lija) los granos ruedan por la superficie del talud. Pero
cuando se construye sobre la superficie lisa (vidrio) se presentan dos fenómenos, primero
las partı́culas empiezan a rodar sobre la superficie del talud y posteriormente (posiblemente
debido al peso y a la baja o nula rugosidad del vidrio) las partı́culas que están en la base
de la pila se deslizan sobre el vidrio. Posiblemente debido a estos dos fenómenos se debe
la diferencia de la forma de las pilas construidas sobre bases con diferente rugosidad.

El material empleado en estos ensayos fue arena del Guamo lavada sobre tamiz #200 (sin
material fino). Los ensayos se realizaron a temperatura y humedad controlada.

94
CAPÍTULO 5. ENSAYOS DE LABORATORIO MIC 2008-I-34

5.1.1. Métodos seleccionados para determinar el ángulo de reposo

[Link]. Método: embudo

El embudo se apoya sobre una superficie plana y luego el material granular completamente
seco se vierte dentro del embudo. Se debe garantizar que el diámetro de salida del embudo
sea mayor o igual que cinco veces el diámetro de la partı́cula (Dembudo ≥ 5dparticula ) de
acuerdo con Mehmet, 2006 [10]. El embudo se levanta lentamente para ir construyendo la
pila por pluviación del material, pero siempre debe estar en contacto con la parte superior
(cima) de la pila (Figura 5.2). Ésto con el fin de evitar que el impacto de las partı́culas
afecte la estructura interna de la pila y el material quede en el estado más suelto posible. Se
debe garantizar durante el ensayo que el embudo esté completamente vertical y el material
caiga en el mismo punto, para que la forma de la pila quede perfectamente cónica.

Este método se realiza de 2 formas: 1. se construye la pila levantando lentamente el embudo,

se toma una foto y posteriormente se mide el ángulo de reposo en la zona lineal del talud
(parte media del talud). 2. una vez construida la pila (de la forma anterior) se hace una
pequeña excavación en la pata del talud como lo recomiendan Herle & Gudehus, 1999 [11].
Esta excavación hace que el ángulo de los taludes de la pila disminuyan.

Figura 5.2. Método Embudo

Para este método se emplearon tres embudos de plástico de diferentes tamaños:

1. Embudo pequeño: diámetro interno del orificio de salida 11.4 cm.

2. Embudo mediano: diámetro interno de salida 14.7 cm.

3. Embudo grande: diámetro interno de salida 19.0 cm.

95
CAPÍTULO 5. ENSAYOS DE LABORATORIO MIC 2008-I-34

Variación de la altura del material dentro del embudo

Para realizar el ensayo por éste método primero fue necesario determinar si la altura del
material dentro del embudo influye en los resultados. Para ésto se realizó un análisis del
flujo del material de dos formas:

1. Altura variable del material dentro del embudo (hvariab ): se llenó el embudo con la
cantidad de material deseada y se dejó fluir el material completamente (sin depositar
más material dentro del embudo durante el ensayo) midiendo el tiempo de salida de
todo el material y la masa de material que salió.

2. Altura constante del materia dentro del embudo (hconst ): se dejaron fluir las partı́culas
manteniendo la altura del material dentro del embudo constante y se midió de igual
forma el tiempo y la masa.

Este análisis se realizó con partı́culas de tamaños variables entre 1.18 mm y 0.60 mm
(material pasa tamiz #8, retiene tamiz #16). En la Tabla 5.1 se presentan los resultados
de este análisis. De estos resultados se concluye que el material fluye más rápido cuando se
deja salir todo (hvariab ), que cuando se mantiene constante la altura del material (hconst ).
La medición del tiempo se realizó observando y cronometrando la salida del material, la
precisión al medir el tiempo no es la mejor y ésto influye en los resultados. Esta influencia
es más significativa en los embudos de mayor diámetro, debido a que por el tamaño del
embudo salen más partı́culas por segundo y el error en la medición del tiempo influye
significativamente en los resultados.

Tabla 5.1. Flujo de las partı́culas dependiendo del embudo y altura del material dentro
del embudo
Embudo Tamaño partı́culas Flujo [gr/s] Variación
Pasa Tamiz # Retiene Tamiz # hvariab hconst [ %]
Grande 126 114 9.8
Mediano 16 30 98 89 8.5
Pequeño 53 51 3.9

Este análisis de flujo se hizo con el fin de determinar si la altura del material dentro del
embudo influye significativamente en el ángulo de reposo. La variación del flujo fue menor
al 10 % y estos ensayos se realizaron dejando caer el material libremente (Figura 5.3). En los
ensayos para determinar el ángulo de reposo las pilas se construyen con el embudo siempre
en contacto con la parte superior de la pila (Figura 5.2), ésto hace que la influencia se

96
CAPÍTULO 5. ENSAYOS DE LABORATORIO MIC 2008-I-34

todavı́a menor. Concluyendo, la influencia de la variación de la altura del material dentro

del embudo en el flujo de las partı́culas se puede despreciar.

Figura 5.3. Flujo

Variación de la Masa dentro del embudo

El segundo análisis realizado antes de medir el ángulo de reposos por este método (embudo),
fue determinar si la cantidad de material dentro del embudo influye en el flujo de las
partı́culas. Ésto con el fin de analizar si el peso de todo el material sobre las partı́culas
saliendo por el orificio del embudo influye.

Este análisis se realizó con altura variable del material dentro del embudo (hvariab ). Los
embudos y masas empleados fueron:

embudo grande, masas de 250, 500 y 2900 gr.

embudo mediano, masas de 250, 500 y 750 gr.

El flujo promedio obtenido variando las masas fue de 126 gr/s con el embudo grande y
de 97.7 gr/s con el embudo mediano (valor obtenido del promedio de 3 valores por cada
masa). La variación del flujo promedio con los valores del flujo medido es menor al 1.5 %.
De acuerdo con estos resultados se concluye que la cantidad del material dentro del embudo
no influye en el flujo de las partı́culas.

[Link]. Método: Cilindro “Santamarina”

En este trabajo para determinar el ángulo de reposo por este método, se empleó un cilindro
de vidrio de 9.17 cm de diámetro y se siguió la metodologı́a planteada por Santamarina &

97
CAPÍTULO 5. ENSAYOS DE LABORATORIO MIC 2008-I-34

Cho [25], descrita en el Capı́tulo 1.

Dentro del cilindro se deposita el material y luego se inclina a un ángulo mayor de 60o
(en este trabajó el cilindro se inclinó 80o aproximadamente). Luego se devuelve lentamente
evitando golpes, para que el material se deposite lo más suelto posible. El ángulo de reposo
(ϕrep ) se mide con un transportador en la parte media del talud, debido a las curvaturas
que se presentan en la parte superior e inferior de la sección (Figura 5.4). El resultado de
ϕc debe ser el promedio de 5 mediciones.

Figura 5.4. Método: cilindro - Santamarina

[Link]. Método: Cilindro “Train”

La metodologı́a seguida en este trabajo para determinar el ángulo de reposo por el método:
cilindro “Train”, fue tomada de los métodos descritos por Train, 1858 [28], presentados en
el capı́tulo 1.

Se deposita el material dentro del cilindro (9.17 cm de diámetro), se tapa y luego se rota
horizontalmente sobre una superficie plana. Cuando el cilindro esta rotando lentamente se
observa que el material toma primero un ángulo de reposo, luego se presenta una avalancha
o “falla” en el talud y este ángulo decrece. Es decir, al rotar lentamente el cilindro el talud
que se forma al inicio no cambia dentro del cilindro, pero con respecto a la horizontal este
ángulo se va incrementando. Luego, si se sigue rotando el cilindro llega un punto en el
que el talud formado dentro del cilindro “falla” y se forma un nuevo talud con un ángulo
menor. En este trabajo se tomó el primer ángulo que se forma tan pronto se gira el cilindro
o el que se presenta tan pronto “falla” el talud (Figura 5.5). No se tomó el mayor ángulo

98
CAPÍTULO 5. ENSAYOS DE LABORATORIO MIC 2008-I-34

como recomienda Train, porque es dificil saber el momento justo antes de presentarse la
avalancha o “falla”, por lo tanto no hay certeza en que el ángulo que se mide sea el mayor.
El ángulo de reposo por este método se midió directamente con un transportador y fue el
promedio de 5 mediciones.

Figura 5.5. Método: cilindro “Train”

[Link]. Método: Cilindro abierto

En este trabajo se realizó este ensayo con un cilindro metálico (diámetro = 10.15 cm y
altura = 11.59 cm) abierto en los dos extremos. El método seguido fue el siguiente:
El material se deposita dentro del cilindro hasta llenarlo completamente. El llenado debe
ser lento y sin vibración o compactación para que el material quede lo más suelto posible.
Luego se levanta lentamente el cilindro y la pila se va formando a medida que el material va
saliendo (Figura 5.6). El ángulo de reposo se determina midiendo el ángulo de inclinación
en la parte media del talud. Al igual que en los otros métodos, ϕrep debe ser el promedio
de 5 mediciones.

99
CAPÍTULO 5. ENSAYOS DE LABORATORIO MIC 2008-I-34

(a)

(b)

Figura 5.6. Método: cilindro abierto

5.1.2. Influencia del método de ensayo y la cantidad de material (masa)

en el ángulo de reposo

[Link]. Método: embudo

El ángulo de reposo determinado por este método fue medido sobre la misma pila de tres
formas diferentes:

1. Midiendo diámetro (D) y altura (H) de la pila y suponiendo la forma de ésta per-
fectamente cónica: ϕc = tan−1 (2H/D)

2. Tomando fotos (desde diferentes ángulos) de la pila construida y por medio de auto-
cad se mide el ángulo de reposo en la parte media del talud (parte lineal).

3. Igual que el numeral anterior pero a la pila se le realiza una pequeña excavación en
la pata del talud, siguiendo la recomendación de Herle & Gudehus [11].

100
CAPÍTULO 5. ENSAYOS DE LABORATORIO MIC 2008-I-34

Tabla 5.2. Influencia de la masa en el ángulo de reposo, método: embudo - Material pasa
tamiz #16, retiene tamiz #30

Método Masa ϕrep promedio [o ]

[gr] diámetro y altura sin excavación con excavación
250 29.1 37.4 37.0
500 27.6 37.0 36.6
Embudo 750 24.9 36.1 35.6
Grande 1000 25.9 35.8 35.5
1500 24.6 35.7 35.3
2000 26.0 35.5 34.8

En la Figura 5.7 se presentan los resultados de éste análisis y se observa la tendencia

decreciente de ϕc con el aumento de la masa. ϕc calculado midiendo el diámetro y altura
de la pila decrece 4.5o , medido de la pila sin excavar decrece 1.9o y medido de la pila con
excavación decrece 2.3o .
40

Variación de ϕrep con la masa

Método: embudo
36
Ángulo de reposo ϕrep [o]

32

28

24
Pila sin excavación
Pila con excavación
Diametro y Altura

20
0 500 1000 1500 2000
Masa [gr ]

Figura 5.7. Influencia de la masa en el ángulo de reposo, método: embudo - Material pasa
tamiz #16, retiene tamiz #30

[Link]. Método: cilindro “Santamarina”

Estos ensayos se realizaron siguiendo el método de Santamarina & Cho [25] pero con el
material seco. Los ensayos se realizaron con el material granular seco y el cilindro de vidrio
con el fin de tener las mismas condiciones en los ensayos realizados. En la Tabla 5.3 y
Figura 5.8 se presentan los resultados de este análisis, y se observa una tendencia creciente

101
CAPÍTULO 5. ENSAYOS DE LABORATORIO MIC 2008-I-34

del ángulo de reposo al aumentar la masa (variación de 3.2o ). Posiblemente para este
método sea necesario determinar alguna relación entre la cantidad de material, tamaño de
partı́culas y diámetro del cilindro.

Tabla 5.3. Influencia de la masa en el ángulo de reposo, método: cilindro “Santamarina”

(material seco) - Material pasa tamiz #16, retiene tamiz #30

Método Masa [gr] ϕrep promedio [o ]

250 34.1
500 34.9
Santamarina 750 35.0
(material seco) 1000 36.1
1500 36.2
2000 37.3

40

Variación de ϕrep con la masa

Método: Santamarina (material seco)

38
Ángulo de reposo ϕrep [o]

36

34

32
0 500 1000 1500 2000
Masa [gr ]

Figura 5.8. Influencia de la masa en el ángulo de reposo, método: cilindro “Santamarina”

(material seco) - Material pasa tamiz #16, retiene tamiz #30

[Link]. Método: cilindro “Train”

Los resultados de la variación de ϕc con la masa, por este método se presentan en la

Tabla 5.4 y Figura 5.9. En este método no se observa ninguna influencia de la masa en los
resultados.

102
CAPÍTULO 5. ENSAYOS DE LABORATORIO MIC 2008-I-34

Tabla 5.4. Influencia de la masa en el ángulo de reposo, método: Train - Material pasa
tamiz #16, retiene tamiz #30

Método Masa [gr] ϕrep promedio [o ]

500 34.8
Train 1000 34.6
1500 34.5
2000 34.6

40

Variación de ϕrep con la Masa

Método: Train
38
Ángulo de reposo ϕrep [o]

36

34

32

30
0 500 1000 1500 2000
Masa [gr ]

Figura 5.9. Influencia de la masa en el ángulo de reposo, método: Train - Material pasa
tamiz #16, retiene tamiz #30

[Link]. Método: cilindro abierto

De acuerdo con la metodologı́a descrita por Pilpel, 1966 [23] este ensayo se realiza llenando
completamente el cilindro, por lo tanto, en este método no se varió la cantidad de material
(1200 gr). El material empleado fue el mismo que en los métodos anteriores (arena del
Guamo pasa tamiz #16, retiene tamiz #30), con el fin de comparar resultados. El ángulo
de reposo promedio determinado por este método fue 33.9o . En los ensayos realizados por
este método se observó que ϕc es función de la velocidad con que se levante el cilindro, si
se levanta muy rápido ϕc es menor.

103
CAPÍTULO 5. ENSAYOS DE LABORATORIO MIC 2008-I-34

[Link]. Cuadro Resumen - influencia del método y la cantidad del material

en el ángulo de reposo

Para responder esta pregunta se determinó el ángulo de reposo por los diferentes métodos
variando la masa desde 250 gr hasta 2000 gr. El material empleado en todos los ensayos fue
arena del Guamo de tamaño variable entre 1.18 mm y 0.60 mm (pasa tamiz #16, retiene
tamiz #30), material completamente seco. En la Tabla 5.5 y Figura 5.10 se resumen los
resultados obtenidos en cada método.

Tabla 5.5. Influencia del método y la masa en el ángulo de reposo

Método - ϕrep promedio

Masa [gr] Embudo Embudo Embudo Cilindro Cilindro Cilindro
diámetro y altura sin excav. con excav. Santamarina Train abierto
250 29.1 37.4 37.0 34.1 -
500 27.6 37.0 36.6 34.9 34.8
750 24.9 36.1 35.6 35.0 -
1000 25.9 35.8 35.5 36.1 34.6 33.9
1500 24.6 35.7 35.3 36.5 34.5
2000 26.0 35.5 34.8 37.3 34.6

40

Variación de ϕrep con la masa y el método

36
Ángulo de reposo ϕrep [o]

32
Pila sin excavación
Pila con excavación
Santamarina
Train
Diametro - altura

28

24
0 500 1000 1500 2000
Masa [gr ]

Figura 5.10. Influencia del método y la masa en el ángulo de reposo

El ángulo de reposo calculado midiendo diámetro y altura de la pila es una aproximación,

104
CAPÍTULO 5. ENSAYOS DE LABORATORIO MIC 2008-I-34

porque en este cálculo no se tiene en cuenta la forma aplanada de la parte superior de la

pila, ni la forma de la pata del talud (forma logarı́tmica [21]). Los resultados obtenidos por
este método son los que presentan mayor diferencia entre ellos y con respecto a los otros
métodos.

De los resultados obtenidos por medio de la pila (método: embudo), medido de las 3
formas, se puede concluir que a mayor masa, menor ángulo de reposo. La excavación que
recomiendan Herle & Gudehus [11] en la pata del talud de la pila, hace que el ángulo de
reposo disminuya menos de 1o con respecto al talud inicial.

El método de Santamarina es el único método en el que la tendencia de ϕrep es creciente,

a mayor masa mayor ángulo de reposo. Para este método deberı́a existir una relación entre
la cantidad de material, tamaño de las partı́culas y diámetro del cilindro.

En el método que menos influencia tiene la masa en el ángulo de reposo es en el cilindro

”Train”(0.25 grados) y en el método que más influye es midiendo altura y diámetro de la
pila (4.5o ). Claro está que como se dijo anteriormente este último método es una aproxi-
mación de ϕrep . Por el método del embudo (pila sin excavación) ϕc decrece 1.9o y por el
método: cilindro “Santamarina” ϕc aumenta 3.2o .

[Link]. Relación lineal entre el ángulo de reposo y la redondez de la partı́cula

(Santamarina & Cho [26])

El ángulo de reposo de la arena del Guamo se determinó con la relación lineal entre ϕrep
y la redondez de la partı́cula, propuesta por Santamarina & Cho:

ϕrep = 42 − 17R (5.1)

En la Tabla 5.6 se presentan los ángulo de reposo calculados a partir de la redondez

determinada para las 9 partı́culas de la arena del Guamo (Capı́tulo 4 - Tabla 3.5)

105
CAPÍTULO 5. ENSAYOS DE LABORATORIO MIC 2008-I-34

Tabla 5.6. Ángulo de reposo de la arena del Guamo determinado a partir de la relación
lineal entre ϕrep y la redondez de la partı́cula

Partı́cula # Redondez Clasificación ϕrep

1 0.55 Redondeado 32.7
2 0.77 Bien Redondeado 28.9
3 0.52 Redondeado 33.2
4 0.74 Bien Redondeado 29.4
5 0.53 Redondeado 33.0
7 0.42 Sub Redondeado 34.9
8 0.42 Sub Redondeado 34.9
9 0.43 Sub Redondeado 34.7
10 0.56 Redondeado 32.5

De este análisis para la arena del Guamo se obtuvo que las partı́culas sub redondeadas
tienen ϕrep promedio de 34.8o ; las redondeadas 32.9o y las bien redondeadas 29.2o . Es difı́cil
comparar estos resultados con los ángulos de reposos obtenidos para la arena del Guamo
por los otros métodos, porque el material con el que se realizaron los otros análisis no fue
clasificado por redondez de la partı́cula. por lo tanto, el ángulo de reposo determinado por
los otros métodos es un promedio de toda la masa de material (partı́culas de diferentes
formas).

5.1.3. Influencia de la velocidad del desplazamiento del embudo en el

ángulo de reposo

Para este análisis se construyeron pilas con 1000 gr de arena del Guamo (material com-
pletamente seco, pasa tamiz #16, retiene tamiz #30) levantando el embudo a diferentes
velocidades (embudo grande). El embudo siempre debe estar en contacto con la parte su-
perior de la pila para evitar que el impacto de las partı́culas influya en la estructura interna
de la pila y el material quede lo más suelto posible. Las velocidades a las que se realizaron
los ensayos fueron:

1. “Rápido”: formando la pila de 1000 gr de material en 20 segundos (50 gr/s).

2. “Lento”: formando la pila de 1000 gr de material en 60 segundos (16.7 gr/s).

106
CAPÍTULO 5. ENSAYOS DE LABORATORIO MIC 2008-I-34

Tabla 5.7. Variación de la velocidad con que se levanta el embudo - Material pasa T16
retiene T30
Método Velocidad ϕrep promedio [o ]
Sin excavación Con excavación
Embudo Rápido 35.8 35.5
(grande) Lento 37.9 37.6

A partir de estos resultados se puede concluir que a mayor velocidad con que se levante el
embudo menor ángulo de reposo. La variación en este análisis fue de 2,1o .

5.1.4. Influencia del diámetro de salida del embudo en el ángulo de re-

poso

El análisis de esta variable se realizó empleando 3 embudos de tamaños diferentes (diámetro

de salida: 11.4, 14.7 y 19 cm) y una masa de 1000 gr. Se seleccionó esta cantidad de material
porque en los resultados obtenidos variando la masa se observó que después de 1000 gr la
variación de ϕrep es pequeña (por el método del embudo aproximadamente 0.5o ). En la
construcción de las pilas el embudo se levantó muy lento y se mantuvo siempre en contacto
con la parte superior de la pila para evitar el impacto de las partı́culas.

El la Tabla 5.8 y Figura 5.11 se presentan los resultados de este análisis. De estos resultados
se puede concluir que la influencia del diámetro de salida del embudo en ϕrep es mı́nima.
La variación observada en este análisis fue de 0.5o , siendo mayor el ángulo de reposo
determinado con el embudo de menor diámetro.

Tabla 5.8. Influencia del diámetro de salida del embudo en ϕrep - Material pasa tamiz
#16, retiene tamiz #30

Embudo Masa [gr] ϕrep promedio [o ]

Grande 1000 35.8
Mediano 1000 36.0
Pequeño 1000 36.3

107
CAPÍTULO 5. ENSAYOS DE LABORATORIO MIC 2008-I-34

37

Variación de ϕrep con el tamaño del embudo

36.5

Ángulo de reposo ϕrep [o] 36

35.5
Embudo grande
Embudo mediano
Embudo pequeño

35
0 400 800 1200 1600 2000
Masa [gr ]

Figura 5.11. Influencia de diámetro de salida del embudo en ϕrep

5.1.5. Influencia del tamaño y distribución de las partı́culas en el ángulo

de reposo

El análisis de esta variable se realizó construyendo una pila (embudo grande) con 1000 gr
de material granular. Las resultados del ángulo de reposo que se presentan en este análisis
fueron mediciones realizadas en las pila sin excavación en la pata del talud. El ϕrep se
determinó con la misma arena pero con 7 gradaciones diferentes. En la Tabla 5.10 y Figura
5.12 se presenta el ángulo de reposo determinado para cada tipo de material.

El material 1 y 2 que se presentan en la Tabla 5.10 están conformados por la siguiente

granulometrı́a:

Tabla 5.9. Descripción de los materiales 1 y 2

Material Cantidad Tamaño

[gr] pasa tamiz # retiene tamiz #
250 8 16
1 250 16 30
250 30 50
250 50 100
2 500 8 16
500 50 100

108
CAPÍTULO 5. ENSAYOS DE LABORATORIO MIC 2008-I-34

Tabla 5.10. Variación de ϕrep con el tamaño y distribución de las partı́culas

Método Material ϕrep promedio

pasa Tamiz retiene Tamiz [o ]
4 8 34.2
8 16 36.1
Embudo 16 30 35.8
(Grande) 30 50 35.5
50 100 34.8
Material 1 35.2
Material 2 35.2

39

Variación de ϕrep con el tamaño de las partículas

38
Ángulo de reposo ϕrep [o]

37

36

35

34

33
0 20 40 60 80 100 120
Material retenido en Tamiz #

Figura 5.12. Variación de ϕrep con el tamaño y distribución de las partı́culas

En estos resultados se observa que ϕrep decrece a medida que el tamaño de las partı́culas
disminuye. Estos resultados coinciden con el análisis que Herle & Gudehus, 1999 presentan
en su investigación: “Al rodar las partı́culas en el proceso de formación de la pila, los granos
grandes tienden a superar esfuerzos más altos. Esto puede explicar un incremento de ϕrep
al aumentar el tamaños de los granos”.

El ϕrep determinado para el material pasa tamiz #4, retiene tamiz #8 no se ajusta a la
tendencia de los otros datos. Es necesario aclarar que este material contiene partı́culas
grandes que no cumplen la condición del diámetro mı́nimo del embudo (demb ≥ 5dpart ) y
durante los ensayos el flujo de este material no fue constante (se tapó el embudo). Con
base en ésto se podrı́a decir que indirectamente se comprobó la relación entre el diámetro
del embudo y el diámetro de la partı́cula de mayor tamaño sugerido por Mehemet, 2006
[10].

109
CAPÍTULO 5. ENSAYOS DE LABORATORIO MIC 2008-I-34

El ángulo de reposo obtenido para los materiales 1 y 2 (mezcla de diferentes tamaños)

dio el mismo, posiblemente ésto se debe a que las 2 mezclas contienen el mismo tamaño
máximo y mı́nimo de las partı́culas. Al contener el mismo tamaño máximo y mı́nimo de
las partı́culas en la mezcla se presenta el mismo relleno de los vacı́os por las partı́culas de
menor tamaño.

En conclusión ϕrep varı́a 1.5o con el tamaño de las partı́culas del material.

5.1.6. Variación del ángulo de reposo con la forma de las partı́culas

El análisis de la influencia de la forma de las partı́culas en el ángulo de reposo se realizó con

3 tipos de materiales granulares diferentes: redondeado (arena del Guamo), alargado (arroz)
y angular (azucar) (Figura 5.13).

El ángulo de reposo para este análisis se determinó por el método del embudo (sin exca-
vación en la pata del talud), con el embudo grande y 1000 gr de material. En la Tabla
5.11 se presentan los resultados de este análisis. Las partı́culas de arena (pasa tamiz #8,
retiene #16) y las partı́culas de azúcar son aproximadamente del mismo tamaño, mientras
que los granos de arroz son de mayor tamaño. Como se analizó anteriormente el tamaño
de las partı́culas también influye en el ángulo de reposo (a mayor tamaño de las partı́culas,
mayor ϕrep ).

Realizar una comparación entre estos resultados es difı́cil debido al fenómenos observado
durante los ensayos realizados con azúcar. Durante los ensayos se realizaron las siguientes
observaciones:

El azúcar es un material susceptible a la humedad del ambiente, por este motivo

determinar el ϕrep de este material fue complicado. Durante la construcción de las 5
pilas para determinar ϕrep , se observó que en las primeras pilas el flujo del material
es continuo. Después, el azúcar de humedece y se adhieren unas partı́culas con otras,
ésto hace que el flujo no sea continuo y la pila quede de forma irregular.

Las pilas construidas con arroz y azúcar son más susceptibles a la excavación en la
pata del talud que las pilas construidas con arena, es decir disminuyen en mayor
proporción el ϕrep con la excavación realizada.

Al comparar el ϕrep del arroz con cualquiera de los otros dos materiales se observa que éste
es mayor 0.5o . Posiblemente esta pequeña diferencia se deba al tamaño de los granos y a
la forma.

110
CAPÍTULO 5. ENSAYOS DE LABORATORIO MIC 2008-I-34

(a) Partı́culas redondeadas: arena del Guamo (b) Partı́culas elongadas: arroz

(c) Partı́culas angulares: azucar

Figura 5.13. Imágenes tomadas con estereoscopio con 10 aumentos. Escala: 1/10 mm

Tabla 5.11. Influencia de la forma de las partı́culas

Material ϕrep promedio [o ]

Azúcar 35.5
Arena (pasa T16 retiene T30) 35.8
Arroz 36.3

111
CAPÍTULO 5. ENSAYOS DE LABORATORIO MIC 2008-I-34

5.1.7. Influencia de la saturación del material en el ángulo de reposo

Este análisis se realizó con el fin de determinar si el ángulo de reposo varia con la satura-
ción del material. Siguiendo la metodologı́a descrita por Santamarina (con material seco y
sumergido) y con material pasa tamiz # 16, retiene tamiz #30 se realizaron lo ensayos del
ángulo de reposo. En la Tabla 5.12 y Figura 5.14 se presentan estos resultados.

Tabla 5.12. Variación de ϕrep con la saturación del material

Masa ϕrep promedio [o ] - Método Santamarina

[gr] Material seco Material sumergido
500 34.9 33.8
750 35.0 34.2
1000 36.1 36.0
1500 36.5 36.4

40

38
Ángulo de reposo ϕrep [o]

36

34

32
Material seco
Material sumergido

30
0 400 800 1200 1600
Masa [gr ]

Figura 5.14. Material seco y sumergido

Comparando estos resultados se puede concluir que el ángulo de reposo no varı́a con la
saturación del material.

5.2. Ángulo de Fricción Crı́tico

El ángulo de fricción crı́tico ϕc se determinó mediante un ensayo triaxial CU. El ensayo

triaxial se realizó con el mismo material empleado en los ensayos del ángulo de reposo
(arena de Guamo pasa tamiz #16, retiene #30) con el fin de comparar estos dos ángulos
(ϕc y ϕrep ).

112
CAPÍTULO 5. ENSAYOS DE LABORATORIO MIC 2008-I-34

(a) (b)

(c) (d)

(e) (f)

(g) (h)

Figura 5.15. Ensayo triaxial arena del Guamo

113
CAPÍTULO 5. ENSAYOS DE LABORATORIO MIC 2008-I-34

La relación de vacı́os inicial de las muestras (eo = 0,6), los parámetros de la arena del
Guamo determinados por Pablo Arias [2] y la definición de densidad relativa, clasifican las
muestras ensayadas como material granular denso (Dr > 2/3):

emax − e
Dr =
emax − emin
1,15 − 0,6
Dr =
1,15 − 0,52
Dr = 0,87

Las muestras se fallaron hasta una deformación axial del 16 % para llevarlas todas hasta
el estado crı́tico (“el estado crı́tico se alcanza al 10 % de deformación” [3]). En la Figura
5.15a, b, c y d se presenta de forma resumida el proceso constructivo de la muestra; en
la Figura 5.15(f) la cámara triaxial durante el proceso de saturación y compresión; en la
Figura 5.15(g) la cámara triaxial en la prensa para fallar la muestra, y en la Figura 5.15(h)
un especimen fallado a un alto porcentaje de deformación (1 ≈ 25 %).

En las Figuras 5.16, 5.17 y 5.18 se presentan los resultados de este ensayo.
2500

2000

1500
q [kPa]

1000
Ensayo Triaxial CU
Arena del Guamo
eo = 0.60
σ3 = 75 kPa
500
σ3 = 150 kPa
σ3 = 300 kPa

0
0 0.05 0.1 0.15 0.2
ε1 [-]

Figura 5.16. Ensayo triaxial CU arena del Guamo, plano q − 1

114
CAPÍTULO 5. ENSAYOS DE LABORATORIO MIC 2008-I-34

250

Ensayo Triaxial CU
Arena del Guamo
eo = 0.60
σ3 = 75 kPa
σ3 = 150 kPa
0
σ3 = 300 kPa
Δu [kPa]

-250

-500
0 0.05 0.1 0.15 0.2
ε1 [-]

Figura 5.17. Ensayo triaxial CU arena del Guamo, plano ∆u − 1

2500

2000

M = 1.64 ϕc = 40o
1500
q [kPa]

1000
Ensayo Triaxial CU
Arena del Guamo
eo = 0.60
σ3 = 75 kPa
500
σ3 = 150 kPa
σ3 = 300 kPa

0
0 500 1000 1500 2000 2500
p' [kPa ]

0
Figura 5.18. Ensayo triaxial CU arena del Guamo, plano q − p

El ángulo de fricción crı́tico de la arena del Guamo determinado a partir de la pendiente

115
CAPÍTULO 5. ENSAYOS DE LABORATORIO MIC 2008-I-34

0
M en el plano q − p fue de 40o

5.3. Celda Hele Shaw

Teniendo en cuenta los ensayos realizados en esta investigación por los diferentes métodos
y la literatura respecto al ángulo de reposo, se decidió diseñar y construir Celdas Hele
Shaw para la determinación del ángulo de reposo de material granular.

El diseño de las celdas Hele Shaw se realizó con base en la descripción de la celda hecha
por Grasselli and Herrmann (1997) [9]. En este artı́culo analizan el efecto de las paredes
de la celda en el ángulo de reposo con partı́culas de diámetro promedio igual a 250µm.
Los Autores variaron el espesor de las paredes desde 1mm hasta 10mm, cada 1mm. En
este análisis obtuvieron que después de 5mm (20φpart ) no varı́a mucho el ángulo de reposo
y después de 10mm (40φpart ) no varı́a nada el ángulo de reposo. Las medidas de la celda
Hele Shaw vertical son 20 x 30 cm y el espesor de las placas depende del tamaño de las
partı́culas. Para que las paredes rı́gidas no juegan un papel importante en el ángulo de
reposo el espesor de la celda debe ser mayor a 20 veces el diámetro máximo de la partı́cula.
Este espesor hace que la movilidad local del grano sea esencialmente tridimensional [21],
aunque el sistema es bidimensional porque no hay grado de libertad en la forma de la
superficie en la dirección perpendicular a la pared.

Debido a la clasificación por tamaños de los materiales granulares, un arena contiene

partı́culas que varı́an desde 4.74 mm hasta 0.08 mm (pasa tamiz #4, retiene tamiz #200),
por lo tanto para cumplir con el espesor mı́nimo de la celda (20 veces el diámetro máximo
de la partı́cula) se diseñaron 2 celdas con espesores diferentes:

Tabla 5.13. Espesor mı́nimo de la celda Hele Shaw (20φpart ) teniendo en cuenta el diáme-
tro máximo de la partı́cula

Material φpart Mayor - φpart Menor 20 φpart Mayor

Pasa Tamiz # Retiene Tamiz # [mm] [mm]
4 8 4.74 - 2.36 94.8
8 16 2.35 - 1.18 47.0
16 30 1.17 - 0.60 23.4
30 50 0.59 - 0.30 11.8
50 100 0.29 - 0.15 5.8
100 200 0.14 - 0.08 2.8

En la Figura 5.19 se presentan las celdas diseñada para determinar el ángulo de reposo de
materiales granulares. La celda de la Figura 5.19(a) se diseñó para material granular con

116
CAPÍTULO 5. ENSAYOS DE LABORATORIO MIC 2008-I-34

partı́culas de diámetro mayor o igual que 2.35mm, es decir material para tamiz #8. La
celda de la Figura 5.19(b) sirve para determinar el ángulo de reposo de partı́culas menores
a 0.59mm (pasa tamiz #30). Estos espesores cumplen con el requerimiento de espesor de
la celda ≥ 20φpart .

(a) (b)

Figura 5.19. Celda Hele Shaw - medidas en cm.

Una de las variables analizadas en este trabajo es la influencia de la saturación del material
en el ángulo de reposo. Debido a ésto se diseñó una celda (Fig 5.20) en la cual se puede
construir la pila de material granular vertiéndolo por pluviación dentro de la celda llena
de agua. En esta celda se puede determinar el ángulo de reposo por el método descrito por
Atkinson, 1993 [3] (resumido en el capı́tulo 1).

Figura 5.20. Ángulo de reposo con material saturado - medidas en cm.

117
CAPÍTULO 5. ENSAYOS DE LABORATORIO MIC 2008-I-34

5.4. Equipo para medir el ángulo de reposo

El diseño del equipo para determinar el ángulo de reposo (Figura 5.21) se realizó tomando
como base la norma DIN 18126 y partiendo del principio que para determinar el ángulo
de reposo el material debe quedar en el estado más suelto posible. Esta norma indica
la metodologı́a para determinar la relación de vacı́os máxima y mı́nima de materiales
granulares. Este equipo se diseñó a partir de la Figura 5 de esta norma [1] donde describen
el aparato para determinar el arreglo más suelto. El mecanismo es prácticamente el mismo
con la diferencia de la falta del cilindro y la mesa vibratoria, que no son necesarios para
determinar el ángulo de reposo.

Figura 5.21. Ángulo de reposo - adaptado Norma Din 18126. 1. base, 2. guı́a para man-
tener centrado el embudo, 3. embudo, 4. sistema de poleas, 5. manivela.

Requerimientos estipulados en la Norma que se adaptan también para este ensayo:

La muestra no debe tener grumos, el material debe estar completamente seco (secado
al horno a 105 ◦C), no debe contener finos y evitar humedad en la muestra antes y

118
CAPÍTULO 5. ENSAYOS DE LABORATORIO MIC 2008-I-34

durante el ensayo.

Si el tamaño de la partı́cula más grande es 2mm, entonces se debe cumplir que el

diámetro del embudo sea 6 veces el diámetro de la partı́cula.

Si la partı́cula de mayor tamaño es 25 mm, el diámetro del embudo debe ser 5 veces
el diámetro de la partı́cula.

El procedimiento consiste en adicionar el material en el embudo antes de iniciar el

levantamiento del embudo (no adicionar material después), luego se levanta el embudo
por medio del sistema de poleas manteniendo el embudo centrado. La velocidad con
que se realice este ensayo debe garantizar que el embudo este en contacto con la pila
durante todo el ensayo.

Se debe hacer el ensayo 3 veces y promediar para dar el ángulo de reposo.

En la Figura 5.22 se muestra el equipo construido para determinar el ángulo de reposo de

materiales granulares.

Figura 5.22. Ángulo de reposo - adaptado Norma Din 18126.

119
Capı́tulo 6

Conclusiones

El modelo hipoplástico presenta una variación aproximada del 10 % en el valor del

esfuerzo desviador máximo, cuando el ángulo de fricción crı́tico varı́a un 7 % (2o en
este trabajo). Una variación del 10 % en la resistencia del suelo genera un aumento
del 5 % en el tamaño de una cimentación (análisis de una zapata cuadrada) y del
10 % en los costos.

La arena del Guamo se clasifica por su granulometrı́a como arena bien gradada y por
la forma de las partı́culas como material redondeado.

La forma de la pila construida por pluviación con un embudo depende de la rugosidad

de la superficie sobre la cual se construya. Si la superficie es rugosa los taludes de la
pila son de forma convexa. Si la superficie es lisa la forma de la pila es: curva en la
parte superior, lineal en la parte media del talud (donde se determina ϕrep ) y curva
en la pata del talud.

En las pilas construidas sobre una base lisa (vidrio) ocurren dos fenómenos en el
proceso de formación: 1. las partı́culas ruedan por la superficie de la pila y 2. las
partı́culas se deslizan en la base. En las pilas construidas sobre una base rugosa las
partı́culas ruedan por la superficie de la pila.

El embudo con el que se construya la pila para determinar el ángulo de reposo, debe
tener una diámetro de salida mayor o igual a 5 veces el tamaño máximo de la partı́cula
(φembudo ≥ 5φpartı́cula ). Ésto para que no haya ninguna influencia del embudo en
el ángulo de reposo y el flujo del material sea continuo.

La pila debe construirse depositando el material de forma lenta y con el embudo

siempre en contacto con la parte superior de la pila. Ésto para evitar que el impacto

120
CAPÍTULO 6. CONCLUSIONES MIC 2008-I-34

de las partı́culas afecte la estructura interna de la pila y permitir que el material

quede de la forma más suelta posible en su interior.

Los ángulos de reposo reportados en este trabajo corresponden al promedio de 5

mediciones. La diferencia entre el valor promedio y cada uno de los datos es ≤ 1o .

La influencia de la cantidad de material en el ángulo de reposo es diferente en cada

uno de los métodos empleado. En el método del cilindro “Train” la masa no influye
en ϕrep . En el método del embudo (pila) la masa influye en los resultados, el ϕrep
decrece con el aumento de la masa (en este trabajo la variación fue de 1.9o ). El
método de Santamarina fue que el que presentó mayor influencia (variación de 3.2o )
y una tendencia creciente del ϕrep con el aumento de la masa.

La velocidad de desplazamiento del embudo durante el proceso de formación de la

pila influye en el ángulo de reposo aproximadamente en 2o .

El diámetro de salida del embudo no influye en el ángulo de reposo (en este trabajo
la variación fue de 0.5o ), siempre y cuando se cumpla que φembudo ≥ 5φpartı́cula .

El tamaño de las partı́culas influye en el ángulo de reposo. La tendencia de ϕrep es

decreciente a medida que el tamaño de las partı́culas disminuye. En este trabajo la
variación de ϕrep con el tamaño de las partı́culas fue de 1.5o .

Determinar el ángulo de reposo con el material seco o sumergido no tiene ninguna

influencia.
0
La lı́nea del estado crı́tico en el plano q − p determinada a partir de los ensayos
realizados en este trabajo fue única. Los ensayos realizados fueron triaxiales CU con
arena del Guamo, estado denso. Estos ensayos se realizaron con tres muestras con la
misma relación de vacı́os inicial y presión de confinamiento variable. La trayectoria
de estas muestras durante la fase de corte del ensayo convergen todas en la misma
lı́nea del estado crı́tico.

El ángulo de reposo determinado por medio de la pila (método: embudo) y con 1000
gr de material fue 35.8o (ϕrep no varia significativamente para masas mayores a 1000
gr). El ángulo de fricción crı́tico calculado a partir del ensayo triaxial fue de 40 o . La
diferencia entre estos dos ángulos (ϕc y ϕrep ) es de 4o aproximadamente.

Teniendo en cuenta los resultados del ángulo de reposo (por diferentes métodos y
analizando diferentes variables) y del ángulo de fricción crı́tico determinados en este
trabajo, se concluye que para continuar empleando la simplificación ϕc = ϕrep se debe
tener en cuenta lo siguiente: 1. emplear una cantidad de material ≥ 1000 gr (arena);

121
CAPÍTULO 6. CONCLUSIONES MIC 2008-I-34

2. construir la pila con un desplazamiento lento del embudo y manteniéndolo siempre

en contacto con la parte superior de la pila; 3. tomar una muestra representativa del
material (todos los tamaños), ya que el tamaño de las partı́culas influye en el ángulo
de reposo (a > tamaño de las partı́culas, > ángulo de reposo); y 4. construir la pila
sobre una superficie rugosa.

122
Bibliografı́a

[1] DIN 18126. Bestimmung der dichte nichtbindiger böden bei lockerster und dichtester
lagerung. 1996.

[2] Pablo Andrés Arias. Modelos de comportamiento de suelos granulares: estudio y

determinación de sus parámetros. Master’s thesis, 2006.

[3] John Atkinson. An introduction to the mechanics of soils and foundations. 1993.

[4] J.H. Atkinson & P.L. Bransby. The mechanics of soil an introduction to critical state
soil mechanics. 1978.

[5] Muni Budhu. Soils mechanics and foundations. 2 edition, 2007.

[6] Roberto Cudmani. Modelación del comportamiento de suelos granulares a través de

una relación constitutiva hipoplástica. Technical report, 2004.

[7] Jake Dodds. Particle shape and stiffness - effects on soil behavior-. Master’s thesis,
Georgia Institute of Technology, 2003.

[8] American Society for testing and materials ASTM. Standard test method for bulk
solids characterization by carr indices. d 6393-99. 1999.

[9] Yan Graselli and Hans Herrmann. On the angles of dry granular heaps. Physica A,
247:301 – 312, 1997.

[10] Mehmet Yasar Gundogdu. Discharge characteristics of polydisperse powders through

conical hoppers. part 2: Predictions for coarse, granular, free flowing powders. Taylor
& Francis, (24):23 – 37, 2006.

123
BIBLIOGRAFÍA MIC 2008-I-34

[11] I. Herle and G. Gudehus. Determination the parameters of a hypoplastic constitutive

model from properties of grain assemblies. Mechanics of cohesive-Frictional Materials,
4:461–486, 1999.

[12] I. Herle and D. Kolymbas. Hypoplasticity for soils with low friction angles. Computers
and Geotechnics, 31(5):365–373, 2004.

[13] Dimitros Kolymbas. A rate-dependent constitutive equation for soils. Mech. Res.
Comm, 4:367–372, 1977.

[14] Dimitros Kolymbas. Introduction to hypoplasticity. Advances in Geotechnical Engi-

neering and Tunelling. Balkema, 2000.

[15] William Lambe and Robert Whitman. Mecánica de suelos. segunda edition, 1996.

[16] Kalman H.; Goder D.; Rivkin M. and Ben-Dor G. The effect of the particle-surface
friction coefficient on the angle of repose. Bulk Solids Handling, 13:123 – 128, 1993.

[17] Mooney M.A.; Finno R.J.; Viggiani M.G.. A unique critical state for sand. pages 1100
– 1108, Noviembre 1998.

[18] James Mitchell. Fundamentals of soil behavior. John Wiley, New York, 3rd edition,
2005.

[19] Britto A.M.; Gunn M.J. Critical state soil mechanics via finite elements. 1987.

[20] Andrzej Niemunis. Extended hypoplastic models for soils. Technical Report Heft 34,
Institut für Grundbau und Bodenmechanik der Ruhr Universität Bochum, Bochum,
Germany, 2003.

[21] Yan Graselli; Hans Herrmann; Gadi Oron and Zapperi Stefano. Shape of heaps and
in silos. The European Physical Journal B, 10:673 – 679, 1999.

[22] Yan Graselli; Hans Herrmann; Gadi Oron and Zapperi Stefano. Effect of impact
energy on the shape of granular heaps. Springer - Velag, 2:97 – 100, 2000.

[23] N. Pilpel. The flow of powders and granular solids. II(7):699 – 702, 1966.

[24] J.C. Santamarina. Soils and waves. John Wiley, 2001.

[25] J.C. Santamarina and G.C. Cho. Determination of critical state parameter in sandy
soils - simple procedure. Geotechnical Testing Journal, 24:185 – 192, 2001.

[26] J.C. Santamarina and G.C. Cho. Soil behaviour: The role of particle shape. The
Skemton conference, 1:604 – 617, 2004.

124
BIBLIOGRAFÍA MIC 2008-I-34

[27] Donald W. Taylor. Principios fundamentales de mecánica de suelos. 1962.

[28] David Train. Some aspects of the property of ancle of repose of powders. pages 127
– 135, 1958.

[29] L. E. Vallejo. Fractal analysis of granular materials. Géotechnique, 45(1):159 – 163,

1995. Technical note.

[30] Z. Chik & L. E. Vallejo. Characterization of the angle of repose of binary granular
materials. Can. Geotech. J., 42:683 – 692, 2005. Note.

[31] P.A. Wolffersdorff. A hypoplastic relation for granular materials with a predefined
limit state surface. Mechanics of Cohesive-Frictional Materials, 1:251–271, 1996.

[32] Wei Wu and Dimitros Kolymbas. Numerical testing of the stability criterion for
hypoplastic constitutive equations. Mechanics of Materials, 9:245–243, 1990.

125