Tesis Javier Guzmán Nápoles Leal PDF
Tesis Javier Guzmán Nápoles Leal PDF
Tesis Javier Guzmán Nápoles Leal PDF
Facultad de Construcciones
Departamento de Ingeniería Civil
TRABAJO DE DIPLOMA
Santa Clara
2016
"Año 57 de la Revolución"
i
PENSAMIENTO
(Carl Sagan).
ii
DEDICATORIA
A mi madre Ivis, que volcó su atención a mí desde que nací, que siempre estará conmigo
en cada paso. A mi segunda madre, mi tía Yixy, por creer siempre en mí, a mi abuela
Iris, con su apoyo incondicional de abuela. A mis tíos Ale y Darién, que me aconsejaron
como si de padres se trataran. A mis dos hermanos Henry y Tony que tanta fe tienen en
mí. Y a mi novia Claudia, la razón de que mis frustraciones sean más llevaderas. Y a
todos los amigos que empujan para llegar a la meta. Gracias por acompañarme.
iii
AGRADECIMIENTOS
Quiero agradecer primero a mi tutor, considerado por mi parte un amigo, que siempre
sacó tiempo de socorro en los momentos en que me aparecían los obstáculos más difíciles
de vencer. A mi abuela, mi novia y mi suegra por preocuparse tanto de mis necesidades,
mantenerme decidido y dispuesto. A mi tía por su inagotable fuente de confianza y
seguridad, que inculcó en mí el deseo de superarme y dar lo mejor. Y a todas las
personas que disfrutan o sufren con cada uno de mis logros o caídas, en estos cinco
años.
iv
RESUMEN
Se presenta un software programado en MATLAB para el diseño automatizado geotécnico
y estructural de cimentaciones superficiales rectangulares aisladas, teniendo en cuenta la
aplicación del método de los Estados Límites para evaluar la seguridad de los mismos,
determinando el costo de la solución. En el trabajo además se incluyen ejemplos de diseño
que contribuyen a entender la forma de emplear el sistema y se comparan con los resultados
de diseño obtenidos por otras hojas de cálculo presentadas por otros autores, como forma
de validarlos. Este sistema constituye un elemento dentro de un sistema más amplio
dirigido a la optimización de conjuntos estructurales, que utiliza las potencialidades de
MATLAB para ello, de ahí su pertinencia. En la concepción del software se utilizan los
últimos avances del método de los Estados Límites aplicado al diseño geotécnico y
estructural de las cimentaciones. El resultado global será una herramienta que asista al
ingeniero o estudiante a resolver problemas en el campo del diseño óptimo de las
cimentaciones o estructuras en general.
v
ABSTRACT
Software is presented programmed in MATLAB for the design automated geotechnical and
structural of isolated rectangular superficial foundations, keeping in mind the application of
the method of the States Limits to evaluate the security of the same ones, determining the
cost of the solution. In the work design examples are also included that they contribute to
understand the form of using the system and they are compared with the design results
obtained by other calculation leaves presented by other authors, like form of validating
them. This system constitutes an element inside a wider system directed to the optimization
of structural groups that uses the potentialities of MATLAB for it, of there its relevancy. In
the conception of the software the last advances of the method of the States Limits are used
applied to the geotechnical and structural design of the foundations. The global result will
be a tool that attends the engineer or student to solve problems in the field of the optimum
design of the foundations or structures in general.
vi
Tabla Contenidos
PENSAMIENTO ..................................................................................................................... i
DEDICATORIA ..................................................................................................................... ii
RESUMEN............................................................................................................................. iv
ABSTRACT............................................................................................................................ v
INTRODUCCIÓN ..................................................................................................................9
Campo de Acción..................................................................................................................10
Hipótesis................................................................................................................................11
Tareas científicas...................................................................................................................12
Valor metodológico...............................................................................................................13
Valor práctico........................................................................................................................13
1.1.6 Metodología para Método del Factor de Seguridad Global. (MFSG). ...........19
Diseño Estructural.........................................................................................................30
MATLAB VS EXCEL..................................................................................................36
1.6. Cálculo de los costos directos a partir de los reglones variantes del PRECONS...41
Referencias............................................................................................................................46
2.5. Determinación del costo del tipo de cimentación de estudio en el país. ................81
Conclusiones: ........................................................................................................................93
Recomendaciones..................................................................................................................94
Bibliografía ...........................................................................................................................95
Anexos ..................................................................................................................................97
Diseño automatizado de cimentaciones superficiales con MATLAB
9
INTRODUCCIÓN
En nuestros tiempos la ingeniería acude a varias ayudas de cálculo para lograr viabilizar los
proyectos estructurales que se presentan. Específicamente en el campo de las estructuras se
logran proyectos cada vez más racionales y a la vez con una evaluación más objetiva de la
seguridad y de sus costos.
Con un uso crítico y creativo del software, confeccionado con los conocimientos de sus
bases, principios y procedimientos, entrada de datos y análisis e interpretación de los
resultados, obtendremos resultados fidedignos. Pero nada sustituye el proceso creativo, el
razonamiento lógico, la interpretación física del fenómeno, el análisis de invariantes y el
examen crítico del problema que un ingeniero pueda hacer.
Diseño automatizado de cimentaciones superficiales con MATLAB
10
Problema científico
Objeto de estudio
Campo de Acción
Objetivo general
Objetivos específicos
estados límites y que incluya los costos directos de la solución, que abarque el caso de base
homogénea y sin presencia de nivel freático.
Elaborar y resolver con dicho sistema programado en MATLAB, un grupo de problemas
tipo en las temáticas que abarca el trabajo, con el fin de validar el programa,
comparándolos con los resultados que se obtienen mediante las hojas de cálculo elaboradas
por otros autores y que permitan verificar los resultados obtenidos.
Justificación de la investigación
Con el desarrollo de una herramienta computacional nueva capaz del diseño de las
cimentaciones superficiales sobre un programa (MATLAB) se tiene otro punto de
comparación, además que por las características del MATLAB puede usarse para optimizar
los resultados en futuros proyectos.
Hipótesis
Tareas científicas
Novedad científica
Realmente el trabajo no presenta una novedad científica notable, ya que existen otros
sistemas computacionales para resolver la tarea que se propone. Tal vez, no se ha
encontrado en la literatura ni en nuestro país un sistema similar programado en MATLAB,
que sea susceptible de ser integrado a un sistema para optimizar conjuntos estructurales que
es lo que se pretende, radicando en este punto el mayor aspecto novedoso del trabajo.
Valor metodológico
Se integran los criterios de diseño de las cimentaciones superficiales, según las últimas
actualizaciones que se dispone en la aplicación del método de los estados límites en el
diseño geotécnico en el país.
Valor práctico
Relevancia Social
Resumen
Introducción
Diseño automatizado de cimentaciones superficiales con MATLAB
14
En este capítulo encontrará una síntesis coherente de los aspectos que abarca el diseño
geotécnico y estructural de cimentaciones superficiales según aparece en la bibliografía
disponible, pero tratando de enfatizar en el empleo de herramientas de software para ello.
En esa dirección se incluye una caracterización consciente del MATLAB, en la que se
resumen las características y potencialidades fundamentales para acometer diseños de
elementos estructurales como este.
a) b)
c)
d) e)
Figura 1.1 Cimentaciones superficiales. Fuente: (Rodriguez, et al., 1989)
Por su morfología (Figura 1.2)
a) recta
b) escalonada
Diseño automatizado de cimentaciones superficiales con MATLAB
17
c) ataluzada
d) aligeradas o nervadas
a) b) c) d) e)
Figura 1.2 Cimentaciones superficiales. Fuente: (Rodriguez, et al., 1989)
Por su forma en planta
a) rectangular
b) cuadrada
c) circular
d) anular
e) poligonal (octogonal, hexagonal)
Vale mencionar que entre las cimentaciones superficiales también se encuentran las zapatas
corridas y las losas o placas. Las zapatas corridas son cimentaciones de gran longitud
comparada con su sección transversal, se utilizan como base de muros o alineación de
pilares. Las losas o placas son elementos de cimentación cuyas dimensiones en planta son
muy grandes comparadas con su espesor y que en general, definen un plano sobre el que
apoyan los pilares o muros de un edificio. (Rodriguez, et al., 1989)
Otro tipo comprenden a los cimientos semiprofundos cuando el terreno apto para cimentar
se encuentra a una profundidad entre 3m y 6m y existe la cimentación por pozos y
cimentación con pedestales intermedios.
Y los profundos que se emplean cuando el suelo o estrato que se encuentra bajo la
estructura no son capaces de soportar la carga, con la adecuada seguridad o con un
Diseño automatizado de cimentaciones superficiales con MATLAB
18
Pilares: Presentan mayor diámetro y permiten una inspección ocular del suelo o roca donde
se apoyarán.
Por características de nuestro país como: los tipos de suelos, la necesidad imperante de
vivienda familiar, las restricciones o limitaciones en la construcción por razones
económicas y climatológicas, es preferible la construcción de cimentaciones superficiales
sobre las semiprofundas y profundas. Podemos notar que es el tipo de cimentación que más
se usa, en comparación con las demás mencionadas y es por tanto nuestro objetivo a
desarrollar.
a) Debe ejecutarse a una profundidad adecuada para impedir los levantamientos las
socavaciones o los daños que puedan producir futuras construcciones cercanas
c) Que no se produzcan asientos superiores a los tolerables por la estructura o por sus
condiciones de servicio.
Con el fin de determinar cuál usar en nuestro proyecto, se propone a continuación una
caracterización breve de cada uno.
Este método se procura asegurar un trabajo de relación tensión vs deformación dentro del
rango de comportamiento lineal del material, y por eso se considera que es suficientemente
seguro. (Hernández Santana & Hernández Caneiro, 2010)
En este método, se asume que el valor de Tensión Admisible definido por el Informe
Ingeniero-Geológico es tal, que al poner a trabajar la base de la cimentación a valores como
máximo igual a su valor, para el caso de combinaciones de carga donde no exista
excentricidad, y en los casos de combinaciones donde exista excentricidad no permitir que
la tensión máxima actuante sea mayor que 1,2 veces la tensión admisible, en dicha base se
garantizan tanto las condiciones de diseño relacionada con el criterio de estabilidad, como
las relacionadas con el criterio de deformación, y por tanto el diseño solo consiste en
comprobar que las tensiones actuantes no sobrepasen los valores establecidos con respecto
a la tensión admisible. (Quevedo Sotolongo, 2010)
En este método se parte de recoger en un único coeficiente todas las incertidumbres que
están presentes en la obtención de los distintos parámetros que intervienen en el diseño, y
que determinan cuánto se aleja el diseño de la zona de falla. (Hernández Santana &
Hernández Caneiro, 2010)
Diseño automatizado de cimentaciones superficiales con MATLAB
20
En este método para determinar el área de la base según el criterio de estabilidad se usan las
solicitaciones y los valores de las propiedades físico-mecánicas del suelo con sus valores
normativos.
Para el diseño o chequeo del criterio de deformación en el método del factor de seguridad
en este método se asume que con los factores de seguridad global aplicados, se garantiza
que el estado tensional en la base de la cimentación , para las combinaciones de carga por
deformación, tenga un comportamiento lineal, y por tanto no hay que chequear la condición
de linealidad.: Sin embargo para suelo como las arenas sueltas o las arcillas blandas si se
recomienda chequear la posible falla local, según lo establecido, lo que ya fue demostrado
que es prácticamente chequear la condición de linealidad. (Quevedo Sotolongo, 2010)
En diseños con bases en la Teoría de Seguridad, los parámetros que intervienen son
tratados como variables aleatorias. En esta teoría el nivel de seguridad es definido en
Diseño automatizado de cimentaciones superficiales con MATLAB
21
El MEL, es en el que se basan nuestras normas vigentes y propuestas por aprobar referentes
al diseño de cimentaciones. Se encuentran los primeros trabajos realizados en Cuba, a
finales de la década de los ´80. (Quevedo Sotolongo, 1994). Llevándose a cabo
posteriormente toda una serie de investigaciones en la misma dirección, que han permitido
tener la base teórica para la introducción general de los estados límites en el campo de la
geotecnia.
Los estados límites se suelen dividir en dos tipos: “estados límites últimos” y “estados
límites de servicio”. Se entiende que se ha producido un estado límite último en una
Diseño automatizado de cimentaciones superficiales con MATLAB
22
El método de estados límites realiza el diseño como si fuera a ocurrir la falla y se garantiza
su no existencia con la introducción de los coeficientes parciales de seguridad, los cuales se
presentan a continuación:
Este método establece dos condiciones límites de diseño, las cuales abordaremos a
continuación.
Este estado incluye varias condiciones que deben ser verificadas. Entre ellas
podemos señalar el chequeo al vuelco (Quevedo Sotolongo, 2007):
En el diseño o revisión de las áreas de las bases de las cimentaciones, según el estado límite
de estabilidad debe garantizarse que dicha cimentación sea segura al posible vuelco,
analizando todos los factores que inciden en él (Figura 1.4):
dónde:
Diseño automatizado de cimentaciones superficiales con MATLAB
23
M estabilizantes – es la sumatoria de los momentos de todas las fuerzas, con sus valores
característicos, que se oponen al vuelco de la cimentación con respecto a la esquina del
mismo.
M desestabilizantes – es la sumatoria de los momentos de todas las fuerzas, con sus valores
característicos, que provocan el vuelco de la cimentación con respecto a la esquina del
mismo.
Figura 1.4 Parámetros necesarios para chequear el factor de seguridad al vuelco. (Quevedo
Sotolongo, 2007)
Para los casos de cimientos rectangulares con pedestal céntrico, la condición (1) se cumple
automáticamente cuando la excentricidad es menor o igual que un tercio del ancho o largo
de la base por lo que para esto bastaría con chequear la siguiente condición:
Diseño automatizado de cimentaciones superficiales con MATLAB
24
𝑙 𝑏
𝑒𝑙 ≤ ó 𝑒𝑏 ≤
3 3
𝐻 ∗ ≤ 𝑁 ∗ ∗ 𝑡𝑔𝜑 ∗ + 0.75 ∗ 𝑏΄ ∗ 𝑙 ΄ ∗ 𝑐 ∗
Con la cual se garantiza el equilibrio entre la acción de las fuerzas cortantes externas y la
resistencia al deslizamiento desarrollada en la superficie de contacto entre la cimentación y
el suelo de la base. Esta función de las características físico-mecánicas del suelo,
determinadas para una probabilidad de diseño del 95%.
𝑁 ∗ ≤ 𝑄𝑏𝑡 ∗
𝑞𝑏𝑟∗ −𝑞 ΄∗
𝑄𝑏𝑡 ∗ = 𝑏΄ ∗ 𝑙 ΄ ∗ ( + 𝑞 ΄∗ )
𝛾𝑠
Suelo y C - .
Diseño automatizado de cimentaciones superficiales con MATLAB
25
𝑞𝑏𝑟∗ = 0.5 ∗ 𝛾2 ∗ ∗ 𝐵 ΄ ∗ 𝑁𝛾 ∗ 𝑠𝛾 ∗ 𝑖 𝛾 ∗ 𝑑𝛾 ∗ 𝑔𝛾 + 𝑐 ∗ ∗ 𝑁𝑐 ∗ 𝑠𝑐 ∗ 𝑖 𝑐 ∗ 𝑑𝑐 ∗ 𝑔𝑐 + 𝑞 ΄∗ ∗ 𝑁𝑞 ∗ 𝑠𝑞
∗ 𝑖 𝑞 ∗ 𝑑𝑞 ∗ 𝑔𝑞
Suelo C ( = o).
𝑞𝑏𝑟∗ = 5.14 ∗ 𝑐 ∗ ∗ (1 + 𝑠𝑐 ΄ + 𝑑𝑐 ΄ − 𝑖 𝑐 ΄ − 𝑔𝑐 ΄ ) + 𝑞 ΄∗
𝑞 ΄∗ = 𝛾1 ∗ ∗ 𝑑 + 𝑞𝑠𝑐 ∗
De la misma manera que la flecha de una viga puede ser un factor de limitación en el
proyecto estructural, el asentamiento del suelo producido por las cargas es con frecuencia el
factor predominante en el proyecto de las cimentaciones. (Sowers & Sowers, 1978).
𝑁
𝑝 ≤ 𝑅΄ 𝑑𝑜𝑛𝑑𝑒: 𝑝 =
𝑏 ∗𝑙
Las presiones actuantes se determinan teniendo siempre presente que las mismas están en
función de la combinación de carga para el diseño por deformación con sus valores
característicos.
La presión del Límite de Linealidad del suelo R´ está en función de las características
físico-mecánicas del suelo, determinadas para una probabilidad de diseño del 85%.
Se debe de garantizar que la base del cimiento este comprimida y no haya peligro al vuelco.
M estabiliza ntes
F .Svuelco 3
M desestabil izantes
n
Hi
SC i 1 6
( iS 4 iC iI )
dónde:
n – cantidad de estratos por debajo del nivel de cimentación hasta una profundidad igual a
la Potencia Activa.
Diseño automatizado de cimentaciones superficiales con MATLAB
27
Hi – espesor del estrato (i) existente por debajo del nivel de cimentación hasta una
profundidad igual a la Potencia Activa.
Vale aclarar que para nuestro proyecto solo trataremos un solo estrato, por tanto la
ecuación se simplifica al ser n=1.
Se recomienda el cálculo de los asentamientos por métodos No Lineales para los siguientes
casos.
siendo:
q1, q2- componentes de la tensión de confinamiento lateral del suelo: q1 que representa la
tensión mínima y q2 como la tensión máxima.
𝜇 1 − 𝑠𝑒𝑛𝜑 2 ∗ 𝑐 ∗ 𝑐𝑜𝑠𝜑
𝑞1 = ∗ 𝑅΄, 𝑞2 = ∗ 𝑞𝑏𝑟 −
1−𝜇 1 + 𝑠𝑒𝑛𝜑 1 + 𝑠𝑒𝑛𝜑
p - presión media real actuante en el suelo, para la cual se van a determinar los
asentamientos No Lineales, tomando la N actuante definida para la combinación de carga
del 2do Estado Límite, y el área de la base obtenidas para las condiciones del 1er Estado
Límite .
Diseño automatizado de cimentaciones superficiales con MATLAB
29
La Presión bruta media actuante en el centro del cimiento, para la determinación del
Slineal, se tomará igual al valor de la Tensión Límite de Linealidad (R´), determinada para
una probabilidad de diseño del 85%. (p = R´) y en función de las correcciones
establecidas por excentricidad.
Para el cálculo del asentamiento Lineal bajo el Punto Característico (PC máx ) es necesario
definir la presión bruta que se considerará actúa en este punto, según la Figura 1.6 y que
puede ser determinada por la expresión.
12 ∗ 𝑀
𝑃𝑧 = 𝑅 ΄ + ∗ 0.37 ∗ 𝑙
𝑏 ∗ 𝑙3
Figura 1.6 - Asentamiento Lineal bajo el Punto Característico de máxima tensión, PC máx .
Fuente: (Quevedo Sotolongo, 2007)
El Asiento No Lineal se calculará bajo el Punto Característico Máximo (PC máx ), y por tanto
la presión bruta actuante considerada para el cálculo será la Presión actuante bajo este
punto del cimiento, debido a las solicitaciones reales del 2 do Estado Límite.
𝑁
𝑃𝑧 𝑛𝑜 𝑙𝑖𝑛𝑒𝑎𝑙 =
𝑏∗𝑙
Diseño automatizado de cimentaciones superficiales con MATLAB
30
Este es el método de seguridad que vamos a incluir en el diseño, ya que además de der el
usado en nuestro país, proporciona un diseño más racional y no necesita de un complejo
análisis estadístico como el de la teoría de la seguridad.
Diseño Estructural.
Para el diseño estructural nos proponemos solo el diseño del plato de la cimentación, aquí
se determinan el peralte, la cuantía, el diámetro, la cantidad de barras, y la posición del
refuerzo, en caso que este sea necesario. (Quevedo Sotolongo, 2010)
El diseño estructural se realizara tomando como base las cargas en sus combinaciones
pésimas para cada aspecto a diseñar, al evaluar los estados limites relativos a capacidad
portante, mientras que la comprobación del estado límite de fisuración, se tomaran las
cargas de servicio, en sus combinaciones más desfavorables todo ello de acuerdo con los
requisitos apropiados de diseño y conforme a lo estipulado en la (Chagoyén & Broche,
2002)
Diseño plato.
Elemento flexible: se entiende cuando existe compatibilidad entre las deformaciones del
plato de la cimentación y de la base de la misma. Para obtener los peraltos en cimientos
flexibles, partimos de los peraltos mínimos para no incumplir esta condición, después
tenemos en cuenta el criterio de punzonamiento, criterio cortante y el criterio de flexión
negativa, y si se trata de una solución masiva de hormigón hay que determinar el peralto
por flexión positiva y decantarnos por el mayor de estos.
Este se calcula a partir del momento último que se calcula en el punto por flexión positiva,
con el cual se obtienen aéreas de acero refuerzo en ambas direcciones y deberán siempre
ser comparadas con un área de acero mínima Aa min 0.002Lh donde h d rec , esto en la
dirección de L, de forma homóloga en la dirección de B. También se considera un refuerzo
superior que dependerá si lo necesita por flexión negativa, la solución para este problema
puede ser tratada también con un aumento en el peralto.
Deberá tenerse en cuenta la necesidad de acero refuerzo en las dos direcciones, o sea se
calculará un área de acero para la dirección de L y una para la dirección de B. Esto se hará
a partir de una fuerza máxima de tracción provocada por las cargas externas actuantes:
FL max y FB max .
En nuestro país se han interesado en el tema de las cimentaciones los investigadores con el
fin de concebir un programa que resuma el diseño de las mismas.
Diseño automatizado de cimentaciones superficiales con MATLAB
32
El programa de computación Galileo fue realizado por (Gómez & Morales, 2000) el cual
aborda el tema de la optimización, elaboraron un modelo matemático cuyo núcleo consistió
en un algoritmo genético para resolver la optimización de vigas de Hormigón armado y
cimentaciones aisladas superficiales. Confeccionan el software de diseño y optimización de
vigas de secciones ¨T¨, rectangulares y de cimentaciones aisladas rectangulares, la función
objetivo fue la de costo, pero aun el sistema de restricciones y ecuaciones de estado no
refleja un nivel actual de los conocimientos sobre diseño geotécnico y estructural de estos
elementos.
Existen dos programas computacionales en lenguaje MSdos los cuales fueron trabajos de
diploma realizados por la Dra. Ing. Ana Virginia Gonzáles Cueto y el Dr. Ing. Luis Orlando
Ibáñez los cuales se encuentran adecuados a la norma cubana y se utilizan como método
práctico en la docencia; el DGCim para desarrollar el diseño geotécnico teniendo en cuenta
los dos estados límites de resistencia y deformación y el DECim con el cual se realiza el
diseño estructural de la cimentación. Estos programas permiten entrada de datos y
devuelven resultados pero en ellos no es posible visualizar todo el proceso de cálculo
necesario.
El GEO es un programa para diseño y análisis Geotécnico, dentro de sus diversas funciones
ofrece programas para diseñar cimentaciones de los tipos superficiales y profundos, estos
determinan la capacidad de carga y diseñan el refuerzo para las secciones transversales.
El programa que tiene relación con los cimientos diseña y verifica cimentaciones aisladas
ya sea de forma rectangular o circular descansando sobre un subsuelo a capas. El programa
realiza la verificación de vuelco y deslizamiento de la cimentación y determina
Diseño automatizado de cimentaciones superficiales con MATLAB
34
obtener la respuesta requerida, además los resultados pueden modificarse en función del
código de diseño utilizado y del tipo de cimentación, lo que los hace más apropiados y
exactos para situaciones específicas de diseño. Esta flexibilidad que permite el programa
hace que se emplee frecuentemente en el campo de la ingeniería para la realización de
cálculos especializados.
Los parámetros de entrada del programa en cuanto a los datos del suelo están determinados
por los resultados que se obtienen en los ensayos del cono de penetración. Este software
puede usarse para diseñar cimentaciones profundas o superficiales para este último caso las
dimensiones del cimiento pueden ser racionalizadas, se requiere de entrada el ancho; la
estabilidad y la capacidad de carga son chequeadas como parte del proceso de diseño. Las
regularidades del diseño están dadas para el código base NEN 6744 sobre cimentaciones
superficiales.
Es otra de estas ayudas que según la promoción hecha en su página de ventas, analiza,
diseña, y dibuja zapatas aisladas y sus pedestales, llegando a manejar situaciones complejas
como flexión biaxial, levantamiento y propiedades de frontera.
electrónicas a las que puedan acceder con facilidad y brinde mayores dividendos a los
estudiantes e ingenieros.
MATLAB VS EXCEL
Con el fin de determinar el programa óptimo para realizar el diseño de las cimentaciones
superficiales pretendemos comparar MATLAB (moderno software con las potencialidades
para la optimización como el que estamos buscando) con otros dos conocidos software con
capacidad de programar.
Podemos comparar MATLAB con EXCEL que es uno de los programas más conocidos y
aplicados. Mientras que MATLAB es más rápido, poderoso, comprensivo, tiene una
programación más fácil, y un ambiente más flexible que las hojas cuadriculadas de
EXCEL. EXCEL es más familiar a los ingenieros en su uso de manera general, su forma
principal son hojas de cálculo, la optimización encogida es más fácil, pero la programación
en él es mucho más compleja.
MATLAB VS MATHCAD
Una comparación más compleja sería con Mathcad, ya que este ha demostrado ser una
herramienta muy útil a los ingenieros, con un diseño eficaz de su hoja de cálculo, donde las
modificaciones actualizan toda la hoja al instante.
Diseño automatizado de cimentaciones superficiales con MATLAB
37
Mathcad:
Permite crear documentos que combinan texto y matemática (significando las ecuaciones,
normalmente), y los gráficos se visualizan en la misma área de trabajo.
MATLAB:
Tiene varias alternativas para crear curvas y plotear variables, muy útiles para graficar
resultados hasta de tres dimensiones para representar.
Caracterizar un programa que abarca tantos aspectos con MATLAB no es nuestro objetivo
principal, mientras que si lo son las definiciones principales del mismo, su modo de trabajo
y entorno de trabajo. Según (Wilson, 1999) como cada programador, MATLAB empezó
Diseño automatizado de cimentaciones superficiales con MATLAB
38
como una manera simple de hacer las cosas poderosas, y se ha vuelto una manera no tan
simple de hacer las cosas muy poderosas.
En (Smith, 2010) se dice que MATLAB se desarrolló para ingenieros para crear,
manipular, y visualizar las series de matrices. En un nivel básico puede realizar las mismas
funciones de una calculadora científica, pero MATLAB se ha extendido más allá de sus
capacidades originales y ahora ha proporcionado un sistema interactivo y programando, con
muchas aplicaciones, incluyendo el análisis financiero así como el cómputo científico y
técnico. Es un programa interactivo para la computación numérica y la visualización de
datos. Permite resolver complicados problemas numéricos sin necesidad de escribir un
programa donde se encuentren implicados elevados cálculos matemáticos y la visualización
gráfica de los mismos.
Los aspectos fundamentales que es necesario dominar para poder programar en MATLAB
son el ambiente o entorno de trabajo con sus diferentes ventanas que serán detalladas
posteriormente; así como sus comandos principales, el concepto de funciones, como se
maneja en MATLAB, su estructura de programación, las estructuras básicas de los
algoritmos lineales y/o no lineales, las sentencias que permiten hacer ciclos y muy
importante el trabajo con arreglos y otros tipos de variables. (Alfio Quarteroni, 2010).
Una ventaja que nos proponemos explotar es que dispone en la actualidad de un amplio
abanico de programas de apoyo especializados, denominados Toolboxes (cajas de
herramientas), que extienden significativamente el número de funciones incorporadas en el
programa principal.
Manipulación de matrices.
Encontrando las raíces de polinomios.
Aplicaciones de procesos digitales (la caja de herramientas).
Trazando: los x-y y los gráficos polares 3D.
Desarrollar algoritmos.
Modelación, simulación y diseño.
Adquisición de datos.
Análisis de datos, exploración, y visualización,
Desarrollo de aplicaciones, incluyendo la construcción de interfaces gráficas de
usuario.
MATLAB ofrece varias ventajas a los usuarios encima de los lenguajes convencionales:
objetos y rasgos. Los gráficos, desplegando vectores y matrices como gráficos, así como
anotando e imprimiendo estos gráficos. Y por último la Interfaz de Programas de
Aplicación (el API).
Hoy hay abundancia de programas con un alto nivel de programación, que son diseñados
para diferentes tipos de aplicaciones. MATLAB es diferente a otros programas en la
esencia de que en su forma básica es un programa no binario, creado por un recopilador.
Cada comando es interpretado y se lleva a cabo por el mismo en el orden en que fue escrito.
(Gustafsson, 2011).
Para trabajar en modo programado hay que utilizar los llamados m–ficheros.
Los scripts están formados simplemente por una serie de instrucciones en MATLAB que se
ejecutan como si estuviéramos en modo interactivo; Un script es un m–fichero que agrupa
Diseño automatizado de cimentaciones superficiales con MATLAB
41
Dentro de estas hay que prestar especial cuidado al vectorizar las variables que como
resultado sean matrices o como las combinaciones de carga que tengan que ser tratadas en
la misma operación varios valores diferentes por estar comprendida dentro de un arreglo.
En estos casos debemos vectorizar el producto (.*), la división (. /), y la exponenciación (.
^), pues la suma o la resta son operaciones que se realizan de forma vectorizada.
Las variables empleadas en una sesión están almacenadas en la memoria a la que se accede
desde la consola (en el workspace) están son las variables que vamos a necesitar pues son el
argumento de entrada y salida de las funciones, todas las variables globales están
almacenadas en una memoria global a la que se puede acceder parcialmente a través del
comando global.
1.6. Cálculo de los costos directos a partir de los reglones variantes del PRECONS.
Uno de los objetivos que tenemos propuestos es que además de los resultados del diseño de
las cimentaciones, obtener dentro del mismo programa los costos directos del servicio de la
construcción para esa solución obtenida; estos se calcularán sobre la base de los
subsistemas de costos y precios, los procedimientos y regulaciones establecidas en el
PRECONS II.
Diseño automatizado de cimentaciones superficiales con MATLAB
42
Con vista a desarrollar una ficha de costos directos que cumpla con lo estipulado en nuestro
país tenemos que hacer un estudio de la principal herramienta, el PRECONS.
Los gastos directos de producción que tienen lugar en la ejecución de una obra se clasifican
en las siguientes partidas de costo:
Materiales
Mano de Obra
Uso de Equipos
Otros Gastos Directos de Obra
Gastos Generales Directos de Obra (Centro de información de la construcción.,
2005)
El contenido de las partidas de costos directos, que componen los costos de los renglones
variantes del PRECONS II se basan en las premisas siguientes:
Para obtener una metodología que determine los costos directos, elaboramos una secuencia
de pasos para obtener las actividades necesarias, los volúmenes y los renglones variantes:
1- Definir los objetos de obra o partes componentes de los trabajos a realizar: Se divide
la obra general en objetos de obra y estos se dividen en etapas, ya sea; Movimiento de
Diseño automatizado de cimentaciones superficiales con MATLAB
43
3- Localizar y organizar los reglones variantes con sus respectivos costos unitarios, estos
deben dar cumplimiento a la secuencia de actividades.
4- Multiplicar cada uno de los costos unitarios por el volumen de trabajo que
corresponda y organizarlos.
5- Obtener los costos directos de materiales, mano de obra y uso de equipos, resultantes
de la suma de estos en cada uno de los renglones variantes.
6- Calcular “Otros Gastos Directos de Obra” 12.5% de los costos directos totales.
7- Calcular “Otros Gastos Directos de Empresa” 10% de los costos directos totales.
Tenemos que señalar que estos costos no incluyen factores como son, gastos de
administración, dirección, aseguramiento, gastos de comunicaciones, tramitaciones legales,
mantenimiento de equipos, reposición de herramientas, utensilios y otros como talleres a
pie de obra, que si se consideran en gastos indirectos y presupuestos independientes.
Diseño automatizado de cimentaciones superficiales con MATLAB
44
El método de los Estados Límites aplicado en Cuba por la norma cubana es el más racional
de todos los métodos de introducir la seguridad en el diseño de cimentaciones y en la
actualidad prácticamente es el método de diseño que con más frecuencia se incorpora a las
actualizaciones de las normas de diseño utilizadas en el mundo. Se considera que es un
diseño donde las cargas y las tensiones a las que está sometido el material que se emplee en
el elemento a diseñar, así como las deformaciones y desplazamientos que en ella se
originan, tanto en el período de construcción como durante su vida útil, estén cerca de los
límites permisibles para cada caso, sin llegar a sobrepasarlos.
Existen varios trabajos y programas tanto en el ámbito nacional como internacional, que
tratan el diseño geotécnico y estructural de cimentaciones rectangulares aisladas, teniendo
como base la entrada de datos y obteniendo una respuesta, pero no existe uno con soporte
en un programa que permita la optimización de conjuntos estructurales como puede
realizarse desde MATLAB o EXCEL y permitan combinar sus resultados con los
resultados de problemas de análisis y diseño de los elementos de la superestructura como
por ejemplo SAP2000, ANSYS, etc.
Aunque los costos directo no representan la totalidad del costo de los servicios de la
construcción, estos si representan los gastos más inmediatos e invariantes como los
materiales, mano de obra y equipos. Otros costos indirectos no están dentro de nuestro
Diseño automatizado de cimentaciones superficiales con MATLAB
45
objetivo pues estos se nos escapan al depender de factores como la empresa constructora
que ejecute la obra y los trámites bancarios que se realicen.
Diseño automatizado de cimentaciones superficiales con MATLAB
46
Referencias
Alfio Quarteroni, F. S. P. G., 2010. Scientific Computing with MATLAB and Octave.
Berlín: Springer.
Ballate Reina, D. R., 2016. Costo de los servicios de construcción y montaje. Ejemplo
Resuelto.. Santa Clara, UCLV.
Quevedo Sotolongo, G., 2007. Propuesta de norma no publicada para el diseño geotécnico
de cimentaciones superficiales por estados límites.. UCLV ed. Santa Clara(Ciudad de La
Habana): Ministerio de la Construcción.
Quevedo Sotolongo, G., 2010. Conferencia 14 Met Generales de diseño. Santa Clara,
UCLV, pp. 7-10.
Quevedo Sotolongo, G., 2010. Introducción al diseño estructural de cimientos. Santa Clara,
UCLV.
Rodriguez, J. M. O., Serra, J. G. & Oteo, C. M., 1989. Curso aplicado de cimentaciones,
cuarta edición. Cuarta Edición ed. Barcelona: GRAFICINCO,S.A.
Wilson, W. &., 1999. MATLAB. Dr. Dobb’s Journal, 7 enero, pp. 50-62.
Diseño automatizado de cimentaciones superficiales con MATLAB
48
El diseño geotécnico consiste en obtener los lados de la base que satisfagan las
restricciones por estabilidad: entiéndase chequeo al vuelco, al deslizamiento y capacidad de
carga y las restricciones por utilización, o sea, deformación, donde el asentamiento debe ser
menor que el permisible según sea el caso de edificación a la que va a pertenecer el
cimiento.
Manualmente hay que introducir los datos del problema como profundidad de cimentación,
la rectangularidad que tendrán los lados de la base (b/l), la distancia entre la base de
Diseño automatizado de cimentaciones superficiales con MATLAB
49
cimentación y el cortante que llega al pedestal o vaso (Hc), si existe o no profundidad del
cimiento dentro del estrato resistente (D), la inclinación del terreno, la existencia de
sobrecarga, las condiciones de trabajo de la cimentación y el tipo de fallo.
También el programa pedirá que especifique las propiedades físico-mecánico del suelo del
estrato por debajo del nivel de la cimentación sobre el cual apoya la estructura, dígase
cohesión, ángulo de fricción interna, y peso específico del suelo. De las propiedades del
suelo por encima del nivel de cimentación se definirá el peso específico del relleno.
Como se mencionó anteriormente una parte de los datos se introduce por una hoja Excel
(Figura 2.1), donde se especifican las cargas y combinaciones que quiera el usuario utilizar,
esta interfaz es necesaria para desde aquí obtener los resultados de otro software con los
resultados de problemas de análisis y diseño de los elementos de la superestructura como
por ejemplo SAP2000 o Staad.pro:
Figura 2.1
Fuente: elaboración propia.
Dentro del scripts principal estos datos se adquieren por la función xlsread, donde hay que
incluir el nombre del documento Excel, el número de la hoja dentro del mismo y el rango
de filas o columnas que se quiere leer Valores_normativos=xlsread('Datos de entrada
DiGEC',1,'J2:N100');. Están propuestas hasta un máximo de cien combinaciones, pero con
Diseño automatizado de cimentaciones superficiales con MATLAB
50
Figura 2.2
Fuente: elaboración propia.
Nótese que la entrada como tal es numérica, porque aunque MATLAB si puede tener
variables definidas como texto con la función string, pero el trabajo posterior no procede,
las condiciones tienen que hacerse con valores numéricos.
Figura 2.3
Fuente: elaboración propia.
Para detallar en las funciones que determinan los coeficientes de minoración y las
propiedades físico- mecánicas del suelo, remitirse al anexo 1.
Figura 2.4
Fuente: elaboración propia.
Figura 2.5
Fuente: elaboración propia.
Se incrementa el valor del lado b y por consiguiente el lado l, hasta que la capacidad del
suelo de soportar las fuerzas horizontales superen a las actuantes, nótese que se usan los
valores de cálculo de las cargas, e igualmente el programa comprueba todas las
combinaciones. Dentro del programa, la función para lograr esto usa la variable
Cantidad_combinaciones que se repetirá en las demás, esta variable está definida como la
longitud de la matriz que recoge las combinaciones. Figura 2.6.
Figura 2.6
Fuente: elaboración propia.
Para el chequeo final del primer estado límite, la capacidad de carga del suelo de la base de
la cimentación se usa el método de Brinch Hansen que establece la norma, teniendo en
cuenta que la capacidad de carga de la base de la cimentación qbr cambia según los
parámetros C-Fi del suelo. Figura 2.7
Figura 2.7
Fuente: elaboración propia.
dganma, dc, dq- Factores que valoran el efecto de la profundidad dentro del estrato
resistente.
gganma gc, gq- Factores de inclinación del terreno.
Todos esto factores se calculan dentro de la misma función, pues algunos se modifican al
aumentar los lados para salir de la condición. El desarrollo de estas se encuentran en el
anexo 2.
qbr * q´*
Qbt * = b´·l´· q´*
s
La condición que hace aumentar los lados del cimiento, hasta que no se cumpla en este caso
es:
𝑁 ∗ > 𝑄𝑏𝑡 ∗
Mientras esta afinación sea cierta se incrementarán las dimensiones, Figura 2.8.
Figura 2.7
Fuente: elaboración propia.
Figura 2.6
Fuente: elaboración propia.
Para el segundo estado límite el programa necesita que el usuario le ofrezca nuevos datos.
Para determinar los coeficientes de las condiciones de trabajo del suelo y el tipo de
estructura, ganmac1 y ganmac2, tiene que especificar el tipo de suelo y le relación longitud,
altura de la estructura de la misma. Otro dato será el coeficiente de confiabilidad de los
datos. Como se propone el cálculo de los asentamientos con métodos que tienen al módulo
de deformación del suelo como parámetro deformación, se pedirá dicho módulo de
deformación y el coeficiente de Possion del suelo. Figura 2.7
Figura 2.7
Fuente: elaboración propia.
Figura 2.8
Fuente: elaboración propia.
En el estado límite de deformación se cheque nuevamente el vuelco partiendo del lado que
se obtuvo en el diseño por primer estado límite, Figura 2.9, de no cumplir se aumentan las
Diseño automatizado de cimentaciones superficiales con MATLAB
59
dimensiones que obtuvimos del primer estado límite. Luego se procede a la comprobación
del estado tenso deformacional lineal del suelo (Condición de Linealidad) , esta condición
se puede considerar cumplida cuando las tensiones actuantes no sobrepasen el valor de la
presión del límite de linealidad del suelo(R’) es decir:
p R´
Figura 2.9
Fuente: elaboración propia.
R´ y p, tienen indicado su modo de obtención en el Capítulo I de este trabajo, si la presión
actuante es menor que la presión límite de linealidad, el cálculo de los asentamientos se
trataran considerando un comportamiento lineal del suelo, de lo contrario se determinaran
los asentamientos para una etapa no lineal del suelo. Figura 2.10
Figura 2.10
Fuente: elaboración propia.
Para cualquier comportamiento del suelo, lo primero será determinar la profundidad por
debajo del nivel de solera donde los desplazamientos o deformaciones sean apreciables, a
esta profundidad se le llamará potencia activa, (Ha). Esta depende del suelo que defina cuál
es la potencia activa final, donde se cumple la siguiente condición:
Diseño automatizado de cimentaciones superficiales con MATLAB
60
Si en el límite hallado por el criterio anterior existen suelos muy compresibles (E0 5000
kPa) entonces se tomará como potencia activa aquella profundidad donde se cumple que:
Si a una profundidad inferior a la determinada por uno de los criterios anteriores existieran
suelos, semi-rocas o rocas poco compresibles (Eo 100,000 KPa) se tomará esta
profundidad como potencia activa.
siendo:
Figura 2.11
Fuente: elaboración propia.
Diseño automatizado de cimentaciones superficiales con MATLAB
61
El factor de influencia Jz, coeficiente que se determina en función de la relación l/b y Ha/b
en los cimientos rectangulares, siendo Ha la distancia vertical entre el punto característico y
el punto donde se desea calcular para el punto característico. El factor Jz de un área
rectangular uniformemente cargada se determinó a partir del punto esquina por la
superposición de los esfuerzos en los puntos esquinas de las cuatro áreas elementales en
que queda dividido el cimiento por el punto característico.
Los asentamiento lineales para cuando M=0, se calculan por el método de la sumatoria de
capas, aunque dicha sumatoria como sólo tratamos un estrato, se simplifica a sólo un
promedio pesado de los asentamientos superior, 4 veces el central y el inferior del estrato,
Figura 2.12
Figura 2.11
Fuente: elaboración propia.
Cuando M≠0 además de los asentamientos producidos por los esfuerzos verticales se le
suman los producidos por el giro, aplicando el principio de superposición de efectos. Donde
los asentamientos absolutos producidos por el momento se calcularan por el método de
solución cerrada. Aquí se suman además del asentamiento lineal, el asentamiento por el
giro Sgiro:
dónde:
(1 − 𝑣 2 ) 𝑀𝑙 (1 − 𝑣 2 ) 𝑀𝑏
𝑡𝑎𝑛 𝛼𝑙 = 8 ∗ 𝑘1 ∗ ∗ 3 , 𝑡𝑎𝑛 𝛼𝑏 = 8 ∗ 𝑘2 ∗ ∗ 3
𝐸 𝑙 𝐸 𝑏
dónde:
Diseño automatizado de cimentaciones superficiales con MATLAB
62
k1, k2 – coeficiente que depende de la relación 1/b y la potencia activa (HA). Estos
coeficientes se determinan a partir de una interpolación lineal.
v – coeficiente de Poisson.
La función para interpolar los valores de k1 y k2, dependen de introducir una matriz y los
intervalos en que varían por cada eje, Figura 2.12
Figura 2.12
Fuente: elaboración propia.
Donde V1, X1, Y1, V2, X2, Y2, están definidas en las tablas:
21
P P1
S nolineal hm
3E
P 2 P
P2 P1q2 q1
1 2
P2 P1
P2 P1 2q2 q1 S lineal
[1]
P1- valor de la presión crítica inferior, puede estar dada por R, 1.2 R o por otro valor
predeterminado, que prefije un límite de comportamiento lineal del suelo.
P2- valor de la presión crítica superior, y va a estar dada por la expresión de capacidad de
carga (qbr) evaluada para los valores medios de las características físicas mecánicas del
suelo.
q1, q2- componentes de la tensión de confinamiento lateral del suelo. q1 que representa la
tensión mínima y q2 como la tensión máxima actuante lateralmente.
P- presión media real actuante en el suelo, para la cual se van a determinar los
asentamientos No Lineales por la expresión
La función para encontrar los asentamientos no lineales también tiene que considerar la
presencia o no de momento, Figura 2.13
Figura 2.13
Fuente: elaboración propia.
12M
pz R´ 0.37l
bl 3
Mientras que la presión que se usa para el cálculo de los asentamientos absolutos no
lineales es:
N
pz NoLineal
bl
Diseño automatizado de cimentaciones superficiales con MATLAB
65
El algoritmo para encontrar las dimensiones para el segundo estado límite, se basa en los
lados obtenidos en el primero:
1. Determinar las combinaciones de carga para el diseño para el segundo estado límite
según lo establecido, utilizando las solicitaciones normativas, y las características
físico-mecánicas del suelo para una probabilidad del 85%.
2. Se comprueba que para las combinaciones del segundo estado límite sea seguro al
vuelco. De no cumplir, se aumenta en 0.05 metros el lado b y por consiguiente hay
un aumento en el lado l al ser este último l=b/rectangularidad.
3. Calcular la distribución neta de contacto (p) para las combinaciones de carga
correspondientes.
4. Definir el comportamiento tenso-deformacional del suelo (Lineal o no lineal), en
función de la tensión límite de linealidad (R’).
5. Utilizando las características físico-mecánicas del suelo con sus valores
característicos, las solicitaciones actuantes normativas correspondientes a este
estado límite y las dimensiones del cimiento calculadas según el criterio de
estabilidad, se procede al cálculo de las deformaciones o desplazamientos de la base
del cimiento (asentamiento o giro). En función también de la existencia de momento
o no.
6. Determinar la potencia activa y Jz en función de la rectangularidad y la
profundidad.
7. Calcular las deformaciones y desplazamientos de la base de los cimiento, en su
etapa lineal o no lineal, según corresponda.
8. Una vez calculados los asentamientos absolutos, se comparan con los admisibles y
en caso de ser mayores se aumentan las dimensiones de la base, sumándoles 0.05
metros al lado b hasta que el asentamiento absoluto sea menor que el admisible.
Diseño automatizado de cimentaciones superficiales con MATLAB
66
Así como con el diseño geotécnico se establecen las dimensiones definitivas de los lados
para que la base de la cimentación cumpla con los criterios de rotura y de utilización, con el
diseño estructural nos proponemos encontrar el peralto del plato y el refuerzo si es el caso
Diseño automatizado de cimentaciones superficiales con MATLAB
67
para que el cimiento resista las solicitaciones de carga y cumpla efectivamente su función
de transmitir las cargas de la superestructura.
Este diseño comienza con la entrada de los datos necesario para esta etapa. Figura 2.14,
como las dimensiones del pedestal, los recubrimientos, la resistencia del hormigón, las
tensión de fluencia del acero y el tipo de barra que se quiere usar, dígase diámetro y
especificar si es lisa o corrugada.
Figura 2.14
Fuente: elaboración propia.
Figura 2.14
Fuente: (Chagoyén & Broche, 2002)
Las funciones para determinar el esfuerzo actuante y los peralto mínimos pueden
encontrarse en el anexo 4.
Criterio de punzonamiento:
Al escoger esta variante el programa solo ejecuta las funciones necesarias para este tipo de
cimiento, obtiene el peralto mínimo para que desde un inicio esta condición esté cumplida.
A continuación se procede a chequear el punzonamiento, mientras los esfuerzos cortantes
provocados por el punzonamiento (Taopz) sea mayor a 0.75 veces el esfuerzo resistente
(Rpz), se le incrementa 0.01 metro al peralto, Figura 2.15
Figura 2.15
Fuente: Elaboración propia.
Figura 2.16
Fuente: (Chagoyén & Broche, 2002)
2
0.16 1 f 'c
d
R pz 0.16 1 f ' c ,
2b0
0.34 f ' c
dónde:
- Relación del lado mayor al lado menor del área cargada, en este caso, del pedestal,
lc
columna o vaso bc .
N pz L ML B MB
aopz x y
A pz JL JB
dónde:
( b d ) (l c d )
N pz N' 1 ( c )
BL
ML, MB -momentos resultantes provocados por la carga externa, en las direcciones de los
x, y-posición del punto de la sección crítica donde pz es máximo, expresadas en el sistema
de coordenadas de ejes ortogonales entre si y paralelos a las direcciones en que actúan las
solicitaciones de flexión y por tanto a los lados de la sección crítica, en metros.
L, B- fracción del momento desequilibrante trasmitido por excentricidad de cortante en la
unión columna, pedestal o vaso-plato en la sección crítica supuesta. Para las columnas,
pedestales o vasos circulares se sustituirá la misma por una cuadrada de igual área.
1 1
L 1 B 1
2 lc d 2 bc d
1 1
3 bc d 3 lc d
,
Diseño automatizado de cimentaciones superficiales con MATLAB
71
JL
d lc d d 2 lc d 2
lc d bc d
2 3 ,
JB
d bc d d 2 bc d 2
bc d lc d
2 3 ,
x
lc d y
bc d
2 , 2
Figura 2.17
V U
l V c l
V U
b V c b
dónde:
Ø=0.75,
f 'c
Vcl = 160 f ' c B d; Vcb = 160 L d,
d - peralto efectivo.
B L
spL 20 df spB 20 df
b' l'
B C L d L C B d
spL 20 df spB 20 df
b' l' l' b'
Vu l spL sL ; Vu b spB sB
El peralto efectivo que se usa para la comprobación inicial de cortante es el obtenido del
criterio del punzonamiento, si este peralto no cumple, se hacen aumentos de 0.01 metros al
mismo hasta que cumpla por este criterio. Figura 2.18
Diseño automatizado de cimentaciones superficiales con MATLAB
73
Figura 2.18
Fuente: Elaboración propia.
De este criterio sale una dimensión de peralto, hv.
El peralto resultante obligatoriamente cumple con los dos criterios anteriores, ahora se
compara con la variable hrigido, esta variable está definida en el programa como la mitad
del vuelo mayor, y es la frontera que usa la norma cubana para delimitar si el cimiento tiene
un comportamiento rígido o flexible. Para pedestal centrado el vuelo es:
l lc b bc
Cl , Cb ,
2 2
Figura 2.20
Fuente: (Chagoyén & Broche, 2002)
Si por el contrario resulta flexible es necesario encontrar el peralto necesario por flexión
positiva, peralto que le permite al hormigón tomar los esfuerzos a flexión.
Los cimiento masivos requieren que las tracciones originadas por las flexiones que se
producen en la base, sean resistidas por el hormigón solamente, es decir a cuenta de su
momento de fisuración de la sección sin agrietarse, sin fisurarse.
Figura 2.21
Fuente: (Chagoyén & Broche, 2002)
6M uL 6M uB
hL hB
B fr L fr
dónde:
Ø=0.9,
f r - Módulo de Ruptura del hormigón a tracción por flexión puede determinarse a partir de
la resistencia especificada a la compresión f' c , por la siguiente expresión:
f r = 0,62 f ' c
Diseño automatizado de cimentaciones superficiales con MATLAB
75
Y MuL y MuB se encuentran a través de las expresiones a partir de comparar también los
lados efectivos l´ y b´ con la posición de la sección crítica a flexión positiva :
Si l’ xS3L ML
x S23 L
B act
b
20 df
2 B
Si b’ CB MB
xS23 B l
L 20 df
2 L
Si b’< CB MB
L l b
2 b x S3 B 20 df x S23 B
2 L 2
De estos dos peraltos hL y hB, se queda con el mayor y a este le denomina, hmasivo.
Pero como el cimiento tiene comportamiento flexible, en algunos casos puede necesitar
aumentar su peralto por flexión negativa.
Dónde actúen cargas de grandes excentricidades con sobrepesos apreciables dentro del área
de cimentación o ubicados a profundidades considerables, pueden ocurrir flexiones
negativas en el vuelo descargado de la misma. Como todavía estamos en la solución de
cimiento masivo, la respuesta a esta posible necesidad será aumentar el peralto. El
programa calcula el peralto para que cumpla con este criterio, si resulta mayor que el
peralto considerado por los criterios anteriores da como respuesta este peralto, de lo
contrario la respuesta final será el resultado mayor de los criterios antes evaluados.
Posición se la sección critica a este tipo de flexión. Figura 2.22. La sección crítica a flexión
negativa S4, se tomará en el punto donde el esfuerzo cortante sea nulo, pero nunca más allá
de la sección crítica a flexión positiva S3.
Diseño automatizado de cimentaciones superficiales con MATLAB
76
Figura 2.22
Fuente: (Chagoyén & Broche, 2002)
Para la distribución uniforme dentro del área efectiva, el cálculo del momento estará
dependiendo de las excentricidades:
Si 2eL xS3 x x S3
Si 2eL< xS3 b
2 eL
x B
b
qd
B
M uneg L
B
2
b
q prom x B
2
x 2eL 2
Si 2eB xS3 x x S3
l
Si 2eB< xS3 2 eB
L
x
l
qprom
L
M uneg B
L
2
l
q prom x 2
L
x 2eB 2
MunegL MunegB
hL hB
0.1667 B fr 0.1667 L fr
dónde:
Ø=0.9,
fr - Módulo de Ruptura del hormigón a tracción por flexión puede determinarse a partir de
la resistencia especificada a la compresión f' c que ya se determinó anteriormente para la
flexión positiva.
El mayor peralto de entre los todos los diseños será reconocido dentro del programa como
ht, y esta será la respuesta 2. Figura 2.23
Figura 2.22
Fuente: Elaboración propia.
Para este diseño se calcula igualmente los esfuerzos, lo lados efectivos de la cimentación, y
el peralto mínimo, funciones y ecuaciones detalladas todas en el anexo 4.
Inicialmente se encontraran los peraltos necesarios para cumplir con los criterios de
punzonamiento y cortante, especificados en el epígrafe anterior.
Diseño automatizado de cimentaciones superficiales con MATLAB
78
Al peralto resultante de estos análisis se le compara igualmente con el peralto para que se
comporte rígido (hrigido), si resulta mayor que el rígido, se calculará el área de acero para
cimientos de hormigón armado con comportamiento rígido:
FLmáx FBmáx
A sL A sB
fy fy
, , siendo:
FLmáx , FBmáx : Fuerza máxima que tracciona al refuerzo en la dirección del lado L y B
respectivamente, provocada por las cargas externas actuantes y que se determina por las
siguientes expresiones:
CL N ' M L H L H c CB N ' M B H B H c
FLmáx
*
FBmáx
d 4 L , d 4 B
Hc distancia desde el nivel de solera hasta el punto de actúan las fuerzas cortantes.
N, ML, MB, HL, HB solicitaciones de cálculo, expresadas en kN, kN-m, kN-m, kN, y kN
respectivamente.
A esta respuesta se le llamara 3, con el peralto que resulto y con el refuerzo inferior
correspondiente a los cimientos rígidos de hormigón armado.
Diseño automatizado de cimentaciones superficiales con MATLAB
79
𝑀𝑢𝐿
(0.85 − √0.7225 − 1.7 ) 𝐵 × ℎ × 𝑓´𝑐
∅ × 𝐵 × ℎ2 𝑓´𝑐
𝐴𝑎𝑐𝑒𝑟𝑜𝑙 =
𝑓𝑦
Dónde:
h- peralto efectivo.
MuL- momento actuante que responde a las ecuaciones antes vistas para el caso de cimiento
masivo, para encontrarle el peralto por flexión positiva, solo que para el caso de hormigón
armado los esfuerzos de flexión lo soportan las armaduras de refuerzo.
El área de acero que resulte se compara con una cuantía mínima de acero, en caso de ser
maro se queda entonces con el área mínima.
𝐴𝑎𝑐𝑒𝑟𝑜
𝑁𝑜.𝑏𝑎𝑟𝑟𝑎𝑠 =
𝜋(𝑑𝑖𝑎𝑚) 2⁄
4
𝐵 − 2 𝑟𝑒𝑐𝑢𝑏. 𝑙𝑎𝑡𝑒𝑟𝑎𝑙
𝑒𝑠𝑝𝑎𝑐𝑖𝑚𝑖𝑒𝑛𝑡𝑜 =
𝑁𝑜.𝑏𝑎𝑟𝑟𝑎𝑠 − 1
Diseño automatizado de cimentaciones superficiales con MATLAB
80
Dónde:
Figura 2.23
Fuente: Elaboración propia.
Nótese que para redondear, se usa la función ceil, que redondea por exceso, y la misma
función, restringe la respuesta con los espaciados máximos y mínimos.
En estos cimientos con comportamiento flexible hay que considerar la flexión negativa
también, el programa calcula el peralto necesario para que cumpla por flexión negativa,
como ya se demostró en el capítulo anterior. Seguidamente compara ese resultado (h_fn),
con el peralto del cimiento flexible, si este criterio no decide, el programa da la respuesta 4.
Pero en el caso que si influya, se le pedirá al usuario que indique como quiere tratar esta
situación. Tiene dos opciones, aumentar el peralto a la dimensión que calculó y que es
mayor a la obtenida por los criterios de punzonamiento y de cortante o calcular un área de
refuerzo que soporte dichas flexiones.
Decantarse por la primera opción, sería la solución 5, donde solo se tiene refuerzo inferior
por flexión positiva y un peralto que soporte la flexión negativa.
Diseño automatizado de cimentaciones superficiales con MATLAB
81
La otra opción implica calcular un refuerzo superior con los momentos por flexión negativa
de manera similar a la forma de calcular el refuerzo por flexión positiva, solo cambie el
momento actuante, por los obtenidos por las ecuaciones explicadas en el epígrafe anterior
para flexión negativa, o sea M uneg L y M uneg B . Y esta será la respuesta 6 en nuestra
Figura 2.24
Fuente: Elaboración propia.
Antes de poder terminar con el diseño estructural, se hacen según la respuesta obtenida, el
chequeo de las longitudes de anclaje, la condición de adherencia por resistencia. Explicadas
ambas en el anexo 5.
Para la determinación de los costos, según la solución se calculan los volúmenes de trabajo.
En el caso más amplio esto comprenden, volumen de excavación, de hormigonado, de
rehincho y peso del acero a utilizar. Las funcione que calculan dichos volúmenes se
encuentran en el anexo 6.
Con dicho volúmenes se proceden a encontrarse los costos por materiales, mano de obra y
uso de equipos, determinando los renglones variantes del PRECONS II que se van a usar.
Estos quedan resumidos en la siguiente tabla:
Diseño automatizado de cimentaciones superficiales con MATLAB
82
12 Excavaciones
3053 Colocación y
ensamblaje de barras
rectas o dobladas
41 Hormigonado
calidad de hormigón.
22 Productos
semielaborados
14 Rehincho o rellenos
No se abarcan todas las variantes de problemas pues solo se considera un solo estrato y sin
presencia de nivel freático. Pero a manera general se logra el diseño geotécnico y
estructural, al encontrar las dimensiones definitiva B y L, para el primer y segundo estado
límite, considerando tanto comportamiento lineal como no lineal del suelo, además de
encontrar el peralto definitivo del plato y el refuerzo que necesite para que resista como
estructura.
Otra limitación del trabajo es que solo se considera el diseño del plato de la cimentación,
sin el pedestal. Además de que la respuesta del mismo no tiene múltiples formas, el usuario
tiene que especificar qué tipo de cimiento quiere que el software diseñe.
El suelo que se trabaja primeramente, será uno con C=35 kPa, y Fi=15 grados, a una
profundidad de 2 metros. Condiciones de trabajo normales y el fallo grave.
Las diferencias son despreciables entre los resultados obtenidos, por tanto quedan
comprobados los tras criterios del primer estado límite.
CM N’ = 1000 kN H= 0 M= 0
CVld N’ = 300 kN H= 0 M= 0
CVcd N’ = 300 kN H= 0 M= 0
= 40o
1 = 16 kn/m3
2 = 18 kn/m3
Resultados qbr (kN/m2 ) Rprima Pactuante Asentamiento Ha (metros) Carga Tensión peso b=l (metros)
(kN) (kN) (metros) impuesta propio *0.2
(kN/m2 ) (kN/m2 )
Mathcad 946.15466 394.8603 371.88034 0.03295 lineal 4.73 23.18314 23.697 1.95
suelo 1
Mathcad 322.13663 331.24595 161.01537 0.01609 lineal 4.46 20.73446 21.156 3.15
suelo 2
El suelo que trataremos será uno con C=35 kPa, y Fi=15 grados, a una profundidad de 1.5 metros.
Con 0.3 metros por ambos lados del pedestal, f´c=25 MPa, y fy=300 MPa.
N (kN) Hl Ml Hb Mb
(kN) (kN-m) (kN) (kN-m)
1200 23 12 21 21
1400 23 2 2 12
1500 2 2 1 12
1500 2 12 12 12
1875 23 12 21 2
1435 23 12 12 12
1245 21 12 21 12
La respuesta estructural que decida el usuario definirán los resultados, el DiGEC no obtiene más de
una respuesta. Tiene que correr el programa de nuevo para que tome otro camino, ya sea solución
masiva o de hormigón armado. Con el ejemplo anterior se realizó un diseño de hormigón armado y
la única incertidumbre fue el caso de los espaciamientos, aunque seguramente se usó la misma
forma de determinación en ambos casos, el manejo de las funciones para aproximar los resultados
no fue el mismo. De esta forma comprobamos el diseño estructural.
92
Diseño automatizado de cimentaciones superficiales con MATLAB
Conclusiones
Los costos directos no se incluyen en las comprobaciones finales, la variante de renglones
variantes que se usó, no está implementada en ningún otro programa.
Nos ajustamos a los casos donde la base es homogénea y sin presencia de nivel freático, que no es
la solución más genérica pero es un primer acercamiento a un programa en MATLAB que diseñe
geotécnicamente con la totalidad de las variantes del problema.
El diseño estructural no incluye el diseño del pedestal, que también es parte estructural del cimento
y es necesario para una solución total.
Diseñar una cimentación incluye muchos factores determinar como la profundidad del cimiento y
el tipo de base en que se va a apoyar, no sólo las dimensiones para que la base no falle, ni el
peralto y refuerzo para que como estructura el cimiento soporte y transmita las cargas de la
superestructura, que es lo que puede diseñarse matemáticamente.
Como concepto el programa DiGEC hecho en MATLAB, no representa nada nuevo en términos de
solución, los trabajos previos, de una forma u otra incluyen todos los aspectos del mismo, pero al
ser MATLAB un programa con las funciones para la optimización ya programadas, el DiGEC abre
las puertas a la opción de optimización de conjuntos estructurales.
Se logran resultados fiables para un número grande de combinaciones de carga, en realidad para
todas las que necesite el usuario, pero no identifica que combinación predomina como crítica sobre
las demás.
93
Diseño automatizado de cimentaciones superficiales con MATLAB
Recomendaciones
Ampliar los casos a problemas que incluyan suelos estratificados y con presencia del nivel freático
para tener así una solución más cercana a los casos reales que pueden aparecerle al ingeniero que
use el programa.
Completar el diseño estructural al incluir al pedestal, tratándose como una columna corta. La
complejidad radica en programarlo, porque una vez hecho basta con copiarlo al script principal, e
integrarlo al programa completo.
Recomendamos que se empleado los usuarios de MATLAB como otra forma rápida y segura de
obtener resultados del diseño geotécnico y estructural de cimentaciones rectangulares superficiales.
94
Diseño automatizado de cimentaciones superficiales con MATLAB
Bibliografía
Alfio Quarteroni, F. S. P. G., 2010. Scientific Computing with MATLAB and Octave. Berlín:
Springer.
Ballate Reina, D. R., 2016. Costo de los servicios de construcción y montaje. Ejemplo Resuelto..
Santa Clara, UCLV.
CEN Comité Europeo de Normalización, 2007. Eurocódigo EC7 Geotechnical design. Bruselas:
BSI.
Gonzáles Cueto, A. V. & Quevedo Sotolongo, G., 2007. Aplicación de la teoría de seguridad al
diseño de cimentaciones en arenas. Chequeo de linealidad, Santa Clara: UCLV.
Hernández Santana, J. J. & Hernández Caneiro, J. A., 2010. Capítulo 3 La seguridad en el diseño
estructural.. En: Hormigón Estructural. Diseño por estados límites.. La Habana: UCLV, pp. 202-
251.
Hernández, J. J., 2012. Capítulo 10: Columnas. Estudio de la flexión compuesta en secciones de
hormigón armado.. En: s.l.:s.n.
95
Diseño automatizado de cimentaciones superficiales con MATLAB
Quevedo Sotolongo, G., 1994. Diseño de cimentaciones superficiales. Manual del Proyectista.
Primera ed. Santa Clara: UCLV.
Quevedo Sotolongo, G., 2007. Propuesta de norma no publicada para el diseño geotécnico de
cimentaciones superficiales por estados límites.. UCLV ed. Santa Clara(Ciudad de La Habana):
Ministerio de la Construcción.
Quevedo Sotolongo, G., 2010. Conferencia 14 Met Generales de diseño. Santa Clara, UCLV, pp.
7-10.
Quevedo Sotolongo, G., 2010. Introducción al diseño estructural de cimientos. Santa Clara,
UCLV.
Quevedo, G. S., 1987. Optimización del proyecto de cimentaciones de edificios industriales en las
condiciones de Cuba.. MICI. Moscú: Instituto de la Ingeniería en Construcción de Moscú..
Rodriguez, J. M. O., Serra, J. G. & Oteo, C. M., 1989. Curso aplicado de cimentaciones, cuarta
edición. Cuarta Edición ed. Barcelona: GRAFICINCO,S.A.
Smith, F., 2012. Aplicación de Visual Basic en el diseño geotecnico de cimentaciones suprficiales.
Internacional Journal of Scientific & Engeneer Research volumen 3.
Wilson, W. &., 1999. MATLAB. Dr. Dobb’s Journal, 7 enero, pp. 50-62.
96
Diseño automatizado de cimentaciones superficiales con MATLAB
Anexos
Anexo 1: Funciones que determinan los coeficientes de minoración y las propiedades físico-
mecánicas del suelo.
Función para encontrar los coeficientes de minoración para una probabilidad de diseño del
95%:
(c 0, y 25º ).
Para > 30o ----- gtg = 1.1, gc = 1.4, g = 1.05
Función para encontrar los coeficientes de minoración para una probabilidad de diseño del
85%:
97
Diseño automatizado de cimentaciones superficiales con MATLAB
Para encontrar las propiedades físico-mecánicas del suelo minorada para el primer estado
límite:
Para encontrar las propiedades físico-mecánicas del suelo minorada para el segundo estado
límite:
98
Diseño automatizado de cimentaciones superficiales con MATLAB
c tan
* c* * tan 1
g gc
g tan
Anexo 2:
N = 2 (Nq - 1) tan *
Nc = (Nq - 1) cot *
Suelo y C -
s = 1 – 0.4 B´ / L´
Suelo C ( = 0)
Suelo y C –
5
0.7 H *
1
i =
N * b´ l´ c * cot *
99
Diseño automatizado de cimentaciones superficiales con MATLAB
5
0.5H *
1
iq =
N * b´ l´ c * cot *
(1 iq )
ic = iq - ( Nq 1)
Suelo C ( = 0)
H*
1
ic´ = 0.5 - 0.5 b ´l´ c *
d, dc, dq. - factores que valoran el efecto de la profundidad del cimiento dentro del estrato
resistente D.
Suelo y C -
D b
d = 1
dc = 1 + 0.4 (D/b)
Db
d =1
100
Diseño automatizado de cimentaciones superficiales con MATLAB
Suelo C ( = 0)
D b
D b
101
Diseño automatizado de cimentaciones superficiales con MATLAB
Anexo 3:
Caso Valor
12
12
102
Diseño automatizado de cimentaciones superficiales con MATLAB
10
40
103
Diseño automatizado de cimentaciones superficiales con MATLAB
d) H 300 S. medio 10
𝑁
𝜎𝑎𝑐𝑡 = 𝑏´×𝑙´ 𝑏´ = 𝐵 − 2𝑒𝑏 𝑙´ = 𝐿 − 2𝑒𝑙
El peralto efectivo mínimo para hacer el cimiento de hormigón armado: (para determinar
el peralto efectivo mínimo solo cambia hmin1 de 0.15m a 0.2m)
104
Diseño automatizado de cimentaciones superficiales con MATLAB
Condición de resistencia
VuL VuB
dl y dl
0.9 n L p L d L 0.9 n B p B d B , donde
dl 0.3 f´ c
- Barras lisas: (MPa)
dl 0.433 f´ c
2
- Barras corrugadas: (MPa)
y debiendo reducirse a las 3/4 partes (75 %) en las zonas de baja adherencia (zonas )
VuL
es la fuerza cortante en la sección crítica, que se determinará según la dirección
considerada, ya fue calculada para el cortante.
Longitud de anclaje:
105
Diseño automatizado de cimentaciones superficiales con MATLAB
Longitud de anclaje
106
Diseño automatizado de cimentaciones superficiales con MATLAB
En cimientos rígidos proporcionar una parte recta (ld - l d ), cuyo valor podrá tomarse
como el mayor de los valores siguientes:
1/3 ld
10
15 cm
En este caso de cimientos flexibles la armadura puede disponerse de lado a lado, sin
ganchos y sin subir por las caras laterales, siempre que los vuelos cumplan con:
ldL y ldB las longitudes de anclaje recto calculadas en las direcciones y para los diámetros
correspondientes, para zona de alta adherencia.
107