Universidad Austral de Chile

Facultad de ciencias de la Ingeniera

Escuela de Ingeniera Civil en Obras Civiles

RECUPERACIN DE MTODOS ANALITCOS Y COMPUTACIONALES EN INGENIERA DE FUNDACIONES ESFUERZOS EN PILOTES

Tesis para optar al ttulo de: Ingeniero Civil en Obras Civiles Profesor Patrocinante: Sr. Lus Collarte Concha. Ingeniero Civil. M. Sc. en Ingeniera Civil. Especialidad Hidrulica Mecnica de Suelos.

ANIBAL DANIEL MEZA PEREZ VALDIVIA CHILE 2007

A mis Padres, pilares fundamentales de mi vida. A mi hermano Alfonso, apoyo incondicional para el logro de mis metas Y a toda mi familia, por su entrega y dedicacin.

INDICE TEMATICO
Pg. CAPITULO I: INTRODUCCIN Y OBJETIVO 1

1.1 1.2 1.3

Introduccin Objetivos Metodologa

2 3 4

CAPITULO II:

ASPECTOS TEORICOS

2.1 2.1.1 2.1.2 2.1.3 2.1.4 2.1.5 2.1.6 2.1.6.1 2.1.6.2 2.1.7 2.2 2.3 2.4 2.4.1 2.4.2

Cimentaciones Profundas Pilotes Uso de los Pilotes Instalacin de Pilotes Comportamiento del Pilote durante la Hinca Efectos del Pilote en el Suelo Capacidad de Carga de los Pilotes Anlisis Esttico Anlisis Dinmico Cargas Laterales Mtodo de los Elementos Finitos FORTRAN MATHCAD13 Caractersticas Principales reas de Aplicacin

6 6 7 8 9 9 11 11 14 15 17 18 19 20 20

CAPITULO III:

ESTUDIO Y ANALISIS DE LOS PROGRAMAS EN FORTRAN IV

21

3.1 3.1.1 3.1.2 3.1.3 3.1.4

Programacin en FORTRAN IV Identificacin de las Variables Lectura de las Variables Clasificacin de los Parmetros Generacin de Parmetros de Clculo

22 22 22 23 24

3.1.5 3.2 3.3

Obtencin de Resultados Problemas de Ejecucin Criterio de Programacin

25 26 28

CAPITULO IV:

PROGRAMACION

29

4.1 4.1.1 4.1.1.1 4.1.1.2

Estudio y Anlisis de la Programacin Lateral Piles Introduccin al Problema Estimacin del Mdulo de Reaccin Horizontal del Suelo, kh

30 30 30 31

4.1.1.3 4.1.1.4

Mtodo de Anlisis, Matriz de Rigidez Matriz Global de Rigidez para un Elemento Dividido en N nodos.

34 35

4.1.1.5 4.1.2 4.1.2.1

Anlisis de Esfuerzos Static Piles Estimacin de la Capacidad de Soporte Por Punta

39 41 42

4.1.2.2

Estimacin de la Capacidad de Soporte Por Friccin Lateral

44

4.2 4.2.1 4.2.2

Programacin Final en Detalle LATERAL PILES STATIC PILES

47 47 53

CAPITULO V:

ANALISIS DE RESULTADOS

61

5.1 5.1.1 5.1.1.1 5.1.1.2 5.1.1.3 5.1.1.4 5.1.2 5.1.2.1 5.1.2.2

Comprobacin de los Resultados LATERAL PILES Problema 1 Problema 1 aplicado a LATERAL PILES Problema 2 Problema 2 aplicado a LATERAL PILES STATIC PILES Problema 1 Problema 1 aplicado a STATIC PILES

62 62 62 64 66 70 74 74 77

CAPITULO VI:

CONCLUSIONES

80

6.1.- Conclusiones 6.2.- Caractersticas de uso y recomendaciones para los programas 6.2.1.- LATERAL PILES 6.2.2.- STATIC PILES

81 82 82 83

BIBLIOGRAFIA

85

ANEXOS

ANEXO A ANEXO B ANEXO C ANEXO D

Programas originales en FORTRAN IV Matriz de Rigidez Global del Sistema Test de Penetracin Standard Programas Finales LATERAL PILES STATIC PILES

87 100 102 104 105 107 109 110 111

ANEXO E

Diagramas de flujo de programacin LATERAL PILES STATIC PILES

INDICE DE FIGURAS Y TABLAS

Pg. Figura 1: Distribucin de esfuerzos y reacciones entre un pilote sometido a carga lateral, y el suelo Figura 2.1: Divisin de un pilote en subelementos. Estado de Reacciones Figura 2.2: Elemento base de diseo Figura 3: Elemento 1 de la subdivisin de un pilote Figura 4: Elemento 2 de la subdivisin de un pilote Figura 5: Elemento 3 de la subdivisin de un pilote Figura 6: Elemento 4 de la subdivisin de un pilote Figura 7: Elemento particionado en 4 subelementos Figura 8: Pilote solicitado por una carga axial esttica Figura 9: Elemento de un pilote para la determinacin de la capacidad De soporte por friccin lateral Figura 10: Distribucin de esfuerzos y desplazamientos de un pilote Cargado lateralmente Figura 11: Distribucin de esfuerzos y rigideces de un pilote Cargado lateralmente (Problema 1) Figura 12: Distribucin de esfuerzos y rigideces de un pilote Cargado lateralmente (Programa) Figura 13: Distribucin de esfuerzos y rigideces de un pilote Cargado lateralmente (Problema 2) Figura 13: Distribucin de esfuerzos y rigideces de un pilote Cargado lateralmente (Programa) Tabla 1: Factores de Capacidad Nc, Nq y N, para un pilote solicitado lateralmente Tabla 2: Factor de Capacidad Nc* dado por Vesic para la estimacin de La capacidad de soporte por punta de un pilote Tabla 3: Factor de Capacidad Nq* dado por Vesic para la estimacin de La capacidad de soporte por punta de un pilote Tabla 4: Factores de Capacidad de resistencia dados por Terzagui y Vesic. (LATERAL PILES). Tabla 5: Factores de Capacidad de resistencia dados por Vesic (STATIC PILES). 55 49 43 43 73 66 65 63 52 34 35 36 36 37 37 39 41 44 33

32

Tabla 6: Valores entregados por el estudio de mecnica de suelos realizada al ex Hotel Pedro de Valdivia para el emplazamiento del nuevo casino. Tabla 7: Clculo de la capacidad total de soporte de un pilote realizado por el profesor gua. Tabla 8: Interpolacin realizada para encontrar los valores exactos correspondientes a la profundidad de 13,4 m Tabla 9: Tabla comparativa de los valores finales de la capacidad de soporte total de un pilote.

74

75

76

79

RESUMEN

En el presente trabajo se realiza la programacin en MATHCAD13 del diseo y verificacin de esfuerzos de pilotes para dos tipos de casos; el primero para un pilote con una carga solicitante lateral, aplicada en su cabeza superior, la cual puede ser tomada como una carga ssmica, y el segundo caso, para un pilote con carga axial esttica, solicitacin proveniente de una estructura.

Para ello se analizaron y estudiaron metodologas aplicadas en programas antiguos, creados en FORTRAN IV, los cuales debido a su extensa metodologa de programacin, se reprogramaron en un software mas actual.

Los datos y resultados obtenidos en los programas, fueron revisados y comparados con problemas propuestos de pilotes para ambos casos.

SUMMARY

This project presents two computer programs designed with the software MATHCAD13. These programs design and verify piles stress applied for two kinds of cases. The first one, was designed for a pile with a lateral loading. It was applied in the head of the pile. Which can be taken as a earthquake stress. The second one, was designed for a pile with a static loading which was originated in the main structure.

In this way were studied and analyzed two applicable methodologies with old programs created in FORTRAN IV.

The data and results obtained in the programs were reviewed and compared using real examples of piles for both cases.

CAPITULO I
INTRODUCCIN Y OBJETIVO

1.1.- INTRODUCCION En Ingeniera de Fundaciones, existen diferentes soluciones ante problemas que surgen del diseo y clculo de estructuras, en donde se debe asegurar un sistema de fundacin adecuado y funcional que sea capaz de responder antes las solicitaciones internas y externas a las cuales ser sometida una estructura o elemento dado.

Una de las soluciones frente a determinados tipos de suelos son las cimentaciones profundas, en donde uno de los principales elementos son los Pilotes, los cuales deben ser previamente diseados, confeccionados y ensayados antes de su utilizacin final en un sistema de fundacin. Existen diversos mtodos de clculo de esfuerzos y reacciones para un pilote, los que han sido usados durante aos, as como tambin, existen programas computacionales que se encargan de disear pilotes, utilizando, mtodos de clculo como el de elementos finitos, mecnica de suelos y teora sobre fundaciones. Existen programas creados en software de lenguajes antiguos de programacin, como los creados por Bowles (1974) que son la base de programas actuales utilizados por ingenieros para el clculo y diseo de fundaciones, y cuya programacin ocupa metodologa de anlisis estructural y clculo matricial. Este trabajo de tesis consta en recuperar mtodos analticos y computacionales en ingeniera de fundaciones, basndose en dos programas que calculan esfuerzos y reacciones en pilotes con distintas solicitaciones, creados en FORTRAN IV por Joseph Bowles en el ao 1974 en su libro Analytical and Computer Methods in Foundation Engineering. Las metodologas recuperadas se programaran en el software de diseo ingenieril MATHCAD13, en el cul se analizaran dos casos de diseo de pilotes; en el primero, se modelara el pilote como un elemento tipo barra con dos grados de libertad por nudo (traslacin y rotacin), y apoyado resortes que representan las caractersticas del suelo, solicitado con una carga lateral aplicada en la cabeza del pilote; el segundo, tomando el modelo clsico dado por autores tales como Vesic o Meyerhoff, para la determinacin de la Capacidad de Resistencia Ultima de un pilote introducido en un suelo de distintas estratigrafas.

1.2.- OBJETIVOS

Objetivo general: Recuperar dos programas de metodologa Analtica y computacional en Ingeniera de Fundaciones, escritos en FORTRAN IV, que calculan los esfuerzos y reacciones en pilotes.

Objetivos especficos:

Estudiar y revisar los procedimientos empleados, revisar la teora aplicada, as como tambin, los mtodos de clculo utilizados en la programacin.

Revisar el lenguaje de programacin en un software adecuado, realizando las correcciones necesarias para el funcionamiento de los programas.

Compilar y ejecutar los programas, realizando las correcciones necesarias para su funcionamiento.

En el caso de no hacer funcionar uno o los dos programas elegidos, programar en un software de Ingeniera adecuado, la metodologa empleada en los programas originales.

Analizar y verificar diferentes problemas de pilotes, revisando los resultados entregados con mtodos aproximados en ingeniera de fundaciones.

1.3.- METODOLOGA

Se comenzar con una fase de anlisis exploratorio al texto que contiene la escritura de programacin, en la cual se estudiaran y elegirn los programas que sern recuperados, para luego buscar apoyo de la literatura existente en relacin al tema de este proyecto, en donde se extraer y recopilar el material que sea necesario para la realizacin de ste. A continuacin se realizara un estudio del lenguaje ocupado y un anlisis de la metodologa empleada en los programas, corroborando las frmulas empleadas, as como tambin los mtodos de clculo aplicados para el funcionamiento del programa. Luego se proceder a la escritura de los programas, en un software adecuado de lenguaje similar, el cual permitir compilar, realizar las correcciones necesarias y ejecutar los programas. En el caso de no poder ejecutar y recuperar los programas, se extraer la metodologa y procedimientos aplicados en el texto original, apoyndose en otros textos, materia de asignaturas correspondientes y en el profesor gua, para reprogramarlos en un software ms actual de mejores caractersticas y propiedades. Una vez ejecutados los programas, se rectificaran sus resultados con problemas vistos en Fundaciones, para realizar un anlisis comparativo y as finalmente, poder validar su funcionamiento. Los programas recuperados, o los programas creados, en el caso de reprogramarlos, sern entregados al Instituto de Obras civiles de la universidad austral de Chile, quedando a disposicin de acadmicos y estudiantes de las carreras de Ingeniera en Construccin y de Ingeniera Civil en Obras Civiles.

CAPITULO II
ASPECTOS TEORICOS

2.1. - CIMENTACIONES PROFUNDAS

El cimiento es aquella parte de la estructura encargada de transmitir las cargas al terreno. Las cimentaciones de la mayor parte de las estructuras se encuentran debajo de la superficie del terreno.

Cuando el suelo situado al nivel en que se desplantara normalmente una zapata o una losa de fundacin, es demasiado dbil o compresible para proporcionar un soporte adecuado, las cargas se deben transmitir a un material mas adecuado a mayor profundidad lo que se denomina como cimentacin profunda. (Terzagui, 1973).

Las cimentaciones profundas se emplean cuando los estratos de suelo o de roca situados inmediatamente debajo de la estructura no son capaces de soportar la carga, con la adecuada seguridad o con un asentamiento tolerable. Como la cimentacin superficial, tambin la cimentacin profunda, incluyendo los estratos de suelo o roca situados debajo, deben ofrecer seguridad y no asentarse excesivamente por efecto de las cargas de la estructura que soportan. (Terzagui, 1973).

Existen dos tipos de cimentaciones profundas generalmente aceptadas: pilotes y pilares. En este proyecto nos evocaremos nicamente al uso de pilotes.

2.1.1.- PILOTES

Los pilotes son fustes relativamente largos y esbeltos que se introducen en el terreno. Son miembros estructurales con un rea de seccin transversal pequea, comparada con su longitud, y usualmente se instalan usando una piloteadora que tiene un martinete o un vibrador. A menudo se hincan en grupos o en filas, conteniendo cada uno suficientes pilotes para soportar la carga de una sola columna o muro.

Segn Bolognesi (2002) las columnas con poca carga, pueden en algunos casos, necesitar un solo pilote. Sin embargo, ya que en las condiciones del trabajo de campo, la posicin real de un pilote puede quedar a varios centmetros de la posicin proyectada, difcilmente pueden evitarse las cargas excntricas. En consecuencia las cabezas de los pilotes aislados usualmente se arriostran en dos direcciones por medio de contratrabes.

Si solo se necesitan dos pilotes las cabezas se unen con un cabezal de hormign, siendo arriostradas solamente en una direccin, perpendicular a la lnea que une los dos pilotes. Los grupos que contienen tres o ms pilotes, estn provistos de cabezales de hormign reforzado y se consideran estables sin apoyarlos con contratrabes. (Bolognesi, 2002).

2.1.2.- USO DE LOS PILOTES

Los pilotes se usan de muchas maneras. Los pilotes de carga que soportan las cimentaciones son los ms comunes. Estos pilotes transmiten la carga de la estructura a travs de estratos blandos a suelos ms fuertes e incompresibles o a la roca que se encuentre debajo, o distribuyen la carga a travs de los estratos blandos que no son capaces de resistir la concentracin de la carga de un cimiento poco profundo. Los pilotes de carga se usan cuando hay peligro de que los estratos superiores del suelo puedan ser socavados por la accin de las corrientes o las olas, o en los muelles y puentes que se construyen en el agua. (Pilotaje, 2003).

Segn Bowles (1996) los pilotes de traccin se usan para resistir fuerzas hacia arriba, como en las estructuras sometidas a supresin; tales son los edificios cuyos basamentos estn situados por debajo del nivel fretico, las obras de proteccin de presas o los tanques sorterrados. Tambin se emplean para resistir el vuelco en muros y presas, y como anclaje de los cables que sirven de contravientos en las torres o retenidas en los muros anclados y en las torres.

Los

pilotes

cargados

lateralmente

soportan

las

cargas

aplicadas

perpendicularmente al eje del pilote y se usan en cimentaciones sometidas a fuerzas horizontales, como son los muros de sostenimiento de tierras, los puentes, las presas y los muelles, y como defensas y duques de alba en las obras de los puertos. Si las cargas laterales aplicadas en la cabeza superior del pilote son grandes, los pilotes inclinados pueden resistirlas ms eficazmente. Estos son pilotes que se hincan con un cierto ngulo y frecuentemente se usa una combinacin de pilotes verticales e inclinados. (Gunaratne, 2006).

Los pilotes se usan algunas veces para compactar el suelo o como drenes verticales en estratos de baja permeabilidad. Los pilotes colocados muy prximos unos de otros y las tablestacas anchas y delgadas unidas entre si, se usan como muros de sostenimiento de tierras, presas temporales o mamparos contra filtraciones.

2.1.3.- INSTALACIN DE PILOTES

Existen dos formas de usar un pilote en una construccin, una es la de introducir un pilote fabricado en hormign, de tubos o perfiles de acero o de madera, y la otra es el uso de pilotes de concreto fabricados In-Situ.

La operacin de introducir el pilote en el terreno se llama hinca del pilote, lo que se realiza mediante el uso de un martinete, en donde el pilote se coloca, entre guas que mantienen la verticalidad del pilote, y debajo del martillo, el cual funciona generalmente con energa mecnica.

La caracterstica ms importante de un equipo para la hinca de pilotes, es su capacidad para guiar el pilote exactamente. Debe ser lo suficientemente fuerte y rgido para mantener el pilote y el martillo en su posicin, y con la inclinacin fijada, a pesar del viento, las obstrucciones bajo el terreno y el movimiento del martillo. (Das, 2001).

Segn bolognesi (2002) los pilotes de hormign fabricados "in situ" son los que ms se usan para cargas entre 30 y 60 toneladas. Estos tipos de pilotes se pueden dividir en dos grupos: pilotes con camisa o tubo de entibacion en los que un tubo de metal de paredes delgadas se hinca en el terreno y sirve de molde, y pilotes sin tuvo de estibacin, en los que el hormign se coloca en un agujero hecho previamente en el suelo, quedando el hormign finalmente en contacto directo con el suelo. Para el desarrollo de este proyecto, y para los efectos del clculo de esfuerzos y reacciones, se considera el uso de pilotes hincados.

2.1.4.- COMPORTAMIENTO DEL PILOTE DURANTE LA HINCA

En suelos muy blandos los primeros golpes del martillo pueden hincar el pilote varios metros; de hecho el pilote puede introducirse en el terreno bajo el peso del martillo solamente; sin embargo en los suelos duros cada golpe del martillo esta acompaado por una deformacin del pilote y la consiguiente prdida de energa.

El golpe del martillo produce inicialmente un movimiento del pilote hacia abajo, pero este es seguido por un rebote que representa la compresin elstica temporal del pilote y del suelo que lo circunda. La penetracin promedio para varios golpes se puede hallar de la resistencia a la hinca, que es l nmero de golpes necesario para hincar el pilote una distancia determinada, generalmente 2.5, 15 o 30 cm. (Pilotaje, 2003).

Cuando el pilote es muy largo y la hinca difcil, el comportamiento del pilote es ms complejo. En el momento del impacto la parte superior del pilote se mueve hacia abajo; la parte inmediatamente debajo se comprime elsticamente y la punta del pilote permanece momentneamente fija. La zona de compresin se mueve rpidamente hacia abajo y alcanza la punta del pilote una fraccin de segundo despus de producirse el impacto. Como resultado de esta onda de compresin, la totalidad del pilote no se mueve hacia abajo en un instante, sino que lo hace en segmentos ms cortos. (Gunaratne, 2006).

2.1.5.- EFECTOS DEL PILOTE EN EL SUELO

La forma de distribucin del esfuerzo, el asentamiento y capacidad mxima de una cimentacin por pilotaje, depende del efecto del pilote en el suelo. El pilote, representado por un cilindro de longitud L y dimetro D, es una discontinuidad en la masa de suelo, que reemplaza el suelo, segn sea instalado por excavacin, como un pilar, o por hinca. La excavacin altera el suelo cambiando la forma de distribucin del esfuerzo; el suelo puede ser comprimido hacia adentro, desorganizndose la estructura de las arcillas y reducindose la compacidad de las arenas. Al forzar un pilote dentro del agujero o al colocar concreto fresco, puede que se fuerce parcialmente el suelo hacia afuera, originndose ms alteracin. (Das, 2001).

La hinca del pilote origina aun mayor alteracin. El desplazamiento producido por la hinca de los pilotes tiene dos efectos. Primero, se produce un levantamiento del terreno en los suelos de arcilla saturada y en los no cohesivos compactos; Segundo, se establece una fuerte presin lateral en el suelo. (Pilotaje 2003).

