PDF

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 118

ANÁLISIS COMPARATIVO DE LAS DEFORMACIONES DE LOS SUELOS

MEDIANTE TEORÍAS TRADICIONALES Y SOFTWARE DE ELEMENTOS


FINITOS

JENNIFER ANDREA REMOLINA MANOSALVA

UNIVERSIDAD INDUSTRIAL DE SANTANDER

FACULTAD DE INGENIERÍAS FISICOMECÁNICAS

ESCUELA DE INGENIERÍA CIVIL

MAESTRIA EN GEOTECNIA

BUCARAMANGA

2016
ANÁLISIS COMPARATIVO DE LAS DEFORMACIONES DE LOS SUELOS
MEDIANTE TEORÍAS TRADICIONALES Y SOFTWARE DE ELEMENTOS
FINITOS

JENNIFER ANDREA REMOLINA MANOSALVA

Trabajo de aplicación modalidad profundización para optar al título de

Magíster en Geotecnia

Directora

Ing. Msc. Luz Marina Torrado Gómez

UNIVERSIDAD INDUSTRIAL DE SANTANDER

FACULTAD DE INGENIERÍAS FISICOMECÁNICAS

ESCUELA DE INGENIERÍA CIVIL

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

Yo, JENNIFER ANDREA REMOLINA MANOSALVA, mayor de edad, vecino de


Bucaramanga, identificado con la Cédula de Ciudadanía No. 63.560.799 de
Bucaramanga, actuando en nombre propio, en mi calidad de autor del trabajo de grado,
del trabajo de profundización, o de la tesis denominada: ANÁLISIS COMPARATIVO DE
LAS DEFORMACIONES DE LOS SUELOS MEDIANTE TEORÍAS TRADICIONALES Y
SOFTWARE DE ELEMENTOS FINITOS, hago entrega del ejemplar respectivo y de sus
anexos de ser el caso, en formato digital o electrónico (CD o DVD) y autorizo a LA
UNIVERSIDAD INDUSTRIAL DE SANTANDER, para que en los términos establecidos en
la Ley 23 de 1982, Ley 44 de 1993, decisión Andina 351 de 1993, Decreto 460 de 1995 y
demás normas generales sobre la materia, utilice y use en todas sus formas, los derechos
patrimoniales de reproducción, comunicación pública, transformación y distribución
(alquiler, préstamo público e importación) que me corresponden como creador de la obra
objeto del presente documento. PARÁGRAFO: La presente autorización se hace
extensiva no sólo a las facultades y derechos de uso sobre la obra en formato o soporte
material, sino también para formato virtual, electrónico, digital, óptico, uso en red,
Internet, extranet, intranet, etc., y en general para cualquier formato conocido o por
conocer.
EL AUTOR – ESTUDIANTE, manifiesta que la obra objeto de la presente autorización es
original y la realizó sin violar o usurpar derechos de autor de terceros, por lo tanto la obra
es de su exclusiva autoría y detenta la titularidad sobre la misma. PARÁGRAFO: En caso
de presentarse cualquier reclamación o acción por parte de un tercero en cuanto a los
derechos de autor sobre la obra en cuestión, EL AUTOR / ESTUDIANTE, asumirá toda la
responsabilidad, y saldrá en defensa de los derechos aquí autorizados; para todos los
efectos la Universidad actúa como un tercero de buena fe.

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:

Jennifer Andrea Remolina Manosalva.

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 hermana por su voz de aliento y cariño en los momentos alegres y difíciles.

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.

Agradezco sinceramente a todos los docentes que estuvieron presentes a lo largo


de esta maestría, por su tiempo, enseñanzas y por compartir sus conocimientos y
experiencias desde el ámbito ingenieril y personal.

A la Escuela de Ingeniería Civíl de la Universidad Industrial de Santander, por ser


haber sido un ente facilitador de herramientas y procesos académicos.

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

Figura 1. Procedimiento para un estudio de suelos en fundaciones superficiales.


............................................................................................................................... 19
Figura 2. Clasificación de las unidades de construcción por categorías. .............. 19
Figura 3. Número mínimo de sondeos y profundidad por cada unidad de
construcción. .......................................................................................................... 20
Figura 4. Curva de capacidad de carga. ............................................................... 21
Figura 5. Superficies de ruptura. ........................................................................... 24
Figura 6. Solución de Skempton arcillas no drenadas (φ = 0). ............................. 36
Figura 7. Efecto del nivel freático Caso I. .............................................................. 37
Figura 8. Efecto del nivel freático Caso II. ............................................................. 37
Figura 9. Muestreador de tubo partido. ................................................................. 39
Figura 10. Montaje del equipo de perforación6...................................................... 40
Figura 11. Deformación de un cimiento rígido o flexible. ...................................... 44
Figura 12. Asentamiento inmediato o elástico3. .................................................... 45
Figura 13. Asentamiento por consolidación primaria y secundaria3. ..................... 45
Figura 14. Asentamiento inmediato en arcillas en un espacio semi-infinito. ......... 47
Figura 15. Asentamiento inmediato en arcillas en una capa limitada por base
rígida. ..................................................................................................................... 49
Figura 16. Asentamiento inmediato en dos o más capas limitadas por una base
rígida. ..................................................................................................................... 50
Figura 17. Asentamiento por consolidación. ......................................................... 52
Figura 18. Variación volumétrica en el suelo. ....................................................... 52
Figura 19. Curva de consolidación. ....................................................................... 53
Figura 20. Curva de consolidación para una arcilla normalmente consolidada. ... 54
Figura 21. Asentamiento para una arcilla normalmente consolidada para un
número arbitrario de capas. ................................................................................... 56
Figura 22. Curva de consolidación para una arcilla pre-consolidada caso ∆σ ˂ σ’p
– σ’0 ....................................................................................................................... 57
Figura 23. Curva de consolidación para una arcilla pre-consolidada caso ∆σ > σ’p
– σ’0 ....................................................................................................................... 58
Figura 24. Asentamiento por consolidación para arcillas pre-consolidada con un
número arbitrario de capas. ................................................................................... 59
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
Figura 29. Ventana de estratigrafía. ...................................................................... 67

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

Tabla 1. Condiciones de forma según Terzaghi. ................................................... 25


Tabla 2. Factores de capacidad de carga para las ecuaciones de Terzaghi. ........ 25
Tabla 3. Factores de capacidad de carga para las ecuaciones de Meyerhof,
Hansen y Vesic. ..................................................................................................... 27
Tabla 4. Factores de forma, profundidad e inclinación de las ecuaciones de
capacidad de carga de Meyerhof (1951). .............................................................. 28
Tabla 5. Factores de profundidad de las ecuaciones de capacidad de carga de
Hansen. ................................................................................................................. 30
Tabla 6. Factores de forma y profundidad de las ecuaciones de capacidad de
carga de Hansen.................................................................................................... 30
Tabla 7. Factores de inclinación de las ecuaciones de capacidad de carga de
Hansen y Vesic. ..................................................................................................... 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. .................................. 32
Tabla 9. Factores de forma de la ecuación de capacidad de carga de Brinch –
Hansen. ................................................................................................................. 34
Tabla 10. Factores de inclinación de la ecuación de capacidad de carga de Brinch
– Hansen. .............................................................................................................. 35
Tabla 11. Resumen de relaciones de energía8...................................................... 42
Tabla 12. Correlaciones para estimar el ángulo de fricción a partir del ensayo de
SPT8....................................................................................................................... 43
Tabla 13. Correlaciones empíricas para estimar el módulo de Young o Elasticidad
(E) a partir del ensayo de SPT8. ............................................................................ 43
Siendo φ el ángulo de fricción................................................................................ 48
Tabla 14. Factores de forma y rigidez en un espacio semi-infinito. ....................... 48
Tabla 15. Factores de forma y rigidez en una capa limitada por base rígida. ....... 49
Tabla 16. Valores máximos de deformaciones o asentamientos diferenciales. .... 61
Tabla 17. Parámetros para cálculo de la capacidad portante. ............................... 73
Tabla 18. Cargas de proyecto actuantes en cimentación. ..................................... 73
Tabla 19. Parámetros del suelo para el cálculo de asentamientos elásticos. ........ 80
Tabla 20. Deformaciones elásticas en los diferentes estratos generados por el
muro. ...................................................................................................................... 81
Tabla 21. Deformaciones elásticas en los diferentes estratos generados por la
carga estructural. ................................................................................................... 83
Tabla 22. Resumen de parámetros para el material de relleno. ............................ 88
Tabla 23. Resumen de parámetros para el material de base granular. ................. 89
Tabla 24. Resumen de parámetros para el material arcilla 1. ............................... 90

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*

AUTOR: JENNIFER ANDREA REMOLINA MANOSALVA**

PALABRAS CLAVES: Deformaciones o asentamientos, capacidad portante,


propiedades geomecánicas, Software de elementos finitos, compresibilidad,
plasticidad, consolidación.

RESUMEN:

En este trabajo de aplicación se realiza un análisis comparativo de las


deformaciones mediante teorías tradicionales y software de elementos finitos.
Para la primera alternativa se tienen las teorías clásicas como las de Terzaghi,
Meyerhof, Vesic, Hansen y Brinch – Hansen, para la segunda alternativa el uso de
software de elementos finitos por medio de Plaxis 2D y 3D.

Además de las deformaciones y la capacidad de carga a la que está sometido el


suelo, y así poder comparar estos dos métodos para saber qué tan desfasados se
encuentra el uno del otro.

Este análisis se realizará para un caso en específico: ‘’ESTACIÓN 3’’, que


consiste en la construcción de un muro en tierra armada, el cual presenta suelos
arcillosos de compresibilidad media, plasticidad baja y presencia de nivel freático.

Para la construcción de este muro en tierra armada se hizo control y registro de


las deformaciones inmediatas por medio de topografía.

Con el presente proyecto se busca comparar las deformaciones obtenidas


mediante las ecuaciones teóricas y elementos finitos contra las deformaciones
producidas por este tipo de suelos en la realidad.

Finalmente, se realizó el análisis comparativo de las teorías y que diferencias


