MAESTRIA EN GEOTECNIA
BUCARAMANGA
2016
ANÁLISIS COMPARATIVO DE LAS DEFORMACIONES DE LOS SUELOS
MEDIANTE TEORÍAS TRADICIONALES Y SOFTWARE DE ELEMENTOS
FINITOS
Magíster en Geotecnia
Directora
MAESTRIA EN GEOTECNIA
BUCARAMANGA
2016
2
3
ENTREGA DE TRABAJOS DE GRADO, TRABAJOS DE
INVESTIGACION O TESIS Y AUTORIZACIÓN DE SU USO A
FAVOR DE LA UIS
Para constancia se firma el presente documento en dos (02) ejemplares del mismo valor y
tenor, en Bucaramanga, a los 23 días del mes de mayo de Dos Mil dieciséis 2016.
EL AUTOR / ESTUDIANTE:
4
DEDICATORIA
Dedico a Dios este logro, que gracias a los talentos recibidos por El, fortaleza,
constancia, disciplina y demás virtudes, pude superarme a mí misma y ser más
fuerte que los obstáculos.
A mis Padres, por brindarme su apoyo, amor infinito e incondicionalidad, por creer
siempre en mí, porque sin ellos esto no tendría el mismo valor.
A mi familia que de una u otra manera han sido participes de mi proceso y que
tuvieron una palabra de apoyo para mí durante mis estudios.
5
AGRADECIMIENTOS
Agradezco a Dios por colocar personas valiosas en este camino, con grandes
carismas y de gran calidad humana.
A mi directora de proyecto Luz Marina Torrado y a los docentes Wilfredo del Toro
y Jose Alberto Rondón, que me guiaron en este camino de aprendizaje, por su
ejemplo de superación y buen profesionalismo.
6
TABLA DE CONTENIDO
1. OBJETIVOS ............................................................................................... 18
1.1 OBJETIVO GENERAL: .............................................................................. 18
1.2 OBJETIVOS ESPECÍFICOS: ..................................................................... 18
2. MARCO TEÓRICO..................................................................................... 19
2.1 FUNDACIONES SUPERFICIALES ........................................................... 19
2.1.1 Número y profundidad de las perforaciones .............................................. 19
2.2 CAPACIDAD DE CARGA (CAPACIDAD PORTANTE / CAPACIDAD DE
SOPORTE) ............................................................................................................. 20
2.2.1 Factor de seguridad ................................................................................... 21
2.3 PARÁMETROS GEOMECÁNICOS PARA EL CÁLCULO DE LA
CAPACIDAD DE CARGA ....................................................................................... 22
2.3.1 Relación de vacíos (e): .............................................................................. 22
2.3.2 Peso saturado (үsat): ................................................................................. 22
2.3.3 Ángulo de fricción (φ):................................................................................ 23
2.4 ECUACIONES DE CAPACIDAD DE CARGA (CAPACIDAD PORTANTE) 23
2.4.1 Terzaghi (1943) ......................................................................................... 23
2.4.2 Meyerhof (1963) ........................................................................................ 26
2.4.3 Hansen (1970) ........................................................................................... 29
2.4.4 Vesic (1975) ............................................................................................... 32
2.4.5 Brinch – hansen ......................................................................................... 33
2.4.6 Capacidad de carga en condiciones no drenadas – (arcilla φ = 0) ............ 35
2.4.6.1 Efecto de la posición del nivel freático ....................................................... 37
2.4.6.1.1 Caso I: Nivel freático por encima de la profundidad de cimentación. ..... 37
2.4.6.1.2 Caso II: Nivel freático por debajo de la profundidad de cimentación y d ≤
B 37
2.5 CAPACIDAD DE CARGA A PARTIR DE ENSAYOS DE CAMPO ............ 38
2.5.1 Ensayo de penetración estándar (spt) ....................................................... 38
2.5.1.1 Factor de corrección por energia ............................................................... 41
7
2.5.1.2 Factor de corrección por confinamiento ..................................................... 42
2.5.1.3 Correlaciones con parámetros geotécnicos ............................................... 43
2.6 DEFORMACIONES EN LOS SUELOS...................................................... 43
2.6.1 Deformaciones elásticas o inmediatas....................................................... 46
2.6.1.1 En arcillas .................................................................................................. 46
2.6.1.1.1 Espacio semi-infinito .............................................................................. 47
2.6.1.1.2 Capa limitada por base rígida ................................................................ 49
2.6.1.1.3 Dos o más capas limitadas por una base rígida .................................... 50
2.6.2 Asentamientos por consolidación .............................................................. 51
2.6.2.1 Asentamiento en arcillas normalmente consolidadas (schmertmann 1970)54
2.6.2.2 Asentamiento en arcillas pre-consolidadas (schmertmann 1970) .............. 57
2.6.2.2.1 Caso ∆σ ˂ σ’p – σ’0 ................................................................................ 57
2.6.2.2.2 Caso ∆σ > σ’p – σ’0 ................................................................................ 58
2.6.3 Deformaciones por consolidación primaria ................................................ 59
2.6.4 Deformaciones por consolidación secundaria ........................................... 60
2.6.5 Asentamiento tolerable de las fundaciones de estructuras ........................ 61
3. METODOLOGÍA ........................................................................................ 62
3.1 CALCULO DE CAPACIDADES PORTANTES .......................................... 62
3.1.1 Datos generales ......................................................................................... 62
3.1.1.1 Database terrenos ..................................................................................... 63
3.1.1.2 Datos sistema de cimentación ................................................................... 63
3.1.1.3 Estratigrafía ............................................................................................... 63
3.1.1.4 Input gráfico ............................................................................................... 63
3.1.1.5 Cargas ....................................................................................................... 68
3.1.2 Carga última .............................................................................................. 69
3.1.3 Gráficos ..................................................................................................... 70
3.1.3.1 Diagrama de tensiones .............................................................................. 70
3.1.3.2 Parámetros gráficos ................................................................................... 70
3.1.3.3 Informe diagramas ..................................................................................... 71
3.1.3.4 Gráfico carga Última – Base (Q-B) ............................................................ 71
8
3.2 DEFORMACIONES O ASENTAMIENTOS ................................................ 79
3.2.1 Deformaciones elásticas o inmediatos mediante la teoria clásica ............. 79
3.2.2 Deformaciones elásticas o inmediatos mediante elementos finitos (plaxis
3d) 85
3.2.2.1 Introducción de datos y pre-procesador 3d................................................ 85
3.2.2.1.1 Geometría de entrada ............................................................................ 85
3.2.2.1.2 Geometría del contorno ......................................................................... 86
3.2.2.1.3 Planos de trabajo ................................................................................... 87
3.2.2.1.4 Propiedades del material ....................................................................... 87
3.2.2.1.5 Borehole................................................................................................. 95
3.2.2.1.6 Aplicación de carga ................................................................................ 95
3.2.2.1.7 Generación de la malla .......................................................................... 96
3.2.2.2 Cálculos ..................................................................................................... 97
3.2.2.2.1 Puntos de control ................................................................................... 98
3.2.2.3 Post-procesos de los resultados ................................................................ 99
3.2.3 Deformaciones o asentamientos por consolidación unidimensional ........ 101
3.2.3.1 Incremento del esfuerzo bajo un terraplén............................................... 101
3.2.3.2 Incremento del esfuerzo debajo de una cimentación por el método 2:1 .. 103
3.2.4 Deformaciones o asentamientos por consolidación unidimensional
mediante elementos finitos (plaxis 2d) .................................................................. 105
3.2.5 Deformaciones o asentamientos diferenciales ........................................ 107
4. ANÁLISIS COMPARATIVO...................................................................... 108
4.1 CAPACIDADES PORTANTES ................................................................ 108
4.2 DEFORMACIONES ELÁSTICAS ............................................................ 108
4.3 DEFORMACIONES POR CONSOLIDACIÓN ......................................... 109
4.4 DEFORMACIONES TOTALES POR LOS METODOS CLÁSICOS Y DE
ELEMENTOS FINITOS ......................................................................................... 110
4.5 MEJORAMIENTO AL MODELO DE PLAXIS 2D ..................................... 111
5. CONCLUSIONES..................................................................................... 116
BIBLIOGRAFIA ..................................................................................................... 117
9
ANEXOS .................................................................. ¡Error! Marcador no definido.8
10
LISTA DE FIGURAS
11
Figura 30. Convención del signo de las solicitaciones. ......................................... 68
Figura 31. Ventana selección del tipo de carga. ................................................... 68
Figura 32. Ventana de cargas para selección de opciones de análisis. ................ 69
Figura 33. Gráfico de tensiones según Boussinesq. ............................................. 70
Figura 34. Tabla de parámetros gráficos. ............................................................. 71
Figura 35. Gráfico Carga última – Profundidad (B-Qlim). ...................................... 72
Figura 36. Perfil de suelo con la carga del muro. .................................................. 80
Figura 37. Perfil de suelo con la carga de la estructura. ....................................... 83
Figura 38. Ventana de configuración general (General Settings).......................... 86
Figura 39. Ventana de workplanes........................................................................ 87
Figura 40. Ventana de propiedades del material arcilloso. ................................... 94
Figura 41. Ventana de Borehole. .......................................................................... 95
Figura 42. Ventana de introducción de cargas. ..................................................... 96
Figura 43. Ventana de generación de la malla 2D. ............................................... 96
Figura 44. Ventana de generación de la malla 3D. ............................................... 97
Figura 45. Ventana de cálculos (definición de fases). ........................................... 98
Figura 46. Ventana de selección de puntos de control en el nivel -3.00 m. .......... 99
Figura 47. Gráfico de deformaciones elásticas generadas por PLAXIS 3D. ....... 100
Figura 48. Gráfico de curva de desplazamientos de la fase carga de la estructura.
............................................................................................................................. 100
Figura 49. Incremento de esfuerzo bajo un terraplén. ......................................... 102
Figura 50. Método 2:1 para encontrar el incremento de esfuerzo debajo de una
cimentación. ......................................................................................................... 104
Figura 51. Modelo generado en PLAXIS 2D. ...................................................... 106
Figura 52. Cimentación con las distancias entre columnas. ............................... 107
Figura 52. Capacidades de carga a un 1.0 m de profundidad. ........................... 108
Figura 53. Comparativo de deformaciones elásticas. ......................................... 109
Figura 54. Comparativo de deformaciones por consolidación. ........................... 110
Figura 55. Comparativo de deformaciones totales. ............................................. 111
Figura 56. Modelo de la malla deformada con refuerzo de geotextil. .................. 112
Figura 57. Deformaciones totales con refuerzo de geotextil. .............................. 112
Figura 58. Fuerzas axiales generadas por el geotextil. ....................................... 113
Figura 59. Deformaciones totales incluyendo el geotextil. .................................. 114
Figura 60. Comparativo de deformaciones por consolidación teniendo en cuenta el
geotextil en PLAXIS 2D. ...................................................................................... 115
Figura 61. Comparativo de deformaciones totales teniendo en cuenta el geotextil
en PLAXIS 2D. ..................................................................................................... 115
12
LISTA DE TABLAS
13
Tabla 25. Resumen de parámetros para el material arcilla 2. ............................... 91
Tabla 26. Resumen de parámetros para el material limo 1. .................................. 92
Tabla 27. Resumen de parámetros para el material limo 2. .................................. 93
Tabla 28. Asentamientos por consolidación generados por la carga del muro. .. 102
Tabla 29. Asentamientos por consolidación generados por la carga estructural. 104
Tabla 30. Rango de la permeabilidad hidráulica parra varios suelos. ................. 106
Tabla 31. Cálculos de asentamientos diferenciales según la NSR-10. ............... 107
14
TÍTULO: ANÁLISIS COMPARATIVO DE LAS DEFORMACIONES DE LOS
SUELOS MEDIANTE TEORÍAS TRADICIONALES Y SOFTWARE DE
ELEMENTOS FINITOS*
RESUMEN:
__________________
*Trabajo de aplicación
**Facultad de Ingenierías Fisico-mecánicas. Escuela de Ingeniería Civíl. Maestría en Geotecnia.
Directora Luz Marina Torrado Gómez.
15
TITLE: COMPARATIVE ANALYSIS OF SOIL DEFORMATION BY
TRADITIONAL THEORIES AND FINITE ELEMENT SOFTWARE*
ABSTRACT:
In addition to the deformations and the load capacity which is subject to the soil, so
you can compare these two methods to find out how outdated is each other.
This analysis is performed for a specific case: ‘’ESTACIÓN 3’’, which involves the
construction of a wall on reinforced earth, which has clay soils, soils of medium
compressibility, low plasticity and the presence of water table.
For the construction of this wall on reinforced earth, control and registration of the
immediate deformation made by topography.
This project also it seeks to compare the strains obtained by the theoretical
equations and finite element against deformations caused by this type of soil in
reality.
Finally, the comparative analysis was performed theories and what quantitative
differences can be obtained; insomuch that theoretically calculated is always an
approximate and not an exact value, especially when you have conditions of
consolidation.
__________________
*Aplication work
**School of Civil Engineering. Faculty of Physico-Mechanical Engineering. Master in Geotechnics.
Director Luz Marina Torrado Gómez.
16
INTRODUCCIÓN
Como bien dice la frase, una buena construcción siempre dependerá como
primera medida de buenos cimientos, de ahí la importancia de conocer todo sobre
la cimentación.
Los cimientos son la base de una construcción, una cimentación mala o pobre es
un augurio de un futuro poco promisorio para toda estructura1.
Otros factores influyentes que garantizan un buen estudio geotécnico son los
ensayos de campo, laboratorios y la buena elección y uso de herramientas de
cálculos que nos ofrece la actualidad.
Por ello, en este trabajo aplicativo tiene por objetivo realizar un análisis
comparativo mediante las ecuaciones tradicionales y el software de elementos
finitos. Esto en pro de la geotecnia y en la búsqueda de nuevas alternativas.
1
IDEAS PARA CONSTRUIR. ‘’Funciones e importancia de la cimentación’’. Internet :
(˂http://ideasparaconstruir.com/n/4015/funciones-e-importancia-de-la-cimientacion.html>).
17
1. OBJETIVOS
18
2. MARCO TEÓRICO
ESTUDIO DE
SUELOS
• Perfil estratigráfico
• Espesores de los estratos
• Características de los suelos
• Posición del nivel freático • Número de
sondeos
• Profundida
• Profundidad de cimentación
d de las
• Capacidad de carga perforacion
• Asentamientos o deformaciones es
Una capacidad de carga segura es aquella que el suelo puede cargar de forma
segura y sin deformarse demasiado; se conoce como capacidad de carga
admisible y se determina por2:
qult
qadm = → 2 ≤ F.S. ≤ 3
F.S.
2
DEL TORO, WILFREDO. ‘’Capacidad de carga’’. Geotecnia II. Universidad Industrial de
Santander. Marzo de 2015.
20
Figura 4. Curva de capacidad de carga.
q
qult
qadm
δ
qult
qadm =
F.S.
qneta = qult − q q = γD f
qult − q
qadm( neta) =
F.S.
21
2.3 PARÁMETROS GEOMECÁNICOS PARA EL CÁLCULO DE LA
CAPACIDAD DE CARGA
G *γ
e = s w − 1
γ unsat
γ w * ( e + Gs )
γ sat =
1+ e
22
2.3.3 Ángulo de fricción (φ):
23
Para los terrenos muy sueltos, Terzaghi aconseja tener en consideración la carga
última; el valor que se calcula con la fórmula anterior pero introduciendo valores
reducidos de las características mecánicas del terreno y precisamente3:
2 2
cred = * c
Φ red = arctan * tan Φ
3 3
Donde:
3
Geostruct Software. Manual de loadcap. Italia. Loadcap versión 8.0.0.
24
Donde:
Nq =
a2 a = e(0.75π −φ / 2) tanφ
2 cos2 (45 + φ / 2)
tan φ K pγ
N c = ( N q − 1) cot φ Nγ = − 1
cos φ
2
2
Φ deg Nc Nq Kpγ
Nγ
10.8
0 5.7 1.0 0.0
14.7
10 9.6 2.7 1.2
25.0
20 17.7 7.4 5.0
35.0
25 25.1 12.7 9.7
25
35 57.8 41.6 42.4 82.0
141.0
40 95.7 81.3 100.4
800.0
50 347.5 415.1 1153.2
• Carga Vertical
• Carga Inclinada
Donde:
φ N γ = ( N q − 1) tan( 1 . 4φ )
N q = eπ tanφ tan2 45 +
2
Φ
N c = ( N q − 1) cot φ K p = tan 2 45 +
2
26
Tabla 3. Factores de capacidad de carga para las ecuaciones de Meyerhof,
Hansen y Vesic.
2tanΦ(1-
∅ Nc Nq Nγ(M) Nγ(H) Nγ(V) Nq/Nc
senΦ)2
133.7
45 134.9 262.7 200.8 271.3 1.009 0.172
3
27
50 266.5 319.1 873.9 568.6 761.3 1.197 0.130
B
s c = 1 + 0.2 K p Cualquier φ
L
B
Forma sq = sγ = 1 + 0.1K p φ>10°
L
sq = sγ = 1 φ=0
D
dc = 1 + 0.2 K p Cualquier φ
B
Profundidad D
dq = dγ =1+ 0.1 Kp φ>10°
B
d q = dγ = 1 φ=0
Cualquier φ
θ°
2
θ = inclinación de la
ic = iq = 1− resultante en la
90° vertical
2
θ°
iγ 1 − φ>0
Inclinación φ °
iγ = 0 φ=0
28
2.4.3 Hansen (1970)
Cuando φ = 0
qult = 5.14Cu(1 + S 'c +d 'c −ic − bc − g c ) + γD
----------------------------------------------------------------------------------
----------------------------------------------------------------------------------
----------------------------------------------------------------------------------
Donde:
N γ = 1 .5 ( N q − 1) tan φ φ
N q = eπ tanφ tan2 45 +
2
N c = ( N q − 1) cot φ
29
Tabla 5. Factores de profundidad de las ecuaciones de capacidad de carga de
Hansen.
FACTORES DE PROFUNDIDAD
D −1 D
d c = 1 + 0.4 d c = 1 + 0.4 tan
B B
2 D 2 −1 D
d q = 1 + 2 tan ϕ (1 − sin ϕ ) d q = 1 + 2 tan ϕ (1 − sin ϕ ) tan
B B
En los factores siguientes las expresiones con ápices (') valen cuando φ=0.
FACTORES DE
FACTORES DE FORMA
PROFUNDIDAD
'
' B d 'c = 0.4 k
s 'c = 0.2
L
d c = 1 + 0.4 k
Nq B
sc = 1 + d q = 1 + 2 tan ϕ (1 − sin φ ) k
Nc L
sc = 1 para cimentaciones continuas
d γ = 1 para cualquier φ
B D D
sq = 1 + tan φ k= si ≤1
L B B
B −1 D D
sγ = 1 − 0.4 k = tan si >1
L B B
30
Tabla 7. Factores de inclinación de las ecuaciones de capacidad de carga de
Hansen y Vesic.
mH
ic' ( H ) = 0.5 − 0.5 1 −
H ic' (V ) = 1 −
Af ca Af ca N c
1− iq
ic = iq − _(Hansen _ y _Vesic)
Nq −1
5 m
0.5 H H
= 1 − = 1 −
V + A c cot φ V + A c cot φ
iq ( H ) iq (V )
f a f a
5
0.7 H
V + Af ca cot φ (
'
iγ ( H ) = 1 − _ η = 0)
m +1
( 0.7H −η° / 450) H
5
H
_ (η > 0) iγ' (V ) = 1 −
iγ ( H ) = 1−
+ φ V + A f ca cot φ
V Af ca cot
2+ B / L 2+ L / B
m = mB = __ H _ paralelo _ a _ B m = mL = __ H _ paralelo _ a _ L
1+ B / L 1+ L / B
Nota: iq , iγ > 0
31
Tabla 8. Factores de inclinación de la superficie y factores de base inclinada de
las ecuaciones de capacidad de carga de Hansen y Vesic.
β° η° η°
g c' = bc' = bc = 1 −
147° 147° 147°
β°
gc = 1 − bq(V ) = bγ (V ) = (1−η tanφ)2
147°
β + η ≤ 90°
g q ( H ) = gγ ( H ) = (1 − 0.5 tan β )
5
β ≤φ
g q (V ) = gγ (V ) = (1 − tan β )
2
N γ = 2 ( N q + 1) tan φ
Los factores de forma y de profundidad que aparecen en las fórmulas del cálculo
de la capacidad portante son iguales a los propuestos por Hansen, como se
muestran en la figura 5.; mientras se dan algunas diferencias en los factores de
32
inclinación de la carga, del terreno (cimentación en talud) y del plano de
cimentación (base inclinada), ver tablas 6 y 7.
Para que una cimentación pueda resistir la carga de proyecto en seguridad con
respecto a la rotura general, para todas las combinaciones de carga relativas al
ELU (Estado límite último), se debe dar la siguiente desigualdad3:
Vd ˂ Rd
Donde:
R
= ( 2 + π ) c u s c ic + q
A'
Donde:
cu Es la cohesión no drenada.
33
(
i c = 0,5 1 + 1 − H / A' c u ) Factor de corrección de la inclinación de la carga
debida a una carga H.
R
= c' N c Sc ic + q' N q S q d q iq + 0.5γ ' B' Nγ Sγ dγ iγ
A'
Donde:
N γ = 1 .5 ( N q − 1) tan φ ' φ
N q = eπ tanφ ' tan2 45 +
2
N c = ( N q − 1) cot φ '
Factores de forma
FACTORES DE FORMA
B'
sq = 1 + senφ ' Para forma rectangular
L'
s q = 1 + senφ ' Para forma cuadrada o circular
sγ = 1 − 0.3
B' Para forma rectangular
L'
sγ = 0.7 Para forma cuadrada o circular
s N −1
q q
sc = Para forma rentangular, cuadrada o circular
qN − 1
34
Tabla 10. Factores de inclinación de la ecuación de capacidad de carga de Brinch
– Hansen.
FACTORES DE Factores
FACTORES DE INCLINACIÓN
inclinación resultante debido a
DEBIDO A UNA CARGA
una carga horizontal H paralela a
HORIZONTAL H PARALELA A L’
B’
De forma rectangular
De forma rectangular
3
1− H iq =
1 − 0 .7 H
iq = i = V − A' c ' cot φ '
g V + A' c ' cot φ '
i N −1
q q
i =
c N −1
De forma cuadrada o circular q
35
Las ecuaciones quedan:
qult = cu * Nc * Sc * γD f
36
2.4.6.1 Efecto de la posición del nivel freático
0 ≤ D1 ≤ D f
0 ≤ D1 ≤ D f
37
El factor γ en el último término queda:
γ = γ '+
d
(γ − γ ')
B
Este método permite obtener parámetros de resistencia del suelo tales como, el
ángulo de fricción (φ), la cohesión (c’), la resistencia al corte no drenado (Su) y
algunos parámetros de deformabilidad como el módulo de Young (E), a partir de
correlaciones propuestas por diferentes autores y que son de amplio uso en el
diseño de diferentes estructuras geotécnicas dada la facilidad de su uso44.
4
SORIANO, Cristian Yair. Aplicación del ensayo de penetración estándar en la determinación de
parámetros geotécnicos de suelos granulares. Ingeniero Civíl. Universidad Nacional de Colombia.
38
Figura 9. Muestreador de tubo partido5.
5
INSTITUTO NACIONAL DE VÍAS. Norma INV E-111-07. Ensayo de Penetración normal (SPT) y
muestreo de suelos con tubo partido. Bogotá D.C., 1996. p.4.
39
Figura 10. Montaje del equipo de perforación66.
El valor que se reporta en este caso es la suma de los dos últimos valores N2 + N3.
6
INSTITUTO NACIONAL DE VÍAS. Norma INV E-111-07. Ensayo de Penetración normal (SPT) y
muestreo de suelos con tubo partido. Bogotá D.C., 1996. p.6.
40
El ensayo presenta unas condiciones de rechazo si:
Los resultados del ensayo deben corregirse a partir de consideraciones que tienen
en cuenta la energía, por lo tanto, el valor que se obtiene del Ncampo debe incluir
los siguientes factores7:
Ncorr = CN * Ncampo * η
η = η1 * η2 * η1 * η3 * η4
Donde:
Este factor representa el rendimiento del impacto del martillo sobre el cabezal de
golpeo, se define como4:
N1 = Er/E60
Donde:
7
BOWLES J.E, Foundation analysis and design. 5th Ed. McGraw-Hill Book Co, 1997.
41
Er es el valor de la energía de referencia de la barra.
95.76
CN =
σ 'vo
8
LOPEZ M.R. Determinación in situ de propiedades ingenieriles de los suelos y su relación con el
ensayo de penetración; Buenos Aires: Universidad Tecnológica Nacional, 2003. p. 23.
9
GONZALEZ A.J. Estimativos de parámetros efectivos de resistencia con el SPT. Bogotá D.C: X
jornadas geotécnicas de ingeniería colombiana. 1999. p.23.
42
2.5.1.3 Correlaciones con parámetros geotécnicos
Tabla 12. Correlaciones para estimar el ángulo de fricción a partir del ensayo de
SPT8.
43
deformaciones pueden clasificarse en asentamientos inmediatos, por
consolidación o por compresión secundaria10.
10
ANTOLINEZ, Lizzet y AVILA, Jefferson. Análisis comparativos entre teorías usadas en el cálculo
de asentamientos y un software basado en elementos finitos. Trabajo de grado pregrado en
Ingeniería Civíl. Bucaramanga: Universidad Industrial de Santander. 2013.
44
Figura 12. Asentamiento inmediato o elástico3.
45
Este proceso se conoce como consolidación primaria y el análisis se lleva a cabo
con los diferentes modelos de la teoría de la consolidación. El asiento que deriva
de este proceso de expulsión del agua de los vacíos intersticiales se denomina
deformación por consolidación.
2.6.1.1 En arcillas
46
2.6.1.1.1 Espacio semi-infinito
(1 − υ 2 )
ρ = qc * B Cd
Es
Siendo:
ν: Coeficiente de Poisson
E: Módulo de elasticidad
47
Para el cálculo de los parámetros geomecánicos del coeficiente de Poisson y
módulo de elasticidad (ver tabla 13), se hizo uso de las siguientes expresiones:
K0 K 0 = 1 − senφ 1 − senφ
ν= µ=
1 + K0 2 − senφ
RECTANGULAR
Largo/Ancho
48
2.6.1.1.2 Capa limitada por base rígida
Figura 15. Asentamiento inmediato en arcillas en una capa limitada por base
rígida.
qc * B * (1 −υ 2 )
ρ= C' d
E
Siendo:
Tabla 15. Factores de forma y rigidez en una capa limitada por base rígida.
CIRCULAR RECTANGULAR
49
0.10 0.09 0.09 0.09 0.09 0.09 0.09 0.09
Figura 16. Asentamiento inmediato en dos o más capas limitadas por una base
rígida.
qc * B * (1 −υ 2 )
ρ= C'd
E
50
Siendo:
ρT = ρ1 + ρ2
Se calcula el asentamiento total de la capa H1+ H2, con las características del
suelo de la capa 2.
Se calcula el asentamiento de la capa 1 con las características del suelo de la
capa 2.
Este último valor es restado del calculado para H1+ H2.
51
Figura 17. Asentamiento por consolidación.
Así como en las consideraciones del suelo en sus fases, la variación volumétrica
es:
∆e
∆dV = dz
1+ e
∆e
ρ= H
1+ e
53
2.6.2.1 Asentamiento en arcillas normalmente consolidadas
(schmertmann 1970)
H σ 0' + ∆σ
ρ= CcLog
1 + e0 σ0
'
Donde:
54
Y según Nagarja y Murty (1985) propusieron:
LL %
Cc = 0.2343 Gs
100
Siendo:
G *γ
e0 = s w − 1
γ unsat
Donde:
Gs : Gravedad específica
55
Para un número arbitrario de capas:
σ0 σ0+Δσ
Hi σ 0' + ∆σ
ρ = ∑ Cci Log
1 + e0i σ 0 i
'
56
2.6.2.2 Asentamiento en arcillas pre-consolidadas (schmertmann 1970)
Figura 22. Curva de consolidación para una arcilla pre-consolidada caso ∆σ ˂ σ’p
– σ’0
Cs
H σ 0' + ∆σ
ρ= CsLog
1 + e0 σ0
'
Donde:
57
LL %
Cs = 0.0463 Gs
100
Figura 23. Curva de consolidación para una arcilla pre-consolidada caso ∆σ > σ’p
– σ’0
Cs
H σ 'p σ ' + ∆σ
ρ= CsLog + CcLog 0 '
σ
1 + e0 σ '0 p
Donde:
58
Figura 24. Asentamiento por consolidación para arcillas pre-consolidada con un
número arbitrario de capas.
H σ ´ σ ´+ ∆σ Hi σ ´+ ∆σ
ρ = Σ i Csi Log p + Cci Log o ρ = Σ * Csi Log o
σ ´
1 + eoi σ o ´ i p 1 + eoi σ o´ i
59
2.6.4 Deformaciones por consolidación secundaria
T
∆H s = HCs Log
T
100
Donde:
Cs es el índice de expansión
11
Reglamento colombiano de construcción sismo resistente, NSR 10 (2010). Bogotá. Título H, y
Título G.
61
3. METODOLOGÍA
Para los cálculos de capacidades portantes se hizo uso del software LOADCAP,
el cual facilita el análisis por las diferentes ecuaciones tradicionales de Terzaghi,
Meyerhoft, Hansen, Vesic y Brinch – Hansen.
Esto con el fin de comparar los resultados obtenidos por los dos métodos.
Permite incluir los datos del proyecto tales como: descripción, proyectista, Cliente,
fecha, tipo del terreno, corrección de parámetros según Terzaghi y la acción
sísmica. Ver figura 25.
62
3.1.1.1 Database terrenos
Este comando permite administrar una base de datos de terrenos a través de una
ventana de diálogo subdividida en tres columnas: en la columna izquierda se
presenta la nomenclatura, en la central los parámetros geotécnicos y la derecha
las texturas. Ver figura 26.
Con este comando LoadCap permite definir todos los datos geométricos relativos
al tipo de cimentación seleccionada bien sea zapata corrida, aislada, losa,
cimentación circular; a la presencia de una subcimentación; del perfil del terreno y
de la profundidad del nivel freático. Ver figura 27.
3.1.1.3 Estratigrafía
Muestra una tabla donde permite observar, revisar y/o modificar los parámetros
geotécnicos (introducidos inicialmente en los materiales de Database terrenos) de
cada estrato y la textura a asociada al mismo: número de estratos, database de
terrenos, altura del estrato, peso específico del suelo, peso específico del suelo
saturado, cohesión, modulo elástico, módulo edométrico, datos adicionales
(coeficientes de consolidación primaria y secundaria), texturas y descripciones del
suelo. Ver figura 29.
Esta opción nos facilita cuatro comandos: introducir estrato, modificar espesor,
asignar terreno y eliminar estrato.
63
Figura 25. Datos generales de entrada al software.
64
Figura 26. Ventana de la base de datos de materiales.
65
Figura 27. Ventana de los datos del sistema de cimentación.
66
Figura 28. Ventana del perfil del suelo.
67
3.1.1.5 Cargas
Las cargas de proyecto se usan para estimar la carga última. Los valores a
introducir son de proyecto por lo tanto tendrán que ser asignados en caso de que
se realice una comprobación al estado límite último y comprenden los coeficientes
de amplificación de las acciones.
Las cargas de servicio se usan para el cálculo de los asientos. Se pueden asignar
al mismo tiempo varias combinaciones de carga de proyecto y de servicio para
efectuar el cálculo según la normativa.
68
3.1.2 Carga última
Para el cálculo de la carga última se puede elegir entre varios autores: Terzaghi,
Meyerhof, Hansen, Vesic y Brinch – Hansen.
69
Para el cálculo del coeficiente de balasto ks el programa propone la fórmula
adoptada por Bowles.
3.1.3 Gráficos
70
Figura 34. Tabla de parámetros gráficos.
Muestra una tabla con los diferentes valores de la carga admisible en función de la
profundidad del plano de cimentación D y de la anchura B, para cada longitud L
seleccionada de entre las que se proponen. La tabla se reconstruye con base en
los parámetros generales de los gráficos asignados en el menú.
71
Figura 35. Gráfico Carga última – Profundidad (B-Qlim).
Ver Anexo.
DATOS GENERALES
======================================================
=======================================================
72
ESTRATIGRAFÍA TERRENO
DH: Espesor del estrato; γunsat: Peso específico seco; γsat: Peso específico
saturado; φ: Ángulo de rozamiento interno; Cohesión; E: Módulo elástico; ν:
Relación de Poisson.
DH γunsat γsat φ c E
ν
[m] [kN/m³] [kN/m³] [°] [kN/m²] [kN/m²]
Presión
normal
Nombre N Mx My Hx Hy
de
Nr. Tipo
combinación proyecto
[kN] [kN·m] [kN·m] [kN] [kN]
[kN/m²]
73
CARGA ÚLTIMA CIMENTACIÓN COMBINACIÓN…Carga última
• Carga última
======================================================
74
Factor inclinación cargas [Iq] 1.0
======================================================
======================================================
======================================================
75
Factor corrección sísmico inercial [zg] 1.0
======================================================
======================================================
======================================================
76
======================================================
======================================================
======================================================
77
Factor inclinación base [Bg] 1.0
======================================================
======================================================
======================================================
78
Factor forma [Sg] 0.7
======================================================
======================================================
qc * B * (1 − υ 2 )
ρ= C 'd
E
Para ello, se tiene el siguiente perfil de suelo junto con sus respectivos
parámetros:
79
Figura 36. Perfil de suelo con la carga del muro.
γunsat E
TIPO DE γsat
ESPESOR ν
MATERIAL ESTRATOS (KN/m3) (KN/m3) (KN/m2)
(m)
80
Inicialmente se calculan los asentamientos elásticos generados por la carga del
muro.
H (m) 0.90
H/B 0.10
B (m) 8.7
ρ1 (m) 0.008
ρ2 (m) 0.027
81
ESTRATO No. 3 (Limo 1) ESTRATO No. 3 (Limo 1)
ρ3 (m) 0.018
ρ4 (m) 0.020
82
Figura 37. Perfil de suelo con la carga de la estructura.
qc = 7.1KN / m 2
H (m) 0.90
H/B 0.90
B (m) 1
ρ1 (m) 0.001
83
ESTRATO No. 2 (Arcilla 2) ESTRATO No. 2 (Arcilla 2)
ρ2 (m) 0.0004
ρ3 (m) 0.00006
ρ4 (m) 0.0001
84
ρ elásticos estructura (m) 0.001
12
R.B.J. Brinkgreve. Manual de plaxis. Países Bajos. Plaxis b.v. 2004.
85
Figura 38. Ventana de configuración general (General Settings).
86
3.2.2.1.3 Planos de trabajo
En Plaxis 3D las propiedades del suelo y de los materiales que constituyen las
estructuras son almacenadas en conjuntos de datos de los materiales.
87
• La pestaña Interfaces contiene los parámetros que relacionan las
propiedades de las interfaces con las propiedades del suelo. En este caso será
Rigido porque se quiere que la interfaz no altere la resistencia del suelo que le
rodea.
GENERAL
Mohr-
Modelo de material Model Coulomb -
PARAMETERS
INTERFACES
88
Tabla 23. Resumen de parámetros para el material de base granular.
GENERAL
Base
Nombre del material Identification granular -
Mohr-
Modelo de material Model Coulomb -
PARAMETERS
INTERFACES
89
Tabla 24. Resumen de parámetros para el material arcilla 1.
GENERAL
Mohr-
Modelo de material Model Coulomb -
PARAMETERS
INTERFACES
90
Tabla 25. Resumen de parámetros para el material arcilla 2.
GENERAL
Mohr-
Modelo de material Model Coulomb -
PARAMETERS
INTERFACES
91
Tabla 26. Resumen de parámetros para el material limo 1.
GENERAL
Mohr-
Modelo de material Model Coulomb -
PARAMETERS
INTERFACES
92
Tabla 27. Resumen de parámetros para el material limo 2.
GENERAL
Mohr-
Modelo de material Model Coulomb -
PARAMETERS
INTERFACES
93
Figura 40. Ventana de propiedades del material arcilloso.
94
3.2.2.1.5 Borehole
Al igual que en la opción workplane permite insertar los estratos de suelo para la
estructura pero adicional se va insertando el tipo de material en la columna
estratigráfica. Dando clic en Hydrostatic se incluye el nivel freático.
Tener en cuenta que el modelo de construyó a partir del nivel cero hacia abajo.
Ver figura 48.
95
este caso en el eje Y, es decir perpendicular al plano de trabajo) se inserta el
valor.
96
Figura 44. Ventana de generación de la malla 3D.
3.2.2.2 Cálculos
Para empezar la fase de cálculo se hace necesario definir qué tipos de cálculos se
han de llevar a cabo y qué tipo de cargas o de etapas de construcción se han de
activar durante los cálculos. Todo ello puede hacerse dentro del programa de
cálculo Calculations.
Dando clic en Phases se insertan las fases en que se dividirá el proyecto, para
este caso se tendrán en cuenta tress fases, la fase inicial y la fase I y la fase II.
97
Figura 45. Ventana de cálculos (definición de fases).
• Fase I (Carga estructura): Con esta etapa se finalizan las fases del
modelo, y es donde se aplica la carga al suelo con el peso de la estructura (7.1
KN/m2) a construir Cuarto de control Estación 3 (adicional a la carga de la Fase I);
para ello se da clic sobre el modelo nuevamente y en el recuadro que aparece se
activa la opción Distributed load.
Luego de insertar las fases se da clic en el icono select points for curves y apare
una ventana con el modelo en 3D en el cual en cada plano de trabajo (niveles 0.0,
-3.0, -3.6, -4.5, -6.75, -8.6 y -10.6 m) se seleccionan los puntos de interés a
analizar para las deformaciones y esfuerzos que calcula el software.
98
Figura 46. Ventana de selección de puntos de control en el nivel -3.00 m.
Para finalizar se da clic en el icono Calculate para que el software realice los
cálculos y se generen las curvas de deformaciones y esfuerzo.
Las deformaciones calculadas por este software son de ρ = 6.805 cm, así como se
muestra a continuación.
99
Figura 47. Gráfico de deformaciones elásticas generadas por PLAXIS 3D.
100
3.2.3 Deformaciones o asentamientos por consolidación unidimensional
H σ ' + ∆σ
ρ= CcLog 0 '
1 + e0 σ0
qc B1 + B2
(α 1 + α 2 ) − 1 (α 2 )
B
∆σ =
π B2 B2
Donde:
qc = γH
B1 + B2 B
α 1 (radianes) = tan −1 − tan −1 2
B2 z
B1
α 2 (radianes) = tan −1
z
101
Figura 49. Incremento de esfuerzo bajo un terraplén.
Teniendo en cuenta el mismo perfil del suelo con la carga del muro (67.50 KN/m2),
se tiene lo siguiente:
Tabla 28. Asentamientos por consolidación generados por la carga del muro.
DATOS GENERALES
qc (KN/m2) 67.5
B1 (m) 8.7
B2 (m) 6.0
102
Profundidad estrato Peso Esfue
Tip Presi
Esf. Esf. rzo
o de ón de
Total
Esp α1 α2 Δσ Cc e0 ρ
Sue Inici Prom Unita estrat Poro
efecti
Final eso o mitad vo
lo al edio Z rio u
rH σ'0
CL 0.00 0.90 0.45 0.90 12.0 10.80 5.40 0.00 5.40 0.02 1.52 33.76 0.30 1.25 0.10
CL 0.90 3.15 2.03 2.25 18.5 19.55 20.58 5.25 15.33 0.09 1.34 33.62 0.25 1.08 0.14
ML 3.15 5.00 4.08 1.85 19.0 17.00 38.85 25.75 13.10 0.17 1.13 33.10 0.25 0.93 0.13
ML 5.00 7.00 6.00 2.00 20.0 20.38 57.54 45.00 12.54 0.27 0.97 35.10 0.20 0.80 0.13
qc * B * L
∆σ =
(B + Z ) + (L + Z )
103
Figura 50. Método 2:1 para encontrar el incremento de esfuerzo debajo de una
cimentación.
Aplicando el mismo perfil del suelo con la carga estructural de 7.10 KN/m2, se
tiene lo siguiente:
DATOS GENERALES
qc (KN/m2) 7.1
B (m) 1.0
L (m) 1.0
104
Tip Profundidad estrato
Esfuerz
o Peso Esf. Esf. Presió
o
de Total n de
∆σ Cc ρ
Inici Promedi Espesor Unitar estrat
Final efectiv
Sue al oZ H io o mitad Poro u
o σ'0
lo
SU 3 KN/m 2 2 2 KN/m
m m m m KN/m 2 KN/m KN/m KN/m 2 (m)
CS
CL 0.00 0.90 0.45 0.90 12.00 10.80 5.40 0.00 5.40 3.38 0.30 0.03
CL 0.90 3.15 2.03 2.25 18.50 19.55 20.58 5.25 15.33 0.78 0.25 0.01
ML 3.15 5.00 4.08 1.85 19.00 17.00 38.85 25.75 13.10 0.28 0.25 0.00
ML 5.00 7.00 6.00 2.00 20.00 20.38 57.54 45.00 12.54 0.14 0.20 0.00
105
Tabla 30. Rango de la permeabilidad hidráulica parra varios suelos.
Tomando para las arcillas CL 1 y CL 2 valor de 8.64E-5 m/día y para los limos ML
1 y ML 2 valor de 8.64 E-3 m/día.
106
3.2.5 Deformaciones o asentamientos diferenciales
B
3.55
A
2.40 4.65
1 2 3
2.40 0.80 cm
4.65 1.55 cm
3.55 1.18 cm
107
4. ANÁLISIS COMPARATIVO
Referente a los resultados obtenidos por las diferentes teorías clásicas, se entra a
evaluar la capacidad que correspondería para la construcción del cuarto de control
de categoría baja. En la tabla 32 se pueden observar los diferentes valores
obtenidos.
500 462.38
CAPACIDAD DE CARGA (KN/m²)
Teniendo en cuenta que las capacidades obtenidas son altas debido al muro de
contención construido, y sabiendo que la carga factorizada a recibir por parte de la
estructura (cuarto de control) es baja del orden de 7.1 KN/m2, se elige la menor
capacidad que en este caso corresponde a la de Terzaghi. De igual manera
siempre se elige la menor capacidad para no exponer el suelo a un esfuerzo
mayor del que pueda llegar a resistir.
108
Figura 53. Comparativo de deformaciones elásticas.
8.000 7.436
6.805
7.000
Deformaciones elásticas ρ (cm)
6.000
5.000
5.000
4.000
3.000
2.000
1.000
0.000
Timoshenko y Goodier Plaxis 3D Práctica
Métodos
109
Figura 54. Comparativo de deformaciones por consolidación.
55.020
55.000
54.500
54.000
53.500 53.397
53.000
52.500
Consolidación unidimensional Plaxis 2D
Métodos
110
Figura 55. Comparativo de deformaciones totales.
62.000
61.825
61.800
Deformaciones totales ρ (cm)
61.600
61.400
61.200
61.000
60.833
60.800
60.600
60.400
60.200
TEORIA CLÁSICA ELEMENTOS FINITOS
Métodos
Dado que el software ofrece una variedad de alternativas para el análisis que
permiten el mejoramiento del modelo acercándolo más a la realidad, se plantea
calibrar nuevamente el modelo incluyendo el refuerzo del geotextil al muro de
contención en tierra armada. Todo esto con el fin de optimizar los diseños y
resultados.
111
Figura 56. Modelo de la malla deformada con refuerzo de geotextil.
Las deformaciones totales incluyendo el refuerzo del geotextil son de 45.08 cm.
112
Las fuerzas axiales que ejercen el geotextil son del orden de 76.62 KN/m; con ello
se va a las tablas de los catálogos y se selecciona un tipo de malla comercial que
cumpla con la especificación de dicha fuerza.
113
Figura 59. Deformaciones totales incluyendo el geotextil.
Lo cual indica que se pasa de tener una deformación de 55.02 cm a 44.89 cm,
minimizando las deformaciones y obteniendo ahora las siguientes deformaciones
totales con una diferencia porcentual del 18.5 % que representa 8.31 cm. La
diferencia entre los dos métodos cambia al 16 %. Ver figura 60.
Finalmente para la deformación total se obtiene una diferencia entre los dos
métodos del 15 % lo cual corresponde a 9.14 cm. Ver figura 61.
114
Figura 60. Comparativo de deformaciones por consolidación teniendo en cuenta el
geotextil en PLAXIS 2D.
56.000
Deformaciones por consolidación ρ (cm)
54.000 53.397
52.000
50.000
48.000
46.000 44.890
44.000
42.000
40.000
Consolidación unidimensional Plaxis 2D
Métodos
62.000 60.833
60.000
Deformaciones totales ρ (cm)
58.000
56.000
54.000
51.695
52.000
50.000
48.000
46.000
TEORIA CLÁSICA ELEMENTOS FINITOS
Métodos
115
5. CONCLUSIONES
Se puede apreciar que las diferencias entre los dos métodos no exceden de
consolidación inicialmente no exceden del 1.63 %, pero la diferencia se encuentra
cuando por medio del método de elementos finitos se calibran los modelos
incorporando y probando posibles alternativas constructivas lo cual adicionan
nuevos parámetros y que permiten aproximar los resultados cada vez más a la
realidad, lo cual indica que este tipo de herramienta sabiéndose usar e interpretar
facilita el análisis y por ende va generado un avance en la geotecnia. Lo cual a la
final marcó una diferencia en un 15 % (9.14 cm).
116
BIBLIOGRAFIA
BOWLES J.E. (1988). Foundation Analysis and Design. Mc Graw Hill Book Co.
117
ANEXOS
118