Segn Bowles (1996) la hinca de pilotes con martillo produce choque y vibracin que se transmite, a travs del terreno, a las estructuras contiguas. Si el suelo es arena muy suelta, fina y saturada, las vibraciones pueden causar una licuefaccin temporal de la misma, con la correspondiente perdida de capacidad de carga, producindose graves daos; aunque esto raramente ocurre. El pilote transfiere la carga al suelo de dos maneras. Primero, por la punta, en compresin, que se llama resistencia por la punta, y segundo, por esfuerzo cortante a lo largo de su superficie lateral, llamado comnmente friccin lateral (aunque una verdadera friccin no se desarrolla en todos los casos). Los pilotes hincados a travs de estratos dbiles hasta que su punta descanse en un estrato duro, transfieren la mayor parte de su carga por la punta y algunas veces se les llama pilotes resistentes por la punta. En suelos homogneos los pilotes transfieren la mayor parte de su carga por friccin lateral y se les llama pilotes de

10

friccin o pilotes flotantes; sin embargo, la mayora de los pilotes desarrollan ambas resistencias. (Pilotaje, 2003).

2.1.6.- CAPACIDAD DE CARGA DE PILOTES

2.1.6.1.- ANLISIS ESTTICO

La mxima capacidad de carga de un pilote o pilar es la suma de la resistencia por la punta y por la friccin lateral en el instante de la carga mxima.

En el pilote resistente por la punta, esta se asemeja a una cimentacin por superficie enterrada profundamente. Segn Bowles (1996) cuando se carga el pilote se forma un cono de suelo no alterado que se adhiere a la punta. Como la punta va penetrando mas profundamente conforme aumenta la carga, el cono fuerza el suelo hacia los lados cortando la masa a lo largo de una superficie curva. Si el suelo es blando, compresible o tiene un modulo de elasticidad bajo, la masa situada mas all de la zona de esfuerzo cortante se comprime o deforma, permitiendo que el cono penetre mas.

Esta es una forma de esfuerzo cortante local similar al descrito para las cimentaciones poco profundas. Si el suelo o la roca son muy rgidos, la zona de esfuerzo cortante se extiende hasta que el desplazamiento total permita al cono perforar el suelo hacia abajo.

Segn Bowles (1996) existen diferentes factores de capacidad de carga que se usan para resolver determinados problemas. Los factores para cimentaciones poco profundas se aplican a los pilotes resistentes por la punta, cuando descansan en estratos duros, y sobre los cuales se encuentran formaciones dbiles. Tambin se aplican a los pilotes embebidos en arcillas blandas y arenas sueltas. Los factores ms altos se aplican solamente a las arcillas mas duras y a las arenas muy compactas, en

11

las cuales la punta del pilote queda embebida a una profundidad de 10D. Los factores para condiciones intermedias entre estos lmites, se pueden hallar por interpolacin.

Si los pilotes son hincados en el suelo, el ngulo de friccin que debe usarse es el que se obtiene despus de la hinca. Segn Terzagui (1973) en las arenas se produce un aumento de 2 a 5 grados sobre el valor obtenido antes de la hinca. Si la colocacin del pilote se hace con chifln de agua o con perforacin previa, el ngulo no cambia prcticamente. La friccin lateral que acta a lo largo del fuste del pilote es igual a la suma de la fraccin de adherencia en la superficie del pilote, ms la resistencia al esfuerzo cortante del suelo inmediatamente adyacente al pilote. Si f es la friccin lateral: f = c' + 'h tan f = ca + 'h tan Donde ca es la adherencia y tan el ngulo de rozamiento entre el suelo y la superficie del pilote.

Los valores de ca y de tan se pueden determinar por una prueba directa de esfuerzo cortante, substituyendo una mitad de la caja del aparato para el ensayo por el material de la superficie del pilote. Los ensayos hechos con pilotes de tamao natural indican los siguientes valores de ca son relacin a los de c en arcillas saturadas, obtenidos en pruebas de esfuerzo cortante sin drenaje: ca = 0.9c c < 0.5 Kg/cm2 (6b) ca = 0.9 + 0.6(0.49c -1) c > 0.5 Kg/cm2 (6c)

La capacidad total de carga de un pilote es nominalmente, la suma de la resistencia por la punta que se ha movilizado y el producto de la friccin unitaria lateral movilizada por el rea de la superficie lateral del pilote.

12

Sin embargo, la carga mxima o de falla, no es necesariamente igual a la suma de la resistencia mxima por la punta y la friccin lateral mxima. Primero, porque puede que no se movilicen simultneamente las resistencias por la punta y por friccin lateral en las diferentes secciones del fuste del pilote. (Pilotaje, 2003).

El pilote hincado tiene, generalmente, una capacidad de carga mxima mayor que la del pilote colocado con excavacin previa o con chifln de agua, porque los valores mximos, tanto el de resistencia por la punta como el de friccin, se alcanzan durante la hinca.

Una segunda causa de la diferencia entre la capacidad de carga calculada y la real de los pilotes, proviene de la friccin negativa. Segn Bowles (1996) los esfuerzos que se desarrollan en el suelo por el pilote y por cualquier carga superficial, como el relleno, no soportada directamente por los pilotes, hace que el suelo no consolide. Si hay algn estrato muy compresible a algn nivel por arriba de la punta del pilote, la consolidacin har que el suelo de arriba se mueva hacia abajo con respecto al pilote.

Estos estratos en vez de soportar el pilote, debido a su movimiento descendente, aaden carga al pilote. Esta friccin negativa ha sido tan grande en algunos casos, que ha producido la falla de la cimentacin por pilotaje y se debe por lo tanto, considerar en el proyecto al momento de anlisis de los programas.

El mtodo mas seguro para determinar la capacidad de carga de un pilote, para la mayora de los lugares, es la prueba de carga. Los ensayos de carga se hacen para determinar la carga mxima de falla de un pilote o grupo de pilotes o para determinar si un pilote o grupo de pilotes es capaz de soportar una carga sin asentamiento excesivo o continuo. (Bolognesi, 2002).

Los ensayos de carga se pueden hacer construyendo una plataforma o cajn en la cabeza del pilote o grupo de pilotes, en la cual se coloca la carga, que puede ser arena, hierro, bloques de concreto o agua. Para hacer un ensayo ms seguro y ms fcilmente controlable, se usan, para aplicar la carga, gatos hidrulicos de gran

13

capacidad cuidadosamente calibrados. La reaccin del gato ser tomada por una plataforma cargada o por una viga conectada a pilotes que trabajaran a traccin.

2.1.6.2.- ANLISIS DINMICO

Tericamente, es posible establecer alguna relacin entre la capacidad de carga del pilote y la resistencia que ofrece a la hinca con un martillo.

Este anlisis dinmico de capacidad de carga del pilote, se ha usado anteriormente en reiteradas ocasiones. En algunos casos, estas formulas han permitido predecir con exactitud la capacidad de carga del pilote, pero en otros no, lo que ha provocado el colapso de algunas estructuras. Segn Gunaratne (2006) la carga que recibe el pilote, y la "falla" producida por la hinca con martillo, ocurren en una pequea fraccin de segundo, mientras que en la estructura la carga se aplica en un lapso que varia de horas a aos. Una relacin fija entre la capacidad de carga de un pilote obtenido dinmicamente y la obtenida a largo plazo, solo puede existir en un suelo en el que la resistencia a esfuerzo cortante sea independiente de la velocidad de aplicacin de la carga. Esto es aproximadamente cierto en un suelo no cohesivo seco y en suelos no cohesivos hmedos de compacidad intermedia o de granos tan gruesos, que el esfuerzo cortante no produzca un esfuerzo neutro apreciable.

En las arcillas y en los suelos no cohesivos, ya sean muy sueltos o compactos, de granos finos y saturados, la resistencia depende de la velocidad con que se desarrolla el esfuerzo cortante; en estos suelos el anlisis dinmico puede no tener validez.

14

2.1.7.- CARGAS LATERALES

Un pilote vertical cargado lateralmente se flexa como una viga en voladizo parcialmente empotrada. Si las cargas son pequeas, la resistencia del suelo es bastante elstica.

Segn

Gunaratne

(2006),

esta

resistencia

se

puede

representar

aproximadamente, suponiendo que el suelo reacciona como una serie de resortes horizontales, cuya rigidez se puede expresar como un modulo de reaccin k. La ecuacin diferencial de la flexin de la viga se puede resolver para flechas y momentos, as como la presin del suelo, por aproximaciones sucesivas o por aproximacin de elementos finitos. Tales soluciones se encuentran en forma grfica dimensional para supuestas variaciones de k. Las curvas desarrolladas por J. E. Bowles, dan los valores de las flechas, momentos flexionantes y presiones del suelo para un pilote de rigidez constante y para un k aumentado linealmente con la profundidad: k = k'z. La mejor manera de determinar los valores de k o k' de la formacin del suelo, es a travs de un ensayo de carga de un pilote de tamao natural.