cuantitativas se pueden obtener, dado que lo calculado teóricamente siempre es
un aproximado y no un valor exacto, sobre todo cuando se tienen condiciones de
consolidación.

__________________
*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*

AUTHOR: JENNIFER ANDREA REMOLINA MANOSALVA

KEYWORDS: Deformations or settlements, bearing capacity, geochemical


properties, finite element software, compressibility, plasticity, consolidation.

ABSTRACT:

In this application work, a comparative analysis of deformations using traditional


theories and finite element software is done. For the first alternative they have
classical theories such as Terzaghi, Meyerhof, Vesic, Hansen and Brinch –
Hansen, for the second alternative the use of finite element software by Plaxis 2D y
3D.

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.

Principalmente la función de la cimentación es asegurar una distribución uniforme


de las cargas y fuerzas en toda su área de trabajo, para evitar futuros
asentamientos diferenciales.

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.

Para este proyecto, y en el tipo de estructura y suelo que se estudiará, se tendrán


en cuenta ciertas condiciones y parámetros debido a la geología típica de la zona,
como lo son suelos arcillosos, de compresibilidad media, plasticidad baja y registro
de nivel freático superficial.

Para el análisis comparativo se realizaran cálculos de capacidad portante del


suelo y deformaciones para los cuales nos apoyaremos por una parte con las
teorías clásicas existentes que son producto de resultados experimentales y
analíticos; y por otra parte con el uso del software de elementos finitos que para
uso de geotecnia se hará a través de Plaxis 2D y 3D, el cual permitirá simular la
estructura. Los resultados que se obtengan mediante estos dos métodos se
confrontaran y se analizaran sus diferencias o que tan desfasados está el uno del
otro, de igual manera sus ventajas y desventajas con respecto a los datos que se
registran en la realidad con relación a las deformaciones elásticas.

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

1.1 OBJETIVO GENERAL:

Realizar un análisis comparativo de deformaciones de suelos en vivienda de


categoría baja mediante diferentes teorías tradicionales y software de elementos
finitos.

1.2 OBJETIVOS ESPECÍFICOS:

• Determinar las propiedades geomecánicas del suelo.

• Determinar la capacidad de carga a la cual está sometido un suelo en


vivienda de categoría baja.

• Determinar los esfuerzos y asentamientos del suelo encontrado por


métodos empíricos y mediante software de elementos finitos.

• Comparación y conclusiones de los resultados – Análisis caso Estación 3.

18
2. MARCO TEÓRICO

2.1 FUNDACIONES SUPERFICIALES

Un análisis de cimentación supone las siguientes fases:

• Debe empezar por conocer el tipo y función de la estructura a cimentar. Por


lo tanto debe tenerse información sobre las cargas o hacer un cálculo aproximado
de ellas.
• Debe conocerse el perfil del suelo y sus características geotécnicas, es
decir, se debe realizar un estudio de suelo.

Figura 1. Procedimiento para un estudio de suelos en fundaciones superficiales.

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

2.1.1 Número y profundidad de las perforaciones

Figura 2. Clasificación de las unidades de construcción por categorías.

Fuente: Norma Sismo Resistente NSR-10.


19
Figura 3. Número mínimo de sondeos y profundidad por cada unidad de
construcción.

Fuente: Norma Sismo Resistente NSR-10

2.2 CAPACIDAD DE CARGA (CAPACIDAD PORTANTE / CAPACIDAD


DE SOPORTE)

La carga última de una cimentación superficial se puede definir como el valor


máximo de la carga con el cual en ningún punto del subsuelo se alcanza la
condición de rotura, o también como aquel valor de carga, mayor que el anterior,
con el cual el fenómeno de rotura se extiende a un amplio volumen del suelo2.

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
δ

Fuente: Norma Sismo Resistente NSR-10.

2.2.1 Factor de seguridad

• Capacidad portante admisible

qult
qadm =
F.S.

• Capacidad portante última neta

qneta = qult − q q = γD f

• Capacidad admisible neta

qult − q
qadm( neta) =
F.S.

21
2.3 PARÁMETROS GEOMECÁNICOS PARA EL CÁLCULO DE LA
CAPACIDAD DE CARGA

Para el cálculo de las capacidades de carga se requieren de ciertos parámetros


geomecánicos tales como peso específico seco, peso saturado, cohesión y angulo
de fricción.

Para el cálculos de estos parámetros geomecánicos se parte asumiendo valores


estimados para el peso específico seco, gravedad específica y cohesión
propuestos por Grundbau Taschenbuch (1980) Tomado de Oteo (1989).

A partir de dichos parámetros y de la clasificación de suelos realizada a las


muestras obtenidas mediante el ensayo SPT, se obtienen los demás parámetros
requeridos.

2.3.1 Relación de vacíos (e):

 G *γ 
e =  s w  − 1
 γ unsat 

2.3.2 Peso saturado (үsat):

γ w * ( e + Gs )
γ sat =
1+ e

22
2.3.3 Ángulo de fricción (φ):

2.4 ECUACIONES DE CAPACIDAD DE CARGA (CAPACIDAD


PORTANTE)

2.4.1 Terzaghi (1943)

Terzaghi pasa del problema plano al problema espacial introduciendo algunos


factores de forma.

Una sucesiva contribución sobre el efectivo comportamiento del terreno ha sido


aportada por Terzaghi.

Terzaghi admite este comportamiento en los terrenos muy compactos.

En éstos, de hecho, la curva cargas-asentamientos presenta un primer tracto


rectilíneo, seguido por un breve tracto curvilíneo (comportamiento elástico-
plástico); la rotura es instantánea y el valor de la carga límite resulta claramente
individuado (rotura general).

En un terreno muy suelto en cambio la relación cargas-asentamientos presenta un


tramo curvilíneo acentuado desde las cargas más bajas por efecto de una rotura
progresiva del terreno (rotura local). Como consecuencia la individualización de la
carga límite no es tan claro y evidente como en el caso de los terrenos
compactos3. Ver figura 5.

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:

φred = Ángulo de fricción reducido.

Cred = Cohesión reducida.

Figura 5. Superficies de ruptura.

Fuente: Tomado y adoptado de DAS Braja M. (1999). Fundamentos de Ingeniería


geotécnica.

La fórmula de Terzaghi se puede escribir así:

qult = cNc Sc + γDNq + 0.5γBNγ Sγ

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 

Tabla 1. Condiciones de forma según Terzaghi.

Largo Redondo Cuadrado

Sc = 1.0 1.3 1.3

Sү = 1.0 0.6 0.8

Tabla 2. Factores de capacidad de carga para las ecuaciones de Terzaghi.

Φ deg Nc Nq Kpγ

10.8
0 5.7 1.0 0.0

5 7.3 1.6 0.5 12.9

14.7
10 9.6 2.7 1.2

15 12.9 4.4 2.5 18.6

25.0
20 17.7 7.4 5.0

35.0
25 25.1 12.7 9.7

30 37.2 22.5 19.7 52.0

25
35 57.8 41.6 42.4 82.0

141.0
40 95.7 81.3 100.4

45 172.3 173.3 297.5 298.0

800.0
50 347.5 415.1 1153.2

2.4.2 Meyerhof (1963)

Meyerhof propuso una fórmula para calcular la carga última parecida a la de


Terzaghi. Las diferencias consisten en la introducción de nuevos coeficientes de
forma3.

Introdujo un coeficiente sq que multiplica el factor Nq, factores de profundidad di y


de pendencia ii para el caso en que la carga trasmitida a la cimentación sea
inclinada en la vertical.

• Carga Vertical

qult = cNc Sc d c + γDNq S q d q + 0.5γBNγ Sγ dγ

• Carga Inclinada

qult = cNc d c ic + γDNq d q iq + 0.5γBNγ dγ iγ

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

0 5.14 1.0 0.0 0.0 0.0 0.195 0.000

5 6.49 1.6 0.1 0.1 0.4 0.242 0.146

10 8.34 2.5 0.4 0.4 1.2 0.296 0.241

15 10.97 3.9 1.1 1.2 2.6 0.359 0.294

20 14.83 6.4 2.9 2.9 5.4 0.432 0.315

25 20.71 10.7 6.8 6.8 10.9 0.515 0.311

26 22.25 11.9 8.0 7.9 12.5 0.533 0.308

28 25.79 14.7 11.2 10.9 16.7 0.571 0.299

30 30.13 18.4 15.7 15.1 22.4 0.611 0.289

32 35.47 23.2 22.0 20.8 30.2 0.653 0.276

34 42.14 29.4 31.1 28.8 41.0 0.699 0.262

36 50.55 37.8 44.4 40.1 56.2 0.747 0.247

38 61.31 48.9 64.1 56.2 77.9 0.798 0.231

40 75.25 64.2 93.7 79.5 109.3 0.853 0.214

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

Tabla 4. Factores de forma, profundidad e inclinación de las ecuaciones de


capacidad de carga de Meyerhof (1951).

FACTORES VALOR PARA

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)

Es una extensión de la fórmula de Meyerhof; las extensiones consisten en la


introducción de bi que tiene en cuenta la eventual inclinación en la horizontal del
nivel de cimentación y un factor gi para terreno en pendiente.

La fórmula de Hansen vale para cualquier relación D/B, ya sean cimentaciones


superficiales o profundas; sin embargo el mismo autor introdujo algunos
coeficientes para poder interpretar mejor el comportamiento real de la cimentación;
sin éstos, de hecho, se tendría un aumento demasiado fuerte de la carga última
con la profundidad3.

La fórmula de Terzaghi se puede escribir así:

qult = cNc Sc d c ic g c bc + γDNq S q d q iq g q bq + 0.5γBNγ Sγ dγ iγ gγ bγ

Cuando φ = 0
qult = 5.14Cu(1 + S 'c +d 'c −ic − bc − g c ) + γD
----------------------------------------------------------------------------------

D/B 0 1 1.1 2 5 10 20 100

----------------------------------------------------------------------------------

d'c 0 0.40 0.33 0.44 0.55 0.59 0.61 0.62

----------------------------------------------------------------------------------

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