Las curvas se expresan en trminos de la rigidez relativa, T en cm.: T=(E*I)*0.2 (k') E e I se refieren a la seccin transversal del pilote. La profundidad se expresa por el coeficiente adimensional Z. Z = L/T

Si la carga lateral es suficientemente grande, la presin en el suelo exceder la resistencia del mismo y el pilote fallara. Algunas veces se calcula la falla de resistencia como falla del empuje pasivo de tierra en la parte superior del pilote.

15

Lo indicado anteriormente no es real, porque en teora de empuje pasivo, la deformacin por cortante se supone bidimensional o plana mientras que el pilote cargado lateralmente fallara por esfuerzo cortante tridimensional a una presin que excede el empuje pasivo generado. (Gunaratne, 2006; Pilotaje 2003).

Los resultados de ensayos de carga tpicos indican que cuando los pilotes verticales estn embebidos totalmente en el suelo, pueden soportar cargas laterales de solamente 1/10 a 1/5 de su capacidad de carga vertical, sin una flecha excesiva (menos de 1.25 cm.). Si es necesaria una mayor rigidez o resistencia lateral, se deben emplear pilotes inclinados.

16

2.2.- MTODO DE ELEMENTOS FINITOS

El mtodo de Elementos Finitos (MEF) es un mtodo numrico avanzado que permite obtener una aproximacin de la solucin de un problema de contorno, asociado a una ecuacin diferencial, ordinaria o en derivadas parciales, bajo ciertas condiciones de frontera. Este mtodo consiste bsicamente, en aproximar la solucin de un problema de frontera de clase C2, por la solucin del problema equivalente planteado sobre un subespacio de dimensin finita. Usualmente la base de este espacio es generado por funciones lineales, que en el caso de mejorar la precisin de la solucin se tendra que realizar un refinamiento de malla, lo que conduce a la bsqueda de algoritmos de convergencia rpida para la solucin de grandes sistemas de ecuaciones lineales. (Obando, 2004).

Este mtodo constituye un mtodo numrico destinado a resolver mediante ecuaciones matriciales las ecuaciones diferenciales que se plantean en sistemas discretos para estructuras, o continuos para campos.

En

general,

el

mtodo

de

los

Elementos

Finitos,

se

programa

computacionalmente para calcular el campo de desplazamientos y posteriormente a travs de relaciones cinemticas y constitutivas, las deformaciones y tensiones, respectivamente, cuando se trata de un problema de mecnica de slidos deformables o mas generalmente un problema de mecnica de medios continuos.

Una importante propiedad del mtodo es la convergencia, si se consideran particiones de elementos finitos sucesivamente mas finas, la solucin numrica calculada converge rpidamente hacia la solucin exacta del sistema de ecuaciones planteada, lo cual se puede calcular rpidamente a travs de un software computacional. (Fras, 2004).

17

2.3.- FORTRAN

El FORTRAN (FORmula TRANslator) es el primer lenguaje cientfico de programacin, siendo desarrollado a fines de la dcada de los 50, pero es an ampliamente utilizado en aplicaciones cientficas y de ingeniera. La ltima versin de FORTRAN actualmente es la versin FORTRAN 90, que reemplaz a la versin FORTRAN 77 a comienzos de la dcada de los 90. El FORTRAN 90 incorpora nuevas estructuras y caractersticas en su sintaxis, si bien an es compatible con la anterior versin.

En los primeros tiempos de la informtica cualquier clculo que implicara la evaluacin de frmulas matemticas haba de hacerse mediante complicados programas, que traducan esas frmulas al lenguaje del ordenador, muy primitivo. EL FORTRAN es un lenguaje de alto nivel, orientado a facilitar las cosas al usuario, y que permite escribir el programa de manera casi idntica a como se escriben las frmulas en un papel. (Velasco, 2001).

En FORTRAN las instrucciones que el usuario transmite al ordenador se encuentran en un fichero de texto (en cdigo ASCII), que es necesario escribir (utilizando un editor de texto cualquiera). El nombre del fichero sigue las convenciones del sistema operativo del ordenador. En DOS un programa FORTRAN suele tener la extensin '.FOR', por ejemplo, 'CALCULO.FOR'; en Unix es tradicional que la extensin sea '.f', por ejemplo, 'calculo.f'. Los nombres de los programas son libres (salvo ciertas limitaciones de longitud y tipos especiales de caracteres que no pueden aparecer en el nombre del fichero).

Al contrario que el BASIC, que es un lenguaje interpretado, el FORTRAN es un lenguaje compilado. Esto significa que, una vez escrito el programa, ste ha de ser traducido en bloque al lenguaje mquina, o sea, el lenguaje que entiende el procesador del ordenador, mediante un proceso llamado compilacin. (Fortran).

18

2.4.- MATHCAD13

Mathcad es un entorno de documentacin tcnica con prestaciones de clculo numrico y simblico, que permite explorar problemas, formular ideas, analizar datos, modelar y chequear escenarios, determinar la mejor solucin... y finalmente documentar, presentar y comunicar los resultados.

Es un programa algebraico de computadora similar a Mathematica, distribuido por Mathsoft. A diferencia de Mathematica, MathCad es ms intuitivo de usar, permite el uso de plantillas de funciones en las que solo es necesario escribir los valores deseados, incluso para graficar funciones.

Algunas de las capacidades matemticas de Mathcad estn basadas en parte del cdigo del programa algebraico Maple (Ncleo MathSoft de Maple o Mathsoft Kernel Maple, MKM). Mathcad se encuentra organizado como una hoja de trabajo, en las que las ecuaciones y expresiones se muestran grficamente, no como simple texto.

Mathcad es una herramienta tcnica de gran flexibilidad que incorpora todas las funcionalidades de interactividad propias de las hojas de clculo, una gran potencia debida a su extensa librera de funciones y una completa integracin con muchas de las herramientas utilizadas en ingeniera (Excel, AutoCAD, MATLAB, Visio, ODBC, etc.). Sus prestaciones de documentacin tcnica, junto con una notacin matemtica clsica, permiten que los clculos sean auto-documentados y ahorrar as mucho tiempo al ingeniero a la hora de crear informes tcnicos sobre los clculos realizados.

Mathcad incorpora adems un conjunto de herramientas software adicionales que le permiten realizar simulacin de sistemas dinmicos, CAD 2D paramtrico y grficos tcnicos. Su precio le convierte en una herramienta muy asequible que puede ser adoptada por las empresas de ingeniera como una herramienta estndar o corporativa.

19

2.4.1.- CARACTERSTICAS PRINCIPALES Mathcad contiene este software, estn: - Calcular, modelar y visualizar las ideas tcnicas. - Resolver clculos y visualizar los resultados. - Documentar el trabajo (incluidos los clculos) con notacin cientfico-matemtica. - Integrar datos entre programas y sistemas: Microsoft Office (Excel, Visio y otros productos de la familia Office), AutoCAD, Axum, SmartSketch, VisSim, MATLAB, bases de datos ODBC (incluidas Microsoft Access y FoxPro). - Distribuir documentos y colaborar a travs de redes. - Crear estndares a travs de toda la organizacin: Mathcad puede ser utilizado por diseadores, tcnicos y directores. herramientas complementarias para desarrollar proyectos

tcnicos desde el concepto hasta su distribucin final. Dentro de las ventajas que tiene

2.4.2.- AREAS DE APLICACIN El amplio abanico de herramientas disponibles en y para Mathcad lo convierten en un entorno de resolucin de problemas en las siguientes reas: Ingeniera mecnica y termodinmica. Ingeniera civil y arquitectura. Ingeniera elctrica y teora de control. Procesado de seales digitales. Ingeniera qumica. Matemtica y estadstica. Finanzas y economa.

20

CAPITULO III
ESTUDIO Y ANALISIS DE LOS PROGRAMAS EN FORTRAN IV

21

3.1.- PROGRAMACION EN FORTRAN Se eligieron dos programas del texo original, el primero es LATERAL PILES, y el segundo, STATIC PILES. Los programas creados por J. E. Bowles, LATERAL PILES y STATIC PILES, se encuentran en el ANEXO A de este proyecto, escritos en FORCE 2.0. El desarrollo de los programas hechos en el software de programacin FORTRAN IV, LATERAL PILES y STATIC PILE, consta del siguiente orden en los procedimientos de ejecucin:

3.1.1.- Identificacin de las variables Se refiere a la clasificacin de los parmetros que se ingresan en los programas, tipo de nmero, si son vectores, si forman matrices, etc., en donde se deben clasificar y definir todas las variables que sern incluidas en el programa como valores reales o enteros. Para la clasificacin de los parmetros utilizados en la programacin, se utilizan los siguientes comandos: DIMENSION: Declaracin que indica el tipo de variable DOUBLE PRECISION: indica el tipo de aproximacin de la variable

3.1.2.- Lectura de las variables En esta etapa, el programa reconoce o lee los valores ingresados al programa para ser ejecutado, asignando a cada valor una letra que la identifique en el desarrollo del programa.

Paralelamente, se muestran los valores ingresados, a travs de un operador de salida, el cual ensea los datos con su leyenda respectiva.

22

La lectura de las variables se realiza a travs del siguiente comando: READ(1,10) donde 1 seala un archivo en donde estn definidas todas las variables que definen un determinado problema a analizar con el programa, las cuales deben ser ingresadas previamente a la ejecucin del programa; 10 seala una posicin en donde existe un determinado formato de escritura o lectura, el cual determina la cantidad de espacios a utilizar por variable, cantidad de decimales, altura de nmero, etc. Cabe sealar que el formato en el cual lee el archivo debe tener el mismo orden para ingresar los datos del problema.

Para visualizar los datos el programa utiliza la siguiente sentencia: WRITE(3,10)

En donde 3 seala la hoja o pantalla en donde sern mostrados los datos ingresados al problema, y 10, indica la posicin en donde se encuentra el formato a utilizar por el programa al mostrar los datos.

3.1.3.- Clasificacin de parmetros

En esta etapa, el programa clasifica los datos entregados, utilizando expresiones lgicas, sentencias, y frmulas predeterminadas, las cuales determinan el uso o no uso de ciertas frmulas para la obtencin de los resultados, o tambin pueden determinar el mtodo de clculo predeterminado en la programacin, siendo algunas de estas las siguientes:

a.LT.b: a es menor que b a.GE.b: a es mayor o igual que b Agregando la sentencia IF: IF(a.GE.b)THEN: si a es mayor o igual que b, entonces

23

En donde THEN nos indica que se debe cumplir con la afirmacin que le sigue, o cierto valor debe ser igual a otro, etc.

3.1.4.- Generacin de parmetros de clculo

Una vez que el programa clasifica los parmetros y determina las frmulas y procedimientos a seguir, se comienza con el mtodo de clculo que esta programado, que para el caso de programas de ingeniera, son los mtodos de clculo a travs de Elementos Finitos, los que generan matrices de nxn, las cuales a travs de iteraciones acotan los valores encontrados, convergindolos a los resultados exactos. A modo de ejemplo: DO 21 I = 1,N IF(a.LT.b)GO TO 21 21 .. Esto quiere decir: el valor de I, se hace correr desde el valor I= 1, hasta el valor I= N, en donde cada iteracin obliga al programa a verificar si el valor de a es menor que b, y si es as, el programa continua ejecutndose en 21, cuya posicin indica el paso siguiente en la ejecucin del programa, si no, el programa asigna un nuevo valor a I, cuantas veces sea necesario hasta que se cumpla la sentencia exigida.

Cabe sealar que las matrices generadas, son mostradas por el programa, seguidas de su respectiva sentencia (leyenda explicativa), a travs de la siguiente sentencia:

WRITE(3,*)(I,J) En donde la matriz con los valores de (I,J) son mostrados en la pantalla 3, sealada anteriormente.

24

3.1.5.- Obtencin de los resultados La obtencin de los resultados, generalmente es a travs de la interaccin de las matrices generadas anteriormente, cuyos resultados son mostrados como se sealo anteriormente. En resumen, se tiene el procedimiento tpico de programacin para resolver un problema de ingeniera a travs del mtodo de elementos finitos, en donde claramente se puede apreciar la declaracin de los parmetros y valores, junto con el ingreso de datos, la generacin de las sentencias y frmulas a emplear, as como tambin las matrices generadas, a travs de los procedimientos de iteracin, y finalmente, se aprecian los clculos y resultados finales.

FORTRAN, como se sealo anteriormente, se basa en una programacin demasiado estricta en cuanto a la declaracin de variables y a su entrega de resultados, en comparacin con otros software mas recientes, en donde la metodologa de programacin es menos estricta, y a la vez, presenta formatos de entrega de datos mucho mas completos que FORTRAN, el cual es completamente operado en el sistema MS-DOS, ya sea para la programacin, como para la ejecucin de los programas.

Otro factor importante que sealar, es que para el ingreso de valores o datos de un determinado problema a analizar en el programa, este, se debe realizar siguiendo en forma muy estricta el formato especificado en el programa, ya que cualquier mnimo detalle conduce a un error en la ejecucin de este.

Segn lo estudiado y modificado para la compilacin de estos programas, nos damos cuenta que el programa de ingreso de datos debe ser programado bajo el mismo formato para cada valor a ingresar; adems se debe hacer mencin a errores de compilacin encontrados al momento de ejecutar los programas, los cuales debieron ser chequeados y mejorados para la ejecucin de estos, pero al no tener el programa de ingreso de datos, la ejecucin de los programas quedaba en carcter de Inactivo.

25

3.2.- PROBLEMAS DE EJECUCIN

A continuacin se sealan los principales problemas encontrados en la programacin original de LATERAL PILES y de STATIC PILES para la compilacin de datos, los cuales fueron arreglados y corroborados con la ayuda de profesores conocedores de FORTRAN, los cuales asistieron en los arreglos:

- Ambos programas no contaban con el manejador o programa principal, que es la sentencia PROGRAM, seguida del nombre del programa, la cual indica el punto de partida del programa principal, diferencindolo de las subrutinas o subprogramas incluidos en el programa principal. - El segundo problema arrojado por el compilador, se refiere a la base de los datos correspondientes a un problema especfico a revisar, cuyo formato de programacin no se encontr en el texo gua, y que bien, pudiendo programar dicha entrada de datos, el reconocimiento de estos valores por parte del programa principal hubiera resultado complicada, debido a los estrictos formatos existentes en la programacin.

A modo de ejemplo, sealamos el siguiente ejemplo, correspondiente a la lectura de datos de STATIC PILES:
7000 READ(1,1000,END=6000)TITLE,UT1,UT2,UT3,UT4,UT5,UT6,UT7,UT8 1000 FORMAT(20A4/4(A4,6X),4(A8,2X)) READ(1,101)NP,NM,NNZF,NC,IPRD,NLC,JTSOIL,NDELT,JJS,LIST,NSTRPT,IU 101 FORMAT(12I5)

Claramente, se observa que la lectura de datos sigue un patrn muy minucioso, que conlleva a programar una entrada de datos demasiado minuciosa, que ocupara valioso tiempo, que puede emplearse en la reprogramacin de estos en un software mucho mas amigable y de mejores caractersticas.

26

- En la formacin de la matriz ASAT, se vio el siguiente problema:

C**** END OF ASAT FORMATION--INVERT ASAT MATRIX--USEE GAUSS-JORDAN METHOD 9168 DO 25 K=1,N DO 20 J=1,N 20 IF(J.NE.K)E(K,J)=E(K,J)/E(K,K) DO 21 I=1,N IF(I.EQ.K) GO TO 21 DO 21 J=1,N IF(J.EQ.K) GO TO 21 E(I,J)=E(I,J)-E(K,J)*E(I,K) 21 CONTINUE

En donde el programa asigna un mismo procedimiento a dos variables distintas, las que varan a travs de iteraciones por separado, lo cual generaba en el programa una iteracin sin avance, solo tomando diversos valores y llevando a una no ejecucin del programa.

La solucin a este problema fue la siguiente:

C**** END OF ASAT FORMATION--INVERT ASAT MATRIX--USEE GAUSS-JORDAN METHOD 9168 DO 25 K=1,N DO 20 J=1,N 20 IF(J.NE.K)E(K,J)=E(K,J)/E(K,K) DO 210 I=1,N IF(I.EQ.K) GO TO 210 DO 21 J=1,N IF(J.EQ.K) GO TO 21 E(I,J)=E(I,J)-E(K,J)*E(I,K) 21 CONTINUE 210 CONTINUE

En donde se cerr la iteracin de cada variable, haciendo posible el desarrollo posterior de la programacin.

27

3.3.- CRITERIO DE PROGRAMAS Por lo sealado anteriormente, se tom la determinacin de reprogramar ambos programas, LATERAL PILES y STATIC PILES, en un software ms avanzado que FORTRAN, que sea de un lenguaje de programacin ms simple, y que entregue resultados en un formato y en una pantalla de mejores caractersticas que FORTRAN, cuya pantalla de ejecucin es en el sistema operativo MS-DOS.

Por lo sealado anteriormente, y dadas sus caractersticas, propiedades de programacin, y entrega de datos, el software que se utilizara para programar nuevamente LATERAL PILES y STATIC PILES, ser el software conocido como MATHCAD13, el cul se menciono en el captulo anterior.

28

CAPITULO IV
PROGRAMACION

29

4.1.- ESTUDIO Y ANALISIS DE LA PROGRAMACION

En esta etapa, se realiza el estudio y anlisis de un pilote para dos casos diferentes de carga, LATERAL PILES para el caso de aplicar una carga lateral en la cabeza superior de un pilote, y STATIC PILES el cual analiza el diseo de un pilote sometido a carga esttica constante.

4.1.1.- LATERAL PILES

4.1.1.1.- Introduccin al Problema

Este programa, calcula la distribucin de esfuerzos y desplazamientos en un pilote que recibe una fuerza de corte externa lateral en su parte superior, la cual puede ser producto de un sismo, o algn tipo de impacto.

El anlisis estructural realizado para obtener los desplazamientos y esfuerzos, se basa en la teora de viga sobre fundacin elstica, que establece la rigidez del pilote como funcin del mdulo del suelo y la rigidez flexional del pilote (EI). Este modelo, permite considerar la variacin de las distintas capas de suelo en el perfil mediante el uso de resortes con diferentes valores de rigidez.

Para obtener la rigidez traslacional, rotacional y traslacional acoplada del pilote, se aplica en la cabeza del pilote, una carga lateral o un momento. La constante de resortes se determina usando un mdulo del suelo Kh que representa la resistencia horizontal del suelo.

Estos resortes aportan una rigidez producida por la reaccin horizontal del suelo, en donde este aporta una distribucin a lo largo del pilote que puede ser constante, lineal, o de carcter parablico. La reaccin horizontal del suelo, se genera mediante el mdulo de reaccin horizontal, el cual, en el caso de no ser constante, es una funcin que varia con la profundidad, generalmente, aumentando su valor directamente proporcional con el aumento de la profundidad Z. 30

Una vez determinadas las rigideces del elemento, se procede a construir la matriz de rigidez, a partir de la particin de un pilote en N-1 elementos, los cuales generan interactuando en conjunto, la matriz global K. Una vez determinada la matriz global K del sistema, se deben determinar, en primer lugar, el vector P, correspondiente a los esfuerzos del pilote, y luego el vector X, correspondiente a las incgnitas de desplazamientos y giros producidos por los esfuerzos, para as, generar la ecuacin general:

[P] = [K ]* [X ]
La cual se debe generar en un programa computacional para encontrar las soluciones del problema.

4.1.1.2.- Estimacin del mdulo de reaccin horizontal del suelo, kh

Como se seal anteriormente, este mdulo de reaccin horizontal, es una funcin que depende de la profundidad Zi, en la cul, se genera un nodo en el elemento. Existen procedimientos dados por diferentes autores, los cuales entregan soluciones aproximadas en cuanto a sus resultados. En este proyecto, se considera el mtodo empleado por el autor J. E. Bowles. Segn Bowles (1996), una manera para estimar el valor de kh para diferentes nodos del elemento, en funcin de la profundidad, es la representada por la ecuacin siguiente:
k h = Ah + Bh * Z n

Donde Ah y Bh son evaluados usando las expresiones de capacidad, obtenindose de la siguiente manera: Ah = Fw1 * C m * C (c * N c + 0.5 * B * N )

Bh = Fw2 * C m * C * * N q
31

En donde Z es la profundidad de evaluacin. Los siguientes son valores sugeridos por Bowles (1996) para las constantes: c B = = = = = = = = Cohesin del suelo en kN/m2 Dimetro del Pilote en m. Densidad del suelo 40 kN/m3 1.5 2.0 0.4 0.6 1, para pilotes de perfiles de acero, y en suelos cohesivos. 1.3 1.7, para pilotes circulares 2.0 4.4, para pilotes circulares.

C Cm n Fw1 Fw2

Fw1, Fw2=

Adems con el valor del ngulo de friccin interna , segn diferentes autores, se pueden determinar las constantes de capacidad, de acuerdo a la tabla 1:

TABLA 1: Factores de capacidad Nc, Nq y N. Fuente: gunaratne (2006)

Con todos estos valores, se puede obtener en funcin de la profundidad, el valor del mdulo de reaccin horizontal del suelo, kh(Z), el cual ser calculado por el programa diseado.

32

En la figura 1, se puede apreciar la distribucin de kh en el suelo, lo cual, genera las reacciones en el pilote:

FIGURA 1: Distribucin de Esfuerzos y reacciones entre un pilote sometido a carga lateral y el suelo; Modelamiento del pilote. Fuente: Elaboracin Propia

En general, se encuentran las rigideces en el elemento para N nodos, y sus valores son incluidos en la generacin de la matriz global de rigidez del sistema.

La resistencia lateral del suelo, para cada elemento, segn Gunaratne (2006), est dada por las siguientes ecuaciones: Para nudos superficiales:

K j = 0.5 * L * B * kh
Para nudos profundos:

K j = L * B * kh
En donde L= B= Longitud tributaria debido a la subdivisin del pilote en m. Ancho o Dimetro del pilote en m.

33

4.1.1.3.- Mtodo de anlisis, matriz de rigidez

Este mtodo conocido como el mtodo de elementos finitos con el cual, se obtendr mediante la particin del elemento en N nodos, la matriz de rigidez general del sistema.

En primer lugar, se discretiza el pilote en un nmero unidimensional de elementos, tal como se muestra en la figura 2:

FIGURA 2.1: Divisin de un pilote en subelemento. Estado de reacciones Fuente: Gunaratne (2006)

En donde se tiene: 1,2,, N P2, P4, P6, P1, P3, P5, X2, X4, X6, X1, X3, X5, K3, K4, K5, Numero de nodo Fuerzas laterales internas en los elementos del pilote, concentradas en los nodos. Momentos internos en los elementos del pilote concentrados en los nodos. deflexin nodal de cada elemento del pilote rotacin nodal de cada elemento del pilote. resistencia lateral del suelo de cada elemento, representada por la rigidez elstica.

34

Basado en las relaciones de rotacin-deflexin del anlisis estructural, se tiene que para el primer elemento a modo de ejemplo, el sistema de ecuaciones matriciales queda definido de la siguiente manera y tomando en cuenta el elemento base de la figura 2.2:

FIGURA 2.2: Elemento base de diseo Fuente: Elaboracin Propia

que es de la forma: P = K*X, donde k es la matriz de rigidez local del elemento.

4.1.1.4.- Matriz global de rigidez para un elemento dividido en N nodos:

Para generar la matriz global de rigidez del pilote, la cual debe ser programada en el software, hay que obtener las matrices de rigidez locales para los elementos, para as, unir los elementos de estas matrices en la matriz global del sistema. Cabe mencionar, que estas matrices de rigidez, fueron estudiadas y obtenidas en los cursos de mtodos matemticos para ingeniera y en el de Computacin para ingeniera estructural.

35

La matriz de rigidez rotacin-deflexin que se ocupa es la siguiente:

4 EI / L 6 EI / L2 k = 2 EI / L 2 6 EI / L

6 EI / L2 12 EI / L3 6 EI / L2 12 EI / L3

2 EI / L 6 EI / L2 4 EI / L 6 EI / L2

6 EI / L2 12 EI / L3 6 EI / L2 12 EI / L3

En donde E es el modulo de elasticidad del pilote, I es su inercia, y L corresponde al largo Li, de cada elemento.

En primer lugar, hay que generar la matriz de rigidez para el primer elemento, lo cual se hace de la siguiente forma: Para el elemento 1:

L1

FIGURA 3: Elemento 1 de la particin de un pilote FUENTE: Gunaratne (2006)

k1,1 k 2 ,1 k1= k 3,1 k 4,1

k1, 2 k 2, 2 k 3, 2 k 4, 2

k1,3 k 2,3 k 3, 3 k 4,3

k1, 4 k 2, 4 k 3, 4 k 4, 4

Para el elemento 2:

K2*X6
P5 L2

P6

K3*X8

P7

P8

FIGURA 4: Elemento 2 de la particin de un pilote FUENTE: Gunaratne (2006)

36

 k 3, 3 k 4,3 k2 = k 5,3 k 6,3

k 3, 4 k 4, 4 + K 2 k 5, 4 k 6, 4

k 3, 5 k 4,5 k 5,5 k 6,5

k 4,6 k 5, 6 k 6, 6 + K 3 k 3, 6

Se debe considerar la rigidez aportada por el suelo a travs de los resortes en cada nodo. para el elemento 3:
P10

K3*X10 K4*X12

P9 P11 L3 P12

FIGURA 5: Elemento 3 de la particin de un pilote FUENTE: Gunaratne (2006)

k 5,5 k 6,5 k3 = k 7 ,5 k 8,5

k 5, 6 k 6,6 + K 3 k 7,6 k 8, 6

k 5, 7 k 6,7 k 7,7 k 8, 7

k 6 ,8 k 7 ,8 k 8,8 + K 4 k 5,8

si por ejemplo, N = 5, se llegara hasta el elemento 4, definido de la siguiente manera:

P13

P14 L4

K4*X14
P15
FUENTE: Gunaratne (2006)

P16

FIGURA 6: Elemento 4 de la particin de un pilote

k 7,7 k 8, 7 k4 = k 9,7 k10,7

k 7 ,8 k 8,8 + K 4 k 9 ,8 k10,8

k 7 ,9 k 8, 9 k 9,9 k10,9

k 7 ,10 k 8,10 k 9,10 k10,10

37

Con estos primeros trminos, para N = 4, nos podemos dar cuenta de la forma que tiene la matriz global K del sistema, para el elemento que se desea evaluar. La correcta ubicacin de cada elemento, es fundamental a la hora de formar la matriz de rigidez global del sistema. De lo anterior, se puede generar una matriz global, a partir de la secuencia encontrada con el anlisis de los elementos anteriores, la cul se seala a continuacin para N nodos, tomando N-1 elementos, y generando una matriz cuadrada del orden de (2N+2)x(2N+2).

La forma de la matriz global de rigidez del sistema para N+2 nodos, se encuentra en el ANEXO B (considerando el elemento de la superficie).

De la matriz de rigidez global del sistema, se puede concluir que para un nmero de nodos N, se generan 2N incgnitas, adems, en otras palabras:

Si: N N-1 2N =5 =4 N de nodos N de elementos

= 10 N de incgnitas

2N+2 = 12 orden de la matriz de rigidez global, o sea K12,12.

Como se defini anteriormente con respecto a los esfuerzos, P1, es el momento generado en la parte superior del pilote; P2, es la fuerza externa horizontal (lateral) solicitante del pilote; P2N-1, corresponde al momento generado en la base inferior del pilote y finalmente P2N corresponde al esfuerzo horizontal en la base inferior del pilote; o sea:

P1 P2

= M1, momento requerido para refrenar la rotacin en la cabeza del pilote. = P, carga axial aplicada en la cabeza superior del pilote (en la superficie)

38

P2N-1 = M2, momento generado en la base inferior del pilote P2N = Pinf, se refiere a la reaccin inferior en la base del pilote

En la figura 7, se aprecia lo descrito anteriormente.

P1 L1

K2 L2 K3 L3

K4 L4

K5

L5

FIGURA 7: Pilote particionado en 4 elementos FUENTE: Gunaratne (2006)

4.1.1.5.- Anlisis de Esfuerzos

Para el anlisis de esfuerzos, se considerara un vector P de orden 2Nx1, el cul, contiene la fuerza externa lateral solicitante, ingresada como dato, adems de estar compuesto por las incgnitas que se definen a continuacin:

39

 M1 P 2 M3 P4 M5 P5 P= M 6 P6 M M M 2 N P2 N

0 2 2 3 3 X = 4 M N 1 N 1 0 0

Para el anlisis de los desplazamientos y giros generados, se considera el vector X de orden 2Nx1, en el cul, se seala el desplazamiento mximo generado directamente por el esfuerzo P2, adems de los giros y desplazamientos generados a medida que se profundiza mas en el pilote. En este vector X, se tienen las siguientes condiciones de borde: X1 = 1 = 0 no existe giro en la cabeza superior

X12 = 6 = 0 se supone el pilote empotrado en la base X11 = 6 = 0 se supone el pilote empotrado en la base En donde la deflexin lateral requerida es .

Finalmente se tiene el sistema sealado anteriormente:

[P] = [K ]* [X ]
El programa que se cree, debe ser capaz de resolver esta ecuacin matricial, ingresando por supuesto las incgnitas y los parmetros que caracterizan un determinado problema de pilotes.

40

4.1.2.- STATIC PILES

Este programa, verifica la capacidad de carga que tiene un pilote. El pilote, transfiere la carga al suelo de dos maneras: Primero, por la punta, en compresin, a lo que se llama resistencia por la punta, y segundo, por esfuerzo cortante a lo largo de su superficie lateral, llamado comnmente friccin lateral. La mxima capacidad ltima de carga Qult de un pilote es la suma de la resistencia ltima por la punta Qult,p, y por la friccin ltima lateral Qult,f, en el instante de carga mxima:

Qult = Qult , p + Qult , f

La cual se observa en la figura 8:

1 1 c1 2 2 c2

h1 h2

5 5

c5

h5

FIGURA 8: Pilote solicitado a una carga axial esttica FUENTE: Elaboracin Propia

Cabe sealar que el peso propio del pilote se debe descontar, debido a que el anlisis de esfuerzos se realiza en la base inferior del pilote, por lo que se tiene la ecuacin de la siguiente forma:

Qult = Qult , p + Qult , f Pp

41

Siendo Pp el peso propio del pilote. Segn Gunaratne (2006) el quiebre de la estructura de un pilote se debe a dos razones; La falla del suelo circundante al pilote, y la excesiva sobrecarga externa del pilote. El diseo de un pilote debe ser econmico para la obra determinada y resistente a las solicitaciones aplicadas, en donde se deben considerar cargas mayoradas para su diseo. Existen dos tipos de mtodos; El primero es el mtodo por tensiones admisibles (ASD), y el segundo es el diseo por resistencia de cargas (LRFD). Segn Gunaratne (2006), para el caso del mtodo de diseo por tensiones admisibles (ASD), se requieren dos condiciones:

1) 2)

Cargas Admisibles Deflexin Admisible

Qapl * Fs Qult
est Tol

En donde Fs, es el factor de seguridad, el cual va entre los valores de 2.5 a 3.5, siendo este intervalo de carcter conservador, Qult es la resistencia ltima en la base inferior del pilote y Qapl, es la solicitacin externa para este. est es la deflexin estimada para un pilote, y Tol, es la deformacin tolerable que puede tener un pilote. La determinacin de la capacidad por punta y por friccin, se realiza a travs del mtodo de Vesic, quien se basa en la teora de expansin de cavidades, basada en parmetros de esfuerzo efectivo.

4.1.2.1.- Estimacin de la capacidad por punta: La capacidad ltima de resistencia por punta, se determina a travs de la siguiente ecuacin:

Qult , p = A p (c * N c + * q * N q )
1 + 2K 0 3

Ap: c:

rea de la seccin del pilote cohesin del suelo

42

Nc y Nq: q: Ko: En donde:

Factores de capacidad de resistencia de fundaciones capacidad de esfuerzo vertical efectivo en Fza/area coeficiente de presin de tierra para asentamiento.

K 0 = 1 sen K 0 = 0.95 sen

para suelos (arenas) granulares para suelos (arcillas) finas.

Los factores de capacidad de resistencia dependen del ngulo de friccin interna , y se determinan segn las tablas 2 Y 3.

TABLA 2: Factor de capacidad de resistencia Nc* dado por Vesic para la estimacin de la capacidad de soporte por punta de un pilote FUENTE: Gunaratne (2006)

TABLA 3: Factores de capacidad de resistencia Nq* dado por Vesic para la estimacin de la capacidad de soporte por punta de un pilote FUENTE: Gunaratne (2006)

43

4.1.2.2.- Estimacin de la capacidad por friccin lateral: La estimacin de soporte por friccin lateral, se determina tomando un elemento del pilote, de altura dh, se realiza el anlisis que se muestra en la Figura 9, y es igual a la suma de la fraccin de adherencia en la superficie del pilote, ms la resistencia al esfuerzo cortante del suelo inmediatamente adyacente al pilote, tal como se menciono anteriormente.

`vi

`Hi
`H
_

`Hi+1 i+1 `vi+1

FIGURA 9: Elemento de un pilote para la determinacin de la capacidad de soporte por friccin lateral FUENTE: Elaboracin Propia

En donde se tienen las siguientes formulas: v: Presin vertical de empuje, en donde v= h : H: densidad del suelo a una capa de altura h, menos 1t/m3 Presin lateral del suelo

H = K 0 * v
K0: coeficiente de presin del suelo, en donde: para suelos (arenas) granulares para suelos (arcillas) finas.

K 0 = 1 sen K 0 = 0.95 sen

44

Presin lateral promedio por elemento:

`H =

`Hi + `Hi +1
2

Friccin entre el suelo y el pilote:

= `H * tan
En donde se tiene que: = para suelos arcillosos

2/3 para suelos granulares La capacidad ltima de resistencia del pilote en cada elemento o por cada altura de capa de suelo, esta dada por la siguiente ecuacin:

qult , fi = h * p* i

En donde la capacidad de soporte por friccin lateral del pilote, se define como:

Qult , f = qult , fi
i =1

Para la determinacin de la capacidad final de soporte lateral se recomienda el uso de una tabla de valores, en donde la primera columna, corresponde a la profundidad en que se encuentra cada capa, seguida por la estimacin de la presin de empuje vertical, seguido por el promedio de esta ltima, luego la determinacin de la presin lateral del suelo, la columna de la friccin entre el pilote y el suelo, para finalizar con la capacidad de soporte por friccin de cada elemento y su correspondiente sumatoria final.

45

Finalmente se tiene:

Qult = Qult , p + Qult , f Pp

En donde la carga aplicada multiplicada por el factor de seguridad FS, debe ser menor que la capacidad de soporte Qult final del pilote, para que se cumpla el diseo por el mtodo ASD o mas bien conocido como el mtodo de tensiones admisibles.

El programa STATIC PILES, nos entrega el valor de esta capacidad mxima de carga, en donde se puede trabajar con la dimensin del pilote hasta que este cumpla con el criterio de diseo establecido.

46

4.2.- PROGRAMACIN EN DETALLE En esta etapa del proyecto, se explica en detalle la programacin realizada para ambos programas. Se debe sealar que en el ANEXO E se encuentra a modo general, el flujo de programacin de LATERAL PILES y de STATIC PILES.

4.2.1.- LATERAL PILES Este programa calcula la deformacin mxima y la distribucin de reacciones en un pilote, producidas por una fuerza solicitante lateral aplicada en su extremo superior. Este programa utiliza el mtodo de elementos finitos para encontrar una matriz de rigidez KT global del sistema generada por la divisin del pilote en 10 nodos.

Propiedades del Pilote

Largo total del elemento Longitud externa superior Modulo de elasticidad Carga lateral

Lt

KPa

P
Desplazamiento

KN

Dimetro del pilote

Longitud de ltimo elemento

Inercia del elemento

mximo del pilote

47

Esta es la entrada de datos del programa correspondiente al pilote.

Largo total del elemento: En esta ventana se debe ingresar el largo total en metros del pilote a analizar; esta medida puede ser modificada cuantas veces sea necesario hasta obtener los resultados que se desean encontrar

Longitud externa superior: En esta ventana, se debe ingresar la longitud del pilote que quedara sobre la superficie del terreno, en caso que el pilote no sea enterrado completamente. Longitud del ltimo elemento: En esta ventana, se debe ingresar la longitud de la ltima porcin del pilote. Esta medida es para efecto de clculos, por lo cual se puede ingresar el valor de cero si se desea.

Dimetro del pilote: En esta ventana se debe ingresar el dimetro o ancho del pilote en metros.

Modulo de elasticidad E: En esta ventana, se ingresa el modulo de elasticidad del elemento en kPa. Inercia del elemento: Se debe ingresar la inercia del pilote en unidades de m4

Carga Lateral: En esta ventana, se debe ingresar el valor de la carga lateral solicitante en kN, la cual esta aplicada en la parte superior del pilote.

Desplazamiento mximo del pilote:

48

En esta ventana se ingresa el valor del desplazamiento mximo admisible que restringe al pilote en su parte superior.

Propiedades del Suelo

Propiedades indice Segn tabla 1 Constantes segn J.E. Bowles Cohesion

N Nq

Nc

Fw1 Fw2 n

Cm C
KN/m3

c
Densidad del suelo

Para modulo de rigidez horizontal del suelo constante ingresar su valor

KN/m3

kh

kN/m3

TABLA 1Factores de capacidad segn Terzagui y Vesic

Esta etapa corresponde al ingreso de datos del suelo:

Propiedades ndice segn tabla 1: En estas ventanas, se deben ingresar estas propiedades, segn el ngulo de friccin interna que se tenga, haciendo doble clic en la flecha correspondiente, aparece la tabla 1, tal como se muestra a continuacin:

Bearing Capacity Factors Terzaghi`s (1943) Expression Nc Nq 5.7 1 0 7.3 1.6 5 9.6 2.7 10 12.9 4.4 15 17.7 7.4 20 25.1 12.7 25 37.2 22.5 30 57.8 41.4 35 95.7 81.3 40 172 173 45

Ny 0 0.5 1.2 2.5 5 9.7 19.7 42.4 100 298

Nc 5.14 6.49 8.34 11 14.8 20.1 30.1 46.4 75.3 134

Vesic`s (1973, 1975) Expression Nq Ny 1 0 1.6 0.4 2.5 1.2 3.9 2.6 6.4 5.4 10.7 12.5 18.4 22.4 33.5 48.1 64.1 109.3 135 271.3

TABLA 4: Factores de capacidad de resistencia dados por Terzagui y Vesic FUENTE: Elaboracin Propia

Si el ngulo de friccin interna se encuentra entre estos valores, se debe tomar el intervalo correspondiente a , interpolando junto con los valores de Nc, Nq o Ny, para encontrar los valores exactos.

49

Constantes segn J. E. Bowles: Estas constantes son recomendadas por J. E. Bowles, para el clculo de pilotes sometidos a una carga lateral y fueron mencionadas anteriormente.

Cohesin: En esta ventana se debe ingresar la cohesin del suelo en unidades de kN/m2.

Densidad del suelo: En esta ventana, se debe ingresar la densidad del suelo en unidades de kN/m3.

Rigidez horizontal kh: Se debe elegir en la ventana desplegable, si esta distribucin de esfuerzos del suelo es constante o variable; si es variable, el programa ejecuta las frmulas correspondientes para determinarlo, pero si es constante, se debe ingresar el valor de esta en la ventana adjunta en unidades de kN/m3.

Procedimiento de clculo del programa En la primera parte, el programa realiza el clculo de la reaccin horizontal del suelo kh, en caso que esta sea variable, utilizando el imput de datos de la siguiente forma:

z :=

for i 1 .. N aux Z
i i

aux

Ah := Fw1 Cm C c Nc + ( 0.5) B ( N )

Bh := Fw2 Cm C Nq

kh :=

for i 1 .. N 1

aux Ah + Bh z i i aux

( n)

50

Obteniendo la rigidez kh del suelo en forma de un vector y correspondiente a cada profundidad zi del suelo

El programa continua realizando el clculo de la distribucin uniforme de rigidez del suelo en el pilote, segn se indica a continuacin:

K :=

for i 1 .. N 1 aux Kh
i

KH :=

if Modulo 1 aux kh aux if Modulo 2 aux K aux

Kt :=

for i 2 .. N 1 aux L
i i+ 1

B KH

aux

for i 1 aux 0.5 L B KH

i i i

aux

Obtenindose as el mdulo de reaccin de rigidez del pilote en cada nodo.

Luego, haciendo doble clic en la pestaa que sigue, se encuentran los valores de z, kh y Kt, en forma de vectores columna

A continuacin, el programa genera la matriz de rigidez global del sistema, tal como se indico anteriormente, uniendo los elementos del pilote, la cual se encuentra en el ANEXO B, debido a su tamao.

En la ventana KT=, del programa, y haciendo doble clic, se puede apreciar la matriz de rigidez global del sistema, la cual es una matriz cuadrada de 22x22.

En la ltima rea del programa, se realiza una subdivisin de la matriz global, para encontrar las soluciones del problema, en forma exacta.

51

El programa una vez que obtiene los valores de estas submatrices, realiza las siguientes operaciones para encontrar los resultados finales:

X2 := K22

P2

( T)

P1 := K12 K22

P2

( T)

P := K22 X2

En donde K22, K12, P2 y X2, corresponden a las submatrices mencionadas anteriormente.

Finalmente, el programa nos muestra los valores del desplazamiento mximo generado en la cabeza del pilote, del generado a nivel del suelo y el momento recuperador que genera la reaccin del suelo, para posteriormente entregar los resultados en una imagen del problema, con los valores encontrados, junto con un grfico que indica el mdulo de reaccin generado para el suelo en donde se analiz el pilote, tal como se indica en la Figura 10.

FIGURA 10: Distribucin de esfuerzos y desplazamientos de un pilote cargado lateralmente FUENTE: Elaboracin Propia

52

4.2.2.- STATIC PILE Este programa verifica el diseo de un pilote a travs de la Capacidad Resistente por punta y por Friccin, para un suelo de hasta 10 capas, y una fuerza F. En la primera parte se encuentra el ingreso de datos, correspondientes al pilote, tal como se indica a continuacin:

1 1 c1 2 2 c2

h1 h2

5 5

c5

h5

Datos del Pilote

Material del pilote Largo total del pilote Solicitacin axial del pilote Si el pilote es de un perfil de acero se debe ingresar el rea y el permetro de la seccin

L
Densidad del material

ton

diametro del pilote

Factor de Seguridad

As p

m2 m

ton/m3

Seccin del pilote Capa de suelo hasta donde llega el pilote

Material de pilote: Esta ventana desplegable nos da tres opciones de seleccin, Hormign, Madera y Acero, los cuales son los materiales que se utilizan para la fabricacin de un pilote.

Densidad del material: En esta ventana, se debe ingresar la densidad del material del cual esta fabricado el pilote, la cul ser utilizada por el programa para calcular el peso propio respectivo del pilote analizado.

53

Largo total del pilote: En esta ventana, se debe ingresar el largo estimado del pilote, tantas veces sea necesario para que se cumplan todos los criterios de diseo que se explicarn mas adelante. Las unidades de longitud deben ser introducidas en metros.

Dimetro del pilote: En esta ventana se debe ingresar el dimetro del pilote en caso que este sea de una seccin circular, o el ancho en el caso que el pilote sea de seccin cuadrada. Esta medida debe ser ingresada en metros.

Seccin del pilote: En esta ventana desplegable, se indican dos opciones, que sealan la seccin del pilote que se esta analizando, ya sea de seccin cuadrada o circular.

Capa de Suelo hasta donde llega el Pilote: En esta ventana, se debe ingresar la capa a la cul se estima que llegar el pilote que se est verificando.

Solicitacin axial del pilote: En esta ventana se ingresa el valor de la fuerza axial solicitante aplicada en la cabeza del pilote, la cul deber ser ingresada en toneladas-fuerza.

Factor de seguridad: Este factor de seguridad es utilizado para mayorar conservadoramente la carga axial solicitante del pilote.

As y P: En el caso que el pilote analizado sea del material de acero, teniendo en cuenta que su seccin puede variar de doble ngulo, doble T, etc., se debe ingresar en la ventana As el rea de la seccin del pilote, y en la ventana P, se debe ingresar su permetro.

54

En la segunda parte se encuentra el ingreso de datos, correspondientes al tipo de suelo en donde se instalaran los pilotes, tal como se indica a continuacin:

Datos del suelo

N de capas del suelo Ver tabla 1 (valores recomendados) Valores correspondientes al angulo de fricin interna de la ltima capa del suelo

Nc

Nq

Tabla 1

N de capas del suelo: En esta ventana desplegable, se debe seleccionar el nmero de capas que entrega la estatigrafa del suelo, las que para efectos de diseo del programa se selecciona de 1 a 10 capas.

Nc y Nq: Estos son propiedades ndice de la ltima capa del suelo, y estas dependen del ngulo de friccin interna de la ltima capa del suelo. Estos valores recomendados por Vesic, se encuentran haciendo doble clic en la pestaa que corresponde, en donde aparece la siguiente tabla:
TABLA 1

N*c y N*q Segn Vesic N*c N*q 0 6.97 1 5 8.99 1.79 10 11.55 3.04 20 18.83 7.85 30 30.03 18.34 35 37.65 27.36 40 47.04 40.47 45 53.66 59.66
TABLA 5: Factores de capacidad de resistencia dados por Vesic FUENTE: Elaboracin Propia

Si el ngulo de friccin interna se encuentra entre estos valores, se debe tomar el intervalo correspondiente a , interpolando junto con los valores de Nc y Nq para encontrar los valores exactos.

El valor de Nc, tambin puede ser calculado en trminos de Nq, segn la siguiente frmula entregada por vesic:
Nq 1

Nc' :=

tan ( 0.01745 ) 2

55

Valor que debe ser ingresado en el casillero correspondiente a Nc, si es que no se cuenta con este valor y se desea aplicar otro criterio.

Ingreso de estatigrafa: Se debe hacer doble clic en la pestaa correspondiente y aparecer el ingreso de datos. En esta seccin del programa, se deben ingresar las caractersticas de cada capa del suelo, su altura hi, su ngulo de friccin interna , su cohesin C en kg/cm2, su densidad en T/m3, y se debe seleccionar el tipo de suelo al que corresponde en la ventana desplegable de la derecha, en donde se elige entre suelo fino (arcilla) y suelo granular (arena), tal como se muestra a continuacin:

Altura de Capas del suelo en m.

Angulo de friccin interna.

Cohesin del suelo en kg/cm2

Densidad del suelo

Tipo de suelo en cada capa

h h h h h h h h h h
10

1 2 3 4 5 6 7 8 9 10

c c c c c c c c c c

1 2 3 4 5 6 7 8 9 10

en T/m3

2 3 4 5 6 7 8 9

2 3 4 5 6 7 8 9

10

Este programa considera cada capa o estrato de suelo como saturada, por lo que automticamente STATIC PILES considera la densidad del suelo menos la densidad del agua. Si la densidad de alguna de las capas ingresadas no esta saturada, se le debe agregar la densidad que se descuenta:

56

Si la densidad de una capa es 1,8 t/m3, y esta no se encuentra saturada, se debe ingresar el valor de 1,8 + 1 t/m3, o sea 2,8 t/m3. A continuacin comienza la etapa de clculo y entrega de datos del programa, en donde se comienza con el clculo de la Capacidad de soporte por punta del pilote:

Clculo de la seccin y del permetro del pilote

Ap :=

if Seccin

1 1
2

P :=

if Seccin

1 1

if Material aux d aux if Material aux d aux if Material aux As aux if Seccin 2 if Material
2

if Material aux

aux 4 d if Material aux if Material 3 aux p aux if Seccin 1 2 if Material 11 + 2 Ko N := d aux 2 3 2 aux if Material 2 2 aux 2 aux if Material aux p aux return aux 3 d 3 2

aux 4 d

aux aux if Material

d2 4

d2 aux 4
aux if Material aux As aux return aux 3

Esfuerzo vertical efectivo

q :=

i=1

1 h ( i ) i

57

Ko :=

for i 1 .. N if S
i

aux 1 sin r i i aux for i 1 .. N if S

i

( )

Qult_p := Ap c Nc + q Nq

(N

aux 0.95 sin r i i aux aux

( )

Esta parte del programa se encuentra bloqueada, para evitar modificaciones o alteraciones por parte de terceros, por lo cul solamente se entrega el valor final del esfuerzo vertical efectivo y la capacidad de soporte por punta: q en T/m2 y Qult_p en Ton. A continuacin se indica la parte del programa que realiza el clculo de la capacidad de soporte por friccin lateral del pilote. Se debe sealar que esta ventana se encuentra bloqueada por las razones explicadas anteriormente:
Clculo de la capacidad de soporte por friccin lateral:

Z := h
1 2 3 4 5 6 7 8 9

1 1 3 4 5 6 7 8 9 2 2 3 4 5 6 7 8 9

Z := h + h

Z := h + Z Z := h + Z Z := h + Z Z := h + Z Z := h + Z Z := h + Z Z := h + Z Z
10

:= h

10

+Z

( ) 1 v 2 := ( 2 1) h + v 1 2 v 3 := ( 3 1) h + v 2 3 v 4 := ( 4 1) h + v 3 4 v 5 := ( 5 1) h + v 4 5 v 6 := ( 6 1) h + v 5 6 v 7 := ( 7 1) h + v 6 7 v 8 := ( 8 1) h + v 7 8 v 9 := ( 9 1) h + v 8 9 v 10 := ( 10 1) h + v 9 10
v 1 := 1 1 h

:=

for i 1 .. N if S
i

1 2 r i 3

aux
i

aux for i 1 .. N if S
i

2
i

aux r i aux return aux

Zt :=

for i 1 .. N aux Z
i i

v' :=

for i 1 .. N aux v i i aux

aux

58

H :=

for i 1 .. N aux Ko Vi
i i

V :=

for i 1 .. N aux
i

v'i + v'i 1

aux

aux
:=

for i 1 .. N aux Hi tan i i aux

Qult_f :=

i=1

( )

Qult :=

for i 1 .. N aux i P h i i aux

Qult

En donde estos son los algoritmos programados para realizar el clculo correspondiente a esta etapa y los cuales fueron explicados anteriormente. Se debe sealar que al hacer doble clic en la pestaa:
Valores finales para la determinacin de la capacidad de soporte por friccin

Se pueden ver los valores generados por los algoritmos a modo de vectores columna, y en donde finalmente se entrega el valor final de la capacidad de soporte por friccin lateral:

Qult_f En la ltima etapa del programa se realiza el clculo del peso propio del pilote y de la resistencia ltima total del pilote Qult_t, la que se compara con la carga axial solicitante multiplicada por el factor de seguridad ingresado, en donde finalmente se debe cumplir el siguiente criterio: Qult_f > Qsol_may Ok Qult_f < Qsol_may Se debe redimensionar el pilote

La gracia del programa es que los datos se pueden ir iterando rpidamente al ingresarlos, debido a que el programa entrega los valores y criterios de diseo en forma inmediata.

59

Los programas finales creados se pueden ver grficamente en el ANEXO D, o se pueden ver directamente al abrir el archivo correspondiente a los programas que se encuentran anexos en el CD de este proyecto. Cabe sealar, que para abrir los programas, es necesario tener en el ordenador instalado MATHCAD13, o alguna versin ms nueva que esta.

60

CAPITULO V
ANALISIS DE RESULTADOS

61

5.1.- COMPROBACIN DE LOS RESULTADOS Para la validacin de los resultados obtenidos, los programas creados se probaron con problemas dados, cuyos resultados pueden ser comparados con los resultados generados, aplicando dichos problemas a LATERAL PILES y A STATIC PILES.

5.1.1.- LATERAL PILES Para la verificacin de este programa, se utilizaron dos problemas encontrados en el libro The Foundation Engineering del autor Manjriker Gunaratne (2006), en el captulo 8, pg. 345, 346, 347 y 348, los cuales se explican a continuacin, seguido de la revisin correspondiente del programa.

5.1.1.1.- Problema 1: Determinar la rigidez elstica del suelo, de un pilote de 300x300 mm. de seccin, de largo 4.75 m. Asumir que la correccin proveniente del test de penetracin estndar (STP), evaluando la sobrecarga en una capa de arena, es 15 y 16,5 kN/m3, respectivamente.

Desarrollo: Para N=15, se determina un ngulo de friccin interna = 34 (La obtencin de este valor, se encuentra explicado en el ANEXO C de este proyecto), y segn la TABLA 3.1 (TABLA 1 de este proyecto), los factores de capacidad de carga son: Nc= 42, Nq= 29 y N= 29. Ocupando los factores definidos por Bowles (1996) y aplicando las ecuaciones 8.28 del texto, se tiene lo siguiente:

Ah = Fw1 * C m * C (c * Nc + 0.5B * N ) = (1.0)(1.5)(40)(0.5) * (16.5)(0.3)(29) = 4306,5 Bh = Fw2 * C m * C * * Nq = 1.0(1.5)(40)(16.5)(29) = 28710

Finalmente, aplicando: k h = Ah + Bh * Z n , se tiene:

k h = 4307 + 28710 Z 0.5

Esta es la rigidez del suelo encontrada para el pilote dado. Su distribucin se muestra en la siguiente figura:

62

1m

1m

12167,5 kN/m 0.75 m

Figura 11: Distribucin de esfuerzos y rigideces de un pilote Cargado lateralmente FUENTE: Gunaratne (2006)

En donde la rigidez aplicada al pilote se obtiene de las siguientes ecuaciones.

K j = LBk h K j = 0.5 LBk h

para nudos profundos

para el nudo mas superficial (nudo 2)

63

5.1.1.2.- Problema 1 aplicado a LATERAL PILES:

Anibal Meza P Ing. Civ il en Obra s Civ iles

LATERAL PILES
Este programa calcul a la defo rmacin mxima y la di stri buci n de reacciones en un pi lote, producida por una fuerza soli citante l ateral apl icada en su extremo supe ri or. Este programa uti l iza el m todo de elementos fini tos para encontrar u na matri z de ri gid ez KT gl obal del si stema gene rada por l a di visi n de l pil ote en 10 n odos.

Tomando en consideracin la entrada de datos correspondiente a las propiedades del suelo en el programa, ocupando los datos entregados en el Problema 1 y seleccionando la distribucin de rigidez kh como variable en la ventana inferior, se tiene:
Propiedades del Pilote
largo total del elemento longitud externa superior modulo de elasticidad Carga lateral

Lt

KPa

KN

Diametro del pilote

longitud de ltimo elemento

lnercia del elemento

Desplazamiento mximo del pilote

Propi edades del Suelo

Propiedades indice Segn tabla 1 Constantes segn J.E. Bowles Cohesion

N Nq

Nc

Fw1 Fw2 n kN/m 3

Cm C
KN/m3

c
Densidad del suelo

Para modulo de rigidez horizontal del suelo constante ingresar su valor

KN/m3

kh

Al hacer doble clic en la flecha que seala Valores Z, kh y Kt, se despliegan en forma de vectores columna los resultados. Cabe indicar que LATERAL PILES divide automticamente el pilote dado en 10 subelementos, por tanto se tiene:
Valores Z, kh y Kt:
Profundidad de cada nodo en m Mdulo de re accin del suelo en cada nodo Mdulo de re accin del Pilote en c ada nodo

0 0.375 0.75 1.125 1.5 z= 1.875 2.25 2.625 3 3.75

4306.5 21887.713 29170.089 34758.054 kh = 39468.925 43619.287 47371.5 50822.016 54033.679

kN/m 3

645.975 2462.368 3281.635 3910.281 kN/m Kt = 4440.254 4907.17 5329.294 5717.477 12157.578

64

En donde se puede ver la distribucin grficamente al final del programa, en donde se visualizan algunos valores de kh:

Figura 12: Distribucin de esfuerzos y rigideces de un pilote Cargado lateralmente FUENTE: Elaboracin propia

Segn los valores obtenidos, se puede apreciar la semejanza del grafico con la distribucin de la rigidez del suelo, adems de los valores de distribucin de rigideces kh y Kt.

65

5.1.1.3.- Problema 2: Un Pilote confeccionado de un perfil de acero de 300 mm de ancho, perteneciente a un grupo de pilotes, es solicitado en su extremo superior, por una carga lateral de 8 kN, generando un momento recuperador en su extremo superior restringiendo la rotacin. El coeficiente de rigidez horizontal del suelo es constante, no varia segn la profundidad, y su valor es de 1000 kN/m3. Determine la rigidez del elemento, su deflexin lateral (desplazamiento) y su esfuerzo recuperador, para 6 nodos. La inercia de la seccin es de 2.2x10-6 m4 y el modulo elstico del acero es de 2.0x106 kPa.

Desarrollo: Aplicando la metodologa de clculo para encontrar kh y Kt, considerando que kh es constante, se obtiene:

Figura 13: Distribucin de esfuerzos y rigideces de un pilote Cargado lateralmente FUENTE: Gunaratne (2006)

En donde se pueden apreciar las rigideces del elemento en cada nodo. Realizando el anlisis de esfuerzos y reacciones, se obtiene el siguiente vector de esfuerzos:
P = [M 1 8 0 0 0 0 0 0 0 0 M 11 P12 ]

66

Y aplicando las restricciones de borde siguientes, en donde el giro en la cabeza superior del pilote es 0, junto con el desplazamiento y giro en su base, se tiene: X1=1 = 0 X12= 6 = 0 X11 = 6 = 0 ,

X = [0 1 2

3 4

0 0]

Las matrices de rigidez correspondiente a los cuatro elemento son de la forma:

En donde se agrega la rigidez de cada nodo en la posicin correspondiente, generando la matriz de rigidez global del elemento:

Para resolver este sistema, aplicando P= K*X, y considerando las incgnitas que se tienen, se puede particionar la matriz K de la siguiente forma:

67

En donde se tiene que K11 es del orden de 3x3, K12, del orden de 3x9, K21 del orden de 9x3 y K22, del orden de 9x9, redefiniendo a su vez los vectores correspondientes a P y X, de la siguiente forma:
P1 = [M 1 P12 ]

M 11

P2 = [8 0 0 0 0 0 0 0 0] X 1 = [0 0 0] = 0

X 2 = [ 2

2 3

3 4

5 ]

En donde se forman subsistemas lineales, cuya solucin se encuentra de la siguiente forma:

[P1 ] = [K ]11 [X ]11 + [K ]12 [X ]2 [P]2 = [K ]21 [X ]1 + [K ]22 [X ]2

En donde el vector X2 queda de la siguiente forma:

[X ]2 = [K ]21 [X ]1 + [K ]22 1 [P]2 = [K ]22 1 [P]2

Y el vector P1 se redefine como:

[P]1 = [K ]12 [K ]22 1 [P]2

En el texto, sealan la utilizacin de un programa computacional para resolver estos sistemas, los cuales arrojan los siguientes resultados:

Del vector X2, se obtiene el desplazamiento lateral superior que sufre el pilote, cuyo valor es = 0.304 m.

68

De la misma manera, del vector P1, se obtiene el momento superior generado por la accin del suelo en el pilote, a travs de la distribucin de resortes, cuyas rigideces se mencionaron anteriormente por el vector Kt. El valor del momento requerido para refrenar el giro es M1 = 16 kNm. Estos son los resultados del ejemplo estudiado en el texo. A continuacin, se realiza la ejecucin de LATERAL PILES para comparar sus resultados.

69

5.1.1.4.- Problema 2 aplicado a LATERAL PILES Introduciendo los datos del Problema 2 a LATERAL PILES, teniendo en cuenta que la longitud del pilote que se encuentra sobre el terreno es de 1 metro, y su ltima longitud es de 0,75 m, adems de definir el mdulo de rigidez del suelo como constante e introduciendo su valor, se tiene:

Anibal Meza P Ing. Civil en Obras Civiles

LATERAL PILES
Este programa calcula la deformacin mxima y la distribucin de reacciones en un pilote, producida por una fuerza solicitante lateral aplicada en su extremo superior. Este programa utiliza el mtodo de elementos finitos para encontrar una matriz de rigidez KT global del sistema generada por la divisin del pilote en 10 nodos. Propiedades del Pilote
largo total del elemento longitud externa superior modulo de elasticidad Carga lateral

Lt

KPa

KN

Diametro del pilote

longitud de ltimo elemento

lnercia del elemento

Desplazamiento mximo del pilote

Propiedades del Suelo

Propiedades ndice Segn tabla 1 Constantes segn J.E. Bowles Cohesion

N Nq

Nc

Fw1 Fw2 n

Cm C
KN/m3

c
Densidad del suelo

Para modulo de rigidez horizontal del suelo constante ingresar su valor

KN/m3

kh

kN/m3

70

En donde las constantes segn J. E. Bowles son ignoradas por el programa debido a que la rigidez del suelo es constante y no se necesita obtener, y se agrega un desplazamiento mximo opcional, que no influye en los clculos, solo en el criterio de desplazamiento. El programa al estar generado para 10 nodos, subdivide automticamente el pilote, entregndonos los siguientes resultados:

Profundidad de cada nodo en m

Mdulo de re accin del suelo en cada nodo

Mdulo de re accin del Pilote en c ada nodo

0 0.375 0.75 1.125 1.5 z= 1.875 2.25 2.625 3 3.75

4306.5 21887.713 29170.089 34758.054 kh = 39468.925 43619.287 47371.5 50822.016 54033.679

kN/m 3

150 112.5 112.5 112.5 Kt = 112.5 112.5 112.5 112.5 225

kN/m

En donde los valores de Kt que son generados por kh, son similares a los entregados en el Problema 2 (5.1.1.3).

La matriz de rigidez global KT generada por el programa para este pilote, considerando 10 nodos por el programa, es la siguiente:

71

 17.6 26.4 8.8 26.4 0 0 0 0 0 0 0 KT = 0 0 0 0 0 0 0 0 0 0 0

26.4 52.8 26.4 52.8 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0

0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 52.8 26.4 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 64.533 161.333 23.467 187.733 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 161.333 1204.044 187.733 1001.244 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 23.467 187.733 0 23.467 187.733 93.867 0 0 0 0 0 0 0 0 0 0 0 0 0 2114.989 187.733 1001.244 0 187.733 1001.244 0 0 0 0 0 0 0 0 0 0 0 0 0 23.467 187.733 0 23.467 187.733 93.867 0 0 0 0 0 0 0 0 0 0 0 0 0 2114.989 187.733 1001.244 0 187.733 1001.244 0 0 0 0 0 0 0 0 0 0 0 0 0 23.467 187.733 0 23.467 187.733 93.867 0 0 0 0 0 0 0 0 0 0 0 0 0 2114.989 187.733 1001.244 0 187.733 1001.244 0 0 0 0 0 0 0 0 0 0 0 0 0 23.467 187.733 0 23.467 187.733 93.867 0 0 0 0 0 0 0 0 0 0 0 0 0 2114.989 187.733 1001.244 0 187.733 1001.244 0 0 0 0 0 0 0 0 0 0 0 0 0 23.467 187.733 0 23.467 187.733 93.867 0 0 0 0 0 0 0 0 0 0 0 0 0 2114.989 187.733 1001.244 0 187.733 1001.244 0 0 0 0 0 0 0 0 0 0 0 0 23.467 187.733 0 0 23.467 187.733 93.867 0 0 0 0 0 0 0 0 0 0 0 0 0 2114.989 187.733 1001.244 0 187.733 1001.244 0 0 0 0 0 0 0 0 0 0 0 0 0 187.733 23.467 0 23.467 187.733 93.867 0 0 0 0 0 0 0 0 0 0 0 0 0 2114.989 187.733 1001.244 0 187.733 1001.244 0 0 0 0 0 0 0 0 0 0 0 0 0 140.8 11.733 46.933 23.467 187.733 70.4 0 0 0 0 0 0 0 0 0 0 0 0 0 0 187.733 1001.244 140.8 1126.4 46.933 125.156 0 0 0 0 0 0 0 0 0 0 0 0 0 0 11.733 46.933 23.467 46.933 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 46.933 125.156 46.933 125.156 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
8.8 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0

26.4

73

En la cul se pueden observar valores similares a los generados por el Problema 2 del texto. A continuacin, el programa realiza los clculos de los vectores, entregndonos los siguientes resultados:

= 0.30303

Desplazamiento generado en la cabeza superior del pilote Desplazamiento del pilote a nivel del terreno Momento recuperador generado en la cabeza superior del pilote

2 = 0.0524 m M1 = 16 kN m

De splaza mi ento = "Ok"

Si el desplazamiento no cumple, se debe redimencionar el pilote hasta que su desplazamiento mx imo generado por la fuerza solicitante F, sea menor que su desplazamiento mximo

Los cuales se pueden apreciar en la siguiente figura que se incluye al final del programa:

Figura 14: Distribucin de esfuerzos y rigideces de un pilote Cargado lateralmente (Programa) Fuente: Elaboracin propia

Cuyo desplazamiento y momento M1 son de igual valor que el ejemplo estudiado en el texto indicado.

73

5.1.2.- STATIC PILES Para la verificacin de los resultados obtenidos por STATIC PILES, se consider la estatigrafa realizada en la ciudad de Valdivia, en la ubicacin del ex Hotel Pedro de Valdivia, en donde se encuentra en marcha la construccin de un Casino de juegos y hotel para la ciudad. Los datos y la informacin, fueron obtenidos de la memoria de clculo y diseo de fundaciones del proyecto, facilitado por el profesor patrocinante. La verificacin de los resultados se realizo directamente con el Profesor patrocinante, la cul ser mencionada como Problema 1.

5.1.2.1.- Problema 1: En la estatigrafa realizada para el proyecto del casino de juegos de la ciudad de Valdivia, se realizaron diversos sondajes de muestreo, en donde se obtuvieron 5 horizontes o capas de suelo, determinndose del estudio de mecnica de suelos, el uso de pilotes de hormign armado preexcabados y adems, que los horizontes H3 y H4 son capaces de entregar la Resistencia por Punta a los pilotes. Debido a la profundidad de excavacin realizada de 14,5 m. para encontrar el sello de fundacin, se consideran el horizonte 2 como primera capa de profundidad para los pilotes. Los datos obtenidos del estudio de mecnica de suelos que se entregan a continuacin, corresponden al primer sondaje realizado, y fueron los siguientes:

Altura de Capa

friccin interna

Cohesin

Densidad

CAPA

h m 13,7 2,2 5,75 13.85

21,8 32 38 24,9

C kg/cm2 0,1 0 0 0,16

T/m3 1,663 1,853 2,085 1,673

Tipo de Suelo

1 2 3 4

Limo, arena limosa Arena fina, limosa Granular (arena) Granular (arena)

Tabla 6: Valores entregados por el estudio de mecnica de suelos realizada al ex Hotel Pedro de Valdivia para el emplazamiento del nuevo casino. Fuente: Elaboracin Propia.

74

Las densidades de las capas se encuentran saturadas, por lo cual se denominan las densidades como sumergidas.

De la memoria de clculo, se concluye que uno de los largos de los pilotes ser de 13,4 m, y se realizarn iteraciones para la determinacin del dimetro de estos, que irn entre los 60, 80, 100 y 110 cm. El anlisis de los datos por parte del Profesor gua, para determinar la capacidad de soporte total de un pilote cargado axialmente entrego la siguiente tabla de resultados. Cabe mencionar que para efectos prcticos de interpretacin de datos, este anlisis se realiza por cada un metro de profundidad:

Clculo de la Capacidad total de soporte de un pilote

Profundidad

Ko 0,5786 0,5786 0,5786 0,5786 0,5786 0,5786 0,5786 0,5786 0,5786 0,5786 0,5786 0,5786 0,5786 0,5786 0,5786

Nsigma 10 10 10 10 10 10 10 10 10 10 10 10 10 30 30

Nc 22,501 22,501 22,501 22,501 22,501 22,501 22,501 22,501 22,501 22,501 22,501 22,501 22,501 46,41 46,41

m 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15

TAU Qult,f Qult,p Qult,f Qult ult i ton ton/m2 ton ton ton 25,15 0,099 0,313 0,312 25,46 28,89 0,198 0,624 0,937 29,83 32,63 0,298 0,937 1,874 34,51 36,38 0,397 1,249 3,122 39,50 40,12 0,497 1,561 4,684 44,80 43,86 0,596 1,874 6,558 50,42 47,60 0,696 2,186 8,744 56,35 51,35 0,795 2,498 11,242 62,59 55,09 0,895 2,811 14,052 69,14 58,83 0,994 3,123 17,175 76,00 62,57 1,094 3,435 20,612 83,18 66,32 1,193 3,747 24,358 90,67 136,83 1,293 4,059 28,417 165,24 149,15 1,554 4,879 33,297 182,44 161,47 1,693 5,319 38,616 200,08

Tabla 7: Clculo de la capacidad total de soporte de un pilote realizado por el profesor gua Fuente: Elaboracin Propia.

De la Tabla 7 correspondiente al clculo de la capacidad de soporte de un pilote realizado por el profesor gua, y tomando en cuenta la profundidad del pilote, el cul se encuentra entre las profundidades de 13, 14 m, y la profundidad de 13,4, que es la profundidad del pilote, interpolando, se tiene:

75

Profundidad

Qult,f 136,83 141,76 149,15

Qult,p 28,42 30,37 33,297

Qult 165,24 172,12 182,44

13 13,4 14

Tabla 8: Interpolacin realizada para encontrar los valores exactos correspondientes a la profundidad de 13,4 m Fuente: Elaboracin Propia.

A continuacin, se realiza la obtencin de los datos obtenidos por STATIC PILES, aplicados a estos datos, para luego realizar la comparacin respectiva.

76

5.1.2.2.- Problema 1 aplicado a STATIC PILES: Se realiza el ingreso de datos correspondientes al punto 5.1.2.1 Problema 1 a STATIC PILES, en donde se tiene:

Anbal Meza Prez

Ing. Civil en Obras Civiles

STATIC PILES
Este programa Verifica el diseo de un pilote a travs de la Capacidad Resistente por punta y por friccin lateral, para un suelo de hasta 10 capas, y una fuerza axial solicitante F, tal como se indica en la figura. Se deben ingresar los datos que se piden a continuacin, seleccionando el nmero de capas para luego agregar la estatigrafa correspondiente a las capas del suelo Datos del Pilote
Material del pilote Largo total del pilote Solicitacin axial del pilote Si el pilote es de un perfil de acero se debe ingresar el rea y el permetro de la seccin

L
Densidad del material

ton

diametro del pilote

Factor de Seguridad

As p

m2 m

ton/m3

Seccin del pilote Capa de suelo hasta donde llega el pilote

Datos del suelo

N de capas del suelo Ver tabla 1 (valores recomendados) Valores correspondientes al angulo de fricin interna de la ltima capa del suelo

Nq Nc
Nc debe ser Segn tabla 1 o calculado en trminos de Nq

El valor de Nc, se puede obtener de la siguiente frmula, para ser ingresado correspondientemente

Nc' = 14.406

Valor obtenido de Nq, el cul debe ser ingresado en la parte superior de este programa

Nc =

Nq 1 tan( n )

Estatigrafa del suelo:

Altura de Capas del suelo en m.

Angulo de friccin interna.

Cohesin del suelo en kg/cm2

Densidad del suelo en T/m3

Tipo de suelo en cada capa

h h h h

1 2 3 4

1 2 3 4

c c c c

1 2 3 4

1 2 3 4

77

Para el clculo de la Capacidad de resistencia ltima por punta:

0 0.579 Ko = 0.47 0.384 0.579

q = 26.349

T/m2 Ton

Esfuerzo vertical efectivo Capacidad ltima de resistencia por punta en el pilote

Qult_p = 133.847

Para el clculo de la Capacidad de resistencia ltima por friccin lateral:

0 0 13.7 9.083 Zt = 15.7 v' = 10.789 21.45 17.028 35.3 26.349

0 4.542 V = 9.936 13.908 21.688

0 2.628 H = 4.672 5.347 12.558

0 1.051 = 1.824 2.531 3.743

0 45.234 Qult = 11.462 45.716 162.865

Qult_f = 56.696

Ton

Capacidad ltima de resistencia por friccin lateral en el pilote

Clculo del peso propio y de la Capacidad de resistencia ltima total:

Pp = 25.258 Qult_t = 165.284

Ton Ton

Peso propio total del pilote Capacidad ltima de resistencia admisible total en el pilote

Considerando que el pilote se encuentra solamente en la primera capa de suelo, se tiene:

0 0.579 Ko = 0.47 0.384 0.579

q = 26.349

T/m2 Ton

Esfuerzo vertical efectivo Capacidad ltima de resistencia por punta en el pilote

Qult_p = 149.951

Qult_f = 45.234

Ton

Capacidad ltima de resistencia por friccin lateral en el pilote Peso propio total del pilote Capacidad ltima de resistencia admisible total en el pilote

Pp = 25.258 Qult_t = 169.927

Ton Ton

78

A continuacin se encuentran los valores finales para el Problema 1 y para el Problema 1 aplicado a STATIC PILES:

Mtodo de anlisis

Capa de ubicacin del pilote

Prof. Gua STATIC PILES STATIC PILES

13.4 1 2

Qult,f ton 30,37 45,234 56,69

Qult,p Qult ton ton 141,76 172,12 149,95 169,93 133,85 165,28

Tabla 9: Tabla comparativa de los valores finales de la capacidad de soporte total de un pilote. Fuente: Elaboracin Propia.

Tal como se aprecia en la Tabla 9, existe una aproximacin de los valores finales de la capacidad de soporte ltima total del pilote, en donde los valores encontrados por STATIC PILES, consideran el peso propio del pilote, a comparacin del mtodo de anlisis utilizado por el Profesor gua. Se debe volver a mencionar, que los valores de la Tabla 9 para la capa de ubicacin de 13,4, fueron obtenidos interpolando los resultados encontrados entre los 13 y 14 m. por parte del mtodo del profesor gua.

79

CAPITULO VI
CONCLUSIONES

80

En primer lugar, se sealan las conclusiones correspondientes a este trabajo de programacin, para luego sealar algunas caractersticas de uso y recomendaciones sobre los programas.

6.1.- Conclusiones:

La metodologa empleada, arrojo resultados con un margen de error muy pequeo, comparando los valores entregados entre LATERAL PILES y los ejemplos estudiados para este caso, y resultados aproximados entre STATIC PILES y el problema revisado.

Para el caso de LATERAL PILES, el margen de error pequeo se debe principalmente a que este programa utiliza 11 nodos para su anlisis, lo que conlleva a un resultado mas exacto comparado con el resultado de los ejemplos, los cuales utilizan 6 nodos. En el caso de STATIC PILES, se utilizo la metodologa propuesta por Vesic, la que es de un carcter ms conservador, comparada con la de Meyerhoff o Terzagui, ya que considerando el efecto SueloPilote por Punzonamiento, realiza el clculo de la capacidad de resistencia en trminos del Esfuerzo efectivo, con base en la teora de Expansin de Cavidades.

La utilizacin del software MATHCAD13 es simple y directa, en cuanto a la programacin de cualquier tipo de elemento estructural, mediante rutinas, arreglos matriciales o mtodos de clculo como el de elementos finitos, los cuales se ejecutan inmediatamente al ingresar los valores correspondientes a un problema de diseo dado, lo que permite realizar iteraciones en forma continua, hasta que se cumplan los criterios de diseo o hasta que se obtengan los valores requeridos por el anlisis.

81

6.2.- Caractersticas de uso y Recomendaciones para los programas:

6.2.1.- LATERAL PILES

Como se sealo en el captulo IV, este programa calcula la deformacin mxima en la cabeza y en la zona de contacto con el suelo, y la distribucin de reacciones en un pilote, producida por una fuerza solicitante lateral de corte que puede ser tomada como una carga ssmica, utilizando el mtodo de modelacin por elementos finitos.

Se debe tener claro la dimensin del pilote que estar en la superficie y la ltima dimensin que corresponde a la porcin del pilote, menor a 0,5 m, la cual puede ser dada aproximadamente. Esta capa se debe ingresar, ya que en su parte superior, el mtodo de diseo determina la ltima rigidez del suelo aplicada al pilote para la determinacin de la matriz de rigidez global del sistema.

Si se conoce el valor constante del mdulo de rigidez del suelo, se debe ingresar en la casilla correspondiente, indicando en la pestaa anexa, su condicin de constante a lo largo de su profundidad. Si se desconoce el valor de este mdulo, se debe ingresar como variable y determinar las constantes indicadas para la determinacin de este.

Existen reas de trabajo en donde se encuentran haciendo doble clic, tablas con valores de las constantes a ocupar, recomendadas por algunos autores, as como tambin reas de trabajo en donde se pueden ver las matrices generadas por los procedimientos de clculo. Existen reas de trabajo bloqueadas con una clave, en donde se encuentran los criterios y frmulas de programacin.

Se deben ingresar las unidades de los valores tal cul se mencionan en el programa, para as obtener resultados con valores conocidos y fciles de interpretar.

82

6.2.2.- STATIC PILES

Como se sealo en el capitulo IV, este programa verifica el diseo de un pilote a travs de la capacidad resistente por punta y por friccin lateral, para una estatigrafa de hasta 10 capas.

Es importante seleccionar el nmero de capas con que se est trabajando, lo cul indica al programa hasta que datos debe interpretar, as como tambin indicar en que capa de la estatigrafa se encuentra el pilote.

Si el material del pilote es acero, se debe seleccionar en la casilla correspondiente y adems agregar el valor de su seccin y de su permetro en las casillas As y p.

El programa realiza el clculo de Nc, dependiendo del valor de Nq, el cul se puede ocupar e ingresar en la casilla correspondiente.

Se puede agregar el esfuerzo mximo axial solicitante, el cul es corroborado por el programa indicando Ok si el valor de este mayorado por el factor de seguridad, es menor que Qult total. Si es incorrecto, se puede variar el dimetro y largo del pilote hasta que se cumpla el criterio de esfuerzos admisibles.

83

BIBLIOGRAFA

BOLOGNESI, A. J. 2002. Ensayos de carga estticos de compresin en pilotes de Hormign. (Disponible en: http://www.saig.org.ar/Bolognesi.pdf. Consultado el 04 de Enero de 2007).

BOWLES, J. E. 1996. Foundation Analysis and Design. 5ta ed. McGraw-Hill, New York.

BOWLES, J. E. 1974. Analytical and Computer Methods in Foundation Engineering; McGraw-Hill, New York.

DAS, B. M. 2001. Fundamentos de Ingeniera Geotecnia. Internacional Thomson Editores, Mxico. Pp. 495-560.

FORCE 2.0. Descarga de Software de Programacin. (Disponible en: http://superdownloads.uol.com.br/download/104/force-fortran-compiler-andeditor/. Consultado el 08 de Enero de 2007).

FORTRAN. Tutorial de FORTRAN. (Disponible en: http://www.fismat.umich.mx/mn1/tutor_fort/index.html. Consultado el 08 de Enero de 2007).

FRIAS, E. V. 2004. El Mtodo de los Elementos Finitos. (Disponible en: http://www.fi.uba.ar/materias/6510/Clase%2003Notas%20sobre%20Elementos%20Finitos.pdf. Consultado el 12 de Marzo de 2007).

GUNARATNE, M. 2006. The Foundation Engineering; Taylor and Francis Group, LLC, Florida.

85

OBANDO, L. L. 2004. Utilizacin de la Interpolacin en el mtodo de elementos finitos. (Disponible en: http://www.fc.uni.edu.pe/publicaciones/rev0901/6_Mantilla_final.pdf. Consultado el 12 de Marzo de 2007).

PILOTAJE.

2003,

Julio.

Pilotaje,

Cimentaciones

Profundas.

(Disponible

en:

http://www.construaprende.com/Trabajos/T3/Pilotajecim.html. Consultado el 04 de Enero de 2007).

TERZAGUI, K.; R. B. PECK. 1973. Mecnica de Suelos en la Ingeniera Practica, 2 ed.; John Wiley and Sons, Inc., New York, NY. Pp. 521-553.

VELASCO, E. 2001. Un breve curso de FORTRAN. (Disponible en: http://www.uam.es/departamentos/ciencias/fisicateoricamateria/especifica/hojas/ kike/FORTRAN/FORTRAN.html. Consultado el 04 de Enero de 2007).

86

ANEXO A
PROGRAMAS ORIGINALES EN FORTRAN IV

87

A continuacin, se entrega el detalle de la programacin original para LATERAL PILES y para STATIC PILES, por J. E. Bowles, en su libro Analytical and Computer Methods in Foundation Engineering. La recuperacin de estos programas se realiz en el software FORCE 2.0.

PROGRAMA ORIGINAL LATERAL PILES

C C C C C C C C C C JE BOWLES MATRIX DISPLACEMENT ANALYSIS OF A LATERALLY LOADED PILE *PILE MAY BE FULLY OR PATIALLY EMBECDDED--JTSS = NODE SOIL BEGINS UNITS = KIPS (KN); KSF (KN/SQ M); KCF (KN/CU M); FT OR M EXCEPT XM XL = PILE LENGTH; BX = WIDTH SQ PILE; DX = DIAM IF ROUND; BX OR DX NI > 0 TO READ MOM OF INERTIA; LIST = A,ASAT IF > 0; JJS = NO NODE S-MATRIX TO CORRECT; JTSS = NODE SOIL STARTS = 1 IF FULLY EMBEDD EXPO = EXPONENT FOR SOIL MOD. VARIATION WITH DEPTH NCL = NC OF LOADING CONDITIONS--ONLY FOR LARGES VALUES OF XMAX IF ANY NODE SPRING = 0 FOR XMAX--ADJACENT NODE SPRING REDUCED 25% (-) P(I) FORCE IS TO RIGH DIMENSION X(47),P(47),F(70),F1(70),SOIR(70),SMOD(70),G(67),PM(67) DIMENSION A(46,67),S(67,2),C(67,46),E(46,46) DIMENSION LSUM(46),V(46),EE(46),SM(70),PC(46),TITLE(20) EQUIVALENCE (A(1,1),E(1,1)) DOUBLE PRECISION UT5,UT6,NCL,NLC 6000 READ (1,1000,END=150) TITLE,UT1,UT2,UT3,UT4,UT5,UT6,FU1,FU2 1000 FORMAT (20A4/4(A4,6X),A8,2X,A8,2X,2F10.2) READ (1,105)KL,JJS,NCL,NI,NP,JTSS,LIST 105 FORMAT(7I5) WRITE (3,1001)TITLE 1001 FORMAT('1', T5,20A4) READ(1,4)XL,BX,DX,ELAS,XMAX 4 FORMAT(6F10.4) 151 WRITE(3,93)XL,UT1,BX,UT1,DX,UT1,ELAS,UT5,JJS,JTSS,NLC 93 FORMAT(//,T5,'******LATERALLY LOADED PILE BY INFINITE ELEMENT METH 1OD',//,T5, 'PILE LENGTH =',F6.2,1X,A2/T5,'PILE WIDTH (IF SQ) =',F6 2.4,1X,A2,5X,'PILE DIAM (IF ROUND) =',F6.4,1X,A2/T5,'PILE MOD OF E 3LAS =',F10.0,1X,A7/T5,'NO OF NODES REQUIRING CORRECT =',12/T5,'NOD 4E SOIL STARTS =',I3,5X,'NO OF LOAD CONDITIONS =',I3//) READ(1,4)AS,BS,EXPO WRITE(3,214)AS,BS,EXPO,XMAX,UT2 214 FORMAT(//,'SORGRADE MODULUS =',1X,F12.2,' + ',F9.2,'*Z**',F5.3,//, 1T10,'MAX LINEAR SOIL DEFORM, XMAX =',F5.2,1X,A2) M = 3*KL + 1 N = 2*KL + 2 KD = M-N KP = KL + 1 KK = 2*KL KM = N-1 K1 = KL + 2 JS = 2 LSUM(JS-1) = 0 INDEX2 = 0 IF(BX.GT.0..AND.NI.LE.0)XI = BX**4/12. IF(DX.GT.0..AND.NI.LE.0)XI = .0491*DX**4 IF(NI.GT.0)READ(1,4)XI IF(DX.GT.0.)BX = DX 206 AKL = KL H = XL/AKL WRITE(3,98)XI,UT1,H,UT1 98 FORMAT(T10,'THE MOMENT OF INERTIA OF THE PILE =',F10.7,1X,A2,'**4' 1//T10,'PILE SEGMENT LENGTHS =',F8.3,1X,A2//)

88

C*** BUILD A-MATRIX ** NON-LINEAR LCCP BEGINS SINCE A e ASAT COMMON DO 308 I = 1,M SMOD(I) = 0. DO 308 J = 1,2 308 S(I,J) = 0. 55 DO 508 I=1,N DO 508 J =1,M 508 A(I,J)=0. A(KP,KK) = 1. A(1,1)=1.0 NN=2 DO 501 J = 2,KL DO 501 I=2,3 A(J,NN)=1.0 501 NN=NN+1 K = KM DO 503 J = K1,N A(J,K)=-1.0 503 K=K+1 NN=1 DO 504 J = K1,KM DO 504 I=2,3 A(J,NN)=1./H 504 NN=NN+1 K3 = KL+3 NN=1 DO 506 J = K3,N DO 506 I=2,3 A(J,NN)=-1./H 506 NN=NN+1 IF(JS.GT.2)GO TO 300 IF(LIST.LE.0)GO TO 9544 WRITE(3,50) 50 FORMAT(////,T30,'THE STATICS MATRIX IS',//) DO 9151 I = 1,N 9151 WRITE(3,38)(I,A(I,J),J=1,M) 38 FORMAT(1X,'ROW',I2,1X,10F9.4,/,(7X,10F9.4)) C***BUILD S-MATRIX IN 2-COLUMNS 9544 DO 107 I = 1,KK DO 107 J = 1,2 S(I,1) = 4.*ELAS*XI/H S(I,2) = 2.*ELAS*XI/H IF(I/2*2.EQ.I)GO TO 103 GO TO 107 103 SAVE = S(I,1) S(I,1) = S(I,2) S(I,2) = SAVE 107 CONTINUE C*** BUIL SOIL K--1ST K IS REDUCED 50% ARBITRARILY WRITE(3,322) 322 FORMAT(//,T5,'THE SOIL MODULUS AT NODE') DO 205 I = JTSS,KP AM = I-JTSS SMOD(I)= AS +BS*(AM*H)**EXPO 205 WRITE(3,323)I, SMOD(I) 323 FORMAT(T8, I2, 3X, F12.5) S(KK+JTSS,1)=(7.*SMOD(JTSS)+6.*SMOD(JTSS+1)-SMOD(JTSS+2))*H*BX/48. S(KK+JTSS,2) = 0. S(M,1) = (7.*SMOD(KP)+6.*SMOD(KP-1)-SMOD(KP-2))*H*BX/24. S(M,2) = 0. KJTS = KK+JTSS+1 MM = M-1 MK = N-3 DO 209 I = KJTS,MM S(I,2) = 0. 209 S(I,1) = H*BX*(SMOD(1-KM)+10.*SMOD(1-KK)+SMOD(1-MK))/12. DO 115 I = KM,M

89

115 SM(I) = S(I,1) C MODIFICATION OF S-MATRIX--READ AS MANY CARDS AS JJS IF(JJS.LE.0)GO TO 99 READ (1,1008) (I,J,S(I,J),II=1,JJS) 1008 FORMAT(2I5,F10.2) C*** FORM P-MATRIX ..EXTERNAL LOADS--NCN-LINEAR SOIL EFFECTS 99 LSUM(JS) = 0 INDEX2 = INDEX2 + 1 C**** EXTERNAL LOADS INTO P-MATRIX READ NON ZERO VALUES TO NP SUML = 0. DO 241 I = 1,N G(I) = 0. 241 P(I) = 0. DO 211 JJ = 1,NP READ (1,212)I,P(I) 212 FORMAT(I5,F10.4) 211 IF(I.GT.KL+1)SUML = SUML + P(I) IF(INDEX2.GT.1.AND.NLC.GT.1)WRITE(3,8490)INDEX2 8490 FORMAT(//,T5,'COMPUTATION FOR LOAD CONDITION',I2,//) C** NON-LINEAR SOIL ENTRIES--G(I) 300 DO 213 I =1,N 213 PM(I) = P(I) + G(I) IF(JS-3)215,216,216 215 WRITE(3,223) 223 FORMAT(//T5,'THE P-MATRIX IS A FOLLOWS', T35, 'THE NON-LINEAR SO 1IL FORCES, G(I), ARE',//) GO TO 225 216 WRITE(3,48) 48 FORMAT(//,T5,'THE MODIFIED P-MATRIX IS', T35, 'THE NON-LINEAR SO 1IL FORCES, G(I), ARE',//) 225 WRITE (3,47)(I,PM(I),G(I),I=1,N) 47 FORMAT(T8,I2,F15.4, T37, F15.4) WRITE(3,226) 226 FORMAT(//,T5, 'THE S-MATRIX IN TWO COLUMNS') WRITE(3,227) (I,S(I,1),S(I,2), I=1,M) 227 FORMAT(T5,I3,2X,F12.2,2X,F12.2) IF(INDEX2.GT.1.AND.NLC.GT.1)GO TO 232 C *** BUILD THE SAT-MATRIX DO 12 I = 1,KK DO 12 J=1,N KA = 1 IF(I/2*2.EQ.1)KA = KA-1 12 C(I,J) = S(I,1)*A(J,KA)+S(I,2)*A(J,KA+1) DO 530 I = KM,M DO 530 J = 1,N 530 C(I,J) = S(I,1)*A(J,I) C *** BUILD THE ASAT-MATRIX DO 10 I = 1,N DO 15 J=1,N EE(J) = 0. DO 15 K=1,M EE(J) = EE(J)+A(I,K)*C(K,J) 15 CONTINUE DO 10 L = 1,N 10 E(I,L) = EE(L) IF(LIST.LE.0)GO TO 9168 WRITE(3,9164) 9164 FORMAT(//,T5,'THE ASAT MATRIX WITH 10**3 FACTORED',//) DO 9165 I = 1,N 9165 WRITE(3,9166)(I,E(I,J),J=1,N) 9166 FORMAT(T2,I2,2X,-3P14F8.3) C**** END OF ASAT FORMATION--INVERT ASAT MATRIX--USEE GAUSS-JORDAN METHOD 9168 DO 25 K=1,N DO 20 J=1,N 20 IF(J.NE.K)E(K,J)=E(K,J)/E(K,K) DO 210 I=1,N IF(I.EQ.K) GO TO 210

90

DO 21 J=1,N IF(J.EQ.K) GO TO 21 E(I,J)=E(I,J)-E(K,J)*E(I,K) 21 CONTINUE 210 CONTINUE DO 22 I=1,N 22 IF(I.NE.K)E(I,K)=-E(I,K)/E(K,K) E(K,K)=1./E(K,K) 25 CONTINUE C*** END OF ASAT INVERSION--COMPUTE X-MATRIX 232 DO 18 I=1,N X(I)=0. DO 18 K=1,N PC(K) = P(K) IF(JS.GT.2)PC(K) = PM(K) 18 X(I) = X(I) + E(I,K)*PC(K) C*****COMPUTE F-MATRIX--F = SAT*X DO 35 I=1,M F(I)=0. DO 35 K=1,N 35 F(I)=F(I)+C(I,K)*X(K) IF(JS.LE.2)GO TO 39 DO 37 I = KM,M 37 F(I) = F(I) + G(I-KD) 39 DO 233 I = K1,N 233 X(I) = X(I)*FU1 WRITE(3,4)UT3,UT4,UT2,UT3,UT4 41 FORMAT(///T11,'THE P-MATRIX',18X,'THE JOIN DEFLECTIONS',17X, 1'THE BEND. MOMENTS'/T8,'(',A4,' CR ',A4,') IS',18X, 2'(RADIANS OR ',A2,') ARE',18X,'(',A4,' OR ',A4,') ARE'/) C NOTE THAT ORIGINAL P-MATRIX IS ALWAYS WRITTEN HERE WRITE(3,44)(I,P(I),I,X(I),I,F(I),I=1,KK) 44 FORMAT(T7,'LOAD DIR.',I3,F10.4,T39,'JOINT DIR.',I3,F13.6,T79,'MOME 1NT',I3,F12.3) DO 49 I=KM,N NUM=I+1 49 WRITE(3,51) I,P(I),I,X(I),I,F(I) 51 FORMAT(T7,'LOAD DIR.',I3,F10.4,T39,'JOINT DIR.',I3,F13.6,T79,'FORC 1E ',I3,F11.4) IF(NUM.EQ.(M+1)) GO TO 43 DO 42 I=NUM,M 42 WRITE(3,52) I,F(I) 52 FORMAT(T79,'FORCE ',I3,F11.4) 43 CONTINUE C*** COMPUTE SEGMENT SHEARS AND MOMENTS WRITE(3,234) 234 FORMAT('1'///,T5,'SHEAR AL EACH',T22,'BEND. MOMENT AT EACH ORDINA 1TE', T55, 'SOIL REACTION AT') WRITE(3,235)UT3,UT4,UT3 235 FORMAT(T6,'SEGMENT,',A4,T31,A4,T58,'EA.ORD.',A4//) VI = 0. K = KP J = KK DO 236 I=1,KL K=K+1 J=J+1 V(I) = VI - P(K) - F(J) 236 VI = V(I) SUM = 0. DO 237 I = 1,KP SOIR(I) = -F(I+KK) 237 SUM = SUM + SOIR(I) L = -1 DO 238 I = 1,KL L=L+2 238 WRITE(3,239)I,V(I), I, F(L), SOIR(I) 239 FORMAT(T5,I2,F10.4, T24, I2, F15.4, T55,F15.4)

91

FX=-F(KK) WRITE(3,242)KP,FX,SOIR(KP) 242 FORMAT(T24,I2,F15.4,T55,F15.4) WRITE(3,243)SUM,SUML 243 FORMAT(T34,'SUM OF SOIL REACTIONS =',T55,F15.4,T72,'(',F10.3,')') C*** TEST NON-LINEAR--ZERO S-MATRIX AND F = K*XMAX = NEG. P-MATRIX ENTR IF(NLC.GT.1)GO TO 149 LSUN = 0 K2 = KP+JTSS DO 245 I = K2,N IF(ABS(X(I)).LT.XMAX)GO TO 245 244 G(I) = SM(I+KD)*XMAX/FU1 S(I+KD,1) = 0. S(I+KL,1) = 0.75*S(I+KL,1) LSUN = LSUN + 1 245 CONTINUE LSUM(JS) = LSUN IF(LSUM(JS)-1)149,246,246 246 IF(JS.EQ.7.OR.LSUM(JS).LE.LSUM(JS-1))GO TO 149 IF(JS.LE.2)GO TO 54 IF(LSUM(JS).GE.KL/2+1)GO TO 148 54 JS = JS+1 INDEX1 = JS-2 WRITE(3,1048)INDEX1,LSUM(JS-1) 1048 FORMAT(//,T10, '****COMPUTATION CYCLE',I3,' XMAX EXCEEDED AT',I3,' 1NODES',//) GO TO 55 148 WRITE(3,1010) 1010 FORMAT(//,T5,'****PILE UNSTABLE--1/2 OR MORE NODES EXCEED XMAX') 149 IF(INDEX2.LT.NLC)GO TO 99 GO TO 6000 150 STOP END

(BOWLES, 1974)

92

PROGRAMA ORIGINAL STATIC PILE

C C C JE BOWLES STATIC PILE ANALYSIS BY FINITE ELEMENT METHOD PILE MAY BE BATTERED; FULLY CR PARTIALLY EMBEDDED AND WITH VARIABL CROSS-SECTION--USE EITHER FPS CR METRIC UNITS DIMENSION ASAT(48,48),EA(6,6),P(48),X(48,2),XIN(20),XMOD(48,2), 1SLIP(20,2),SS(20),DELX(20),SK(20),SUMDX(20),NEC(10),XC(10,10), 2YC(10,10),PF(51),ELEMX(20),XL(20),XI(20),PER(20),BMEM(20),B(20), 3AREA(20),ES(6,6),ESAT(6,6),EASAT(6,6),F(6),FRIC(20),PAXIS(20), 4MND(20),NPE(20,6),EAOL(20),SM(20),RSPR(20,2),TITLE(20) DOUBLE PRECISION UT5,UT6,UT7,UT8 7000 READ(1,1000,END=6000)TITLE,UT1,UT2,UT3,UT4,UT5,UT6,UT7,UT8 1000 FORMAT(20A4/4(A4,6X),4(A8,2X)) READ(1,101)NP,NM,NNZF,NC,IPRD,NLC,JTSOIL,NDELT,JJS,LIST,NSTRPT,IU 101 FORMAT(12I5) WRITE(3,1001)TITLE 1001 FORMAT('1',///,T5,20A4) NMP = NM +1 FU1 = 12. IF(IU.EQ.2)FU2 = 98.07 READ(1,102)E,PIL,POINTK,AREAP,ALPHA,PERPP,POINTX,XMAX 102 FORMAT(8F10.4) WRITE(3,625)E,UT8,PIL,UT1,POINTK,UT7,AREAP,UT2,POINTK,UT1,PERPP,NM 1,JTSOIL,NLC,XMAX,UT2 625 FORMAT(//,T5,'E =',F12.3,1X,A8, 5X,'PILE LENGTH =',F7.3,1 1X,A2/T5, 'POINT SPRING =',F12.3,1X,A4,5X, 'PILE POINT X-SECTION =' 2,F8.3, ' SQ ',A2/ T15,'INITIAL POINT DEFL =',F7.4,1X,A2/ T15,'ASSU 3MED POINT LOAD =',F8.4,' PERCENT',5X,'NO OF PILE SEGMENTS =',I3/T2 40,'JOINT SOIL STARTS =',I3,3X,'NO OF LOAD CONDITIONS =',I3,/,T7, ' 5MAX POINT DEFLTO REDUCE POINTK 25% =',F7.4,1X,A2) IF(NDELT.GT.0)WRITE(3,632) 632 FORMAT(/,T7,'*** P-DELTA EFFECT IS CONSIDERED IN COMPUTATIONS **') READ(1,102)AS,BS,EXPO WRITE(3,1005)AS,BS,EXPO,UT6 1005 FORMAT(//,T5,'SOIL MODULUS =',F8.1,' + ',F7.2,'Z**',F5.3,1X,A8) SAVEK = POINTK E = E*FU2 AREAP =AREAP/(FU1*FU1) DO 103 I=1,NP DO 103 J=1,NP 103 ASAT(I,J)=0.0 JCOUN = 1 7500 II=0 JJ=0 JJJ = 1 SUMH = 0. DO 4 I=1,NP X(I,2) = 0. 4 P(I) = 0. WRITE(3,546) 546 FORMAT(//,T5,'THE INITIAL P-MATRIX') DO 5 IP=1,NNZP READ(1,6)I,P(I) 6 FORMAT(I5,F10.4) 5 WRITE(3,547)I,P(I) 547 FORMAT(T3,I3,3X,F12.3) NMP1 = NM+1 DO 7510 I= 1,NMP1 PAXIS(I) = 0. 7510 FRIC(I) = 0. ALPHAR = ALPHA/57.2957731 XSIN=SIN(ALPHAR) XCOS=COS(ALPHAR) PAXIS(1) = P(1)*XSIN + P(2)*XCOS IPAXIS = PAXIS(1) IF(IPRD.LE.0)PAXIS(1) = 0. IF(IPRD.GT.0)PHOLDV = P(1)

93

IF(IPRD.GT.0)PHOLDH = P(2) IF(NLC.GT.0.AND.JCOUN.GT.1)GO TO 2000 C READ LOAD TRANSFER DATA AT 8 COORDS PER CARD X = SLIP; Y = SHEAR S C NEC(JJ) = NO OF ELEMENTS FOR CORRESPONDING LOAD TRANSFER CURVE DAT DO 1012 JJ = 1,NC READ(1,6)NEC(JJ) READ(1,102)(XC(I,JJ),I=1,NSTRPT) READ(1,102)(YC(I,JJ),I=1,NSTRPT) WRITE(3,548)NEC(JJ),(XC(I,JJ),I=1,NSTRPT),UT2 548 FORMAT(//T10,'NO OF PILE ELEMENTS FOR FOLLOWING CURVE DATA =',I3// 1T2,'XC=',8F10.5,2X,A2) 1012 WRITE(3,549)(YC(I,JJ),I=1,NSTRPT),UT5 549 FORMAT(T2,'YC =',8F10.5,2X,A7) WRITE(3,177)NC,(NEC(JJ),JJ=1,NC) 177 FORMAT(//,T5,'NO OF SHEAR STR CURVES =',I3,3X, 'NO OF ELEMENTS PER 1 CURVE IN SEQUENCE =', 5(I3,2X)) WRITE(3,104)UT1,UT2,UT1,UT1,UT2 104 FORMAT( ///,T7, 'MEMNO',2X,'NP1', 2X,'NP2',2X,'NP3',2X,'NP4',2X 2,'NP5',2X, 'NP6', 3X, 'ALPHA',4X,'L ',A2,2X,'I, ',A2,'**4',2X, 3'WIDTH, ',A2,2X,'PERIM, ',A2,2X,'A SQ',A2//) K=1 DO 520 I = 1,NM MNO(I) = I NPE(I,1) = K NPE(I,2) = K+1 NPE(I,3) = K+2 NPE(I,4) = K+3 NPE(I,5) = K+4 NPE(I,6) = K+5 K = K+3 READ(1,102)XL(I),XIN(I),BMEM(I),PER(I),AREA(I),RSPR(I,1),RSPR(I,2) B(I) = BMEM(I)/FU1 WRITE(3,108)MND(I),(MPE(I,IZ),IZ=1,6),ALPHA,XL(I),X1N(I),B(I),PER( AI),AREA(I) 108 FORMAT(5X,7I5,3X,F7.3,2X,F6.2,2X,F8.2,3(4X,F6.2)) AREA(I) = AREA(I)/(FU1*FU1) EAOL(I) = E*AREA(I)/XL(I)*2. 520 XI(I) = X1N(I)/(FU1**4) AL = 0. DO 522 I = 1,NMP SK(I)= AS + BS*AL**EXPO IF(I.LT.JTSOIL)SK(I) = 0.0 522 IF(I.GE.JTSOIL.AND.I.LE.NM)AL=AL+XL(I) C*****MODIFY SUBGRADE MODULUS IF JJS>0 212 FORMAT(I5,F10.1) IF(JJS.GT.0)READ(1,212)(JTNC,SK(JTND),M=1,JJS) 517 WRITE(3,887) 887 FORMAT(//,T5,'NODE',3X, 'SOIL MODULUS AT NODE',3X,'MEMND',5X,'SOIL 1 SPRING 1',3X, 'SOIL SPRING 2') 8000 DO 205 IP = 1,NM 8050 EIOL = E*XI(IP)/XL(IP) SINOL = XSIN/XL(IP) COSOL = XCOS/XL(IP) SSPR1 = 0.0 SSPR2 = 0.0 IF(IP.LE.JTSOIL-1)GO TO 528 IF(IP.EQ.NM)GO TO 523 SSPR1 = XL(IP)*B(IP)*(7.*SK(IP)+6.*SK(IP+1)-SK(IP+2))/24. SSPR2 = XL(IP)*B(IP)*(3.*SK(IP)+10.*SK(IP+1)-SK(IP+2))/24. IF(RSPR(IP,1).GT.0)SSPR1 = RSPR(IP,1) IF(RSPR(IP,2).GT.0)SSPR2 = RSPR(IP,2) IF(IP.EQ.JTSOIL.AND.RSPR(IP,1).LE.O)SSPR1 = .5*SSPR1 GO TO 528 523 SSPR1 = XL(IP)*B(IP)*(3.*SK(IP+1)+10.*SK(IP)-SK(IP-1))/24. SSPR2 = XL(IP)*B(IP)*(7.*SK(IP+1)+6.*SK(IP)-SK(IP-1))/24. 528 IF(II.LE.0)WRITE(3,888)IP,SK(IP),IP,SSPR1,SSPR2 888 FORMAT(T5,I4,10X,F12.3,T35,I4,5X,F12.3,5X,F12.3)

94

DO 1010 IJ = 1,6 DO 1010 JJ = 1,6 EA(IJ,JJ) = 0.0 1010 ES(IJ,JJ) = 0.0 EA(1,1)=+COSOL EA(1,2)=+COSOL EA(1,3)=+XSIN EA(1,5)=-XCOS EA(2,1)=-SINOL EA(2,2)=-SINOL EA(2,3)=+XCOS EA(2,5)=+XSIN EA(3,1)=1.0 EA(4,1)=-COSOL EA(4,2)=-COSOL EA(4,4)=-XSIN EA(4,6)=-XCOS EA(5,1)=+SINOL EA(5,2)=+SINOL EA(5,4)=-XCOS EA(5,6) = +XSIN EA(6,2) = 1.0 EA(1,1) = 4.*EIOL EA(1,2) = 2.*EIOL EA(2,1) = 2.*EIOL EA(2,2) = 4.*EIOL IF(IP.LT.NM)XLMAX = MAX1(XL(IP),XL(IP+1)) IF(IP.LT.NM)XLMIN = MIN1(XL(IP),XL(IP+1)) IF(IP.GT.1)XLMAX2 = MAX1(XL(IP-1),XL(IP)) IF(IP.GT.1)XLMIN2 = MIN1(XL(IP-1),XL(IP)) RATIO = XLMAX/XLMIN SUMSQ = (1.+RATIO)**2 FAC2 = SUMSQ/(RATIO**2) F1 = SUMSQ IF(IP.EQ.1)GO TO 756 IF(XL(IP).GT.XL(IP+1))F1 = FAC2 RATIO2 = XLMAX2/XLMIN2 SUMSQ2 = (1.+RATIO2)**2 FAC4 = SUMSQ2/(RATIO2**2) F3 = SUMSQ2 IF(XL(IP).GT.XL(IP-1))F3 = FAC4 756 IF(IP.EQ.1)F3 = 1. ES(3,3) = EAOL(IP)/F3 IF(II.GT.0)ES(3,3) = EAOL(IP) IF(IP.EQ.1)ES(3,3) = EAOL(IP) ES(4,4) = EAOL(IP)/F1 IF(II.GT.0)ES(4,4) = EAOL(IP) IF(IP.EQ.NM)ES(4,4) = EAOL(IP) ES(5,5) = SSPR1 ES(6,6) = SSPR2 DO 202 I = 1,6 DO 202 J=1,6 ESAT(I,J)=0.0 DO 202 K = 1,6 202 ESAT(I,J)=ESAT(I,J)+ES(I,K)*EA(J,K) IF(II.GT.0) GO TO 605 203 DO 204 I=1,6 DO 204 J=1,6 EASAT(I,J)=0.0 DO 204 K = 1,6 204 EASAT(I,J)=EASAT(I,J)+EA(I,K)*ESAT(K,J) DO 205 I=1,6 NS1 = NPE(IP,I) DO 205 j=1,6 NS2 = NPE(IP,J) 205 ASAT(NS1,NS2)=ASAT(NS1,NS2)+EASAT(I,J) 521 WRITE(3,889)NMP,SK(NMP)

95

889 FORMAT(T5,I4,10X,F12.3) IF(LIST.LE.0)GO TO 518 WRITE(3,304) 304 FORMAT(//,T5,'THE ASAT MATRIX',/) DO 302 I=1,NP 302 WRITE(3,303)I,(ASAT(I,J),J=1,NP) 303 FORMAT(I3,2X,12F10.2) C ASAT INVERSION 518 M = NP DO 25 K=1,M DO 20 J=1,M 20 IF(J.NE.K)ASAT(K,J)=ASAT(K,J)/ASAT(K,K) DO 210 I=1,M IF(I.EQ.K) GO TO 210 DO 21 J=1,M IF(J.EQ.K) GO TO 21 ASAT(I,J)=ASAT(I,J)-ASAT(K,J)*ASAT(I,K) 21 CONTINUE 210 CONTINUE DO 22 I=1,M 22 IF(I.NE.K)ASAT(I,K)=-ASAT(I,K)/ASAT(K,K) 25 ASAT(K,K)=1.0/ASAT(K,K) C *******END OF ASAT INVERSION--COMPUTE X-MATRIX***** 2000 P11B = P(1) P21B = P(2) P31B = P(3) DO 7 I=1,NP X(I,1)=X(I,1)+ASAT(I,J)*P(J) IF(ABS(X(I,1)).LT..0000001)X(I,1) = 0.0 7 XMOD(I,1) = X(I,1) XPDEL = X(2,2) WRITE(3,890)JJJ 890 FORMAT(//,T5,'THE X-MATRIX FOR CYCLE',I3,'IS') ANP = NP/6 INP = ANP IF(INP.LT.ANP)INP=INP+1 IZ = 1 DO 188 I = 1,INP IZP = IZ+5 IF(IZP.GT.NP)IZP= NP WRITE(3,623)(NN,X(NN,1),NN=IZ,IZP) 188 IZ = IZP+1 623 FORMAT(T3,I3,F10.6,5(I4,F10.6)) C ****** COMPUTE ELEMENT DEFLECTIONS ***** MM = 1 DO 26 I = 1,NM XPER1 = XL(I)/(XL(I)+XL(I+1)) IF(I.GT.1)GO TO 27 ELEMX(I) = X(MM,1)*XSIN + X(MM+1,1)*XCOS + X(MM+3,1)*XPER1*XSIN + AX(MM+4,1)*XPER1*XCOS GO TO 26 27 XPER2 = XL(I)/(XL(I-1)+XL(I)) IF(I.EQ.NM) GO TO 28 ELEMX(I) = X(MM,1)*XPER2*XSIN + X(MM+1,1)*XPER2*XCOS + X(MM+3,1)* 1XPER1*XSIN + X(MM+4,1)*XPER1*XCOS GO TO 26 28 ELEMX(NM) = X(NP-2,1)*XSIN + X(NP-1,1)*XCOS + X(MM,1)+XPER2*XSIN + A X(MM+1,1)*XPER2*XCOS IF(P(NP-2).GT.0.0.OR.POINTK.LE.0.)GO TO 26 POINTF = P(NP-2)/XSIN C NOTE PTDEF1 IS APPROX VALUE OF PRECED CYCLE IF P(NP-2)<0. POINTK = SAVEK 311 PTDEF1 = POINTF/POINTK IF(POINTK.LE..75*SAVEK)GO TO 26 IF(ABS(PTDEF1*FU1).GT.XMAX)POINTK = .75*SAVEK IF(POINTK.LE..75*SAVEK)GO TO 311 26 MM=MM+3

96

IF(JJJ.LT.2.AND.PTDEF1.GT.POINTX)PTDEF1 = POINTX IF(JJJ.GT.2)PTDEF1 = (PTDEF1+PTDEF2)/2 MM=NM DO 9 I=1,NM IF(I.GT.1) GO TO 10 SUMDX(MM) = ELEMX(MM)+PTDEF1 DELX(MM)=ELEMX(MM) SLIP(MM,1) = (DELX(MM) + PTDEF1)*FU1 GO TO 11 10 DELX(MM)=DELX(MM+1)+ELEMX(MM) SUMDX(MM) = SUMDX(MM+1) + ELEMX(MM) SLIP(MM,1) = (DELX(MM) + PTDEF1)*FU1 11 IF(SLIP(MM,1).LT.0)SLIP(MM,1) = -SLIP(MM,1) 9 MM=MM-1 WRITE(3,30)UT1,UT1,UT2,PTDEF1,UT1 30 FORMAT(//,T5,'ELEM NO',2X,'TOT ELEM DISP',A2,3X,'ELEM MOVE,',A2, 15X,'ELEM SLIP',A2,4X,'PTDEF1 =',F10.7,1X,A2) DO 32 I=1,NM 32 WRITE(3,31)I,SUMDX(I),ELEMX(I),SLIP(I,1) 31 FORMAT(T7,I3,5X,F12.6,5X,F12.6,3X,F12.6) IF(LIST.GT.0)WRITE(3,626) 626 FORMAT(//,T10,'PRINT COUNTERS USED TO RECOMPUTE P-MATRIX') 3060 K = 1 KN = 1 I2 = 1 KCC = NEC(KN) DO 501 KZ = 1,NP 501 IF(KZ/3*3.NE.KZ)P(KZ)=0. DO 52 I = 1,NM IF(I.GT.KCC)I2=KN+1 IF(I2.NE.KN)KCC = KCC+NEC(I2) DO 54 JJ = 1,8 IF(XC(JJ,I2).GT.SLIP(I,1))GO TO 43 IF(XC(JJ,I2).EQ.SLIP(I,1))GO TO 44 54 CONTINUE 43 JJ = JJ-1 44 C1 = ABS(YC(JJ+1,I2)-YC(JJ,I2)) DELYC = 0. C2 = XC(JJ+1,I2)-XC(JJ,I2) IF(YC(JJ+1,I2).LT.YC(JJ,I2))C1 = -C1 DELTAX = SLIP(I,1) - XC(JJ,I2) IF(ABS(C2).LE..0000001)C2 = 0 IF(C2.EQ.0.)GO TO 261 DELYC = DELTAX*C1/C2 261 SS(I) = YC(JJ,I2)+DELYC FRIC(I) = XL(I)*PER(I)*SS(I) JPLUS = I+1 IF(IPRD.GT.0)GO TO 45 IF(IPAXIS)234,235,235 234 PAXIS(I+1) = PAXIS(I) - FRIC(I) GO TO 46 235 PAXIS(I+1) = PAXIS(I) + FRIC(I) GO TO 46 45 IF(PAXIS(I))81,80,80 80 PAXIS(I+1) = PAXIS(I)-FRIC(I) IF(PAXIS(I+1).LT.0..AND.FRIC(I).NE.0.)PAXIS(I+1) = 0.0 GO TO 46 81 PAXIS(I+1) = PAXIS(I)+FRIC(I) IF(PAXIS(I+1).GT.0..AND.FRIC(I).NE.0.)PAXIS(I+1) = 0.0 46 IF(LIST.GT.0)WRITE(3,154)NFC(I2),XC(JJ,I2),XC(JJ+1,I2),JJJ,JJ,I2, 1KCC,JPLUS,PAXIS(I+1),UT3,FRIC(I),UT3,PAXIS(I),UT3 154 FORMAT(T15,'NEC(I2)=',I3,3X,'XC(JJ,I2)=',F7.5,3X,'XC(JJ+1,I2) =', 1F7.5/ T10,'JJJ=',I3,3X,'JJ=',I3,3X,'I2=',I3,3X,'KCC=',I3/ T12, 2'CURRENT VALUE PAXIS(',I3,') =',F10.3,1X,A4,5X,'ELEM. FRIC =', 3F10.3,1X,A4/ T14,'CURRENT VALUE PAXIS(I) =',F10.3,1X,A4) 52 KN = I2 C POINT EFFECTS INCLUDED IN P-MATRIX--***--NOTE POINTX OR PERPP ARE

97

READ DATA: PTDEFL IS COMPUTED POINT DEFLECTION POINTF = POINTF1*POINTK IF(PTDEF1.GT.0.)POINTF = 0.0 IF(IPRD.GT.0)GO TO 151 IF(PERPP.GT.0.)GO TO 48 47 PAXIS(NM+1) = PAXIS(NM+1)+POINTF GO TO 151 48 PAXIS(NM+1) = PAXIS(NM+1)*(1.+PERPP) POINTF = PAXIS(NM+1)-(1.-PERPP)*PAXIS(NM+1) PTDEF1 = POINTF/POINTK 151 WRITE(3,208)POINTF,UT3,PTDEF1,UT1 208 FORMAT(//T5,'THE POINT DEFL DUE TO',F12.4,1X,A4,' IS',F12.7,1X, 1A2//) L=1 DO 37 I = 1,NM IF(IPRD.LE.0.AND.IPAXIS.LT.0.)FRIC(I) = -FRIC(I) IF(PHOLDV.EQ.0..AND.PHOLDH.EQ.0..AND.IPRD.GT.0.)GO TO 41 IF(IPRD.GT.0.AND.I.EQ.1)GO TO 36 IF(IPRD.GT.0)GO TO 38 IF(IPRD.LE.0)GO TO 39 36 P(1) = PHOLDV P(2) = PHOLDH 38 P(L+3) = PAXIS(I+1)*XSIN P(L+4) = PAXIS(I+1)*XCOS GO TO 37 39 IF(I.GT.1)GO TO 35 P(1) = PAXIS(NM+1)*XSIN P(2) = PAXIS(NM+1)*XCOS P(L+3) = P(1)-FRIC(I)*XSIN P(L+4) = P(2)-FRIC(I)*XCOS GO TO 37 35 P(L+3) = P(L)-FRIC(I)*XSIN P(L+4) = P(L-1)-FRIC(I)*XCOS GO TO 37 C THIS BRANCH TAKES CARE OF PILE WITH MOMENT ONLY 41 P(L+3) = 0. P(L+4) = 0. 37 L = L+3 C ****** APPLY P-DELTA EFFECT IF(NDELT.LE.0)GO TO 158 LK=0 DO 924 KI=1,JTSOIL P(LK+3) = -(X(2,1) - X(LK+2,1))*P(1) 924 LK=LK+3 158 WRITE(3,906)UT2,UT5 906 FORMAT(///,T4,'NM',2X,'SLIP,',A2,2X,'SHEAR,',A7,6X, 1'P-MATRIX IN 3-COLS') 907 IP1 = 1 DO 74 I = 1,NM IP2 = IP1+1 IP3 = IP1+2 WRITE(3,893)I,SLIP(I,1),SS(I),IP1,P(IP1),IP2,P(IP2),IP3,P(IP3) 893 FORMAT(T3,I3,F10.6,2X,F10.5,3X,I2,F10.3,2X,I2,F10.3,2X,I2,F10.3) 74 IP1 = IP3+1 NMP2 = NP-2 WRITE(3,895)(IM,P(IM),IM=NMP2,NP) 895 FORMAT(T31,I2,F10.3,2X,I2,F10.3,2X,I2,F10.3) JJJ = JJJ+1 IF(JJJ.LE.2)GO TO 1120 C TEST FOR SLIP CONVERGENCE--------KCOUN = 0 DO 58 I = 1,NM SLIPCH = SLIP(I,1)-SLIP(I,2) 58 IF(ABS(SLIPCH).LE.0.008)KCOUN = KCOUN +1 IF(JJJ.GT.7)GO TO 77 1120 DO 896 J = 1,NM X(J,2) = X(J,1)

98

SLIP(J,2) = SLIP(J,1) 896 XMOD(J,2) = XMOD(J,1) PTDEF2 = PTDEF1 IF(JJJ.LE.2.OR.KCOUN.LT.NM)GO TO 2000 IF(ABS(X(2,1)-XPDEL).GT.0.020.AND.NDELT.GT.0)GO TO 2000 77 II = II+1 JJJ = JJJ-1 WRITE(3,170)JJJ 170 FORMAT('1',///,T5,'THE FINAL COMPUTED PILE ELEMENT FORCES AND 1 OTHER DATA AFTER',I2,' ITERATIONS FOLLOWS') WRITE(3,59) UT1 59 FORMAT(//,T5,'MEMNO',5X,'END MOMENTS',10X,'AXIAL FORCE',8X,'SOIL 1REACTIONS',3X,'ELEM.DEFL,',A2,4X,'SOIL SPRONGS') GO TO 8000 C CYCLE TO FIND PILE ELEMENT FORCES AS F = SAT*X 605 DO 851 I = 1,6 F(I) = 0. DO 851 K = 1,6 N = NPE(IP,K) IF(IP.LT.NM)XPER1 = XL(IP)/(XL(IP)+XL(IP+1)) IF(IP.GT.1)XPER2 = XL(IP)/(XL(IP-1)+XL(IP)) XMOD(N,1) = X(N,1) IF(I.LE.2)GO TO 851 IF(I.GE.5)GO TO 851 IF(K/3*3.EQ.K)GO TO 851 849 IF(IP.EQ.1.AND.K.LE.3)GO TO 851 IF(IP.EQ.1.AND.K.GT.3)XMOD(N,1)=X(N,1)*XPER1 IF(IP.GT.1.AND.K.LT.3)XMOD(N,1) = X(N,1)*XPER2 IF(IP.GT.1.AND.IP.NE.NM.AND.K.GT.3)XMOD(N,1)=X(N,1)*XPER1 851 F(I) = F(I) + ESAT(I,K)*XMOD(N,1) SUMH = SUMH + (F(5)+F(6))*XSIN WRITE(3,62)IP,(F(I),I=1,6),SUMDX(IP),SSPR1,SSPR2 62 FORMAT(T5,I5,6F10.3,F12.5,2X,2F12.3) IF(IP.EQ.NM)GO TO 609 GO TO 205 609 WRITE(3,63)SUMH,UT3 63 FORMAT(//T10,'SUM HOR. SOIL REACTIONS =',F10.2,1X,A4) DELV = SUMDX(1)*XSIN*FU1 DELH=X(2,1)*FU1 ROT = X(3,1) WRITE(3,175)DELV,UT2,DELH,UT2,ROT 175 FORMAT(//T15,'PILE TOP MOVEMENTS: VERT =',F10.5,1X,A2/ T36, 1'HORZ=',F10.5,1X,A2/ T37,'ROT =',F10.5,' RAD =') POINTF = F(4)/XSIN WRITE(3,174)P11B,UT3,P21B,UT3,P31B,UT4,POINTF,UT3 174 FORMAT(//T5,'****** THE ULTIMATE PILE RESISTANCE =',F12.4,1X,A4, 1' VERTICAL AND',F9.3,1X,A4,' HORIZ'/ T12,'THE ORIGINAL MOMENT APPL 2IED =',F12.4,1X,A4/ T15,'POINT FORCE = ',F12.3,1X,A4) IF(IPRD.LE.0)WRITE(3,241) 241 FORMAT(T15,'***PILE RESISTANCE IS COMPUTED') JCOUN = JCOUN+1 IF(JCOUN.GT.NLC)GO TO 7000 IF(NLC.GT.1)READ(1,102)POINTX WRITE(3,7701)JCOUN,POINTX,UT1 7701 FORMAT(//T10,'CURRENT PROB NO =',I3,3X,'POINTX =',F10.6,1X,A4) GO TO 7500 6000 STOP END

(BOWLES, 1974)

99

ANEXO B
MATRIZ DE RIGIDEZ GLOBAL DEL SISTEMA

100

MATRIZ DE RIGIDEZ GLOBAL KT

4 EI L 1 6 EI 2 L1 2EI L 1 6 EI L2 1 0 0 0 0 M 0 0 0 0 6 EI L12 12 EI L13 6 EI L12 12 EI L13 0 0 0 0 M 0 2 EI L1 6EI L12 4EI 4EI + L1 L2 6EI 6EI 2 + 2 L1 L2 2 EI L2 6EI 2 L2 0 0 M 0 6 EI L12 12EI L13 + 6EI
2

0 0 2 EI L2 6EI L2 2 4 EI 4EI + L2 L3 6 EI 6 EI 2 + 2 L2 L3 2 EI L3 6 EI 2 L3 M K K

0 0 6 EI L2 2 12 EI L23
3

0 0 0 0
2

0 0 0 0 6EI L32 12 EI L33

2

0 0 0 0 0 0
2

0 0 0 0 0 0 6 EI 6 EI L4 2 12 EI L43 O +

K K K K K K L K 6 EI LN 12 + K N 1 O 2 EI LN 1 6 EI LN 12 4 EI LN 1 6 EI LN 12

6EI L12

L2 12 EI 12EI + + K2 L13 L23 6EI L2 12EI L23 0 0 M 0

2

6 EI L2

L3 12 EI 12 EI + + K3 L23 L33 6 EI L3 12 EI L33 M 0

2

6 EI

2EI L3 6 EI L32 4 EI 4 EI + L3 L4 6 EI 6 EI 2 + 2 L3 L4 O 2 EI LN 2 6 EI LN 2 0 0
2

6 EI L3

L4 12 EI 12 EI + + K4 L33 L43 O 6 EI LN 2 2 12 EI LN 2 0 0
3

6 EI

2 EI L4 6 EI L4 2 O

4 EI 4 EI + LN 2 LN 1 6 EI LN 2
2

LN 2 2
3

6 EI LN 12

12 EI LN 2

12 EI LN 13 6 EI LN 12

2 EI 0 0 0 0 0 0 K K 0 0 LN 12 6 EI LN 12

12 EI LN 13

0 0 0 0 0 0 0 M 6 EI LN 12 12 EI LN 13 6 EI LN 12 12 EI LN 13 0

ANEXO C
TEST DE PENETRACION STANDARD

102

TEST DE PENETRACION STANDARD STP Segn Gunaratne (2006), el autor R. H. Perry en 1977 en su investigacin Estimating bearing capacity of sand from STP values, Journal of the Geotechnical Engineering Division, ASCE, 103(GT9): 1014 1019, determin lo siguiente: Con las expresiones de capacidad de resistencia de una zapata extendida en un suelo cohesivo, se determina:

s q = 30 N 55 25.4
Donde: N55 : s: Correccin de penetracin standard correspondiente a una eficiencia del 55% en el golpe del martillo o martinete Ajuste en milmetros (B/L) Del STP, se obtiene una correlacin estimativa y equivalente al ngulo de friccin interna del suelo, aplicando la siguiente ecuacin:

N 55 = 25 + 28 q
Donde: q: Estrs efectivo de sobrecarga en la evaluacin.

En lo que corresponde al punto 5.1.1.1 Problema 1, la determinacin del ngulo de friccin interna, se realiza de la siguiente forma: De los datos entregados, se utiliza por criterio del autor, el valor de N= 15, y utilizando s= 8,2 mm, se obtiene lo siguiente:

8.2 q = 30 15 = 145 25.4

Luego, aplicando la segunda ecuacin, se tiene:

15 2 = 25 + 28 = 34 145
En donde el ngulo de friccin interna obtenido es = 34, con cuyo valor se ingresa en la tabla 1 de este proyecto para obtener Nq, Nc y N.

103

ANEXO D
PROGRAMAS FINALES

LATERAL PILES y STATIC PILES

104

Anbal Meza P Ing. Civil en Obras Civiles

LATERAL PILES
Este programa calcula la deformacin mxima y la distribucin de reacciones en un pilote, producida por una fuerza solicitante lateral aplicada en su extremo superior. Este programa utiliza el mtodo de elementos finitos para encontrar una matriz de rigidez KT global del sistema generada por la divisin del pilote en N nodos.

Propiedades del Pilote

Largo total del elemento Longitud externa superior Modulo de elasticidad Carga lateral

Lt

E
inercia del elemento

KPa
4

P
N de nodos

KN

Dimetro del pilote

Longitud de ltimo elemento

Propiedades del Suelo

Propiedades indice Segn tabla 1 Constantes segn J.E. Bowles Cohesin

N Nq

Nc

Fw1 Fw2 n

Cm C
KN/m3

c
Densidad del suelo

Para modulo de rigidez horizontal del suelo constante

KN/m3

kh

kN/m3

TABLA 1 Factores de capacidad segn Terzagui y Vesic Determinacin del coeficiente de rigidez del suelo kh Determinacin de la rigidez Kt del elemento Ver valores de Z, kh y Kt Generacin de matriz KT KT=

105

Resultados finales Criterio de verificacin

= 0.303 2 = 0.248
M1 = 16

m m kN m

Desplazamiento generado en la cabeza superior del pilote Desplazamiento del pilote a nivel del terreno Momento recuperador generado en la cabeza superior del pilote

Desplazamiento = "Ok"

Si el desplazamiento no cumple, se debe redimencionar el pilote hasta que su desplazamiento mximo generado por la fuerza solicitante F, sea menor que su desplazamiento mximo .

Si los resultados de desplazamiento son positivos y de valores excesivos, existe una sobredimencin del pilote, ante lo cual se pueden disminuir las dimensiones de este, o bien, puede dejarse tal cual.

106

Anbal Meza Prez

Ing. Civil en Obras Civiles

STATIC PILES
Este programa Verifica el diseo de un pilote a travs de la Capacidad Resistente por punta y por friccin, para un suelo de hasta 10 capas, y una fuerza axial solicitante F, tal como se indica en la figura. Se deben ingresar los datos que se piden a continuacin, seleccionando el nmero de capas para luego agregar la estatigrafa correspondiente a las capas del suelo

1 1 c1 2 2 c2

h1 h2

5 5

c5

h5

Datos del Pilote

Material del pilote Largo total del pilote Solicitacin axial del pilote Si el pilote es de un perfil de acero se debe ingresar el rea y el permetro de la seccin

L
Densidad del material

ton

diametro del pilote

Factor de Seguridad

As p

m2 m

ton/m3

Seccin del pilote

Datos del suelo

N de capas del suelo Ver tabla 1 (valores recomendados) Valores correspondientes al ngulo de friccin interna de la ltima capa del suelo

Nc

Nq

TABLA 1- valores recomendados por Vesic para Nc y Nq

Ingreso de estatigrafa del suelo:

107

Altura de Capas del suelo en m.

Angulo de friccin interna.

Cohesin del suelo en kg/cm2

Densidad del suelo en T/m3

Tipo de suelo en cada capa

h h h h h h h h h h

1 2 3 4 5 6 7 8 9

1 2 3 4 5 6 7 8 9 10 c

c c c c c c c c c

1 2 3 4 5 6 7 8 9

1 2 3 4 5 6 7 8 9 10

10

10

Clculo y desarrollo:

Clculo de Ko y q, capacidad de soporte por punta:

q = Qult_p =

Clculo de la capacidad de soporte por friccin: Valores finales para la determinacin de la capacidad de soporte por friccin

Qult_f =
Criterio final de diseo

Pp

Peso propio total del pilote Capacidad ltima de resistencia admisible total en el pilote Carga axial solicitante Mayorada Criterio final de diseo, si este "No cumple", se debe redimensionar el pilote hasta que el criterio sea "Ok".

Qult_t = Qsol =

Criterio = ___

108

ANEXO E
DIAGRAMAS DE FLUJO DE PROGRAMACION

LATERAL PILES y STATIC PILES

109

DIAGRAMA DE FLUJO DE PROGRAMACIN DE LATERAL PILES

INICIO

INGRESO DE DATOS 1) PROPIEDADES DEL PILOTE 2) PROPIEDADES DEL SUELO

Rigidez Horizontal del Suelo Kh Variable Calculada por el Programa segn las constantes de J. E. Bowles predeterminadas

Constante Ingresar Valor

Calculo y Desarrollo del Programa

Determinacin de la Rigidez "Kt" de los Resortes del pilote a travs de Kh

Modelacin del Pilote por el Mtodo de Elementos Finitos

Generacin de la Matriz Global de Rigidez del Sistema

Ingreso de Carga Lateral Solicitante "P"

Clculo de Esfuerzos y deformaciones

Redimencionar el Pilote

NO

Criterio de Diseo Amax > A

SI ENTREGA FINAL DE RESULTADOS A: Desplazamiento Superior A2: Desplazamiento en la Superficie M: Momento requerido para refrenar la rotacin en la cabeza del pilote

FIN

110

1) Propiedades del Pilote: - Lt: Largo total - B: Dimetro - I: Inercia - E: Mdulo de elasticidad - Amax: Desplazamiento mx. permitido En la cabeza del pilote 2) Propiedades del Suelo: - Propiedades ndice: Nc, Nq, N - : Densidad del suelo - c: Cohesin del suelo Nota: A = ; y =

Superficie

Modelamiento del Pilote Empotrado Fijo Empotrado Movil

111

DIAGRAMA DE FLUJO DE PROGRAMACIN DE STATIC PILES

INICIO

INGRESO DE DATOS 1) PROPIEDADES DEL PILOTE 2) PROPIEDADES DEL SUELO

Circular

Seccin del pilote

Cuadrada

Hormign

Material del pilote

Acero

Madera

Clculo y Desarrollo

Coeficiente de Presin de Tierra Ko Esfuerzo Vertical Efectivo "q" Capacidad Ultima de Resistencia por Punta "Qult,p" Clculo de: Presin Vertical de empuje Presin lateral del suelo Friccin entre suelo y pilote

Capacidad Ultima de Resistencia por Friccin Lateral "Qult,f"

Mayoracin de Carga Axial solicitante "Qsol" Clculo del peso propio del pilote "Pp" Capacidad de Resistencia Ultima del Suelo Qult = Qult,p + Qult,f - Pp

Redimencionar el Pilote

NO

Criterio de Diseo Qult > Qsol SI

Entrega Final de Resultados Pp, Qult, Qsol, Cumplimiento de criteio

FIN

112

1) Propiedades del Pilote: - Lt: Largo total - B: Dimetro - : Densidad del material - P: Solicitacin axial. - FS: Factor de seguridad 2) Propiedades del Suelo (estatigrafa): - N de capas del suelo - h: altura de cada capa - : ngulo de friccin interna - c: cohesin del suelo por capa - : densidad saturada del suelo - Propiedades ndice: Nq, Nc

113