Para valores de D/B ˂ 1.0 Para valores de D/B > 1.0

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.

Tabla 6. Factores de forma y profundidad de las ecuaciones de capacidad de


carga de Hansen.

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.

FACTORES DE INCLINACIÓN (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.

FACTORES DE BASE (BASE


INCLINACIÓN DE LA SUPERFICIE
INCLINADA)

β° η° η°
g c' = bc' = bc = 1 −
147° 147° 147°

Para Vesic usa


bq ( H ) = e −2η tan φ
Nγ =−2senβ __ para _φ = 0 bγ ( H ) = e −2.7η tan φ

β°
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

2.4.4 Vesic (1975)

La fórmula de Vesic es análoga a la fórmula de Hansen, con Nq y Nc como en la


fórmula de Meyerhof y Nγ como se indica a continuación:

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.

2.4.5 Brinch – hansen

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:

Vd es la carga de proyecto al ELU, normal en la base de la cimentación, que


incluye también el peso de la cimentación misma; mientras

Rd es la carga última de proyecto de la cimentación con respecto a las cargas


normales, teniendo en cuenta también el efecto de cargas inclinadas o
excéntricas.

En la evaluación analítica de la carga última de proyecto Rd se deben considerar


las situaciones a corto y a largo plazo en los terrenos de grano fino.

La carga última de proyecto en condiciones no drenadas se calcula como:

R
= ( 2 + π ) c u s c ic + q
A'
Donde:

A’ = B’ L’ Es el área de la cimentación efectiva de proyecto, entendida, en caso


de carga excéntrica, como el área reducida en cuyo centro se aplica la resultante
de la carga3.

cu Es la cohesión no drenada.

q Es la presión litostática total en el plano de cimentación.

sc = 1 + 0,2 (B’/L’) Factor de forma para cimentaciones rectangulares

sc = 1,2 Factor de forma para cimentaciones cuadradas o circulares.

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.

En condiciones drenadas la carga última de proyecto se calcula como sigue3:

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

Tabla 9. Factores de forma de la ecuación de capacidad de carga de Brinch –


Hansen.

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

i N −1 De forma cuadrada o circular


q q
i =
c N −1 1− H
3
q i =
 
γ V + A' c ' cot φ ' 

Además de los factores correctivos de la tabla anterior, se consideran los


complementarios de la profundidad del plano de cimentación y de la inclinación y
de la inclinación del plano de cimentación y del plano terreno (Hansen)3.

2.4.6 Capacidad de carga en condiciones no drenadas – (arcilla φ = 0)

Arcillas saturadas presentan parámetros de resistencia no drenadas Cu y φ = 0,


determinados en un ensayo triaxial UU.

Nc= 5.7 Terzaghi


Cuando Ф=0 Nc= 5,14 Meyerhoff y Hansen

35
Las ecuaciones quedan:

Terzaghi qult= 5.7CuSc+γDf


Meyerhof qult = 5.14CuScdc+γDf
Hansen q = 5.14Cu(1+S’c+d’c-i’c-g’c-b’c)+γDf
ult

Figura 6. Solución de Skempton arcillas no drenadas (φ = 0).

Fuente: DAS Braja M. (1999). Fundamentos de Ingeniería geotécnica.

qult = cu * Nc * Sc * γD f

36
2.4.6.1 Efecto de la posición del nivel freático

2.4.6.1.1 Caso I: Nivel freático por encima de la profundidad de cimentación.

Figura 7. Efecto del nivel freático Caso I.

0 ≤ D1 ≤ D f

Fuente: Material de clase de Geotecnia II UIS. Wilfredo del Toro. (2015).

El factor de sobrecarga queda: q = γD1 + γ ' D2

El factor γ en el último término queda: γ ' = γ sat − γ w

2.4.6.1.2 Caso II: Nivel freático por debajo de la profundidad de cimentación y d


≤B

Figura 8. Efecto del nivel freático Caso II.

0 ≤ D1 ≤ D f

Fuente: Material de clase de Geotecnia II UIS. Wilfredo del Toro. (2015).

El factor de sobrecarga queda:


q = γD f

37
El factor γ en el último término queda:
γ = γ '+
d
(γ − γ ')
B

2.5 CAPACIDAD DE CARGA A PARTIR DE ENSAYOS DE CAMPO

La capacidad de carga de los suelos puede ser determinada a partir de ensayos


de campo, cuyos datos o resultados están relacionados directa e indirectamente
con la capacidad de soporte.

Ensayos de campo comúnmente utilizados son:

• Ensayo de penetración estándar (SPT)


• Ensayo de penetración de cono estático o cono holandés (CPT)
• Ensayo de veleta
• Pruebas de carga
• Prueba de presurímetro
• Prueba de Dilatómetro

En este caso se utilizó el ensayo de penetración estándar (SPT).

2.5.1 Ensayo de penetración estándar (spt)

El ensayo de penetración estándar SPT (Standard Penetration Teste) es


ampliamente utilizado en exploración geotécnica dado que existe una gran
variedad de correlaciones entre el número de golpes del que se obtienen tras la
ejecución del ensayo y algunos parámetros geotécnicos de gran importancia en el
diseño de cimentaciones.

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.

El ensayo consiste en hincar un muestreador de tubo partido el cual debe cumplir


con unas dimensiones establecidas (Ver Figura 9), mediante la caída de un
martillo de 140 lb (63.5 Kg) desde una altura normalizada de 30’’ (76 cm)3, el
montaje del sistema se puede observar en la Figura 10.

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 ensayo S.P.T. (Standard Penetration Teste) consiste básicamente encontrar el


número de golpes (N) que se necesitan para introducir dentro un estrato de suelo,
un toma-muestras (cuchara partida) de 30 cm de largo, a diferentes profundidades
(generalmente con variación de 0.90 m en 0.90 m). El toma-muestras es golpeado
bajo energía constante, con una pesa en caída libre de 140 lb y una altura de
caída libre de 76.2 cm.

Este ensayo se realiza en depósitos de suelo arenoso y de arcilla blanda; no es


recomendable llevarlo a cabo en depósitos de grava, roca o arcilla consolidada,
debido a los daños que podría sufrir el equipo de perforación al introducirlo dentro
de dichos estratos.

El valor que se define como el N de campo es la suma de los golpes necesarios


para avanzar 12’’ (30 cm) de perforación y se expresa en golpes/pie, este valor
corresponde a la suma de los dos últimos valores reportados en el tramo de 18’’
(45 cm)4.

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:

• No se han alcanzado 50 golpes en alguno de los incrementos de 6’’


• Cuando se han acumulado 100 golpes
• Cuando no se observa avance del muestreador después de aplicar 10
77
golpes .

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:

Ncorr = valor del N campo corregido

CN = factor de corrección dado por la sobrecarga efectiva del suelo

η1 = factor de corrección por energía del martillo

η2 = factor de corrección por longitud de la varilla

η3 = factor de corrección por resistencia interna del toma muestras

η4 = factor de corrección por diámetro de la perforación

De acuerdo con la literatura de González (1999) relacionada con el tema, se


considera que los factores de η2, η3 y η4 equivalen a 1.0 y solo se tendrá en
cuenta los factores de CN y η1.

2.5.1.1 Factor de corrección por energia

Este factor representa el rendimiento del impacto del martillo sobre el cabezal de
golpeo, se define como4:

N1 = Er/E60

Donde:

E60 es el valor de la energía de referencia base que usualmente es 60%

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.

En la Tabla 11 se observan los factores de corrección donde se observa que Er


varía de acuerdo con los equipos y su uso en diferentes países, en el caso de
Colombia usualmente se trabaja con una Er de 45%.

Tabla 11. Resumen de relaciones de energía88.

2.5.1.2 Factor de corrección por confinamiento

Existen diferentes propuestas para el cálculo de este factor de corrección, sin


embargo la siguiente expresión propuesta por Liao y Whitman (1986)8, es la que
se utiliza comúnmente, representa aproximadamente el promedio de dichas
expresiones y es la más fácil de utilizar.

 95.76 
CN =  
 σ 'vo 

Donde σ’vo (kPa) es el esfuerzo efectivo vertical a la profundidad donde se desea


obtener el valor del N corregido. En general se recomienda9 CN ˂ 2.0 por lo que
para efectos del presente artículo cuando se obtengan valores superiores a 2.0 en
este factor, se tomará este último valor con el fin de no incurrir en sobre
estimaciones en los parámetros que se van a analizar.

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

A continuación se presentan las siguientes correlaciones para estimar los ángulos


de fricción y módulos de Young:

Tabla 12. Correlaciones para estimar el ángulo de fricción a partir del ensayo de
SPT8.

Tabla 13. Correlaciones empíricas para estimar el módulo de Young o Elasticidad


(E) a partir del ensayo de SPT8.

2.6 DEFORMACIONES EN LOS SUELOS

La deformación vertical de los suelos ocurre debido a las cargas aplicadas


externamente o en consecuencia al propio peso del suelo; los asentamientos o

43
deformaciones pueden clasificarse en asentamientos inmediatos, por
consolidación o por compresión secundaria10.

Figura 11. Deformación de un cimiento rígido o flexible.

Fuente: DAS Braja M. (1999). Fundamentos de Ingeniería geotécnica.

La aplicación de una carga de dimensiones finitas sobre un terreno cohesivo


genera una serie de fenómenos, los cuales se pueden ilustrar como en la figura
12.

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.

Figura 13. Asentamiento por consolidación primaria y secundaria3.

• Durante la fase de carga, en el terreno se producen subpresiones del agua


intersticial ∆u, y dada la baja permeabilidad del terreno, se puede asumir que, en
el ámbito de las usuales velocidades de aplicación de la carga, nos encontremos
en condiciones no drenadas. El estrato de arcilla se deforma a volumen casi
constante y el consiguiente asiento se indica como deformación elástica.

• La instauración del drenaje, con la transferencia progresiva de la carga de la


base fluida a la armazón sólida, comporta ulteriores asientos cuya velocidad en el
tiempo depende principalmente de las condiciones de drenaje.

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.

• Por último, aun cuando las sobrepresiones intersticiales se disipan (∆u=0),


continúan a darse en el tiempo asientos debidos a deformaciones viscosas en
condiciones drenadas y este asiento se conoce como deformación secundaria.

El comportamiento de los terrenos de grano grueso (arenas y gravas) resulta


macroscópicamente diferente a los de grano fino (limos y arcillas) debido a la
diferente permeabilidad.

En este caso se habla de asiento inicial y secundario para indicar las


deformaciones de naturaleza elasto-plástica y los debidos a fenómenos de
naturaleza viscosa3.

2.6.1 Deformaciones elásticas o inmediatas

Es una deformación ocasionada de forma casi instantánea al aplicarse cargas


externas, ésta deformación se debe al cambio en la forma y no al cambio de
volumen. Las deformaciones elásticas dependen de las propiedades elásticas de
los suelos a deformaciones pequeñas y de la rigidez y extensión del cimiento. La
relación esfuerzo-deformación en los suelos cohesivos saturados tiene un carácter
aproximadamente elástico lineal, en los suelos granulares esta relación obedece a
un comportamiento más complejo, posiblemente elasto-plástico o plástico. Estas
deformaciones se calculan de las ecuaciones derivadas de la teoría de
elasticidad10.

Para el cálculo de las deformaciones elásticas se hizo uso de las siguientes


teorías:

2.6.1.1 En arcillas

Los asientos de una cimentación rectangular de dimensiones B'xL' ubicada en la


superficie de un semiespacio elástico, se pueden calcular por medio de una
ecuación basada en la teoría de la elasticidad (Timoshenko e Goodier (1951))3:

46
2.6.1.1.1 Espacio semi-infinito

Figura 14. Asentamiento inmediato en arcillas en un espacio semi-infinito.

Fuente: Tomado y adoptado de DAS Braja M. (1999). Fundamentos de Ingeniería


geotécnica.

(1 − υ 2 )
ρ = qc * B Cd
Es

Siendo:

qc: Presión de contacto neta entre el cimiento y el suelo.

B: Ancho del cimiento

ν: Coeficiente de Poisson

E: Módulo de elasticidad

Cd: Coeficiente o factor de forma que depende de la geometría del problema,


rigidez del plano cargado y forma de la carga.

A continuación se muestra la tabla para el factor de forma y rigidez, Cd para el


cálculo de asentamiento de puntos bajo áreas cargadas en un espacio semi-
infinito.

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:

Coeficiente de Poisson (ν):

K0 K 0 = 1 − senφ 1 − senφ
ν= µ=
1 + K0 2 − senφ

Siendo φ el ángulo de fricción.

Tabla 14. Factores de forma y rigidez en un espacio semi-infinito.

MITAD DEL MITAD DEL


FORMA CENTRO ESQUINA LADO LADO PROMEDIO
CORTO LARGO

Circular flexible 1.00 0.64 0.64 0.64 0.85

circular rígido 0.79 0.79 0.79 0.79 0.79

cuadrado 1.12 0.56 0.76 0.76 0.95

cuadrado rígido 0.99 0.99 0.99 0.99 0.99

RECTANGULAR

Largo/Ancho

1.5 1.36 0.67 0.89 0.97 1.15

2 1.52 0.76 0.98 1.12 1.30

3 1.78 0.88 1.11 1.35 1.52

5 2.10 1.05 1.27 1.68 1.83

10 2.53 1.26 1.49 2.12 2.29

100 4.00 2.00 2.20 3.60 3.70

1000 5.47 2.75 2.94 5.03 5.15

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.

Fuente: Tomado y adoptado de DAS Braja M. (1999). Fundamentos de Ingeniería


geotécnica.

qc * B * (1 −υ 2 )
ρ= C' d
E

Siendo:

C’d: Tiene el mismo sentido de Cd y depende además del espesor de la capa H.

En la siguiente tabla se muestran los valores de C’d para diferentes geometrías.

Tabla 15. Factores de forma y rigidez en una capa limitada por base rígida.

FORMA DEL CIMIENTO

CIRCULAR RECTANGULAR

H/B diam=B L/B=1 L/B=1.5 L/B=2 L/B=3 L/B=5 L/B=10

0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00

49
0.10 0.09 0.09 0.09 0.09 0.09 0.09 0.09

0.25 0.24 0.24 0.23 0.23 0.23 0.23 0.23

0.50 0.48 0.48 0.47 0.47 0.47 0.47 0.47

1.00 0.70 0.75 0.81 0.83 0.83 0.83 0.83

1.50 0.80 0.86 0.97 1.03 1.07 1.08 1.08

2.50 0.88 0.97 1.12 1.22 1.33 1.39 1.40

3.50 0.91 1.01 1.19 1.31 1.45 1.56 1.59

5.00 0.94 1.05 1.24 1.38 1.55 1.72 1.82

∞ 1.00 1.12 1.36 1.52 1.78 2.10 2.53

2.6.1.1.3 Dos o más capas limitadas por una base rígida

Figura 16. Asentamiento inmediato en dos o más capas limitadas por una base
rígida.

Fuente: Tomado y adoptado de DAS Braja M. (1999). Fundamentos de Ingeniería


geotécnica.

qc * B * (1 −υ 2 )
ρ= C'd
E
50
Siendo:

C’d: Tiene el mismo sentido de Cd y depende además del espesor de la capa H.

Teniendo en cuenta las siguientes consideraciones:

• El asentamiento inmediato total es:

ρT = ρ1 + ρ2

• El asentamiento de la capa 1 es calculado como es el caso 2.6.1.1.2,


considerando que la capa rígida se encuentra a una profundidad H1.
• El asentamiento de la capa 2 es determinado así:

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.

2.6.2 Asentamientos por consolidación

Se dá a lo largo del tiempo, y ocurre en suelos arcillosos saturados cuando son


sometidos a una carga creciente causada por la construcción de una cimentación.

51
Figura 17. Asentamiento por consolidación.

Fuente: Tomado y adoptado de DAS Braja M. (1999). Fundamentos de Ingeniería


geotécnica.

Así como en las consideraciones del suelo en sus fases, la variación volumétrica
es:

Figura 18. Variación volumétrica en el suelo.

Fuente: Tomado y adoptado de DAS Braja M. (1999). Fundamentos de Ingeniería


geotécnica.
52
Es decir que la deformación total por compresión se puede calcular por la
expresión:

∆e
∆dV = dz
1+ e

∆e
ρ= H
1+ e

En la cual ∆e es la variación de la relación de vacíos por los esfuerzos colocados,


e es la relación de vacíos donde inicia compresión y H es el espesor de la capa
que se comprime.

A partir de la curva de consolidación, la ecuación se aplica de la siguiente manera:

Figura 19. Curva de consolidación.

Fuente: Tomado y adoptado de DAS Braja M. (1999). Fundamentos de Ingeniería


geotécnica.

53
2.6.2.1 Asentamiento en arcillas normalmente consolidadas
(schmertmann 1970)

Figura 20. Curva de consolidación para una arcilla normalmente consolidada.

Fuente: Tomado y adoptado de DAS Braja M. (1999). Fundamentos de Ingeniería


geotécnica.

H  σ 0' + ∆σ 
ρ= CcLog  

1 + e0  σ0
'

Donde:

H: Espesor del estrato

e0: Relación de vacíos

Cc: Índice de compresión

σ ’0: Presión de sobrecarga efectiva promedio inicial.

∆ σ : Incremento de la presión vertical debido a la aplicación de la carga.

Para el valor de Cc Skempton (1944) dio la siguiente correlación empírica para el


índice de compresión:

Para arcillas inalteradas: Cc = 0.009(LL − 10 )

54
Y según Nagarja y Murty (1985) propusieron:

LL %
Cc = 0.2343 Gs
100

Siendo:

LL : Límite líquido en la muestra de suelo.

Gs: La gravedad específica.

Para la relación de vacíos se tiene la siguiente expresión:

 G *γ 
e0 =  s w  − 1
 γ unsat 

Donde:

Gs : Gravedad específica

үw: Peso del agua


үunsat: Peso específico seco

55
Para un número arbitrario de capas:

Figura 21. Asentamiento para una arcilla normalmente consolidada para un


número arbitrario de capas.

σ0 σ0+Δσ

Figura 22. Asentamiento por consolidación para un número arbitrario de capas.

Fuente: Tomado y adoptado de DAS Braja M. (1999). Fundamentos de Ingeniería


geotécnica.

 Hi  σ 0' + ∆σ  
ρ = ∑ Cci Log  
1 + e0i  σ 0 i 
'

56
2.6.2.2 Asentamiento en arcillas pre-consolidadas (schmertmann 1970)

2.6.2.2.1 Caso ∆σ ˂ σ’p – σ’0

Figura 22. Curva de consolidación para una arcilla pre-consolidada caso ∆σ ˂ σ’p
– σ’0

Cs

Fuente: Tomado y adoptado de DAS Braja M. (1999). Fundamentos de Ingeniería


geotécnica.

H  σ 0' + ∆σ 
ρ= CsLog  
1 + e0  σ0
'

Donde:

Cs: Índice de expansión

Según Nagarja y Murty (1985) propusieron la siguiente correlación empírica:

57
LL %
Cs = 0.0463 Gs
100

2.6.2.2.2 Caso ∆σ > σ’p – σ’0

Figura 23. Curva de consolidación para una arcilla pre-consolidada caso ∆σ > σ’p
– σ’0

Cs

Fuente: Tomado y adoptado de DAS Braja M. (1999). Fundamentos de Ingeniería


geotécnica.

H  σ 'p   σ ' + ∆σ 
ρ= CsLog   + CcLog  0 ' 
 σ 
1 + e0   σ '0   p 

Donde:

Cs: Índice de expansión

Para un número arbitrario de capas se determina la profundidad desde la cual el


incremento de esfuerzo es mayor y menor que el de pre-consolidación.

58
Figura 24. Asentamiento por consolidación para arcillas pre-consolidada con un
número arbitrario de capas.

∆σ > σ 'p − σ o' σ0 σ0+Δσ

σ’p ∆σ < σ 'p − σ o'

Fuente: Tomado y adoptado de DAS Braja M. (1999). Fundamentos de Ingeniería


geotécnica.

 H  σ ´  σ ´+ ∆σ   Hi  σ ´+ ∆σ 
ρ = Σ i  Csi Log  p  + Cci Log  o  ρ = Σ * Csi Log  o  
 σ ´ 
 1 + eoi    σ o ´ i  p  1 + eoi  σ o´ i 

2.6.3 Deformaciones por consolidación primaria

Es la reducción del volumen de la masa de un suelo originada por la aplicación de


una carga permanente y la expulsión del agua de los vacíos, acompañada por una
transferencia de carga del agua a las partículas sólidas del suelo.

En la consolidación primaria el más importante son las arcillas inorgánicas y


limos10.

59
2.6.4 Deformaciones por consolidación secundaria

Reducción en el volumen de la masa del suelo, causada por la aplicación de una


carga permanente y el acomodo de la estructura interna de su masa, luego de que
la mayor parte de la carga ha sido transferida a las partículas sólidas del suelo.

Es la consecuencia del ajuste plástico de la estructura del suelo y es sucedido de


la consolidación primaria, es decir, después que los excesos de presión en el agua
se han disipado y que toda la carga ha sido entregada a la estructura de suelo. Es
un asentamiento que sucede en el tiempo a esfuerzo constante, es consecuencia
del reajuste de la estructura del suelo. Se presume que materiales con alto
contenido orgánico presenten este comportamiento (arcillas orgánicas). Cabe
mencionar que además existen otros tipos de deformaciones debidos a causas en
los procesos constructivos, es tal el caso de los asentamientos generados por las
excavaciones en las cercanías, por vibraciones, por aumentarse la carga en las
zonas aledañas, etc. Además de éstos existen los colapsos o grandes
desplazamiento del suelo de soporte cuando se sobrecarga una fundación o a los
derrumbes debidos a erosión y desplazamientos geológicos, estos dos últimas
deformaciones son prácticamente indescifrables, por lo cual no se tratan de
cuantificar sino que se trata de utilizar medidas preventivas o evitar áreas donde
pueda presentarse10.

El asiento secundario se calcula según la razón:

 T 
∆H s = HCs Log 
T
 100 
Donde:

H es la altura del estrato en fase de consolidación;

Cs es el índice de expansión

T tiempo en que se desea calcular el asiento secundario;

T100 tiempo necesario para terminar el proceso de consolidación primaria.

Los supuestos de este método son:

• La consolidación secundaria inicia después de terminado el proceso de


consolidación primaria;
• El valor de Cs se puede considerar constante mientras se da el asiento
secundario.
60
2.6.5 Asentamiento tolerable de las fundaciones de estructuras

Es de vital importancia establecer cuál es la deformación tolerable de la estructura,


para poder confrontarlo con la deformación calculada, pues se debe garantizar
que las deformaciones no afectaran el comportamiento estructural de la
edificación.

Para establecer los límites de deformaciones totales se aconseja dirigirse al


capítulo H.4.9.2 de la NSR-10:

Las deformaciones totales calculados a 20 años se deben limitar a los siguientes


valores:

• Para construcciones aisladas 30 cm, siempre y cuando no se afecten la


funcionalidad de conducciones de servicios y accesos a la construcción.
• Para construcciones entre medianeros 15 cm, siempre y cuando no se
afecten las construcciones e instalaciones vecinas.

H.4.9.3 Los límites para deformaciones diferenciales: Las deformaciones


diferenciales calculados se deben limitar a los valores fijados en la tabla 16.,
expresados en función de L, distancia entre apoyos o columnas de acuerdo con el
tipo de construcción11.

Tabla 16. Valores máximos de deformaciones o asentamientos diferenciales.

TIPO DE CONSTRUCCIÓN ∆ máx

Edificaciones con muros y acabados susceptibles de dañarse


L/1000
con asentamientos menores

Edificaciones con muros de carga en concreto o en mampostería L/500

Edificaciones con pórticos en concreto, sin acabados


L/300
susceptibles de dañarse con asentamientos menores.

Edificaciones en estructuras metálicas, sin acabados


L/160
susceptibles de dañarse con asentamientos menores.

Fuente: Norma Sismo Resistente NSR-10

11
Reglamento colombiano de construcción sismo resistente, NSR 10 (2010). Bogotá. Título H, y
Título G.

61
3. METODOLOGÍA

A partir de la recopilación de antecedentes y en la consulta de estudios previos


realizados en el área de interés: ‘’CONSTRUCCIÓN DEL CUARTO DE CONTROL
ESTACIÓN 3 ZONA INDUSTRIAL - YONDÓ’’, nos permite la recolección de
ensayos de laboratorio y registros de deformaciones elásticas ‘’reales’’ medidos
en campo por medio de la topografí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.

Para el análisis de los asentamientos inmediatos y por consolidación se realiza por


medio de las teorías tradicionales y mediante el uso del software de elementos
finitos PLAXIS 3D (asentamientos inmediatos) y PLAXIS 2D (asentamientos por
consolidación).

Esto con el fin de comparar los resultados obtenidos por los dos métodos.

Dentro de los parámetros requeridos para la modelación utilizando los programas


se hizo necesario la obtención de otros parámetros geomecánicos tales como:
pesos unitarios, módulos de Young, coeficientes de Poisson, ángulo de fricción y
cohesiones.

Dichos parámetros fueron calculados a partir de correlaciones, debido a la


dificultad para la obtención de muestras inalteradas en campo.

Ver en anexos cuadro de resumen de parámetros geomecánicos.

A continuación se muestra el procedimiento para el uso del software de LOADCAP


para los cálculos de capacidades portantes:

3.1 CALCULO DE CAPACIDADES PORTANTES

3.1.1 Datos generales

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.

3.1.1.2 Datos sistema de cimentación

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.

3.1.1.4 Input gráfico

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.

Figura 29. Ventana de estratigrafía.

67
3.1.1.5 Cargas

Las cargas sobre la cimentación se diferencian entre cargas de proyecto o de


servicio. Ambos tipos se deben introducir respetando las convenciones indicadas.

Figura 30. Convención del signo de las solicitaciones.

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.

El tipo de carga, proyecto o servicio, se selecciona en la última columna (Tipo).

Figura 31. Ventana selección del tipo de carga.

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.

La carga última vertical y horizontal se calcula para cada combinación de proyecto.

Con el menú Opciones de análisis, en la ventana de cálculo de la carga última, se


puede seleccionar el tipo de análisis a efectuar. Para este caso aplica la condición
drenada y el cálculo según el promedio ponderado de las estratificaciones, dado
que esta opción permite estimar la carga última de la cimentación considerando
como parámetros geotécnicos aquellos obtenidos con el promedio ponderado de
los parámetros de los estratos individuales.

Figura 32. Ventana de cargas para selección de opciones de análisis.

Con cada combinación es posible reducir los parámetros característicos del


terreno según los coeficientes de reducción establecidos y efectuar las
Correcciones sísmicas de acuerdo con los enfoques de proyectos seleccionados.

69
Para el cálculo del coeficiente de balasto ks el programa propone la fórmula
adoptada por Bowles.

3.1.3 Gráficos

3.1.3.1 Diagrama de tensiones

Muestra un diagrama donde se ven las variaciones de las tensiones en el


subsuelo, calculadas al centro de la cimentación, en función de la profundidad z,
con la teoría de Boussinesq o de Westergaard (según la selección efectuada en
Bulbo tensiones). En este caso se hizo uso de la teoría de Boussinesq.

Figura 33. Gráfico de tensiones según Boussinesq.

3.1.3.2 Parámetros gráficos

Permite modificar los parámetros de los gráficos de output (carga última-base de


la cimentación, carga última-profundidad del plano de cimentación), seleccionar el
autor para el cálculo de las tensiones en el subsuelo, intervenir, con la opción de
variación de longitud de la cimentación, en la tabla Informe diagramas.

70
Figura 34. Tabla de parámetros gráficos.

3.1.3.3 Informe diagramas

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ú.

3.1.3.4 Gráfico carga Última – Base (Q-B)

Muestra un diagrama donde se ven las variaciones de la carga última calculada


con los métodos elegidos en función de la anchura de la cimentación B. Para cada
conjunto de valores del plano de cimentación D y de la longitud L se construye un
diagrama, que se puede copiar y/o imprimir (menú Edición) o exportarlo (menú
Exportar). Cada uno de estos comandos se puede también seleccionar apretando
la parte derecha del mouse.

71
Figura 35. Gráfico Carga última – Profundidad (B-Qlim).

Ver Anexo.

Finalmente se da clic en la opción Exportar para así obtener el resumen de los


cálculos realizados de las capacidades portantes por los diferentes autores,
obteniendo los siguientes resultados:

DATOS GENERALES

======================================================

Anchura cimentación 1.0 m

Longitud cimentación 1.0 m

Profundidad plano de cimentación 1.0 m

Altura de empotramiento 1.0 m

Inclinación talud 27.0 °

Distancia cimentación talud 1.0 m

Profundidad nivel freático 4.5

=======================================================

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.

Tabla 17. Parámetros para cálculo de la capacidad portante.

DH γunsat γsat φ c E
ν
[m] [kN/m³] [kN/m³] [°] [kN/m²] [kN/m²]

1.5 18.5 20.1 35.0 2.0 13500.0 0.31

1.5 18.5 20.1 35.0 2.0 13500.0 0.31

0.6 20.2 22.5 41.0 1.0 19000.0 0.20

0.9 12.0 17.0 28.0 21.1 5813.0 0.34

1.13 13.0 18.5 26.0 21.1 4721.0 0.37

1.13 13.0 18.5 26.0 21.1 4721.0 0.37

1.85 14.0 19.0 30.0 21.1 5695.0 0.33

2.0 15.0 20.0 30.0 21.1 5565.0 0.33

Tabla 18. Cargas de proyecto actuantes en cimentación.

Presión
normal
Nombre N Mx My Hx Hy
de
Nr. Tipo
combinación proyecto
[kN] [kN·m] [kN·m] [kN] [kN]
[kN/m²]

1 Carga última 306.79 0.00 0.00 0.00 0.00 0.00 Proyecto

73
CARGA ÚLTIMA CIMENTACIÓN COMBINACIÓN…Carga última

Autor: TERZAGHI (1955)

Carga última [Qult] 920.37 kN/m²

Resistencia de proyecto [Rd] 306.79 kN/m²

Tensión [Ed] 306.79 kN/m²

Factor de seguridad [Fs=Qult/Ed] 3.0

Condición de comprobación [Ed<=Rd]Comprobado

COEFICIENTE DE ASENTAMIENTO BOWLES (1982)

Constante di Winkler 36814.91 kN/m³

• Carga última

Autor: HANSEN (1970) (Condición drenada)

======================================================

Factor [Nq] 26.29

Factor [Nc] 46.12

Factor [Ng] 33.92

Factor forma [Sc] 1.72

Factor profundidad [Dc] 1.0

Factor inclinación cargas [Ic] 1.0

Factor inclinación talud [Gc] 1.0

Factor inclinación base [Bc] 1.0

Factor forma [Sq] 1.7

Factor profundidad [Dq] 1.0

74
Factor inclinación cargas [Iq] 1.0

Factor inclinación talud [Gq] 1.0

Factor inclinación base [Bq] 1.0

Factor forma [Sg] 0.6

Factor profundidad [Dg] 1.0

Factor inclinación cargas [Ig] 1.0

Factor inclinación talud [Gg] 1.0

Factor inclinación base [Bg] 1.0

Factor corrección sísmico inercial [zq] 1.0

Factor corrección sísmico inercial [zg] 1.0

Factor corrección sísmico inercial [zc] 1.0

======================================================

Carga última 1174.15 kN/m²

Resistencia de proyecto 391.38 kN/m²

Condición de comprobación [Ed<=Rd]Comprobado

======================================================

Autor: TERZAGHI (1955) (Condición drenada)

======================================================

Factor [Nq] 26.29

Factor [Nc] 46.12

Factor [Ng] 42.43

Factor forma [Sc] 1.3

Factor forma [Sg] 0.8

Factor corrección sísmico inercial [zq] 1.0

75
Factor corrección sísmico inercial [zg] 1.0

Factor corrección sísmico inercial [zc] 1.0

======================================================

Carga última 920.37 kN/m²

Resistencia de proyecto 306.79 kN/m²

Condición de comprobación [Ed<=Rd]Comprobado

======================================================

Autor: MEYERHOF (1963) (Condición drenada)

======================================================

Factor [Nq] 26.29

Factor [Nc] 46.12

Factor [Ng] 37.15

Factor forma [Sc] 1.74

Factor profundidad [Dc] 1.0

Factor inclinación cargas [Ic] 1.0

Factor forma [Sq] 1.37

Factor profundidad [Dq] 1.0

Factor inclinación cargas [Iq] 1.0

Factor forma [Sg] 1.37

Factor profundidad [Dg] 1.19

Factor inclinación cargas [Ig] 1.0

Factor corrección sísmico inercial [zq] 1.0

Factor corrección sísmico inercial [zg] 1.0

Factor corrección sísmico inercial [zc] 1.0

76
======================================================

Carga última 1387.13 kN/m²

Resistencia de proyecto 462.38 kN/m²

Condición de comprobación [Ed<=Rd]Comprobado

======================================================

Autor: VESIC (1975) (Condición drenada)

======================================================

Factor [Nq] 26.29

Factor [Nc] 46.12

Factor [Ng] 48.03

Factor forma [Sc] 1.72

Factor profundidad [Dc] 1.0

Factor inclinación cargas [Ic] 1.0

Factor inclinación talud [Gc] 1.0

Factor inclinación base [Bc] 1.0

Factor forma [Sq] 1.7

Factor profundidad [Dq] 1.0

Factor inclinación cargas [Iq] 1.0

Factor inclinación talud [Gq] 1.0

Factor inclinación base [Bq] 1.0

Factor forma [Sg] 0.6

Factor profundidad [Dg] 1.0

Factor inclinación cargas [Ig] 1.0

Factor inclinación talud [Gg] 1.0

77
Factor inclinación base [Bg] 1.0

Factor corrección sísmico inercial [zq] 1.0

Factor corrección sísmico inercial [zg] 1.0

Factor corrección sísmico inercial [zc] 1.0

======================================================

Carga última 1252.45 kN/m²

Resistencia de proyecto 417.48 kN/m²

Condición de comprobación [Ed<=Rd]Comprobado

======================================================

Autor: Brinch - Hansen 1970 (Condición drenada)

======================================================

Factor [Nq] 26.29

Factor [Nc] 46.12

Factor [Ng] 45.23

Factor forma [Sc] 1.59

Factor profundidad [Dc] 1.0

Factor inclinación cargas [Ic] 1.0

Factor inclinación talud [Gc] 1.0

Factor inclinación base [Bc] 1.0

Factor forma [Sq] 1.57

Factor profundidad [Dq] 1.0

Factor inclinación cargas [Iq] 1.0

Factor inclinación talud [Gq] 1.0

Factor inclinación base [Bq] 1.0

78
Factor forma [Sg] 0.7

Factor profundidad [Dg] 1.0

Factor inclinación cargas [Ig] 1.0

Factor inclinación talud [Gg] 1.0

Factor inclinación base [Bg] 1.0

Factor corrección sísmico inercial [zq] 1.0

Factor corrección sísmico inercial [zg] 1.0

Factor corrección sísmico inercial [zc] 1.0

======================================================

Carga última 1205.1 kN/m²

Resistencia de proyecto 401.7 kN/m²

Condición de comprobación [Ed<=Rd]Comprobado

======================================================

3.2 DEFORMACIONES O ASENTAMIENTOS

3.2.1 Deformaciones elásticas o inmediatos mediante la teoria clásica

Para el cálculo de los asientos elásticos (inmediatos) se calculan con una


ecuación basada en la teoría de la elasticidad de Timoshenko y Goodier (1951).

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.

Tabla 19. Parámetros del suelo para el cálculo de asentamientos elásticos.

γunsat E
TIPO DE γsat
ESPESOR ν
MATERIAL ESTRATOS (KN/m3) (KN/m3) (KN/m2)
(m)

Relleno muro 3.00 18.5 20.10 13500 0.31

Base granular 0.60 20.2 22.5 19000 0.20

Arcilla 1 0.90 12.0 17.0 5813 0.34

Arcilla 2 2.25 13.0 18.5 4721 0.37

Limo 1 1.85 14.0 19.0 5695 0.33

Limo 2 2.00 15.0 20.0 5565 0.33

80
Inicialmente se calculan los asentamientos elásticos generados por la carga del
muro.

La carga del muro se calcula así:

qc = γ unsat * H = 18.5 * 3 + 12 * 0.6 = 67.5 KN / m 2

Luego haciendo uso de la ecuación de la teoría de elasticidad (Timoshenko y


Goodier), se obtienen los siguientes asentamientos elásticos generados por el
muro:

Tabla 20. Deformaciones elásticas en los diferentes estratos generados por el


muro.

ESTRATO No. 1 (Arcilla 1)

H (m) 0.90
H/B 0.10
B (m) 8.7

L (m) 14.7 L/B 1.69

qc (KN/m2) 67.56 C'd 0.09

ρ1 (m) 0.008

ESTRATO No. 2 (Arcilla 2) ESTRATO No. 2 (Arcilla 2)

H (m) 3.15 H (m) 0.90


H/B 0.36 H/B 0.10
B (m) 8.7 B (m) 8.7

L (m) 14.7 L/B 1.69 L (m) 14.7 L/B 1.69

qc (KN/m2) 67.56 C'd 0.34 qc (KN/m2) 67.56 C'd 0.09

ρ2a (m) 0.037 ρ2b (m) 0.010

ρ2 (m) 0.027

81
ESTRATO No. 3 (Limo 1) ESTRATO No. 3 (Limo 1)

H (m) 5.00 H (m) 3.15


H/B 0.57 H/B 0.36
B (m) 8.7 B (m) 8.7

L (m) 14.7 L/B 1.69 L (m) 14.7 L/B 1.69

qc (KN/m2) 67.56 C'd 0.54 qc (KN/m2) 67.56 C'd 0.34

ρ3a (m) 0.050 ρ3b (m) 0.050

ρ3 (m) 0.018

ESTRATO No. 4 (Limo 2) ESTRATO No. 4 (Limo 2)

H (m) 7.00 H (m) 5.00


H/B 0.80 H/B 0.57
B (m) 8.7 B (m) 8.7

L (m) 14.7 L/B 1.69 L (m) 14.7 L/B 1.69

qc (KN/m2) 67.56 C'd 0.75 qc (KN/m2) 67.56 C'd 0.54

ρ4a (m) 0.071 ρ4b (m) 0.051

ρ4 (m) 0.020

ρ elásticos muro (m) 0.073

ρ elásticos muro (cm) 7.307

Posteriormente se realiza el mismo procedimiento para la carga de la estructura.

82
Figura 37. Perfil de suelo con la carga de la estructura.

Carga dada por el ingeniero estructural para el cimiento:

qc = 7.1KN / m 2

Tabla 21. Deformaciones elásticas en los diferentes estratos generados por la


carga estructural.

ESTRATO No. 1 (Arcilla 1)

H (m) 0.90
H/B 0.90
B (m) 1

L (m) 1 L/B 1.00

qc (KN/m2) 7.1 C'd 0.68

ρ1 (m) 0.001

83
ESTRATO No. 2 (Arcilla 2) ESTRATO No. 2 (Arcilla 2)

H (m) 3.15 H (m) 0.90


H/B 3.15 H/B 0.90
B (m) 1 B (m) 1

L (m) 1 L/B 1.00 L (m) 1 L/B 1.00

qc (KN/m2) 7.1 C'd 1 qc (KN/m2) 7.1 C'd 0.68

ρ2a (m) 0.0013 ρ2b (m) 0.0009

ρ2 (m) 0.0004

ESTRATO No. 3 (Limo 1) ESTRATO No. (Limo 1)

H (m) 5.00 H (m) 3.15


H/B 5.00 H/B 3.15
B (m) 1 B (m) 1

L (m) 1 L/B 1.00 L (m) 1 L/B 1.00

qc (KN/m2) 7.1 C'd 1.05 qc (KN/m2) 7.1 C'd 1

ρ3a (m) 0.00117 ρ3b (m) 0.00111

ρ3 (m) 0.00006

ESTRATO No. 4 (Limo 2) ESTRATO No. 4 (Limo 2)

H (m) 7.00 H (m) 5.00


H/B 7.00 H/B 5.00
B (m) 1 B (m) 1

L (m) 1 L/B 1.00 L (m) 1 L/B 1.00

qc (KN/m2) 7.1 C'd 1.12 qc (KN/m2) 7.1 C'd 1.05

ρ4a (m) 0.0013 ρ4b (m) 0.0012

ρ4 (m) 0.0001

84
ρ elásticos estructura (m) 0.001

ρ elásticos estructura (cm) 0.129

Los asentamientos elásticos totales obtenidos fueron de:

ρ elásticos totales (m) 0.074

ρ elásticos totales (cm) 7.436

Seguidamente se procede a analizar el muro por medio de del software de


elementos finitos PLAXIS 3D (versión 8.5)12, para el cálculo de los asentamientos
inmediatos.

3.2.2 Deformaciones elásticas o inmediatos mediante elementos finitos


(plaxis 3d)

El planteamiento que a continuación se va a describir tiene por objeto modelar y/o


simular el posible comportamiento de la estructura (muro en tierra armada), el cual
cuenta con las siguientes fases para el manejo del software: Introducción de datos
y pre-procesador 3D, cálculos y Post-procesos de los resultados. En este análisis
se pretende el cálculo de deformaciones en la estructura.

3.2.2.1 Introducción de datos y pre-procesador 3d

3.2.2.1.1 Geometría de entrada

Inicialmente se establecen los parámetros básicos del modelo de elementos


finitos. Esto se hace en la ventana de File – New. En esta ventana se incluye la
descripción del problema, el tipo de análisis, los elementos de tipo básico, las
unidades básicas y las dimensiones para el dibujo.

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).

3.2.2.1.2 Geometría del contorno

Una vez realizada la configuración general, se empieza a dibujar la geometría con


la línea geométrica (Geometry Line). Al seleccionar este icono, se introducen las
coordenadas escribiéndolas en la línea de instrucciones de la siguiente manera: (x
y).

Coordenadas de la geometría a interés: (-4.35 -7.35), (4.35 -7.35), (-4.35 7.35),


(4.35 7.35).

86
3.2.2.1.3 Planos de trabajo

Esta opción de planos de trabajo (workplanes) permite definir las profundidades


dependiendo del número de estratos elegidos para el diseño del muro.

Figura 39. Ventana de workplanes.

3.2.2.1.4 Propiedades del material

En Plaxis 3D las propiedades del suelo y de los materiales que constituyen las
estructuras son almacenadas en conjuntos de datos de los materiales.

La base de datos de materiales puede ser activada seleccionando la opción


Materials – Soils and interfaces. Si se desea crear una base nueva de materiales
que no existen en la actual base, se da clic en new en la parte inferior de la
ventana material sets. Se pueden crear varios conjuntos de datos para distinguir
entre diferentes capas de suelo. Las propiedades que se especifican en cada
conjunto de datos se clasifican en tres categorías: General, Parameters e
Interfaces. Cada categoría tiene asociada una pestaña dentro de la
correspondiente ventana.

• La pestaña general, contiene el tipo de modelo constitutivo del suelo, el tipo


de comportamiento del suelo (drenado, no drenado o poroso) y las propiedades
generales del suelo tales como el peso específico.
• La pestaña Parameters contiene los parámetros de rigidez y de resistencia
correspondientes al modelo constitutivo seleccionado.

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.

Ver figura 39.

Tabla 22. Resumen de parámetros para el material de relleno.

PARÁMETRO SIMBOLO VALOR UNIDAD

GENERAL

Nombre del material Identification Relleno -

Mohr-
Modelo de material Model Coulomb -

Tipo de material Type Drained -

Peso específico del suelo γ unsat 18.5 KN/m³

Peso específico saturado del


suelo γ sat 20.1 KN/m³

PARAMETERS

Módulo de Young E ref 13500.0 KN/m²

Coeficiente de poisson ν 0.31

Cohesión c ref 2.0 KN/m²

Ángulo de fricción φ 35.0 °

Ángulo de dilatancia Ψ 0.0 °

INTERFACES

Tipo de material Strength Rigid -

88
Tabla 23. Resumen de parámetros para el material de base granular.

PARÁMETRO SIMBOLO VALOR UNIDAD

GENERAL

Base
Nombre del material Identification granular -

Mohr-
Modelo de material Model Coulomb -

Tipo de material Type Drained -

Peso específico del suelo γ unsat 20.2 KN/m³

Peso específico saturado del


suelo γ sat 22.5 KN/m³

PARAMETERS

Módulo de Young E ref 19000.0 KN/m²

Coeficiente de poisson ν 0.20

Cohesión c ref 1.0 KN/m²

Ángulo de fricción φ 41.0 °

Ángulo de dilatancia Ψ 0.0 °

INTERFACES

Tipo de material Strength Rigid -

89
Tabla 24. Resumen de parámetros para el material arcilla 1.

PARÁMETRO SIMBOLO VALOR UNIDAD

GENERAL

Nombre del material Identification Arcilla 1 -

Mohr-
Modelo de material Model Coulomb -

Tipo de material Type Drained -

Peso específico del suelo γ unsat 12 KN/m³

Peso específico saturado del


suelo γ sat 17 KN/m³

PARAMETERS

Módulo de Young E ref 5813.0 KN/m²

Coeficiente de poisson ν 0.34

Cohesión c ref 21.1 KN/m²

Ángulo de fricción φ 28.0 °

Ángulo de dilatancia Ψ 0.0 °

INTERFACES

Tipo de material Strength Rigid -

90
Tabla 25. Resumen de parámetros para el material arcilla 2.

PARÁMETRO SIMBOLO VALOR UNIDAD

GENERAL

Nombre del material Identification Arcilla 2 -

Mohr-
Modelo de material Model Coulomb -

Tipo de material Type Drained -

Peso específico del suelo γ unsat 13 KN/m³

Peso específico saturado del


suelo γ sat 18.5 KN/m³

PARAMETERS

Módulo de Young E ref 4721.0 KN/m²

Coeficiente de poisson ν 0.37

Cohesión c ref 21.1 KN/m²

Ángulo de fricción φ 26.0 °

Ángulo de dilatancia Ψ 0.0 °

INTERFACES

Tipo de material Strength Rigid -

91
Tabla 26. Resumen de parámetros para el material limo 1.

PARÁMETRO SIMBOLO VALOR UNIDAD

GENERAL

Nombre del material Identification Limo 1 -

Mohr-
Modelo de material Model Coulomb -

Tipo de material Type Drained -

Peso específico del suelo γ unsat 14 KN/m³

Peso específico saturado del


suelo γ sat 19 KN/m³

PARAMETERS

Módulo de Young E ref 5695.0 KN/m²

Coeficiente de poisson ν 0.33

Cohesión c ref 21.1 KN/m²

Ángulo de fricción φ 30.0 °

Ángulo de dilatancia Ψ 0.0 °

INTERFACES

Tipo de material Strength Rigid -

92
Tabla 27. Resumen de parámetros para el material limo 2.

PARÁMETRO SIMBOLO VALOR UNIDAD

GENERAL

Nombre del material Identification Limo 2 -

Mohr-
Modelo de material Model Coulomb -

Tipo de material Type Drained -

Peso específico del suelo γ unsat 15 KN/m³

Peso específico saturado del


suelo γ sat 20 KN/m³

PARAMETERS

Módulo de Young E ref 5565.0 KN/m²

Coeficiente de poisson ν 0.33

Cohesión c ref 21.1 KN/m²

Ángulo de fricción φ 30.0 °

Ángulo de dilatancia Ψ 0.0 °

INTERFACES

Tipo de material Strength Rigid -

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.

Luego de haber creado el conjunto de datos (material de relleno, base granular,


arcilloso y limoso) para este proyecto se arrastran de la ventana de materiales a la
columna estratigráfica.

Figura 41. Ventana de Borehole.

3.2.2.1.6 Aplicación de carga

Para la aplicación de la carga se hace uso de la opción Loads – Distributed loads


(horiz. Planes), en seguida de ser activada esta opción se da clic sobre el contorno
de la geometría. Para modificar el valor de la carga se da doble clic en la línea
geométrica correspondiente y seleccionando el sistema de carga de interés (en

95
este caso en el eje Y, es decir perpendicular al plano de trabajo) se inserta el
valor.

Figura 42. Ventana de introducción de cargas.

3.2.2.1.7 Generación de la malla

Una vez el modelo geométrico se encuentre completamente definido y las


propiedades de los materiales han sido asignadas, la geometría ha de ser dividida
en elementos finitos con el fin de llevar a cabo los consiguientes cálculos.

La generación de la malla se hace dando clic en Generate Mesh 2D y luego se da


clic en Generate Mesh 3D que funciona como un refinamiento con el fin de que los
resultados obtenidos en los cálculos sean mucho más precisos.

Figura 43. Ventana de generación de la malla 2D.

96
Figura 44. Ventana de generación de la malla 3D.

3.2.2.2 Cálculos

Para el inicio de cálculos, con la generación de la situación inicial del modelo de


elementos finitos, se hace clic en Calculate en el cual apare un recuadro de
diálogo en el que pide que se guarde los datos.

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.

• Fase Incial (Initial phase): Primeramente se verifica que en los niveles se


encuentre el tipo de material adecuado que fue seleccionado inicialmente en
borehole. Luego ubicándose en el nivel 0.00 y dando clic sobre el modelo se
desactiva Soil below, y Distributed load, es decir que en esta fase no se encuentra
la carga del muro, ni la carga distribuida. En general, esta fase comprende las
condiciones iniciales de nivel freático, la configuración de la geometría inicial y el
estado inicial de tensión efectiva.

97
Figura 45. Ventana de cálculos (definición de fases).

• Fase I (Carga muro): Esta fase adicional a la anterior tiene en cuenta el


peso del muro en tierra armada. Por tanto dado click sobre el modelo solo se
activa la opción Soil below, es decir que ya se tiene en cuenta a partir de ese nivel
el suelo hacia abajo.

• 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.

3.2.2.2.1 Puntos de control

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.

3.2.2.3 Post-procesos de los resultados

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.

Para obtener la gráfica de desplazamiento se da clic en Curves y en seguida


aparece una ventana y se da clic en New Chart y se selecciona el proyecto de
interés, luego se seleccionan el orden de los ejes X y Y para finalmente obtener la
curva de desplazamiento.

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.

Figura 48. Gráfico de curva de desplazamientos de la fase carga de la estructura.

100
3.2.3 Deformaciones o asentamientos por consolidación unidimensional

Según las condiciones de suelo que se tienen se asume el caso de arcillas


normalmente.

H  σ ' + ∆σ 
ρ= CcLog  0 ' 
1 + e0  σ0 

3.2.3.1 Incremento del esfuerzo bajo un terraplén

Para esta condición bidimensional de carga, el incremento del esfuerzo vertical se


expresa como7:

qc  B1 + B2  
(α 1 + α 2 ) − 1 (α 2 )
B
∆σ = 
π  B2  B2 

Donde:

qc = γH

γ = Peso específico del suelo del terraplén

H = Altura del terraplén

 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.

Fuente: DAS Braja M. (1999). Fundamentos de Ingeniería geotécnica.

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

SU KN/m KN/m KN/m KN/m KN/m KN/m


m m m m 3 2 2 2 2 2 m
CS

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

ρ muro (m) 49.96 cm

3.2.3.2 Incremento del esfuerzo debajo de una cimentación por el método


2:1

En cimentaciones es común usar un método aproximado para determinar el


incremento del esfuerzo con la profundidad causado por la construcción de una
cimentación. Según este, el incremento del esfuerzo a la profundidad z es7:

qc * B * L
∆σ =
(B + Z ) + (L + Z )

Esta ecuación se basa en la hipótesis de que el esfuerzo se dispersa desde la


cimentación a lo largo de líneas con una pendiente de 2 vertical a 1 horizontal.

De igual manera que en los asentamientos elásticos, en estos asentamientos


también se asume el análisis generado por la carga del muro y la carga
estructural7.

103
Figura 50. Método 2:1 para encontrar el incremento de esfuerzo debajo de una
cimentación.

Fuente: DAS Braja M. (1999). Fundamentos de Ingeniería geotécnica.

Aplicando el mismo perfil del suelo con la carga estructural de 7.10 KN/m2, se
tiene lo siguiente:

Tabla 29. Asentamientos por consolidación generados por la carga estructural.

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

ρ estructura (cm) 3.44

Obteniendo así un asentamiento total por consolidación de:

ρ consolidación (cm) 53.40 cm

Finalmente para el asentamiento total generado en el suelo, corresponde a la


sumatoria de los asentamientos elásticos totales más los asentamientos totales
por consolidación, dando como resultado:

3.2.4 Deformaciones o asentamientos por consolidación unidimensional


mediante elementos finitos (plaxis 2d)

Para el cálculo de la deformación por consolidación en elementos finitos PLAXIS


2D, el procedimiento para este software es prácticamente el mismo que se llevó a
cabo para la modelación de PLAXIS 3D, solo que se tienen en cuenta otros
parámetros adicionales tales como índice de compresión, índice de expansión,
relación de vacíos y permeabilidad.

Para la permeabilidad hidráulica se utilizó la siguiente relación:

105
Tabla 30. Rango de la permeabilidad hidráulica parra varios suelos.

Fuente: DAS Braja M. (1999). Fundamentos de Ingeniería geotécnica.

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.

Considerando en el mismo tipo de perfil que en la consolidación unidimensional se


obtiene el siguiente modelo en PLAXIS 2D.

Figura 51. Modelo generado en PLAXIS 2D.

Las deformaciones obtenidas en PLAXIS 2D son del orden de ρ = 55.02 cm.

106
3.2.5 Deformaciones o asentamientos diferenciales

Para finalizar se realizan los cálculos de asentamientos diferenciales según la


Norma NSR-10, como se denota en la tabla 18.

Para este tipo de construcción se adopta la expresión: ∆ máx = L/300.

A continuación en el plano de cimentación se pueden ver todas las distancias


entre columnas (L) posibles por analizar.

Figura 52. Cimentación con las distancias entre columnas.

B
3.55

A
2.40 4.65
1 2 3

Tabla 31. Cálculos de asentamientos diferenciales según la NSR-10.

LONGITUD ∆ máx UND

2.40 0.80 cm

4.65 1.55 cm

3.55 1.18 cm

107
4. ANÁLISIS COMPARATIVO

4.1 CAPACIDADES PORTANTES

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.

Figura 52. Capacidades de carga a un 1.0 m de profundidad.

500 462.38
CAPACIDAD DE CARGA (KN/m²)

450 401.70 417.48


391.38
400
350 306.79
300
250
200
150
100
50
0
TERZAGHI MEYERHOF HANSEN BRINCH - VESIC
HANSEN
AUTORES TEORÍAS CLÁSICAS

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.

4.2 DEFORMACIONES ELÁSTICAS

Para la deformaciones inmediatas calculadas y haciendo uso de la teoria clasica


de elasticidad (Timoshenko y Goodier) versus PLAXIS 3D, se obtuvieron los
siguientes resultados:

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

Por el método clásico y el método de elementos finitos se presentan porcentajes


de error del 49 % y 36 % para las deformaciones elásticas.

4.3 DEFORMACIONES POR CONSOLIDACIÓN

Para el análisis comparativo de deformaciones por consolidación se muestran los


resultados obtenidos por la teoría clásica de consolidación unidimensional de
Terzaghi por el método de incrementos de esfuerzos de un terraplén versus la
consolidación de elementos finitos en PLAXIS 2D.

109
Figura 54. Comparativo de deformaciones por consolidación.

Deformaciones por consolidación ρ (cm) 55.500

55.020
55.000

54.500

54.000

53.500 53.397

53.000

52.500
Consolidación unidimensional Plaxis 2D
Métodos

La diferencia porcentual obedece al 3.34 %.

4.4 DEFORMACIONES TOTALES POR LOS METODOS CLÁSICOS Y DE


ELEMENTOS FINITOS

Las deformaciones totales (elásticos + consolidación) obtenidas por las teorías


clásicas y por el método de elementos finitos mostraron las siguientes diferencias:

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

Y una diferencia porcentual de 1.63 %.

4.5 MEJORAMIENTO AL MODELO DE PLAXIS 2D

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.

Para la alternativa del refuerzo de geotextil se tiene el siguiente modelo:

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.

Figura 57. Deformaciones totales con refuerzo de geotextil.

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.

Figura 58. Fuerzas axiales generadas por el geotextil.

La deformación total teniendo en cuenta el geotextil da un valor de 44.89 cm.

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

Figura 61. Comparativo de deformaciones totales teniendo en cuenta el geotextil


en PLAXIS 2D.

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

Para las propiedades geomecánicas necesarias para la modelación de los


programas de PLAXIS 2D y 3D se hizo mediante el uso de correlaciones para la
obtención de pesos saturados, relaciones de vacíos, módulos de elasticidad,
gravedad específica, ángulos de fricción, índice de compresión, expansión y
permeabilidad.

En cuanto al análisis de las capacidades portantes la teoría de Terzaghi que


es la más clásica muestra lo conservadora que es en el momento de calcular y
para este caso es la que se podría elegir para el diseño por lo que se requiere una
carga estructural baja (7.1 KN/m2). Las teorías de Hansen, Brinch – Hansen y
Vesic mantienen valores similares y cercanos entre sí, pero la que más se aleja y
ofrece mayor capacidad portante es la teoría de Meyerhof.

Según las comparaciones realizadas anteriormente en las deformaciones


elásticas el método de elementos finitos se acercó más al valor real controlado en
campo con una diferencia porcentual del 36 %. La diferencia entre los dos
métodos (Plasticidad versus Plaxis 3d) es de 8.5 %.

Aunque no se tienen mediciones reales respecto a las deformaciones por


consolidación los dos métodos presentan diferencias pequeñas con un margen de
diferencia del 3.34 %; solo que el método de elementos finitos PLAXIS 2D permite
llevar el modelo más a realidad.

Para optimizar el modelo de consolidación y lograr disminuir las


deformaciones se propone incluir al modelo capas de cada 0.40 m con refuerzo de
geotextil, lo cual permite disminuir las deformaciones en un 18 % que corresponde
a 8.31cm. La diferencia entre los dos métodos (Consolidación versus Plaxis 2D) es
de 16 %.

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

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.

BOWLES J.E. (1988). Foundation Analysis and Design. Mc Graw Hill Book Co.

CUTA, Viviana y RENGIFO, María Fernanda. Análisis comparativo de capacidad


portante en medios homogéneos y estratificados para cimientos superficiales.
Trabajo de grado pregrado en Ingeniería Civíl. Bucaramanga: Universidad
Industrial de Santander. 2012.

DAS Braja M. (1999). Fundamentos de Ingeniería geotécnica. Thomson Learning.


California State University.

DAS Braja M. (2001). Principios de Ingeniería de cimentaciones. Thomson


Learning. Cuarta edición California Sate University.

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.

INSTITUTO NACIONAL DE VÍAS. Norma INV E-111-13. Ensayo de Penetración


normal (SPT) y muestreo de suelos con tubo partido. Bogotá D.C., 1996.

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.

Reglamento colombiano de construcción sismo resistente, NSR 10 (2010). Bogotá.


Título H, y Título G.

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.

117
ANEXOS

118

También podría gustarte