Trabajo de Titulacion Cargas Ciclicas

UNIVERSIDAD DE CUENCA
FACULTAD DE INGENIERÍA
ESCUELA DE INGENIERÍA CIVIL
ESTUDIO DEL COMPORTAMIENTO E INESTABILIDAD DEL

REFUERZO LONGITUDINAL EN ESTRUCTURAS DE CONCRETO
REFORZADO SOMETIDAS A CARGAS CÍCLICAS
TESIS PREVIA A LA OBTENCIÓN DEL

TÍTULO DE INGENIERO CIVIL
AUTORES:
ANDREA FERNANDA MONCAYO TELLO

C.I. 1400567572
MAURO ANDRÉS VALLEJO BOJORQUE

C.I. 0105736490
DIRECTOR:
INGENIERO HOMERO FERNANDO CARRIÓN CABRERA MsC.

C.I. 0104201629
CUENCA-ECUADOR
2017
Universidad de Cuenca
RESUMEN
Durante sismos, el desprendimiento del recubrimiento del concreto en estructuras de hor-

migón armado deja expuesto el acero del refuerzo longitudinal, pudiendo generar su falla por
inestabilidad. El presente trabajo tiene por objeto proponer una relación entre la deformación
crı́tica de pandeo y la esbeltez de barras del acero longitudinal de estructuras de hormigón
armado bajo cargas cı́clicas, generando información para el diseño sismorresistente eficiente.
Como punto de partida y en base al estudio de ciclos histeréticos y los modelos revisados, se
propone un nuevo modelo histerético del acero bajo cargas cı́clicas de compresión-tracción
antes de la falla por pandeo. En segunda instancia, se desarrolla un modelo que represente
las restricciones laterales y de giro proporcionadas por los estribos y el diámetro del refuer-
zo longitudinal respectivamente, por medio del cual se puede obtener el esfuerzo crı́tico de
pandeo que este último presenta, considerando la existencia de excentricidades y su efecto
causado. Con los modelos propuestos, se obtienen curvas de esfuerzo crı́tico de pandeo vs.
relación de esbeltez y deformación crı́tica de pandeo vs. relación de esbeltez, calibradas y
optimizadas con datos experimentales obtenidos de la revisión bibliográfica. Se compara las
excentricidades calculadas y los desplazamientos laterales que pudiera presentar el hormigón
en muros y columnas de la base de datos, para validar los resultados. Finalmente, se repre-
senta la curva final deformación crı́tica de pandeo vs. relación de esbeltez mediante funciones
matemáticas simples, para facilitar su uso.
Palabras Clave: Cargas cı́clicas, hormigón armado, refuerzo longitudinal, ciclos histeréticos,
pandeo, esbeltez, excentricidad, desplazamientos laterales.
Andrea Moncayo Tello 2

Mauro Vallejo Bojorque
ABSTRACT
During earthquakes, the detachment of the concrete covering in reinforced concrete struc-
tures exposes the steel of the longitudinal reinforcement, which can generate its failure by
instability. The present work aims to propose a relation between the critical strain of buc-
kling and the slenderness of longitudinal steel bars of reinforced concrete structures under
cyclic loads, generating information for efficient seismic resistant design. As a starting point
and based on the study of hysteretic cycles and the revised models, a new hysteretic model
of the steel is proposed under cyclic loads of compression-traction before the buckling fai-
lure. In the second instance, a model is developed that represents the lateral and turning
restrictions provided by the stirrups and the diameter of the longitudinal reinforcement,
respectively, through which the critical buckling stress can be obtained, considering the exis-
tence of eccentricities and its effect caused. With the proposed models, critical buckling stress
vs. slenderness ratio and critical buckling strain vs. slenderness ratio curves are obtained,
calibrated and optimized with experimental data obtained from the bibliographic review. We
compare the calculated eccentricities and the lateral displacements that could present the
concrete in walls and columns of the database, to validate the results. Finally, the buckling
critical strain vs. slenderness ratio curve is represented by simple mathematical functions,
to facilitate its use.
Keywords: Cyclic loads, reinforced concrete, longitudinal reinforcement, hysteretic cycles,

buckling, slenderness, eccentricity, lateral displacements.

Índice general
Resumen 2
Abstract 3
Índice de figuras 6
Índice de cuadros 8
Lista de sı́mbolos 9
1. Introducción 15
1.1. Antecedentes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
1.2. Justificación . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
1.3. Objetivos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
1.4. Metodologı́a . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
2. Estado del arte 18

2.1. Comportamiento del acero bajo carga de tensión . . . . . . . . . . . . . . . . 18
2.1.1. Zonas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
2.1.1.1. Zona elástica lineal . . . . . . . . . . . . . . . . . . . . . . . 19
2.1.1.2. Zona de fluencia . . . . . . . . . . . . . . . . . . . . . . . . 19
2.1.1.3. Zona de endurecimiento por deformación . . . . . . . . . . . 20
2.2. Modelos del Acero . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
2.2.1. Modelo de Mander . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
2.2.2. Modelo de Park y Paulay . . . . . . . . . . . . . . . . . . . . . . . . 21
2.2.3. Modelo de Ahmad y Shah . . . . . . . . . . . . . . . . . . . . . . . . 22
2.3. Comportamiento del acero bajo carga de compresión . . . . . . . . . . . . . 22
2.4. Comportamiento del acero bajo cargas cı́clicas . . . . . . . . . . . . . . . . 23
2.4.1. Introducción a los ciclos histeréticos . . . . . . . . . . . . . . . . . . . 23
2.4.2. Efecto Bauschinger . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
2.4.3. Modelos de curva esfuerzo vs. deformación bajo carga cı́clica . . . . . 26
2.4.3.1. Modelo de Chang y Mander . . . . . . . . . . . . . . . . . . 26
2.4.3.2. Modelo de Kent y Park . . . . . . . . . . . . . . . . . . . . 28
2.5. Pandeo inelástico . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
2.5.1. Teorı́a del Módulo Tangente . . . . . . . . . . . . . . . . . . . . . . . 31
2.5.2. Teorı́a del Doble Módulo Reducido . . . . . . . . . . . . . . . . . . . 32
2.5.3. Pandeo inelástico en elementos cargados excéntricamente . . . . . . . 35
2.6. Métodos de optimización . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36

3. Pandeo en barras de acero bajo cargas cı́clicas 38

3.1. Modelo propuesto del comportamiento del acero bajo cargas cı́clicas . . . . . 38
3.1.1. Curva esfuerzo vs. deformación bajo carga monotónica . . . . . . . . 38
3.1.2. Curva esfuerzo vs. deformación bajo carga cı́clica. . . . . . . . . . . . 40
3.1.3. Modelo propuesto . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41
3.2. Modelo de barra . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
3.3. Relación entre el pandeo y la esbeltez . . . . . . . . . . . . . . . . . . . . . . 54
3.3.1. Optimización y obtención del factor de la longitud efectiva . . . . . . 59
4. Pandeo en barras de acero del refuerzo longitudinal de estructuras de

hormigón armado bajo cargas cı́clicas 65
4.1. Optimización y obtención del factor de la longitud efectiva . . . . . . . . . . 69
4.2. Análisis de excentricidades y validación de las curvas: obtenida y propuesta . 73
5. Conclusiones y Recomendaciones 76
A. Codificación del programa Matlab del modelo propuesto del comporta-

miento del acero bajo cargas cı́clicas 81
B. Codificación del programa Matlab de las curvas esfuerzo crı́tico de pandeo

vs. relación de esbeltez y deformación crı́tica de pandeo vs. relación de
esbeltez 86
C. Codificación del programa Matlab de la obtención de K mediante optimi-

zación 89

Índice de figuras
2.1. Esfuerzo vs. deformación del acero bajo carga de tracción. . . . . . . . . . . 18

2.2. Punto experimental perteneciente a la zona de endurecimiento por deformación 21
2.3. Material inelástico. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
2.4. Efecto de carga y descarga, con inversión del sentido de carga [1]. . . . . . . 25
2.5. Ciclos de histéresis del acero de refuerzo. Obtenido de [2] . . . . . . . . . . . 25
2.6. Modelo de Chang y Mander. . . . . . . . . . . . . . . . . . . . . . . . . . . 28
2.7. Modelo de Kent y Park. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
2.8. comportamiento de la ecuación 2.31 con diferentes valores del parámetro r. . 30
2.9. Módulo Tangente . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
2.10. Comportamiento de la sección transversal del acero . . . . . . . . . . . . . . 32
2.11. Elemento cargado excéntricamente . . . . . . . . . . . . . . . . . . . . . . . 35
3.1. Curva esfuerzo vs. deformación bajo carga a tracción . . . . . . . . . . . . . 39

3.2. Curva esfuerzo vs. deformación bajo carga a compresión . . . . . . . . . . . 40
3.3. Zona elástica y zona de descarga . . . . . . . . . . . . . . . . . . . . . . . . . 42
3.4. Curva histerética en zona de tracción . . . . . . . . . . . . . . . . . . . . . . 43
3.5. Curva histerética en zona de compresión . . . . . . . . . . . . . . . . . . . . 45
3.6. Graficación de ciclos histeréticos. . . . . . . . . . . . . . . . . . . . . . . . . 47
3.7. Modelo propuesto con la curva de esbeltez 2.5 de Botero [3]. . . . . . . . . . 48
3.8. Modelo propuesto con la curva de esbeltez 4 de Botero [3]. . . . . . . . . . . 48
3.9. Modelo propuesto con la curva de esbeltez 6 de Botero [3]. . . . . . . . . . . 49
3.10. Modelo de barra . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
3.11. Esfuerzos en los extremos de la barra . . . . . . . . . . . . . . . . . . . . . . 50
3.12. Corte de longitud x en la barra . . . . . . . . . . . . . . . . . . . . . . . . . 50
3.13. Esfuerzo de pandeo vs. relación de esbeltez . . . . . . . . . . . . . . . . . . . 56
3.14. Deformación de pandeo crı́tico propuesto por Botero [3] . . . . . . . . . . . . 57
3.15. Deformación de pandeo vs. relación de esbeltez . . . . . . . . . . . . . . . . . 58
3.16. Función error vs. número de iteraciones . . . . . . . . . . . . . . . . . . . . . 61
3.18. Esfuerzo crı́tico de pandeo vs. esbeltez . . . . . . . . . . . . . . . . . . . . . 62
3.19. Deformación crı́tica de pandeo vs. esbeltez . . . . . . . . . . . . . . . . . . . 62
4.1. Gráfica deformación crı́tica vs. relación de esbeltez con datos experimentales
obtenidos en [4] y [5] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67
4.2. Gráfica relación de excentricidad vs. relación de esbeltez, con k=0.5 . . . . . 68
4.3. Gráfica relación de excentricidad vs. relación de esbeltez, con k=0.6 . . . . . 68
4.4. Gráfica relación de excentricidad vs. relación de esbeltez, con k=0.75 . . . . 69

4.6. Gráfica relación de excentricidad vs. relación de esbeltez, con k=0.632 . . . . 71

4.7. Gráfica deformación crı́tica de pandeo vs. relación de esbeltez . . . . . . . . 72
4.8. Relación entre los desplazamientos laterales del hormigón y las excentricidades 75
4.9. Curva obtenida y curva propuesta . . . . . . . . . . . . . . . . . . . . . . . . 75

Índice de cuadros
3.1. Parámetros de las probetas ensayadas por Botero [3]. . . . . . . . . . . . . . 39

3.2. Parámetros de barra ensayada por Botero con esbeltez 2.5 . . . . . . . . . . 41
3.3. Factor R. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41
3.4. Esfuerzos crı́ticos de pandeo de probetas ensayadas por Botero [3] . . . . . . 56
3.5. Deformaciones crı́ticas de pandeo de probetas ensayadas por Botero [3] . . . 58
3.6. Esfuerzos y deformaciones crı́ticos experimentales y calculados. . . . . . . . . 63
3.7. Error entre esfuerzos crı́ticas experimentales y calculadas . . . . . . . . . . . 63
3.8. Error entre deformaciones crı́ticas experimentales y calculadas . . . . . . . . 64
4.1. Resultados de ensayos experimentales en muros y columnas bajo cargas cı́cli-

cas obtenidos en [4] y [5]. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 66
4.2. Deformaciones crı́ticas experimentales y calculados. . . . . . . . . . . . . . . 70
4.3. Relación entre deformaciones crı́ticas experimentales y calculadas . . . . . . 70
4.4. Error entre deformaciones crı́ticas experimentales y calculadas . . . . . . . . 73
4.5. Excentricidades calculadas y desplazamientos laterales del hormigón . . . . . 74

Lista de sı́mbolos
En la siguiente lista se encuentra la simbologı́a utilizada en el presente estudio, en el orden

en el que son introducidas.
Sı́mbolos Generales
σ Esfuerzo.
ε Deformación.
Es Módulo Elástico.
fy Esfuerzo de fluencia.
Pcrit Fuerza crı́tica de pandeo.
I Inercia.
L Longitud.
R Radio.
A Área.
D Diámetro.
Z Módulo plástico.
f 0c Esfuerzo del hormigón a compresión.
Estado del Arte
εy Deformación del esfuerzo de fluencia.

εsh Deformación del inicio de la zona de endurecimiento por deforma-
ción.
εsu Deformación última.
fsu Esfuerzo último.
p Parámetro de Mander.
εs Deformación bajo carga de tracción.
fs Esfuerzo bajo carga de tracción.
εsc Deformación bajo carga de compresión.
fsc Esfuerzo bajo carga de compresión.
fch Tensión caracterı́stica de Kent y Park.
R Parámetro que define la curvatura de Mander.
r Parámetro que define la curvatura de Kent y Park.
εp1 Deformación plástica del ciclo anterior de Kent y Park.
Et Módulo tangente.
Er Módulo reducido.

Universidad de Cuenca - Cláusula de derechos de autor
Mmax Momento máximo resistente a flexión.

e Excentricidad.
Le Longitud efectiva.
r Radio de giro.
d Distancia al eje neutro de la sección transversal.
s Radio del núcleo de Timoshenko.
σmax Esfuerzo máximo.
X Espacio de búsqueda.
Ck Restricciones en el espacio de búsqueda.
n Número de parámetros.
f Función objetivo.
Pandeo en barras de acero bajo cargas cı́clicas
α0 Restricción lateral al inicio de la barra.

αL Restricción lateral al final de la barra.
β0 Restricción rotacional al inicio de la barra.
βL Restricción rotacional al final de la barra.
P Carga axial.
Q Cortante.
MA Momento al inicio de la barra.
MB Momento al final de la barra.
K Factor de longitud efectiva.
σcrit Esfuerzo crı́tico de pandeo.
De Diámetro del estribo.
Dl Diámetro del refuerzo longitudinal.
K1 Factor de reducción de la rigidez rotacional.
εp∗ Deformación crı́tica de pandeo.
ε0∗ Punto de cambio de carga de tensión a carga de compresión.
εp Deformación de pandeo.
εexpi Deformaciones experimentales.
εe i Deformaciones calculadas.
Pandeo en barras de acero del refuerzo longitudinal de estructuras

de hormigón armado bajo carcas cı́clicas
K Factor de longitud efectiva optimizado.

µ Módulo de Poisson.
δx Desplazamientos laterales.



Universidad de Cuenca - Cláusula de propiedad intelectual

Universidad de Cuenca - Cláusula de propiedad intelectual

Capı́tulo 1
INTRODUCCIÓN
1.1. Antecedentes
A lo largo de los años, la ingenierı́a civil ha puesto mayor énfasis en el estudio de estructuras
bajo cargas dinámicas como el viento, los sismos, maquinaria rotatoria, olas marinas, entre
otras. Por tanto es preciso desde el enfoque estructural, velar por la integridad de las estruc-
turas bajo dichas cargas. Se ha demostrado que en situaciones donde se presentan cargas
cı́clicas, como en sismos, un daño frecuente es la desintegración del hormigón que trabaja
como recubrimiento, por lo que los refuerzos longitudinales de acero quedan expuestos y da
como resultado el pandeo de los mismos, una de los modos de falla más comunes.
La falla por pandeo se ha presentado en varios terremotos como el de México en 1985,

Northridge en 1994, Kobe en 1995, Manzanillo en 1995, Chile en 2010, por mencionar al-
gunos, en los cuales se observa evidencias de inestabilidad en el refuerzo longitudinal[3].
Varios reglamentos como el del ACI 318-14 de Estados Unidos o el NSR-10 de Colombia,
consideran este problema y establecen separaciones máximas entre estribos de manera de
retrasar el pandeo del acero de refuerzo longitudinal; no obstante, las separaciones máximas
establecidas se basan en contados estudios del problema de pandeo del acero [3].
Autores como Mander [6], Salmon y Jhonson [7], Monti y Nuti [8] entre otros, han estu-
diado el pandeo de barras de acero de refuerzo longitudinal bajo cargas monotónicas. No
obstante, existen pocos estudios del pandeo del acero bajo cargas cı́clicas. Botero [3] estu-
dió el pandeo en barras de acero bajo cargas cı́clicas; sin embargo, la curva esfuerzo vs.
deformación utilizada era bajo cargas monotónicas.
1.2. Justificación
Ecuador se sitúa en una zona de alta sismicidad y se encuentra geológicamente atravesado
por varias fallas. Durante un evento sı́smico, debido a cargas cı́clicas, las barras de acero
pertenecientes al refuerzo longitudinal en estructuras de hormigón armado son sometidas a
cargas de tensión y compresión, esto, junto con la exposición causada por el desprendimiento
del recubrimiento o las fisuras que se ocasionan en el concreto por la tracción del acero en
el borde de muros y columnas, provocan su inestabilidad y susceptibilidad al pandeo [9];
siendo de vital importancia que en el diseño de estructuras de hormigón armado se consi-
dere el comportamiento del refuerzo de acero longitudinal bajo cargas sı́smicas con mayor

aproximación a la realidad.
1.3. Objetivos
Objetivo General
Proponer una relación entre la deformación crı́tica de pandeo y la esbeltez de barras de

acero pertenecientes al refuerzo longitudinal de estructuras de hormigón armado bajo cargas
cı́clicas, misma que pudiera ser fuente de información y ayuda para el diseño sismorresistente
eficiente.
Objetivos Especı́ficos
1. Proponer un modelo que represente el comportamiento del acero bajo cargas cı́clicas
antes de producirse la falla por pandeo, en base a modelos de diversos autores.
2. Proponer un modelo que represente las restricciones laterales y de giro proporcionadas por
los estribos y el diámetro del refuerzo longitudinal respectivamente, que permita calcular el
esfuerzo crı́tico de pandeo del acero longitudinal.
3. Obtener una curva de esfuerzo crı́tico de pandeo vs. relación de esbeltez L/D y una curva
de deformación crı́tica de pandeo vs. relación de esbeltez L/D del acero bajo cargas cı́cli-
cas, optimizarla y validarla con resultados de ensayos obtenidos de la bibliografı́a consultada.
4. Obtener una curva de deformación crı́tica de pandeo vs. relación de esbeltez L/D del
acero como refuerzo longitudinal en estructuras de hormigón armado sometidas a cargas
cı́clicas, considerando la presencia de excentricidades; optimizarla y validarla con resultados
de ensayos obtenidos de la bibliografı́a consultada.
5. Representar la última curva con funciones matemáticas simples, para mejorar el manejo
y la utilidad del resultado del presente trabajo de titulación.
1.4. Metodologı́a
Se parte de una revisión literaria sobre el pandeo del acero bajo cargas cı́clicas y modelos
histeréticos del material. En base a modelos del acero bajo cargas monotónicas propuestos
por Mander [10], Park y Paulay [11], y Ahmad y Shah [12]; y el modelo del acero bajo cargas
cı́clicas de Mander[6], se propone un modelo de curva esfuerzo vs. deformación del acero ante
cargas cı́clicas que depende solamente de las caracterı́sticas del acero y la deformación máxi-
ma de cada ciclo. Se considera que las restricciones laterales y de giro del modelo de estudio
son proporcionadas por los estribos y el diámetro del refuerzo longitudinal respectivamente.
Implementado el modelo propuesto en Matlab, se obtienen curvas esfuerzo crı́tico de pandeo
vs. relación de esbeltez y deformación crı́tica de pandeo vs. relación de esbeltez. Los paráme-
tros desconocidos del modelo se obtienen al ajustar las curvas con datos experimentales.

Posteriormente, se realiza un análisis de la existencia de excentricidades y su efecto en ensa-

yos de muros y columnas bajo cargas cı́clicas. Se obtiene una curva final deformación crı́tica
de pandeo vs. relación de esbeltez que considera el efecto de las excentricidades presentes,
de manera que se aproxime adecuadamente al comportamiento real de las barras de acero
pertenecientes al refuerzo longitudinal en estructuras de hormigón armado. Se compara las
excentricidades calculadas y los desplazamientos laterales que podrı́a presentar el hormigón
en los ensayos de muros y columnas, y se valida la última curva obtenida. Finalmente, el
producto de la investigación es la curva final ajustada mediante funciones matemáticas sim-
ples, con el objetivo de mejorar el manejo y la utilidad de los resultados obtenidos.

Capı́tulo 2
ESTADO DEL ARTE
2.1. Comportamiento del acero bajo carga de tensión

La curva esfuerzo vs. deformación a tracción del acero suele considerarse como elastoplástica
en el análisis sı́smico de estructuras de hormigón armado. Es decir, que el material tiene un
comportamiento elástico hasta llegar a un esfuerzo igual a fy , al cual le corresponde una
deformación εy ; mientras que para deformaciones mayores a εy , el material posee esfuerzos
iguales a fy . Esta curva puede representar adecuadamente el comportamiento de aceros
de baja resistencia. No obstante, Bonett [13] señala: “para el caso de diseño sı́smico, los
requerimientos de ductilidad, pueden implicar evaluar el esfuerzo del acero a deformaciones
mayores que la de fluencia”.
Esfuerzo, fs
1. Zona elástica lineal.
2. Zona de fluencia.
3. Zona de endurecimiento por deformación.
fsu
3
2
fy
y sh su
Deformación, s
Figura 2.1: Esfuerzo vs. deformación del acero bajo carga de tracción.

2.1.1. Zonas
En la literatura revisada se encontraron varios modelos que dividen a la curva esfuerzo vs.
deformación en tres zonas:
1. Zona elástica lineal.
2. Zona de fluencia.
3. Zona de endurecimiento por deformación.
2.1.1.1. Zona elástica lineal

La zona elástica lineal pertenece a deformaciones comprendidas entre 0 y εy , se caracteriza
por cumplir la ley de Hooke, la relación entre esfuerzo y deformación es constante e igual
al módulo de Young (Es ). Se representa por una lı́nea recta que va desde el origen hasta el
punto de fluencia del acero.
El valor del esfuerzo en esta zona se puede obtener utilizando la siguiente expresión [3]:
fs = Es · εs (2.1)
dfs
Es = (2.2)
dεs
Donde Et , fs y εs son el módulo elástico, el esfuerzo y la deformación en cualquier pun-

to de la curva dentro de la zona elástica lineal.
En esta zona el proceso es reversible, pues, si se retira la fuerza aplicada, el acero vuel-
ve a su posición inicial siguiendo el mismo recorrido de carga.
2.1.1.2. Zona de fluencia

La zona de fluencia pertenece a deformaciones comprendidas entre εy y la deformación del
inicio de la zona de endurecimiento por deformación (εsh ). En esta zona las distorsiones de la
estructura cristalográfica se vuelven permanentes y el material se comporta plásticamente,
por lo tanto, no se requiere un aumento de carga para que los cristales que forman la estruc-
tura del material se puedan mover libremente [3]. Esto corresponde a una reacomodación de
los cristales dentro de su estructura.
Esta zona no es totalmente plana, pero debido a sus pequeñas variaciones, su forma puede
suponerse de tal manera [14].
El valor del esfuerzo en esta zona se obtiene utilizando la siguiente expresión:
fs = fy (2.3)
Es ≈ 0 (2.4)

2.1.1.3. Zona de endurecimiento por deformación

La zona de “endurecimiento por deformación”pertenece a deformaciones comprendidas entre
εsh y la deformación correspondiente al esfuerzo último fsu (εsu ). En esta zona los cristales
que forman la estructura interna del material sufren grandes distorsiones con el aumento de
carga, pero a diferencia de la zona elástica, el material no se comporta en forma lineal sino,
debido al deslizamiento de los cristales, el material va perdiendo rigidez [3].
Los autores de los modelos revisados, discrepan entre ellos en la zona de endurecimiento
por deformación. Cabe aclarar que después del esfuerzo último (fsu ), la curva esfuerzo vs.
deformación tiene una zona descendente, sin embargo, ninguno de los modelos consultados
considera de importancia esta zona, ya que se pretende evaluar la capacidad máxima del
acero, la cual se encuentra asociada al esfuerzo último, siempre que el modo de falla se pre-
sente por este efecto [15].
A continuación, se presentan los distintos modelos revisados:
2.2. Modelos del Acero

2.2.1. Modelo de Mander
Según el modelo de Mander [10], los esfuerzos dentro de la zona de endurecimiento por de-
formación se pueden obtener con la siguiente expresión:
p
εsu − εs
fs = fsu + (fy − fsu ) · (2.5)
εsu − εsh
Donde p es el parámetro que permite definir la forma de la curva esfuerzo vs. deformación
dentro de la zona de endurecimiento por deformación.
Para la obtención de p mediante la ecuación 2.5, se puede utilizar la definición de loga-

ritmo natural y haciendo que fs = fs1 y εs = εs1 , donde (εs1 , fs1 ) es un punto experimental
perteneciente a la zona de endurecimiento por deformación (Figura 2.2), p queda definido
por la siguiente expresión:

fsu − fs1
ln
fsu − fy
p= (2.6)
εsu − εs1
ln
εsu − εsh
Mediante un análisis de las expresiones necesarias para obtener la curva esfuerzo vs. de-
formación, se puede observar en la ecuación 2.5 que los parámetros para obtener fs entre
0 < εs < εsu son: fy , fsu , εsh , εsu , Es y p [16]. Derivando dicha ecuación con respecto a εs y
haciendo que εs = εsh , se obtiene:
dfs p · (fsu − fy ) · 1
Esh = = (2.7)
dεs εsu − εsh

Esfuerzo, fs
fsu
fs1
fy
sh s1 su
Deformación, s
Figura 2.2: Punto experimental perteneciente a la zona de endurecimiento por deformación
Por tanto:

εsu − εsh
p = Esh · (2.8)
fsu − fy
Se ha encontrado que para evaluar el valor de p, la ecuación 2.8 tiene buena relación entre
los resultados analı́ticos y experimentales [10].
Una desventaja del uso de Esh como parámetro de referencia radica en que pequeñas varia-
ciones en los valores experimentales pueden conducir a cambios significativos en el valor de
p.
2.2.2. Modelo de Park y Paulay

Los valores de fs dentro de la zona de endurecimiento por deformación según el modelo de
Park y Paulay [11] se obtienen mediante la siguiente expresión:

m · u + 2 (60 − m) + u
fs = fy · + (2.9)
60 · u + 2 3(30 · r + 1)2
Los parámetros m, r y u se pueden obtener mediante pruebas experimentales, o mediante

las siguientes expresiones:

fsu
(30r + 1)2 − 60r − 1
fy
m= (2.10)
15r2
r = εsu − εsh (2.11)

u = εs − εsh (2.12)
Donde fy es el esfuerzo de fluencia nominal en el acero, fsu es el esfuerzo último, εs la defor-

mación unitaria del acero, εsu la deformación unitaria última y εsh la deformación unitaria
del acero en la cual empieza la zona de endurecimiento por deformación [11].
2.2.3. Modelo de Ahmad y Shah

El modelo propuesto por Ahmad y Shah [12], obtiene el fs dentro de la curva de endureci-
miento por deformación mediante las siguientes expresiones:
fs = fy + Y · (fsu + fy ) (2.13)
A · X − (B − 1) · X 2
Y = (2.14)
A + (A − 2) · X + B · X 2
εs − εsh
X= (2.15)
εsu − εsh
εsh = 0.0145 − 1.2801e−6 · f y (2.16)
Donde las constantes A y B tienen los siguientes valores obtenidos experimentalmente [12]:
A = 1.735
B = 3.62
2.3. Comportamiento del acero bajo carga de compre-

sión
Los modelos descritos anteriormente, hacen referencia al comportamiento del acero bajo car-
gas de tracción. Sin embargo, Dodd y Restrepo [14] establecieron un procedimiento para
determinar la curva esfuerzo vs. deformación del acero bajo carga de compresión a partir de

una curva esfuerzo vs. deformación bajo carga de tracción.
Este procedimiento pretende formar una relación entre deformaciones y esfuerzos de ba-
rras de acero bajo carga de tensión y bajo carga de compresión, mediante la teorı́a de que
no existe perdida de volumen de la barra de acero ensayada, lo que evidencia que las di-
mensiones transversales y longitudinales de la barra en estudio sólo varı́an según el efecto
de Poisson [14].
Las expresiones para obtener fsc y εsc que hacen referencia al esfuerzo y deformación bajo
carga de compresión en función de la curva esfuerzo vs. deformación bajo carga de tracción
son las siguientes [14]:
fsc = −fs · (1 + εs )2 (2.17)
−εs
εsc = (2.18)
1 + εs
Combinando las ecuaciones se obtiene el parámetro p que define la forma de la zona de

endurecimiento por deformación del modelo de Mander [16], para compresión:

fsuc − fsc
ln
fsuc − fyc
pc = (2.19)
εsuc − εsc
ln
εsuc − εshc
2.4. Comportamiento del acero bajo cargas cı́clicas

2.4.1. Introducción a los ciclos histeréticos
Como se ha visto anteriormente, el comportamiento elástico de un material permite que el
mismo regrese a su estado original luego de retirar una fuerza impuesta, por lo que no existe
una deformación permanente producto de haber aplicado y retirado dicha fuerza. Este pro-
ceso es el que ocurre con el acero en su parte elástica, es decir hasta que su deformación es
igual a εy .
Sin embargo, cuando la deformación del acero es mayor a εy , éste se encuentra en su zona
inelástica, pues cuando se quita una fuerza aplicada sobre el material, la descarga no ocurre
siguiendo la misma trayectoria de carga, y se presenta una deformación permanente. Como
se puede apreciar en la figura 2.3.
De acuerdo a Popov [17]; la figura 2.4 muestra un material que es cargado desde 0 hasta el
punto A; posteriormente es descargado, y toma la trayectoria AB. Tanto en la parte inicial
de la carga como en la descarga, el material tiene una respuesta esencialmente elástica con el
módulo de elasticidad inicial del material. Debido a que incursiona el rango inelástico antes

Esfuerzo, fs
Carga
Descarga
Deformación, s
Deformación
permanente
Figura 2.3: Material inelástico.
de llegar al punto A, presenta una deformación permanente; parte de la energı́a de defor-

mación necesaria para alcanzar el punto A es acumulada como energı́a potencial y se libera
en la descarga. En la Figura (2.4) puede observarse que al cargar nuevamente el material
a partir del punto B; se comporta como un material elástico, hasta que encuentra la curva
original en el punto C; y que al continuar la carga sigue la curva original. Si en el punto D se
retira nuevamente la carga, el material llega al punto E de esfuerzo aplicado nulo. Cargándo-
lo nuevamente en el sentido contrario a partir del punto E, desde F el material comienza a
comportarse inelásticamente.
Cuando el material experimenta una serie de ciclos de carga, descarga, y carga en el sentido
opuesto; donde los esfuerzos sobrepasan el lı́mite elástico del material, se observa este tipo
de comportamiento conocido como respuesta histerética. La Figura 2.5 muestra el resultado
de dos tipos de ensayos de carga cı́clica de acero de refuerzo. La gráfica (2.4.a) solo presenta
deformaciones positivas, mientras que en la figura (2.4.b), las deformaciones son positivas y
negativas, produciendo ciclos de histéresis simétricos.
El caso uno se presenta en el refuerzo de las vigas en zonas de rótulas plásticas que ex-
perimentan grandes deformaciones inelásticas a flexión y poca carga axial; mientras que el
segundo caso se puede producir durante la respuesta de una columna a carga cı́clica con
niveles de carga axial moderada o alta [2].
Para ambos casos, se ha observado que la curva obtenida bajo cargas monótonicas resul-
ta una envolvente de las curvas bajo cargas cı́clicas [2].

Esfuerzo, fs
D
A
C
Carga
Recarga
fa Descarga Descarga
B E Deformación, s
Deformación
permanente fb
Figura 2.4: Efecto de carga y descarga, con inversión del sentido de carga [1].
Figura 2.5: Ciclos de histéresis del acero de refuerzo. Obtenido de [2]
2.4.2. Efecto Bauschinger

En la figura (2.4), fa y fb representan las magnitudes del esfuerzo fs para las cuales el ma-
terial deja de comportarse elásticamente y pasa a tener un comportamiento plástico, en la
carga inicial y en la carga realizada en sentido contrario respectivamente. Se observa que
fb es menor a fa, lo que implica que el comportamiento no lineal del material inicia con

una magnitud del esfuerzo menor al cargar el material en sentido contrario. Esto se conoce
como efecto de Bauschinger, quien fue el primero en observarlo, y por esto lleva su nombre [1].
Cuando el acero de refuerzo se somete a cargas cı́clicas, dentro del rango inelástico. La
plataforma de fluencia desaparece y la curva esfuerzo-deformación exhibe el efecto “Baus-
chinger”, según el cual, el comportamiento no lineal comienza para un nivel de esfuerzo
mucho más bajo que la resistencia inicial de fluencia. En consecuencia, las curvas esfuerzo
vs. deformación exhiben una forma redondeada durante la fluencia [13].
2.4.3. Modelos de curva esfuerzo vs. deformación bajo carga cı́cli-

ca
2.4.3.1. Modelo de Chang y Mander
El siguiente modelo propuesto por Chang y Mander [6], se basa en etapas de carga y des-
carga, es decir, en fases de tracción y compresión del acero. Mantiene como envolventes las
curvas del comportamiento monotónico del acero uniendo sus puntos finales con una curva
A, basada en las ecuaciones de Menegotto-Pinto [18], definida para cada una de las etapas,
como se muestra en la figura 2.6. La curva A considera el “Efecto Bauschinger” mencionado
en la sección 2.4.2, puesto que presenta transiciones suaves en la capacidad del acero desde
el inicio hasta el final de la curva, y presenta la siguiente formulación donde f0 y ε0 son la
tensión y deformación al inicio de la curva; ff y εf son la tensión y deformación al final de
la curva, R es el parámetro que representa el efecto Bauschinger [19]:
( )
1−Q
fs = f0 + E0 · (εS − ε0 ) Q + (2.20)
{1 + ρ}1/R
Donde:
R
εS − ε0
ρ = E0 · (2.21)
fS − f0
Esec
−a
E0
Q= (2.22)
1−a
Con:
ff − f0
Esec = (2.23)
εf − ε0
" R #−1/R
E0
a= 1+ (2.24)
Esec

Los parámetros E0 y R fueron calibrados mediante ensayos de fatiga de bajo ciclo reali-
zados por Mander [6] para las fases de carga y descarga respectivamente:
Fase de descarga:
E0 = Es · (1 − ∆ε) (2.25)
13
fy
R = 14· · (1 − 14∆ε) (2.26)
Es
Fase de carga:
E0 = Es · (1 − ∆ε) (2.27)
13
fy
R = 20 · · (1 − 18∆ε) (2.28)
Es
Donde:
εf + ε0
∆ε = (2.29)
2
R>0 (2.30)

Esfuerzo,fs
( f ,f f )
( 0,f0)
( 0,f0) Deformación, s
( f ,ff )
Figura 2.6: Modelo de Chang y Mander.
2.4.3.2. Modelo de Kent y Park

El modelo propuesto por Kent y Park [20] se puede observar en la figura 2.7, idealiza la curva
que sigue el comportamiento del acero bajo cargas cı́clicas, utiliza las curvas del compor-
tamiento monotónico como envolventes y basa el comportamiento cı́clico en la formulación
propuesta por Ramberg-Osgood [21], donde f0 y ε0 son la tensión y deformación al inicio
de la curva, ff y εf son la tensión y deformación al final de la curva, fch es la tensión
caracterı́stica y r es un parámetro que define la curvatura y depende del número de ciclos
incursionados anteriormente:
( )
fs − f0 r−1

(fs − f0 )
εs − ε0 = 1 + (2.31)
Es fch − f0
La figura 2.8 muestra el comportamiento de la ecuación 2.31 con diferentes valores del
parámetro r.
Con el propósito de linealizar la ecuación anterior, Kent y Park realizaron un análisis de

regresiones tomando los logaritmos en ambos lados de la ecuación y comparando con prue-
bas experimentales, de esta manera propusieron las siguientes expresiones empı́ricas para
encontrar los valores de la tensión caracterı́stica y el parámetro r de Ramberg-Osgood:
fch = fy · (0.973 − 9.806εp1 ) (2.32)

45.07 9.77
r = 12.33 + − (2.33)
exp(1000εp1 ) loge (1000εp1 + 2)
Donde fy es la tensión de fluencia del material y εp1 es la deformación plástica del ciclo
anterior medida experimentalmente.
Esfuerzo,fs
Envolvente de la ( 0,f0) para

curva esfuerzo vs carga de
deformación para compresión
carga de tensión.
Deformación, s
TanѲ=Es p1
Envolvente de la
curva esfuerzo vs
deformación para
carga de
comprensión.
( 0,f0) para
carga de
tracción
Figura 2.7: Modelo de Kent y Park.

r=1
r=5
r=20
1 r=∞
0.5
TanѲ=Es
0 1 2 3 4
Figura 2.8: comportamiento de la ecuación 2.31 con diferentes valores del parámetro r.
2.5. Pandeo inelástico

Actualmente, el problema del pandeo inelástico presenta considerable importancia debido a
que se han implementado diferentes usos del acero y sus aleaciones de alta resistencia en la
mayorı́a de miembros estructurales pertenecientes a casas, edificios, puentes, vı́as, etc. Por
esta razón, se ha presentado la necesidad de realizar varios estudios e investigaciones sobre
el comportamiento de este material y su utilidad.
Uno de los primeros problemas estudiados fue sobre inestabilidad elástica en relación al
pandeo lateral de miembros a compresión realizado por Leonhard Euler, quien presentó una
solución teórica enfocada en barras doblemente apoyadas, en la cual la carga crı́tica de pan-
deo es representada por la siguiente expresión [22]:
π 2 EI
Pcrit = (2.34)
L2
Donde E es el módulo de elasticidad del material, I la inercia de la sección transversal y L

la longitud de la barra.
Sin embargo, una columna o barra con una esbeltez moderada se pandea cuando parte
de sus fibras han alcanzado la fluencia [23], estudios analı́ticos y ensayos a compresión ex-
perimentales han demostrado que esta solución solo es válida cuando el pandeo se presenta
en la zona elástica del material debido a su comportamiento lineal pero es errónea su con-
sideración después de haber alcanzado el lı́mite de proporcionalidad puesto que el módulo

de elasticidad deja de ser constante, aspecto mencionado por Lamarle en 1845 y debido al
cual se plantean teorı́as sobre variaciones que se deben realizar para utilizar la ecuación de
Euler, en el rango inelástico.
2.5.1. Teorı́a del Módulo Tangente

Engesser propuso una teorı́a para el pandeo inelástico que supone el incremento de la des-
viación lateral simultáneamente con el incremento de la carga axial. En tal caso, la tensión
en el lado convexo de la barra disminuye en las etapas iniciales de flexión pero puede ser
compensada por el aumento de esfuerzo de compresión directa debido a que la fuerza axial
crece continuamente [22].
Para encontrar la carga crı́tica de pandeo tomando en cuenta esta teorı́a, basta con rempla-
σj
σi
εi εj ε
Figura 2.9: Módulo Tangente
zar el módulo de elasticidad Es por el módulo tangente Et en la ecuación de Euler:
π 2 EI
Pcrit = (2.35)
L2

2.5.2. Teorı́a del Doble Módulo Reducido

Considere indicó en 1889, que en una columna sometida a carga de compresión y que ha
sobrepasado el lı́mite de proporcionalidad, los esfuerzos aumentan en el lado cóncavo al inicio
de la flexión regidos por el módulo tangente Et mientras que disminuyen en el lado convexo
siguiendo una pendiente aproximadamente igual al módulo elástico del material Es . Por lo
tanto, la resistencia máxima de la columna no depende del módulo tangente ni del módulo
elástico, sino de un módulo reducido Er comprendido entre estos dos [3].
d1
A1,Et
X
A2,Es
d2
Disminución de deformación
debido a la flexión
Aumento de deformación
debido a la flexión
Figura 2.10: Comportamiento de la sección transversal del acero
En 1985, Engesser consideró los comentarios realizados por Considere y planteo la “Teorı́a
del Doble Módulo Reducido” la cual supone lo siguiente [3]:
Los desplazamientos laterales de los ejes de la columna son pequeños comparados con
las dimensiones de su sección transversal
Las secciones transversales planas permanecen planas y normales al eje deformado

después de la flexión.
El comportamiento de cualquiera de las fibras longitudinales del material se representa

con la curva esfuerzo vs. deformación de este.
El plano de flexión es un plano de simetrı́a de todas las secciones transversales.
A continuación se presenta la demostración de la teorı́a del doble módulo reducido, lo que

permite evaluar el esfuerzo de pandeo una barra de refuerzo bajo carga de compresión [24].
En la figura (2.10) se puede apreciar que el lado cóncavo de la sección presenta un módulo

tangente Et , mientras que el lado convexo presenta un módulo de elasticidad E; debiendo

encontrar un módulo equivalente, es decir, el doble módulo reducido Er . Para obtener el
doble módulo reducido Er , se realiza un promedio de pesos de las áreas 1 y 2 , obteniéndose
[7]:
Er I = Et I1x + Es I2x (2.36)

Donde:
Er es el módulo reducido.
Et es el módulo tangente.
Es es el módulo elástico.
I es el momento de inercia de la sección circular de la barra.
I1x es el momento de inercia del sector circular 1 respecto al eje X.
I2x es el momento de inercia del sector circular 2 respecto al eje X.
Despejando Er de la ecuación 2.36 se obtiene:

Et I1x Es I2x
Er = + (2.37)
I I
El momento de inercia I de una sección circular es:
πR4
I= (2.38)
4
Reemplazando el momento de inercia I en la ecuación 2.37, se obtiene:

4
Er = · [Et I1x + Es I2x ] (2.39)
πR4
Por definición de momento de inercia, de la figura 2.10 se obtiene:

Z d1
I1x = z1 2 dA1 (2.40)
0
Z d2
I2x = z2 2 dA2 (2.41)
0
Sustituyendo las dos expresiones anteriores en la ecuación 2.39:

Z d1 Z d2
4 2 2
Er = · Et z1 dA1 + Es z2 dA2 (2.42)
πR4 0 0

Con el propósito de simplificar, se hacen los siguientes cambios de variables [3]:

Z d1
1
X= 4· z1 2 dA1 (2.43)
R 0
Z d1
1
Y = 4· z2 2 dA2 (2.44)
R 0
Al sustituir las dos expresiones anteriores en la ecuación 2.42, se obtiene:

Er 4 Et
= · X· +Y (2.45)
Es π Es
Como se puede apreciar en la expresión 2.45, el módulo reducido depende de la geometrı́a

de la sección, del módulo elástico E y del módulo tangente Et .
Para obtener el módulo tangente Et , se calcula el primer momento de área respecto al eje
neutro X:
A1 Y1 Et = A2 Y2 Es (2.46)
Donde:
A1 y A2 son las áreas del sector circular 1 y 2 respectivamente, definidas en la figura

2.10.
Y1 y Y2 son los centroides del sector circular 1 y 2 respectivamente.

2 seno2α
A1 = R · α − (2.47)
2

2 seno2α
A2 = R · (π − α) + (2.48)
2
 
 2 seno(α)3 
Y1 = R ·  − cosα  (2.49)
 seno2α 
3 α−
2
 
 
2 · seno(α)3

 

Y2 = R · + cosα (2.50)
seno2α
 3 · (Π − α) +

 


2

Sustituyendo las cuatro expresiones anteriores en la ecuación 2.46 y despejado Es /Et se ob-
tiene:
seno2α
h 2(senoα)3 i
Es α − 2 3(α− seno2α )
− cosα
2
= h 3
i (2.51)
Et (π − α) + seno2α
2(senoα)
seno2α + cosα
2 3((π−α)+ 2
)
La expresión anterior permite obtener el módulo tangente Et necesario para evaluar el módu-
lo reducido equivalente Er .
De manera similar a la “Teorı́a del módulo tangente”, basta con remplazar el módulo de
elasticidad del material Es por el módulo reducido Er en la ecuación del pandeo de Euler:
π 2 Er I
Pcrit = (2.52)
L2
2.5.3. Pandeo inelástico en elementos cargados excéntricamente

El pandeo inelástico ha sido también estudiado tomando en cuenta las imperfecciones que se
presentan en el elemento o estructura, tales como excentricidades debidas a otros elementos
parte de la misma estructura, ası́ como excentricidades involuntarias causadas por la fabri-
cación, armado o ejecución. Tal análisis se basa en la falla por pandeo producida una vez
que la estructura alcanza el máximo momento resistente a flexión, anulando la capacidad de
seguir recibiendo carga.
Para un mejor enfoque se considera un elemento cargado excéntricamente como se mues-

tra en la figura 2.11.
P e e P
x
w max
Figura 2.11: Elemento cargado excéntricamente
Los momentos producidos en los extremos son:
M0 = P · e (2.53)

El máximo momento resistente a flexión, como lo demuestra [22] es:

r !
P Le
Mmax = M0 · sec · (2.54)
EI 2
Donde Le es la longitud efectiva para cualquier condición de apoyo del elemento.
Por tanto, el esfuerzo máximo resistente en el lado cóncavo del elemento se puede repre-
sentar por la siguiente ecuación denominada Ecuación de la Secante:
" r !#
P Mmax · d P ed Le P
σmax = + = 1 + 2 · sec (2.55)
A I A r 2r AE
Donde d es la distancia entre la fibra del extremo y el eje neutro de la sección transversal y
ed
r es su radio de giro. El factor 2 es conocido como la relación de excentricidad.
r
Introduciendo el concepto del radio del núcleo (s), propuesto por Timoshenko [25], se pre-
senta la siguiente igualdad:
Z I r2
s= = = (2.56)
A Ad d
donde Z es el módulo plástico de la sección transversal.
Por tanto se puede escribir una variación de la Ecuación de la Secante de la siguiente forma:
" r !#
P e Le P
σmax = 1 + · sec (2.57)
A s 2r AE
Al usar esta ecuación se puede encontrar el esfuerzo crı́tico de pandeo de un elemento o

estructura cargada excéntricamente.
2.6. Métodos de optimización

La optimización puede plantearse como el ajuste de las entradas o caracterı́sticas de un
proceso industrial, modelo matemático o experimento para lograr que la salida alcance un
punto mı́nimo o máximo dependiendo de los intereses. Si la salida es ganancia, el objetivo
es maximizar; y si son pérdidas, la salida es minimizar. El problema a optimizar se puede
plantear mediante un modelo matemático descrito por f : X ⊆ Rn → R donde X es el
espacio de búsqueda. Ası́ el problema de optimización puede definirse como:
min f (x)X⊆Rn sujeto a Ck ≤ 0 (2.58)
Donde X ⊆ Rn → R es un conjunto en el espacio Euclidiano de dimensión n, f (x) es un

conjunto de números reales, k es el número de restricciones y Ck proporciona la forma de la

restricción, la cual es una función que limita el espacio de búsqueda, a esto se le llama pro-
blema de optimización restringido. Un punto Xi se le llama punto factible si satisface todas
las restricciones y al conjunto de todos los puntos factibles se le llama conjunto factible de
búsqueda. En el caso donde las restricciones no se definen es un problema de optimización
no restringido, por lo mismo, el dominio de la función es el espacio factible de búsqueda [26].
Los diferentes campos de la optimización pueden clasificarse por las propiedades de la función
objetivo: la forma de las restricciones, el dominio, naturaleza de la función objetivo, etc. [26].
Matlab es un software de cálculo numérico que ofrece un ambiente para realizar operaciones
aritméticas y vectoriales en lı́nea de comandos. Además, ofrece un entorno de programación
de scripts que facilitan la implementación de códigos para la solución de problemas de gran
envergadura. Matlab ha sido utilizado en la implementación de algoritmos de optimización,
debido a su facilidad de programación y su poder de cómputo en el manejo de vectores y
matrices; y a sus herramientas para graficar funciones o series de datos en 2 y 3 dimensiones,
lo que permite analizar los resultados de la optimización de forma visual. Cuenta con un
conjunto de funciones predefinidas para realizar operaciones matriciales, para que el usuario
no invierta tiempo en la programación de tales operaciones.
Este software incluye módulos adicionales llamados “toolbox”, los cuales contienen entornos
para el desarrollo de aplicaciones más especı́ficas. Dentro de los toolbox se encuentra el de
optimización, el cual consta de un conjunto de comandos que determinan el mı́nimo de una
función lineal o no lineal que puede ser o no multivariable con restricciones de igualdad y
desigualdad.
De los comandos de optimización de Matlab se escogió fminsearch, una función que resuelve
los problemas de optimización de funciones de una o más de una variable sin restricciones,
que tiene la siguientes sintaxis:
f minsearch (f un, x0)
Donde fun es la función objetivo y x0 es un vector de valores iniciales de iteración. El co-

mando:
[x, f valor] = f minsearch (f un, x0, options)
devuelve el parámetro minimizado (x), y el valor de la función en x. En options se especifica la

estructura con la que el comando fminsearch optimizará, se utiliza optimset para establecer
los valores de options descritos en [26].

Capı́tulo 3
PANDEO EN BARRAS DE ACERO

BAJO CARGAS CÍCLICAS
3.1. Modelo propuesto del comportamiento del acero

bajo cargas cı́clicas
Como punto de partida, se realiza un análisis comparativo de las diferentes curvas esfuerzo
vs. deformación propuestas en la bibliografı́a consultada, tanto para cargas monotónicas
como para cargas cı́clicas.
3.1.1. Curva esfuerzo vs. deformación bajo carga monotónica

Para la obtención de una curva esfuerzo vs. deformación bajo carga monotónica, se requieren
los parámetros fy , fsu , εsh , εsu , Es y p en el caso del modelo de Mander. Estos parámetros
se obtuvieron de las pruebas experimentales realizadas por [3], ellos ensayaron 100 probetas
procedentes de distintos lugares de la ciudad de México en el año 1993. Las probetas en-
sayadas fueron de 600mm de longitud con un espacio libre entre mordazas igual a 80mm.
La velocidad del ensayo fue constante de manera que el mismo duró en promedio tres mi-
nutos. Las probetas ensayadas fueron 50 de diámetro grande dentro del cual se encuentran
diámetros de 5/8 pulg, 3/4 pug y 1 1/2 pulg, y 50 de diámetros de 5/16 pulg, 3/8 pulg y
1/2 pulgada pertenecientes a los diámetros pequeños. Se usaron 50 para la obtención del
parámetro fy y 45 para los demás parámetros debido a que 5 probetas de diámetro pequeño
se detectaron con posibles fallas en medición.
Los valores obtenidos de los distintos parámetros se muestran en la tabla 3.1
Cabe aclarar que el valor de p se obtuvo con ayuda de un punto experimental de la curva
esfuerzo vs. deformación (εsh1 , fsh1 ) en la zona de endurecimiento por deformación, y εsh1 se
obtuvo como un promedio entre εsh y εsu . Y en estos ensayos se consideró Es constante e
igual a 2 · 106 kg/cm2 .
Debido a que en muchos ensayos de probetas con diámetros pequeños, se obtuvo mucha
variabilidad en las estadı́sticas correspondientes, se optó por tomar como parámetros base
los obtenidos de los diámetros grandes.

Media Diámetros Media Diámetros

Media Global
Grandes Pequeños
fy (kg/cm2 ) 4591 4577 4609
fsu (kg/cm2 ) 7465 7491 7436
εsh 0.0075 0.0088 0.006
εsu 0.1289 0.1171 0.142
p 3.421 3.474 3.362
Cuadro 3.1: Parámetros de las probetas ensayadas por Botero [3].
Con los parámetros establecidos, se obtiene la curva de esfuerzo vs. deformación de una
barra de acero bajo carga de tensión con los modelos de Mander [6], Park y Paulay [11] y
Ahmad y Shah [12].
En la gráfica 3.1 se muestra las tres curvas obtenidas con los distintos modelos, con el propósi-
to de compararlas entre sı́.
Curva esfuerzo vs deformación bajo carga de tracción

8000
7000
6000 Modelo de Mander

Modelo de Park y Paulay
Modelo de Ahmad y Shah
Esfuerzo kg/cm 2
5000
4000
3000
2000
1000
0
0 0.02 0.04 0.06 0.08 0.1 0.12
Deformación
Figura 3.1: Curva esfuerzo vs. deformación bajo carga a tracción
Se considera que el modelo de Mander se adapta de mejor manera a las curvas experimenta-
les, debido a que el modelo de Park y Paulay y el modelo de Ahmad y Shah se basan en el
uso de una función cuadrática para representar la zona de endurecimiento por deformación.
Sin embargo, el modelo de Mander considera el parámetro p, el cual depende de un punto

dentro de la curva experimental en la zona de endurecimiento por deformación, por tanto

este modelo cuenta con una mejor relación con la curva experimental.
Posterior a la obtención de las curvas esfuerzo vs. deformación de acero bajo carga de tensión,
se procede a graficar las curvas esfuerzo vs. deformación bajo carga de compresión mediante
las fórmulas 2.17 y 2.18 de Dodd y Restrepo [14]. En la gráfica 3.2 se muestra las tres curvas
obtenidas con los distintos modelos junto a una curva experimental de una barra de acero
bajo carga de compresión, la cual pertenece al mismo lote de barras ensayadas bajo carga
de tensión para obtener los parámetros fsu , fy , εsh y εsu .
Figura 3.2: Curva esfuerzo vs. deformación bajo carga a compresión
Como se puede apreciar en la figura 3.2, el modelo de Mander es el que más se aproxima a
la curva experimental.
3.1.2. Curva esfuerzo vs. deformación bajo carga cı́clica.

Tanto el modelo de curva esfuerzo vs. deformación bajo carga cı́clica propuesto por Mander
[6], como el modelo de Kent y Park [20], utilizan parámetros experimentales que se obtienen
de un ensayo de tensión previo. Por esta razón, se propone un modelo aproximado del
comportamiento del acero bajo cargas cı́clicas que nos permita obtener la curva esfuerzo-
deformación basado en las caracterı́sticas básicas del material, las cuales se especifican en la
mayorı́a de catálogos que disponen las empresas distribuidoras.

3.1.3. Modelo propuesto

El siguiente modelo se basa en las ecuaciones propuestas por Menegotto-Pinto [18] y el mode-
lo propuesto por Chang y Mander [6], y tiene la finalidad de representar el comportamiento
del acero bajo cargas cı́clicas antes de que ocurra la falla del material.
Descripción del Modelo
Datos requeridos:
Los datos necesarios son las caracterı́sticas básicas del material que comúnmente se especifica
en los catálogos de las empresas distribuidoras, éstas son: el esfuerzo de fluencia a tracción
fy , y el módulo de elasticidad del material Es .
Además, es preciso ubicar las deformaciones máximas en cada uno de los ciclos compre-
sión-tracción que se vayan a representar en el modelo.
Datos supuestos:
Los valores iniciales de los esfuerzos y deformaciones últimas, y del factor p correspondientes
a las caracterı́sticas del acero, se obtienen del promedio de los ensayos producidos por Botero
[3] bajo carga de compresión para una esbeltez de 2.5. La tabla 3.2 muestra los valores de
las caracterı́sticas mencionadas
Tracción Compresión
fsut (kg/cm2 ) 10800 fsuc (kg/cm2 ) 11333
εsut 0.1949 εsuc 0.1629
p 2.092 p 3.474
Cuadro 3.2: Parámetros de barra ensayada por Botero con esbeltez 2.5
Para hacer uso de la ecuación 2.20 propuesta en el modelo de Chang y Mander [6], se
supone inicialmente los valores correspondientes al factor R tanto para tracción como para
compresión:
Tracción Compresión
Caracterı́stica Valor Caracterı́stica Valor
Rt 0.547 Rc 0.580
Cuadro 3.3: Factor R.
Datos calculados:
La deformación de la fluencia del material en tracción, se calcula al despejar la ecuación 2.1
y remplazar los datos que se han denominado como requeridos. Se obtiene la deformación del
inicio de la zona de endurecimiento εsh con la ecuación 2.16. Con el uso de las ecuaciones 2.17
y 2.18, se calcula el esfuerzo de fluencia del material Fyc , las deformaciones de fluencia εyc
y de inicio de la zona de endurecimiento por deformación εshc , correspondientes a compresión.
Zona de descarga y Zona elástica lineal

El modelo presentado supone que el material es sometido a compresión, posteriormente,

existe una descarga hasta que el esfuerzo es nulo, y empieza a someterse a tracción. De igual
manera, existe la zona de descarga al final de la tracción y de esta manera se completa un
ciclo, como se puede apreciar en la figura 3.3. La descarga tiene un comportamiento lineal
con pendiente igual al módulo elástico del material Es .
Esfuerzo,fs
Zona de descarga
Deformación, s
Zona de descarga
Zona elástica lineal
Figura 3.3: Zona elástica y zona de descarga
Zona de endurecimiento por deformación
Cuando el material se comprime por primera vez, su comportamiento en la zona de endure-

cimiento por deformación se representa mediante la ecuación 2.5 propuesta en el modelo de
Mander [10] para acero sometido a cargas monotónicas, que se demostró es la que mejor se
aproxima a los resultados experimentales de Botero [3].
Posterior a la primera compresión, el comportamiento del material en la zona de endureci-

miento por deformación se representa con las ecuaciones de Menegotto-Pinto [18] propuestas
en el modelo de Chang y Mander [6] para acero sometido a cargas cı́clicas, con las variaciones
que se describen a continuación:
Tracción

Esfuerzo,fs
( sut ,f tf )
to,0)
Deformación, s
Figura 3.4: Curva histerética en zona de tracción
La deformación al inicio de la curva es εt0 y es igual a la deformación que presenta el

material al final de la zona de descarga ante cargas de la compresión.
La deformación al final de la curva es la deformación última del material en tracción

εsut y su valor es constante.
El esfuerzo al inicio de la curva es nulo.
El esfuerzo al final de la curva es ftf y su valor es igual al del esfuerzo último del
material en tracción fsut multiplicado por un factor de reducción que depende de εt0 .
El parámetro R, que regula la curvatura de la curva en tracción, es igual a Rt del

cuadro 3.3 más un factor que depende de la deformación inicial εt0 .
Después de un análisis de varias curvas histeréticas experimentales, se observa que mientras

mayor es la deformación inicial de la curva εt0 , el esfuerzo final disminuye y la curvatura de
la curva aumenta, por lo tanto el parámetro R. Debido a esto, se establecieron las últimas
dos variaciones mencionadas.
Las ecuaciones propuestas para tracción son:
ftf = fsut · (1 − 2.4207 · εt0 ) (3.1)
ftf
Esec = (3.2)
εsut − εt0

εsut + εt0
∆ε = (3.3)
2
E0 = Es · (1 − ∆ε) (3.4)
R = Rt − (εt0 · 4.4715) (3.5)
" R #−1/R
E0
a= 1+ (3.6)
Esec
Esec
−a
E0
Q= (3.7)
1−a
R
εS − εt0
ρ = E0 · (3.8)
ftf
( )
1−Q
fs = f0 + E0 · (εS − εt0 ) Q + (3.9)
{1 + ρ}1/R
Compresión:

Esfuerzo,fs
c0,0)
Deformación, s
( suc ,fcf )
Figura 3.5: Curva histerética en zona de compresión
La deformación al inicio de la curva es εc0 y es igual a la deformación que presenta el

material al final de la zona de descarga después de la tracción.
La deformación al final de la curva es la deformación última del material en compresión

εsuc y su valor es constante.
El esfuerzo al inicio de la curva es nulo.
El esfuerzo al final de la curva es fcf y su valor es igual al del esfuerzo último del
material en compresión fsuc multiplicado por un factor de reducción que depende de
εc0 .
El parámetro R, que regula la curvatura de la curva en compresión, es igual a Rc del

cuadro 3.3.
Las ecuaciones propuestas bajo carga de compresión son:
fcf = fsuc · (1 − 0.9529 · εc0 ) (3.10)
fcf
Esec = (3.11)
εsuc − εc0

εsuc + εc0
∆ε = (3.12)
2
E0 = Es · (1 − ∆ε) (3.13)
R = Rc (3.14)
" R #−1/R
E0
a= 1+ (3.15)
Esec
Esec
−a
E0
Q= (3.16)
1−a
R
εS − εc0
ρ = E0 · (3.17)
fcf
( )
1−Q
fs = f0 + E0 · (εS − εc0 ) Q + (3.18)
{1 + ρ}1/R
A continuación se indica el proceso que sigue el modelo propuesto para la graficación de los
ciclos histeréticos, mostrado en la figura 3.6.
1. Se calcula la curva AB con el modelo de Mander bajo cargas monotónicas, y se grafica

hasta la deformación deseada ε1 . Posteriormente, se grafica la descarga que tiene pendiente
igual a el módulo elástico.
2. Se calcula la curva CD con el modelo de Mander bajo cargas cı́clicas junto a las varia-
ciones propuestas anteriormente para tracción, y se grafica hasta la deformación deseada ε2 .
Seguidamente, se grafica la descarga.
3. Se calcula la curva EF con el modelo de Mander bajo cargas cı́clicas junto a las variaciones
propuestas anteriormente para compresión, y se grafica hasta la deformación deseada ε3 .
4. Se repite el paso 2 y 3 para la gráfica de más ciclos.

Esfuerzo,fs
fsut
D ftf
suc 3 1 C to E
A co 2 sut
Deformación, s
fcf F
fsuc B
Figura 3.6: Graficación de ciclos histeréticos.
En las figuras 3.7, 3.8 y 3.9, se compara las curvas esfuerzo vs. deformación bajo cargas cı́cli-
cas del modelo propuesto junto con las curvas experimentales obtenidos por Botero [3]. El
modelo propuesto tiene una aproximación aceptable con relación a los resultados experimen-
tales comparados. Este modelo representa adecuadamente el comportamiento del acero bajo
cargas cı́clicas y tiene como ventaja que los datos requeridos son propiedades del material
que se pueden obtener fácilmente y las deformaciones finales de cada ciclo histerético.
El modelo se elaboró en el software matemático Matlab. El código del programa se presenta

en el Apéndice A.

Curva esfuerzo vs deformación bajo cargas cíclicas

14000#
12000
10000
8000
6000
Esfuerzo kg/cm 2
4000
2000
0
-2000
-4000
-6000
-8000
-10000
Curva experimental Botero y Rodríguez
-12000
Modelo propuesto
-14000
-0.15 -0.125 -0.1 -0.075 -0.05 -0.025 0 0.025 0.05 0.075 0.1 #0.125 0.15
Deformación
Figura 3.7: Modelo propuesto con la curva de esbeltez 2.5 de Botero [3].

10000
8000
6000
4000
Esfuerzo kg/cm 2
2000
-2000
-4000
-6000
-8000 Curva experimental Botero y Rodríguez

Modelo propuesto
-1#0000
-0.05 -0.04 -0.03 -0.02 -0.01 0 0.01 0.02 0.03 0.04 0.05
Deformación
Figura 3.8: Modelo propuesto con la curva de esbeltez 4 de Botero [3].


10000
8000
6000
4000
Esfuerzo kg/cm 2
2000
-2000
-4000
-6000
-8000 Curva experimental Botero y Rodríguez

Modelo propuesto
-1#0000
-0.03 -0.02 -0.01 0 0.01 0.02 0.03
Deformación
Figura 3.9: Modelo propuesto con la curva de esbeltez 6 de Botero [3].
3.2. Modelo de barra

A continuación, se establece el esfuerzo crı́tico de pandeo para una barra con cualquier con-
dición de apoyo. Con el fin de representar una barra entre dos estribos.
Las hipótesis son:
Una barra de sección constante, con cualquier condición de apoyo.
Figura 3.10: Modelo de barra
Existe una carga axial P centrada en el eje de la barra, ası́ como un cortante Q y dos mo-
mentos genéricos en los extremos de la misma [27].
Al analizar una longitud x, se puede representar la ecuación de equilibrio de la siguiente

manera [27]:

Figura 3.11: Esfuerzos en los extremos de la barra
Figura 3.12: Corte de longitud x en la barra
M + P · w − MA + Q · x = 0 (3.19)
Al considerar M = EIwxx , la ecuación puede reescribirse:
EIwxx + P · w = M A − Q · x (3.20)
Derivando dos veces la ecuación con respecto a x :
EIwxxxx + P wxx = 0 (3.21)
La ecuación diferencial 3.21 representa el comportamiento de la barra sometida a carga axial

y cualquier condición de apoyo. Las soluciones de la ecuación son:
w=0 (3.22)
Con multiplicidad 2, y
r
P
w=± i (3.23)
EI

La solución general es:

P P
w = A1 seno x + A2 cos x + A3 + A4 x (3.24)
EI EI
Condiciones de frontera:
x=0
−[EIwxxx + P wx ] = α¯0 w (3.25)
EIwxxx = β¯0 w (3.26)
x=L
−[EIwxxx + P wx ] = −α¯L w (3.27)
EIwxxx = −β¯L w (3.28)
Al considerar:
ᾱ0
α0 = (3.29)
EI
ᾱL
αL = (3.30)
EI
β̄0
α0 = (3.31)
EI
β̄L
αL = (3.32)
EI
Se obtiene:
wxxx + K 2 wxx = 0 (3.33)
wxxx (0) + K 2 wx (0) + α0 w(0) = 0 (3.34)
wxx (0) − β0 wx (0) = 0 (3.35)

wxxx (L) + K 2 wx (L) + αL w(L) = 0 (3.36)
wxx (L) − βL wx (L) = 0 (3.37)
Al sustituir las condiciones de frontera:
α0 A2 + K 2 A3 + α0 A4 = 0 (3.38)
β0 KA2 + K 2 A2 + β0 A3 = 0 (3.39)
(αL seno(KL))A1 + (αL cos(KL))A2 − (K 2 − αL L)A3 + αL A4 = 0 (3.40)
(KβL cos(KL) − K 2 seno(KL))A1 − (βL Kseno(KL) + K 2 cos(KL))A2 + βL A3 = 0 (3.41)
Para que exista una solución diferente de la trivial:

2

0 α 0 K α0

2

β 0 K K β0 0
=0
2

αL sin(KL) αL cos(KL) αL L − K αL
(KβL cos(KL) − K 2 sin(KL)) −(βL Ksin(KL) + K 2 cos(KL)) βL 0
Al resolver el determinante se obtiene la siguiente ecuación caracterı́stica:
Si KL = u
u6 u4 u2

−(α0 + αL ) 6 + {β0 βL (α0 + αL ) + α0 αL L} 4 + α0 αL (β0 + βL − β0 βL L) 2 sin(u)
L L L
5 3

u u u
+ (α0 + αL )(β0 + βL ) 5 − α0 αL L(β0 + βL ) 3 − 2α0 αL β0 βL cos(u)
L L L
u
+2α0 αL β0 βL = 0
L
(3.42)
Cuya solución general se puede escribir de la siguiente manera y hace referencia a la ecuación
de la carga de pandeo crı́tico para una columna con cualquier condición de apoyo[28].
u2 EI
P crit = (3.43)
L2

u2 EI
σcrit = (3.44)
AL2
Si u = π/k donde k es el factor de longitud efectiva, y reemplazando la inercia del cı́rculo
π · D4 π · D2
I= y su área A = :
64 4
π2E
σcrit = 2 (3.45)
KL
16 ·
D
A continuación se presentan un ejemplo para mejor entendimiento. Para el caso de que la

barra se encuentre articulada en sus dos extremos, α0 y αL son infinito ya que no existe
un desplazamiento, mientras que β0 y βL son 0 debido a que la articulación acepta el giro.
Reemplazando estos valores de α y β en la ecuación 3.42:
u6 u4

−(αo + αL ) 6 + {αo · αL · L} 4 sin(u) = 0
L L
Dividiendo para ∞2 y aplicando los lı́mites de α0 → ∞ y αL → ∞:
u4
· sen (u) = 0
L3
Por lo tanto, u = π y la fuerza crı́tica de pandeo es igual a la fuerza de Euler:
π 2 EI
P crit =
L2
Debido a que el modelo de la barra presentado previamente representa al comportamiento

de un refuerzo longitudinal entre dos estribos, se puede asumir que α0 y αL representan la
rigidez lateral proporcionada por el estribo inferior y superior respectivamente, por lo que
se encuentran en función del diámetro de los estribos; mientras que β0 y βL representan la
rigidez rotacional del refuerzo longitudinal y dependen del diámetro del mismo.
E·A E · π · De2
α= = (3.46)
L 4·L
4 · Er · I Er · Dl4
β= · k1 = · k1 (3.47)
L 16 · L
Donde:
A = área de la sección transversal del refuerzo longitudinal
I = inercia de la sección transversal del refuerzo longitudinal

De = diámetro de los estribos

Dl = diámetro del refuerzo longitudinal
L = separación entre estribos
E = módulo de elasticidad del acero
Er = doble módulo reducido del acero
k1 = factor de reducción de la rigidez rotacional
Finalmente, con ayuda de la ecuación de la secante 2.55 o su variación 2.57, se determina el

esfuerzo crı́tico de pandeo del refuerzo longitudinal para una excentricidad e o relación de
ed
excentricidad 2 conocida.
r
3.3. Relación entre el pandeo y la esbeltez

Como se ha indicado anteriormente, para la obtención de la carga crı́tica con pandeo inelásti-
co se puede utilizar el módulo tangente o el doble módulo reducido.
Ciertos resultados experimentales obtenidos con especı́menes de distintas secciones trans-

versales, indican que las cargas de pandeo obtenidas se encuentran entre las predichas por
estas dos teorı́as. Estas dudas fueron finalmente aclaradas por Shanley en 1946 [29].
Shanley demostró que cuando se somete a compresión a una columna recta hasta alcanzar
una posición deformada en equilibrio, los esfuerzos del lado cóncavo de la sección aumentan
asociados al módulo tangente Et , pero al aumentar un diferencial de carga, no se presenta
una disminución inmediata de los esfuerzos del lado convexo ya que éstos son compensados
por los esfuerzos de compresión directa. Esto genera que en el instante en el que ocurre el
pandeo, el valor obtenido por el doble módulo reducido disminuya, por lo que el esfuerzo de
pandeo de la columna se encuentra entre los predichos por la teorı́a del módulo tangente y la
teorı́a del doble módulo reducido. Por tanto, se puede indicar que la teorı́a del módulo tan-
gente puede ser considerada como un lı́mite inferior mientras que la teorı́a del doble módulo
reducido puede considerarse como un lı́mite superior [7].
En este estudio, se emplea el doble módulo reducido, ya que considera el aumento y re-
ducción de los esfuerzos en el lado convexo y cóncavo respectivamente, producidos por la
flexión [3]. A más de esto, Mander[10] ha demostrado que la teorı́a del doble módulo reducido
presenta una mejor aproximación que la teorı́a del módulo tangente. Por lo tanto, el esfuerzo
crı́tico de pandeo es:
u2 Er I π 2 Er
σcrit = = 2 (3.48)
AL2 16( kL )
D
Para obtener el esfuerzo crı́tico de pandeo de una barra bajo carga cı́clica, se necesita la
gráfica esfuerzo vs. deformación del acero que represente el comportamiento del mismo para
la obtención del doble módulo reducido. La gráfica a utilizar es la curva experimental es-
fuerzo vs. deformación de una barra de acero de esbeltez 2.5 sometida a carga cı́clica con
ciclos simétricos obtenida por Botero [3]. Se escogió esta gráfica ya que es de interés en este
estudio el comportamiento del acero antes del pandeo puesto que el objetivo es ubicar el
punto donde empieza el mismo. Debido a una esbeltez de 2.5, la barra no llega al pandeo en

ninguno de sus ciclos.
Las caracterı́sticas del acero para la obtención de su curva esfuerzo vs. deformación son
fy de 4836 kg/cm2 y Es de 2000000kg/cm2 . Utilizando estas caracterı́sticas y con ayuda del
modelo propuesto en el inciso 3.1.3, se obtiene la gráfica esfuerzo vs. deformación bajo carga
cı́clica del acero estudiado por Botero [3].
Una vez obtenida la gráfica de esfuerzo vs. deformación deseada, con ayuda de la ecua-
ción 3.44, se procede a graficar el esfuerzo de pandeo σp vs. (kL/D) que hace referencia al
factor de longitud efectiva por la relación de esbeltez, es decir la separación entre estribos
sobre el diámetro de la barra.
El proceso que se realiza para la obtención de la curva esfuerzo vs. relación de esbeltez
es el siguiente:
1. Se escogen las condiciones de apoyo de la barra que se requiere analizar, se obtienen tres
curvas con k de 0.5, 0.75 y 1 respectivamente.
2. Se escoge una relación de esbeltez L/D.
3. Se escoge un esfuerzo de la curva esfuerzo vs. deformación obtenida anteriormente σ1 y

se obtiene su módulo reducido.
4. Con el módulo reducido obtenido y la relación de esbeltez escogida, se utiliza la ecuación

2.55 ó 2.57, para una excentricidad o relación de excentricidad conocida, y se obtiene σ2 .
5. Si σ1 y σ2 son iguales, entonces σ1 es el esfuerzo de pandeo crı́tico de la barra con las

condiciones de apoyo y la relación de esbeltez indicadas.
6. Si σ1 y σ2 no son iguales se debe escoger de la gráfica esfuerzo vs. deformación otro σ1

hasta que éste sea igual a σ2 .
7. Se repite el mismo proceso para una relación de esbeltez desde 2 hasta 13, cada 0.1. No se
tomaron en cuenta relaciones de esbeltez menores a 2 debido a que no se registra pandeo.
8. Como a cada esbeltez le corresponde un esfuerzo de pandeo crı́tico para un mismo k, se

procede a graficar.
La figura 3.13 presenta las gráficas de esfuerzo de pandeo crı́tico vs. relación de esbeltez con
un k de 0.5, 0.75 y 1, junto con los puntos experimentales que hacen referencia al esfuerzo
crı́tico de las barras de distintas relaciones de esbeltez ensayadas por Botero [3]. Debido a
que los datos experimentales se obtienen de ensayos en barras de acero aisladas, se consideró
una relación de excentricidad nula.
Una vez obtenida las curvas esfuerzo de pandeo crı́tico vs. relación de esbeltez, se procede a
la elaboración de curvas deformación de pandeo crı́tico vs. relación de esbeltez.
A cada esfuerzo de pandeo crı́tico le corresponde una esbeltez y una deformación de pandeo
crı́tico debido a la curva esfuerzo vs. deformación, por lo tanto se grafica la deformación de

12000
Ensayos cíclicos reversibles con ciclos simétricos
11000 Ensayos cíclicos reversibles con ciclos asimétricos
10000
9000
Esfuerzo de pandeo
8000
7000
6000 K=0.5
5000
K=0.75
4000
K=1.0
3000
2000
0 2 4 6 8 10 12 14
L/D
Figura 3.13: Esfuerzo de pandeo vs. relación de esbeltez
Esfuerzo de pandeo
Esbeltez
(kg/cm2 )
4 Ciclos simétricos 8486

4 Ciclos asimétricos 7300
Cuadro 3.4: Esfuerzos crı́ticos de pandeo de probetas ensayadas por Botero [3]
pandeo crı́tico junto a la esbeltez correspondiente.
Para obtener la deformación correspondiente al esfuerzo de pandeo crı́tico, se utiliza el mo-

delo planteado por Botero [3]. Donde la deformación εo ∗ corresponde al punto en el cambio
de carga de tensión a carga de compresión, en el ciclo donde ocurre el pandeo para la defor-

mación εp . Por tanto, se considera que la deformación crı́tica de pandeo es igual a:
εp ∗ = ε0 ∗ −εp (3.49)
Esfuerzo
+ +
(fm,εm)
Deformación
(fp,εp)
- -
εp ε+o
(fm,εm)
εp*
Figura 3.14: Deformación de pandeo crı́tico propuesto por Botero [3]
Los puntos experimentales ensayados por Botero [3], mostrados en las figuras 3.13 y 3.15,
siguen la forma de la curva propuesta debido a que pertenecen a barras del mismo lote del
acero de esbeltez 2.5 al que pertenece la curva esfuerzo vs. deformación considerada.
La programación en el programa Matlab para la obtención de las gráficas se muestra en

el Apéndice B del presente trabajo.

0.2
0.18 Ensayos cíclicos reversibles con ciclos asimétricos
0.16
Deformación de pandeo
0.14
0.12
0.1
0.08
0.06
0.04
K=0.5
0.02
K=0.75
K=1.0
0
0 2 4 6 8 10 12 14
L/D
Figura 3.15: Deformación de pandeo vs. relación de esbeltez
Deformación de
Esbeltez
pandeo εp
4 Ciclos simétricos 0.0722

4 Ciclos asimétricos 0.0690
Cuadro 3.5: Deformaciones crı́ticas de pandeo de probetas ensayadas por Botero [3]

3.3.1. OPTIMIZACIÓN Y OBTENCIÓN DEL FACTOR DE LON-

GITUD EFECTIVA
Las gráficas 3.13 y 3.15 muestran que los puntos experimentales de las probetas ensayadas
por Botero se acercan más a la curva con un factor de longitud efectiva k igual a 0.75 que
a las curvas de k igual a 0.5 y 1. Esto sucede porque al momento de realizarse los ensayos
de las probetas, no se puede garantizar un empotramiento perfecto debido a la existencia de
pequeños movimientos entre las probetas y el dispositivo de agarre, lo que produce un giro
en los extremos. Por tal razón, el factor de longitud efectiva de los ensayos es mayor a un
k de 0.5 que pertenece al empotramiento perfecto. Se consideran los valores de la longitud
efectiva k de 0.5 y 1 como los lı́mites superior e inferior que representan las condiciones de
apoyo doblemente articulada y doblemente empotrada respectivamente.
La optimización tiene como objetivo calibrar parámetros determinados de un modelo para

minimizar las diferencias entre los resultados experimentales y los resultados obtenidos por
el modelo. En este caso el parámetro a calibrar será el factor de longitud efectiva k, de ma-
nera que se obtenga las mı́nimas diferencias entre los puntos experimentales de esfuerzo de
pandeo crı́tico y deformación de pandeo crı́tico, y las curvas esfuerzo crı́tico vs. relación de
esbeltez y deformación crı́tica vs. relación de esbeltez.
En el presente trabajo, la calibración de la curva se realiza mediante el programa Matlab.

Ya que representa una buena herramienta de optimización de una o varias variables.
Se define εexpi como el rango de deformaciones de pandeo crı́tico experimentales, x un vector

de un término que representa el parámetro a optimizar, en este caso el factor de longitud
efectiva y εi (x) el rango de valores de las deformaciones de la curva deformación crı́tica vs.
relación de esbeltez propuesta con cada esbeltez correspondiente a las deformaciones expe-
rimentales. El error en cada punto de deformación es:
Ri (x) = εexpi − εi (x) (3.50)
Se define a la función de error igual a la suma de los cuadrados de los errores, por lo tanto,
la función error queda:
X
f (x) = Ri 2 (x) (3.51)
Aplicando la función fminsearch del programa Matlab, se obtiene el parámetro (x ) óptimo

para que la función error sea mı́nima.
Se procede a indicar la obtención de la función error, para la cual se necesitaron una se-
rie de funciones en Matlab que se explican a continuación una a una y se presenta el código
de cada función en anexos.

1. La primera función tiene el nombre de Erm, pide como dato de entrada la deformación
crı́tica experimental de cada ensayo y tiene como finalidad devolver el valor del módulo
reducido perteneciente a cada deformación crı́tica de entrada.
2. La segunda función llamada K tiene como datos de entrada la esbeltez, el diámetro de

la barra, y la deformación crı́tica de pandeo de cada ensayo experimental. Primero se
obtiene la longitud de la barra igual a la multiplicación de la esbeltez por el diámetro de
la barra. Se procede a obtener α y β en función de k1 (factor de reducción de la rigidez
rotacional) con las ecuaciones 3.46 y 3.47. Debido a que los ensayos de Botero [3] tenı́an
barras sujetas en sus extremos, se considera que la rigidez de los estribos es muy grande
y por lo tanto se utiliza un diámetro de estribo muy grande. Mediante la función Erm
se obtiene el módulo reducido para cada deformación crı́tica y con ayuda de la ecuación
3.42, esta función devuelve el factor de longitud efectiva en función de k1 .
3. La función numero tres es Fs, la cual tiene como datos de entrada la esbeltez de cada
ensayo, la deformación crı́tica de pandeo, el diámetro de cada barra y el diámetro de los
estribos. Tiene como objetivo devolver la carga crı́tica de pandeo en función de k1 con
ayuda de la ecuación 3.48, obteniéndose el módulo reducido de la función Erm y el factor
de longitud efectiva de la función K.
4. La cuarta función se llama Df, la cual tiene como objeto devolver la deformación de la
gráfica esfuerzo vs. deformación perteneciente al esfuerzo de entrada obtenido de la fun-
ción Fs.
5. La quinta función Fe hace referencia a la función error, ya que suma las diferencias al
cuadrado entre las deformaciones obtenidas de la función Df y las deformaciones crı́ticas
de cada ensayo experimental.
6. Finalmente, con ayuda de fminsearch se optimiza la función Fe de manera que se obtiene

el k1 que produce el valor mı́nimo de la función error.
Se realiza el proceso de optimización tanto con las deformaciones crı́ticas de pandeo como
con los esfuerzos crı́ticos de pandeo de los datos obtenidos por Botero [16]. En el caso de los
esfuerzos, la función Erm tiene como datos de entradas los esfuerzos crı́ticos experimentales
y el error se calcula como el cuadrado de la diferencia entre los esfuerzos calculados en el
paso 3 y los de entrada.
En la gráfica 3.16 se presenta el valor de la función error vs. el número de iteraciones del
proceso de optimización, para la optimización con las deformaciones.

Current Function Value: 0.00154827

0.02
0.018
0.016
0.014
Function value 0.012
0.01
0.008
0.006
0.004
0.002
0
0 5 10 15 20 25
Stop Pause Iteration
Figura 3.16: Función error vs. número de iteraciones
La gráfica 3.17 muestra el valor de la función error vs. el número de iteraciones del proceso
de optimización, para la optimización con los esfuerzos.
×10 8 Current Function Value: 8.73244e+06

3.5
2.5
Function value
1.5
0.5
0
0 10 20 30 40 50 60 70 80 90 100
El factor de reducción de la rigidez rotacional k1 optimizado es igual a 1.26e-05 y 1.068e-05

para deformaciones y esfuerzos respectivamente, los mismos que corresponden a un factor
de longitud efectiva promedio de 0.755
En la tabla 3.6, se presenta el valor de las cargas crı́ticas de pandeo y las deformacio-
nes crı́ticas de pandeo calculadas con el k de 0.755
La tabla 3.7 presenta las cargas crı́ticas de pandeo experimentales y las cargas calculadas

12000
11000 Ensayos cíclicos reversibles con ciclos asimétricos
10000
9000
Esfuerzo de pandeo
8000
7000
6000
5000
K=0.755
4000
3000
0 2 4 6 8 10 12 14
L/D
Figura 3.18: Esfuerzo crı́tico de pandeo vs. esbeltez
0.2
0.18 Ensayos cíclicos reversibles con ciclos asimétricos
0.16
0.14
0.12
0.1
0.08
0.06
0.04
0.02 K=0.755
0
0 2 4 6 8 10 12 14
L/D
Figura 3.19: Deformación crı́tica de pandeo vs. esbeltez
junto a la relación de estos dos valores. La relación promedio entre los esfuerzos crı́ticos
experimentales y calculados es de 0.98, con un coeficiente de variación del 6 % por lo que es

Esfuerzo Esfuerzo Deformación

Deformación
Esbeltez de de pandeo de pandeo
de pandeo
pandeo calculado calculado
4 Ciclos simétricos 8486 8056 0.072 0.05742
4 Ciclos asimétricos 7300 8056 0.069 0.05742
Cuadro 3.6: Esfuerzos y deformaciones crı́ticos experimentales y calculados.
aceptable.
Esfuerzo de Esf exp

Esfuerzo de
pandeo
pandeo Esf calculado
calculado
8486 8056 1.05
8000 8056 0.99
7930 8056 0.98
7300 8056 0.89
7250 8056 0.88
6900 6304 1.09
6681 6304 1.06
6484 6304 1.03
6030 6304 0.96
5956 6304 0.94
5204 5182 1.00
4750 5182 0.92
5186 5182 0.99
5088 5182 0.97
Cuadro 3.7: Error entre esfuerzos crı́ticas experimentales y calculadas
En la tabla 3.8 se presenta las deformaciones crı́ticas de pandeo experimentales y las defor-
maciones calculadas junto con la relación de estos dos valores. La relación promedio entre las
deformaciones crı́ticas experimentales y calculadas es de 1.15, con un coeficiente de variación
de 17 % por lo que es aceptable.

Deformación Def exp

Deformación de
de pandeo
pandeo Def calculado
calculado
0.0722 0.05742 1.26
0.05 0.05742 0.87
0.05 0.05742 0.87
0.069 0.05742 1.20
0.069 0.05742 1.20
0.039 0.02584 1.53
0.031 0.02584 1.31
0.0211 0.02584 0.77
0.0472 0.02584 1.71
0.0459 0.02584 1.67
0.0123 0.01514 0.63
0.0097 0.01514 0.48
0.0236 0.01514 1.47
0.019 0.01514 1.18
Cuadro 3.8: Error entre deformaciones crı́ticas experimentales y calculadas

Capı́tulo 4
PANDEO EN BARRAS DE ACERO

DEL REFUERZO LONGITUDINAL
DE ESTRUCTURAS DE
HORMIGÓN ARMADO BAJO
CARGAS CÍCLICAS
Como parte de una estructura de hormigón armado, el esfuerzo crı́tico de pandeo en barras
de acero pertenecientes al refuerzo longitudinal se ve afectado por la resistencia del hor-
migón, el armado del refuerzo transversal y las excentricidades que se presentan. La relación
de esbeltez depende directamente de la separación del refuerzo transversal y por lo tanto se
relaciona con la capacidad resistente de las barras de acero del refuerzo longitudinal, mientras
que las excentricidades que se presentan en el momento de aplicar una carga axial, reducen
la resistencia de las barras de acero del refuerzo longitudinal y se deben principalmente a la
expansión ocurrida en el hormigón al ser comprimido, que a su vez depende de la resistencia
del mismo.
En este capı́tulo, se utilizan los resultados obtenidos en ensayos experimentales de muros rec-
tangulares de concreto reforzado sometidos a cargas cı́clicas ensayados por [30], [31] y [32];
y presentados en [4] y ensayos experimentales de columnas circulares de concreto reforzado
sometidas a cargas cı́clicas realizados por Mervyn Kowalsky [5].
Como un primer enfoque, se muestra en la figura 4.1 las gráficas deformación crı́tica de
pandeo vs. relación de esbeltez obtenidas en el Capı́tulo 3, con factores de longitud efectiva
k iguales a 0.5, 0.75 y 1, junto con los datos de elementos experimentales. Se observa que los
puntos de los muros se encuentran entre las curvas con factor de longitud efectiva k igual a
0.5 y 0.75. Además, los puntos de las columnas presentan deformaciones crı́ticas de pandeo
bajas en relación con las curvas mostradas, una de las causas posibles es la presencia de
excentricidades. Razón por la cual, se realiza el siguiente análisis que permite comprobar el
efecto que las excentricidades tienen en el refuerzo longitudinal de estructuras de hormigón
armado.

Relación Deformación
Denominación de esbeltez de pandeo
L/D ε∗
WSH3 6.3 0.032
WSH2 7.50 0.028
WSH5 6.30 0.031
WSH6 4.50 0.039
Muros
RW1 8.00 0.031
RW2 5.30 0.040
W-MC-N 5.00 0.049
TW1 8.00 0.021
N-9 2.67 0.061
N-13 3.66 0.05
N-14 5.33 0.0457
N-15 3.67 0.048
N-16 2.00 0.063
N-19 2.67 0.054
N-20 2.67 0.0514
N-21 2.67 0.059
Columnas N-22 2.67 0.056
N-23 2.67 0.058
N-24 2.67 0.052
N-25 2.29 0.052
N-26 2.29 0.05
N-27 2.67 0.0618
N-28 2.67 0.058
N-29 2.67 0.055
N-30 2.67 0.046
Cuadro 4.1: Resultados de ensayos experimentales en muros y columnas bajo cargas cı́clicas
obtenidos en [4] y [5].

0.2
Ensayos de muros rectangulares bajo cargas cíclicas
0.18 Ensayos de columnas circulares bajo cargas cíclicas
0.16
0.14
0.12
0.1
0.08
0.06
0.04
K=0.5
0.02 K=0.75
K=1.0
0
0 2 4 6 8 10 12 14
L/D
Figura 4.1: Gráfica deformación crı́tica vs. relación de esbeltez con datos experimentales
obtenidos en [4] y [5]
.
Las curvas deformación crı́tica vs. relación de esbeltez mostradas en la figura 4.1 representan
el comportamiento de las barras de acero pertenecientes al refuerzo longitudinal de muros
y columnas, sin presencia de excentricidades. Para analizar el efecto de las mismas, se usa
ed
la ecuación de la secante(Ec.2.55), con el fin de obtener la relación de excentricidad 2 de
r
cada uno de los elementos ensayados. Los demás términos de la ecuación se obtienen de la
siguiente manera:
P
El factor es el esfuerzo crı́tico de pandeo de cada elemento. Debido a que disponemos
A
solamente de las deformaciones crı́ticas de pandeo, se utiliza el modelo del material
propuesto en este trabajo para determinar el esfuerzo correspondiente a cada una de
las deformaciones, con las caracterı́sticas de la barra de acero de esbeltez 2.5, de acuerdo
a lo explicado en el Capı́tulo 5.
El término E es el doble módulo reducido correspondiente a cada una de las defor-

maciones crı́ticas de los elementos y se obtiene de la misma manera que el anterior
factor.
El radio de giro r depende de las caracterı́sticas geométricas de cada barra de acero

obtenidas de [4] y [5].
La longitud efectiva Le = k·L, es la separación del refuerzo transversal de cada elemento

multiplicado por el factor de longitud efectiva que, para este análisis se utilizan valores

de 0.5, 0.6 y 0.75.
Las figuras 4.2, 4.3 y 4.4 presentan la variación de la relación de excentricidad con respecto
a la relación de esbeltez para los diferentes factores de longitud efectiva usados en el análisis:
k=0.5
0.4
Columnas - Kowalsky
Muros - Iñiguez
0.35
Línea de tendencia
0.3
0.25
ed/r2
0.2
0.15
0.1
0.05
0
2 3 4 5 6 7 8
L/D
Figura 4.2: Gráfica relación de excentricidad vs. relación de esbeltez, con k=0.5
k=0.6
0.35
Columnas - Kowalsky
Muros - Iñiguez
0.3 Línea de tendencia
0.25
0.2
ed/r2
0.15
0.1
0.05
0
2 3 4 5 6 7 8
L/D

k=0.75
0.25
Columnas - Kowalsky
Muros - Iñiguez
Línea de tendencia
0.2
ed/r2 0.15
0.1
0.05
0
2 3 4 5 6 7 8
L/D
Como se puede apreciar en las figuras 4.2, 4.3 y 4.4, la relación de excentricidad disminuye
con el aumento de la relación de esbeltez hasta alcanzar valores muy cercanos a cero en
relaciones de esbeltez de 6 en adelante aproximadamente.
Se procede a realizar la optimización y obtención del factor de longitud efectiva de la curva

deformación crı́tica vs. relación de esbeltez obtenida en el capı́tulo anterior con los datos
de deformación crı́tica de los muros de esbelteces mayores e iguales a 6, debido a que su
capacidad de resistencia es afectada de manera mı́nima por la presencia de excentricidades
o éstas son nulas.
4.1. Optimización y obtención del factor de la longitud

efectiva
Siguiendo el mismo procedimiento para la optimización y obtención del factor de longitud
efectiva descrito en el Capı́tulo 3, con la ayuda del software Matlab, se calibra la longitud
efectiva k con la variación del factor que reduce la rigidez rotacional k1.
En la gráfica 4.5 se muestra el valor de la función error vs. el número de iteraciones del
proceso de optimización, generado por Matlab, para la optimización con las deformaciones.
El k1 optimizado es igual a 1.15e-04, el mismo que corresponde a un factor de longitud efec-

tiva promedio de 0.632
En la tabla 4.2, se presenta el valor de las deformaciones crı́ticas de pandeo calculadas

con el el factor de longitud efectiva k igual a 0.632.

×10 -3 Current Function Value: 0.00371055

3.722
3.72
3.718
Function value
3.716
3.714
3.712
3.71
0 10 20 30 40 50 60 70 80
Deformación
Deformación
Esbeltez de pandeo
de pandeo
calculado
6.3 Muro-WSH3 0.032 0.3296
6.3 Muro-WSH5 0.031 0.3296
7.5 Muro-WSH2 0.028 0.0237
8 Muro-RW1 0.031 0.02106
8 Muro-TW1 0.021 0.02106
Cuadro 4.2: Deformaciones crı́ticas experimentales y calculados.
La tabla 4.3 presenta las deformaciones crı́ticas de pandeo experimentales y calculadas junto
a la relación entre los dos valores anteriores. El promedio de la relación entre las deforma-
ciones crı́ticas de pandeo experimentales y calculadas es 1.1, por lo que es aceptable.
Deformación
Deformación de Def exp
de pandeo Def calculado
pandeo
calculado
0.032 0.03296 0.97
0.031 0.03296 0.94
0.028 0.0237 1.18
0.031 0.02106 1.47
0.021 0.02106 1
Cuadro 4.3: Relación entre deformaciones crı́ticas experimentales y calculadas

Una vez obtenido el factor de longitud efectiva óptimo, se realiza nuevamente el análisis del
comportamiento de la relación de excentricidad con respecto a la relación de esbeltez con
k=0.632, mediante la ecuación 2.57, con el fin de encontrar una función que represente dicho
comportamiento y tenerlo en cuenta para la obtención de la curva deformación crı́tica de
pandeo vs. relación de esbeltez de elementos de hormigón armado.
La figura 4.6 muestra las relaciones de excentricidad correspondientes a las relaciones de
esbeltez de los diferentes elementos y la función que representa su comportamiento obtenida
con la herramienta CurveF itting del software Matlab. Se utiliza una ecuación racional cuyo
denominador es un polinomio de tercer grado y su numerador un polinomio de primer grado.
La ecuación obtenida es:

L
−0.04206 · + 0.4084
ed D
= 3 2 (4.1)
r2

L L L
− 6.399 · + 13.81 · − 8.74
D D D
Se considera que para relaciones de esbeltez menores a 2, la relación de excentricidad se

mantiene constante en un valor de 0.25, puesto que no es razonable que este factor crezca
indefinidamente.
k=0.632
0.35
Columnas - Kowalsky
0.3 Muros - Iñiguez
Función Obtenida
0.25
0.2
ed/r2
0.15
0.1
0.05
-0.05
0 2 4 6 8 10 12 14
L/D
Con los resultados obtenidos de los análisis que se realizaron en este capı́tulo, se modifica la
curva deformación de pandeo vs. relación de esbeltez de manera que se considere el efecto
que las excentricidades provocan en la capacidad resistente de las barras de acero pertene-
cientes al refuerzo longitudinal en elementos de hormigón armado y represente de la manera

más aproximada a la realidad los valores de deformación crı́ticos de pandeo de dichas ba-
rras. Por lo tanto, utilizando la ecuación 2.55 y considerando la relación de excentricidad con
respecto a la relación de esbeltez, se obtiene una nueva curva deformación de pandeo vs. rela-
ción de esbeltez y se compara con las deformaciones de muros y columnas para su validación.
0.07
Ensayos de muros rectangulares bajo cargas cíclicas
Ensayos de columnas circulares bajo cargas cíclicas
0.06
Curva Obtenida
0.05
0.04
0.03
0.02
0.01
0
2 4 6 8 10 12 14
L/D
Figura 4.7: Gráfica deformación crı́tica de pandeo vs. relación de esbeltez
Como se puede observar, la curva obtenida 4.7 tiene una buena aproximación con las defor-
maciones crı́ticas de pandeo tanto para las columnas ensayadas por Kowalsky [5] como para
los muros obtenidos en [4]. La tabla 4.4 muestra las relaciones entre las deformaciones de
pandeo experimentales y calculadas, donde la relación promedio es de 1.05 con un coeficiente
de variación igual a 11 %, por lo tanto responde de manera correcta al comportamiento de
las barras de acero pertenecientes al refuerzo longitudinal en elementos de hormigón armado
sometidos a cargas cı́clicas.

Deformación
Deformación de Def
de pandeo Def calculado
pandeo
calculado
0.032 0.031 1.03
0.028 0.024 1.16
0.031 0.031 0.99
0.039 0.046 0.86
Muros
0.031 0.022 1.42
0.04 0.039 1.03
0.049 0.041 1.18
0.021 0.022 0.96
0.061 0.052 1.18
0.05 0.050 1.00
0.045 0.039 1.18
0.048 0.050 0.96
0.063 0.061 1.04
0.054 0.052 1.04
0.051 0.052 0.99
0.059 0.052 1.14
Columnas 0.056 0.052 1.08
0.058 0.052 1.12
0.052 0.052 1.00
0.052 0.056 0.93
0.050 0.056 0.90
0.061 0.052 1.19
0.058 0.052 1.12
0.055 0.052 1.06
0.046 0.052 0.89
Cuadro 4.4: Error entre deformaciones crı́ticas experimentales y calculadas
4.2. Análisis de excentricidades y validación de las cur-

vas: obtenida y propuesta
Una razón de la existencia de excentricidades en elementos de hormigón armado, puede ser
las expansiones laterales del hormigón bajo carga de compresión debido al efecto de Poisson.
El mismo que es la relación entre las deformaciones transversales y las longitudinales. Por
tal razón se realiza un análisis comparativo entre las excentricidades que se obtienen con
las correspondientes relaciones de excentricidad calculadas para los datos experimentales de
muros obtenidos en [4] y columnas ensayadas por Kowalsky [5] respectivamente y los des-
plazamientos producto de la expansión del hormigón en cada elemento.
Con el fin de obtener las excentricidades de cada uno de los elementos, se calcula el ra-
Z 4R
dio del núcleo , que para una sección circular maciza es igual a . Obtenido el radio del
A 3π
núcleo, con ayuda de la ecuación 2.57, se despeja la excentricidad.

Para la obtención de los desplazamientos laterales, el coeficiente de Poisson relativo a las

deformaciones elásticas puede tomarse entre 0.1 y 0.2 [33]. Por lo tanto, se procede a obtener
las deformaciones laterales de las columnas y los muros con un coeficiente de Poisson de 0.15.
La tabla 4.5 contiene las excentricidades calculadas junto a los desplazamientos laterales
del hormigón δx . Como se aprecia en la tabla se han ignorado los datos de los muros y co-
lumnas con esbeltez mayor e igual a 6, debido a que la relación de excentricidad de estos
puntos es cercana a 0, las excentricidades son muy pequeñas, motivo por el cual son ignoradas.
Excentricidad f’c
Esbeltez δx
(cm) (kg/cm2 )
4.5 0.06 465 0.04
Muros 5.3 0.015 349.6 0.014
5 0.009 333.24 0.015
2.67 0.06 478.8 0.1
3.66 0.034 428.9 0.07
3.67 0.036 508.3 0.06
2.67 0.07 445.03 0.08
2.67 0.074 454.9 0.08
2.67 0.061 449 0.08
2.67 0.066 459.1 0.08
Columnas
2.67 0.063 464.7 0.079
2.67 0.073 454.8 0.08
2.29 0.1 442.2 0.1
2.29 0.106 414 0.11
2.67 0.06 432.4 0.08
2.67 0.063 438.7 0.08
2.67 0.068 416 0.08
Cuadro 4.5: Excentricidades calculadas y desplazamientos laterales del hormigón
Como puede apreciarse en la gráfica 4.8 las excentricidades calculadas son semejantes a los
desplazamientos laterales del hormigón de los muros y las columnas, pues la relación entre
ellas se aproxima a 1,verificándose ası́ que las excentricidades tienen una buena relación con
la realidad, por lo tanto se valida la curva deformaciones de pandeo vs. relación de esbeltez,
tanto obtenida con el estudio realizado como la propuesta en el presente trabajo.
Finalmente, con el fin de mejorar el manejo y la utilidad del presente trabajo, se propone una
aproximación de la curva obtenida mediante una función hiperbólica simple y una función
lineal como se muestra en la gráfica 4.9.
La curva propuesta responde a la función lineal y = 0.051 para esbelteces menores a 4 y a

la función parabólica y = 0.191 · x−0.614 − 0.031 desde esbelteces iguales y mayores a 4.

Desplazamiento lateral del hormigón / Excentricidad

3
-1
-2
1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6
Relación de esbeltez L/D
Figura 4.8: Relación entre los desplazamientos laterales del hormigón y las excentricidades
0.07
Curva obtenida
Curva Propuesta
0.06 y= 0.051
0.05
y= 0.191x-0.614 - 0.031
0.04
0.03
0.02
0.01
0
2 4 6 8 10 12 14
L/D
Figura 4.9: Curva obtenida y curva propuesta

Capı́tulo 5
CONCLUSIONES Y
RECOMENDACIONES
CONCLUSIONES
La relación propuesta en el presente trabajo de titulación, entre la deformación crı́ti-

ca de pandeo y la esbeltez de barras de acero pertenecientes al refuerzo longitudinal
de estructuras de hormigón armado, tiene una aproximación aceptable con el com-
portamiento real de dichas barras ante cargas cı́clicas, debido a que la relación entre
deformaciones experimentales y calculadas tiene una media de 1.05 y un coeficiente de
variación de 11 %. La relación antes mencionada se evidencia en la correlación con los
datos experimentales obtenidos en [4] y [5] en los ensayos en muros y columnas como
muestra la gráfica 4.7, representa una buena fuente de información y ayuda para un
diseño sismorresistente eficiente.
Del Capı́tulo 3, el modelo de Mander de la curva esfuerzo vs. deformación del acero
bajo carga de compresión es el que más se acerca a la curva experimental obtenida
por Botero [3], por lo que es el modelo que mejor representa la realidad; razón por la
cual, se escoge los modelos de Mander bajo carga monotónica y bajo cargas cı́clicas
como base para la obtención de un nuevo modelo que se basa solamente en el módulo
elástico, el esfuerzo de fluencia del acero y en las deformaciones máximas de cada ciclo.
Se puede observar en las gráficas 3.7, 3.8 y 3.9 que el modelo propuesto tiene una bue-
na aproximación con las gráficas experimentales obtenidas por Botero [3], por lo que
podemos decir que el modelo da resultados aceptables. Cabe recalcar que el modelo de
Mander de curva de esfuerzo vs. deformación bajo carga monotónica hace referencia
a carga de tracción, y debido a que las caracterı́sticas de las curvas a tracción son
ligeramente diferentes a compresión, se debe usar las ecuaciones 2.17 y 2.18 de Dodd
y Restrepo [14] para relacionarlas.
Los parámetros de rigidez lateral α y rotacional β, pretenden representar condiciones

de apoyo de una barra de refuerzo longitudinal de manera real. Es un modelo de mucha
utilidad, debido a que con el diámetro de la barra de refuerzo longitudinal y el diáme-
tro de los estribos, se obtiene una buena aproximación de la condición de apoyo de la
barra entre estribos. Sin embargo, cabe recalcar que en el desarrollo de este modelo
no se consideró la acción del hormigón, por lo que este trabajo puede ser base para

futuras investigaciones.
Se demostró que la utilización del doble módulo reducido da esfuerzos crı́ticos de pan-
deo que se aproximan a la realidad, debido a que la curva de esfuerzo crı́tico de pandeo
vs. relación de esbeltez L/D se ajustan a los ensayos realizados por Botero.
El factor de longitud efectiva K optimizado para los ensayos de barras aisladas de Bote-
ro, es 0.755, valor que se justifica por las condiciones del ensayo, debido a la existencia
de un pequeño giro entre las probetas y el dispositivo de agarre, lo que produjo un giro
en los extremos de las probetas.
Se considera que los ensayos de los muros obtenidos en [4] con relación de esbeltez
igual o mayor que 6 presentan una excentricidad despreciable, consecuencia de que
las relaciones de excentricidad de estos ensayos son muy cercanas a 0. Al optimizar el
factor de longitud efectiva de la curva deformación de pandeo vs. relación de esbeltez
sin excentricidad con dichos ensayos, se obtiene k de 0.632, esto es producto de que en
ensayos experimentales es difı́cil conseguir un empotramiento perfecto.
Con el fin de validar las excentricidades calculadas, y por ende la curva obtenida, se
calcula la expansión lateral del hormigón bajo carga de compresión debido al efecto
Poisson, puesto que es una de las causas de la presencia de excentricidades en las
estructuras de hormigón armado. Se observa en la gráfica 4.8 que las excentricidades
calculadas y los desplazamientos del hormigón tienen una relación aceptable.
RECOMENDACIONES
Se recomienda calibrar el modelo con ensayos realizados con barras de acero usadas en
Ecuador. Los ensayos utilizados se obtuvieron de la bibliografı́a con material de México
y Alaska, debido a que no se realizaron ensayos propios por la falta de equipo.
Para futuras lı́neas de investigación, se recomienda considerar los esfuerzos que ejerce
el hormigón sobre el refuerzo longitudinal con el objeto de obtener un estudio más
completo y detallado.

Bibliografı́a
[1] Luis Enrique Garcı́a. Dinámica estructural aplicada al diseño sı́smico. Universidad de
los Andes, Colombia, 1998.
[2] T. Paulay and M. J. N Priestly. Seismic design of reinforced concrete and masonry
buildings. 1992.
[3] J. C. Botero and M. E. Rodriguez. Pandeo en barras de refuerzo sometidas a cargas

monotónicas y cı́clicas reversibles. Master’s thesis, Universidad Nacional Autónoma de
México, 1996.
[4] M. A. Iñiguez. Resistencia a flexo-compresión y capacidad de deformación lateral de mu-

ros rectangulares de concreto reforzado en zonas sı́smicas. Master’s thesis, Universidad
Nacional Autónoma de México, 2015.
[5] J. C. Goodnight, Y. Feng, M. J. Kowalsky, and J. M. Nau. The effects of load history and
design variables on performance limit states of circular bridge columns. North Carolina
State University, 2015.
[6] G. A. Chang and B. Mander. Seismic energy based fatigue damage analysis of bridge
columns:Part I- evaluation of seismic capacity, 1994.
[7] C. G. Salmon and J. E. Johnson. Steel Structures: Design and Behavior. Ed. Harper
Collins Publishers, 1990.
[8] G. Monti and C. Nuti. Nonlinear cyclic behavior of reinforcing bars including buckling.
Journal of Structural Engineering,ASCE, 1992.
[9] M. A. Iñiguez and M. E. Rodriguez. Resistencia a flexo-compresión y capacidad de

deformación lateral de muros rectangulares de concreto reforzado en zonas sı́smicas.
Master’s thesis, Universidad Nacional Autónoma de México, 2015.
[10] J. B. Mander, M. J. N. Priestley, and R. Park. Seismic Design of Bridge Piers, 1984.
[11] R. Park and T. Paulay. Estructuras de concreto reforzado. Limusa, México, DF, 1986.
[12] H. Ahmad and S. Shah. Structural properties of high strength concrete and its impli-
cations for precast prestressed concrete. Journal Prestressed Concrete Institute, 1985.
[13] R. L. Bonett. Vulnerabilidad y riesgo sı́smico de edificios. Aplicación a entornos urbanos

en zonas de amenaza alta y moderada. PhD thesis, Universidad Politécnica de Cataluña,
España, Barcelona, 2003.

[14] L. L. Dodd and J. I. Restrepo. Model for predicting cyclic behavior of reinforcing steel.
ASCE, 121(3), 1995.
[15] T. Andriono and R. Park. Seismic design considerations of the properties of new zea-
land manufactured steel reinforcing bars. New Zealand National Society for Earthquake
Engineering, 19, 1986.
[16] J. C. Botero and M. E. Rodriguez. Criterios de diseño sı́smico de estructuras de concreto

reforzado considerando las propiedades mecánicas de aceros de refuerzo producidos en
México, 1994.
[17] E. P. Popov. Introduction to the Mechanics of Solids. Prentice Hall Inc., 1968.
[18] M. Menegotto and P. E. Pinto. Method of Analysis for Cyclically Loaded Reinforced
Concrete Plane Frames Including Changes in Geometry and Non-elastic Behaviour of
Elements Under Combined Normal Force and Bending. 1973.
[19] C. Lacaze. Estudio y Modelamiento del impacto del pandeo en la fatiga de bajos ciclos
en barras longitudinales para hormigón armado, 2009.
[20] D. C. Kent and R. Park. Inelastic behavior of reinforced concrete members with cyclic
loading. Master’s thesis, University of Canterbury, 1969.
[21] W. Ramberg and W. R. Osgood. Description of stress-strain curves by three parameters,

1943.
[22] Timoshenko and Gere. Theory of Elastic Stability. Ed. McGraw Hill International Book
Company, 1970.
[23] O. Ruiz. Diseño y comportamiento de miembros estructurales de acero por el método

LRFD, 2002.
[24] D. Bernal. Ductilidad y rotaciones plásticas en secciones crı́ticas de marcos de concreto

reforzado, 1982.
[25] S. Timoshenko. Strength of materials, volume 210. van Nostrand New York, 1955.
[26] N. Hernández, J. Medina, and J. C. Seck Tuoh. Introducción a Matlab para Resolver
Problemas de Ingenierı́a Aplicando Algoritmos Genéticos, 2012.
[27] A. E. Martı́nez, A. E. Puertas, and R. Gallego. Inestabilidad de barras comprimidas y

pandeo global de estructuras, 2015.
[28] A. E. Martı́nes, M. E. Puertas, and R. Gallego. Inestabilidad de estructuras. Inestabi-

lidad de barras comprimidas y pandeo global de estructuras.
[29] B. G. Johnston. Column buckling theory: Historic highlights. Journal of Structural

Engineering,ASCE, 1983.
[30] Dazio. Alessandro, Beyer. Katrin, and Bachmann. Hugo. Quasi-static cyclic tests and
plastic hinge analysis of rc structural walls. Engineering Structures, 31(7):1556–1571,
2009.

[31] Thomsen. John H and Wallace. John W. Displacement based design of reinforced concre-
te structural walls: an experimental investigation of walls with rectangular and t-shaped
cross-sections: a dissertation. PhD thesis, Clarkson University Postdam, 1995.
[32] Villalobos Fernandez. Enrique Javier. Seismic response of structural walls with geome-
tric and reinforcement discontinuities. PhD thesis, PURDUE UNIVERSITY, 2014.
[33] P. J. Montoya, A. G. Meseguer, and F. M. Cabré. Hormigón Armado. Gustavo Gili,,

1991.

Apéndice A
Codificación del programa Matlab del

modelo propuesto del comportamiento
del acero bajo cargas cı́clicas
%Entrada de datos
%Se pide el ingreso de las caracterı́sticas básicas del material del cual se desea obtener la
curva esfuerzo vs. deformación..
Es=input(’Ingrese el Es: ’);
Fy=input(’Ingrese el Fy: ’);
%Se coloca los datos correspondientes a las caracterı́sticas del material supuestos en el mo-
delo.
fsu=10800;
ey=(Fy/Es);
esh=ey;
esu=0.1946;
F yc = −F y ∗ ((1 + (ey))2 );
fsuc=-11333;
eyc=-((ey)/(1+ey));
eshc=-((esh)/(1+esh));
esuc=-0.1629;
P=3.474;
Pc=2.092;
R=0.58;
Rt=0.5472;
%Se pide el ingreso de las deformaciones que se desean alcanzar en cada ciclo, comenzando
por la compresión.
ec=input(’Ingrese la primera deformación: ’);
et=input(’Ingrese la segunda deformación: ’);
ec1=input(’Ingrese la tercera deformación: ’);
et1=input(’Ingrese la cuarta deformación: ’);
ec2=input(’Ingrese la quinta deformación: ’);
et2=input(’Ingrese la sexta deformación: ’);
ec3=input(’Ingrese la séptima deformación: ’); i=1;
%Se prosigue a calcular el valor de esfuerzos para deformaciones en intervalos de 0.0001

hasta las deformaciones ingresadas por el usuario, para cada una de las zonas.
%Zona elástica lineal.
for a=0:-0.0001:eyc
defc(i)=a;
i=i+1;
end
%Zona de fluencia.
if eshc < ec
for a=eyc:-0.0001:eshc
defc(i)=a;
i=i+1;
else
defc(i)=a;
i=i+1;
end
%Zona de endurecimiento por deformación.
for a=eshc:-0.0001:ec
defc(i)=a;
i=i+1;
end
%Zona de descarga.
ec0=(ec-((fsc(i-1))/(Es)));
for a=ec:0.0001:ec0
defc(i)=a;
i=i+1;
end
%Zona de tracción tomando en cuenta el efecto Bauschinger.
for a=ec0:0.0001:et
dE=abs((esu-ec0))/2;
Eo=Es*(1-dE);
c = (1 + (Eo/Esec)( Rt − (ec0 ∗ 4.4715)))( − 1/(Rt − (ec0 ∗ 4.4715)));
(Rt-(ec0*4.4715))
Q=((Esec/Eo)-c)/(1-c);
f sc(i) = Eo ∗ (a − ec0) ∗ (Q + ((1 − Q)/((1 + (Eo ∗ ((a − ec0)/
(f su + (ec0 ∗ 2.4207 ∗ f su))))( Rt − (ec0 ∗ 4.4715)))( 1/(Rt − (ec0 ∗ 4.4715))))));

defc(i)=a;
i=i+1;
end
%Zona de descarga.
et0=(et-((fsc(i-1))/(Es)));
for a=et:-0.0001:et0
defc(i)=a;
i=i+1;
end
%Zona de compresión tomando en cuenta el efecto Bauschinger.
for a=et0:-0.0001:ec1
dE=abs((esuc+et0))/2;
Eo=Es*(1-dE);
c = (1 + (Eo/Esec)R )( − 1/R);
f sc(i) = Eo ∗ (a − et0) ∗ (Q + ((1 − Q)/((1 + (Eo ∗ ((a − et0)/
((f suc − (f suc ∗ et0 ∗ 0.9529)))))R )( 1/R))));
defc(i)=a;
i=i+1;
end
%Zona de descarga.
ec01=(ec1-((fsc(i-1))/(Es)));
for a=ec1:0.0001:ec01
defc(i)=a;
i=i+1;
end
for a=ec01:0.0001:et1
Eo=Es*(1-dE);
defc(i)=a;
i=i+1;
end

%Zona de descarga.
et01=(et1-((fsc(i-1))/(Es)));
for a=et1:-0.0001:et01
defc(i)=a;
i=i+1;
end
for a=et01:-0.0001:ec2
Eo=Es*(1-dE);
c = (1 + (Eo/Esec)R )( − 1/R);
f sc(i) = Eo ∗ (a − et01) ∗ (Q + ((1 − Q)/((1 + (Eo ∗ ((a − et01)/
((f suc − (f suc ∗ et01 ∗ 0.9529)))))R )( 1/R))));
defc(i)=a;
i=i+1;
end
%Zona de descarga.
ec02=(ec2-((fsc(i-1))/(Es)));
for a=ec2:0.0001:ec02
defc(i)=a;
i=i+1;
end
for a=ec02:0.0001:et2
Eo=Es*(1-dE);
defc(i)=a;
i=i+1;
end
%Zona de descarga.

et02=(et2-((fsc(i-1))/(Es)));
for a=et2:-0.0001:et02
defc(i)=a;
i=i+1;
end
for a=et02:-0.0001:ec3
Eo=Es*(1-dE);
c = (1 + (Eo/Esec)R )( − 1/R);
f sc(i) = Eo ∗ (a − et02) ∗ (Q + ((1 − Q)/
((1 + (Eo ∗ ((a − et02)/((f suc − (f suc ∗ et02 ∗ 0.9529)))))R )( 1/R))));
defc(i)=a;
i=i+1;
end
%Se grafica como resultado la curva obtenida.

plot(defc,fsc,’–r’)
title(’Curva esfuerzo vs deformación bajo cargas cı́clicas’)
xlabel(’Deformación’)
ylabel(’Esf uerzokg/cm2 ’)
grid
legend(’Modelo propuesto’)

Apéndice B
Codificación del programa Matlab de

las curvas esfuerzo crı́tico de pandeo
vs. relación de esbeltez y deformación
crı́tica de pandeo vs. relación de
esbeltez
%Entrada de datos
%Como primer punto se coloca los datos correspondientes a las caracterı́sticas del material
para poder elaborar la curva esfuerzo vs. deformación.
H2(1)=1000000000;
H1(1)=1000000000;
H(1)=1000000000;
Es=2100000;
i=1;
j=1;
l=1;
n=1;
F(1)=1;
Fy=4680;
ey=(Fy/Es);
esh=ey;
esu=0.1649;
fsu=10800;
eshc=-0.0026;
esuc=-0.1629;
fsuc=-11333;
F yc = −F y ∗ ((1 + (ey))2 );
eyc=((ey)/(1+ey));
ec=esuc;
fsc=fsuc;
e0=0.0334;
Gcr2(1)=0;

Gcr1(1)=0;
Gcr(1)=0;
%Se procede a graficar la curva esfuerzo vs. deformación para intervalos de deformación
de 0.0001.
for e=(e0-0.0001):-0.0001:ec
s=0.0001:0.0001:e0+abs(ec);
Esec= (fsc)/(ec-e0);
dE=abs((ec-e0))/2;
R=0.58;
Eo=Es*(1-3*dE);
c = (1 + (Eo/Esec)R )( − 1/R);
f s(i) = Eo ∗ (e − e0) ∗ (Q + ((1 − Q)/((1 + (Eo ∗ ((e − e0)/(f sc)))R )( 1/R))));
a(i)=e;
ac(j)=e;
Gc(j)=fs(i);
j=j+1;
i=i+1;
end
% A continuación se obtiene el Gcr que hace referencia a los esfuerzos crı́ticos correspon-
dientes a un rango de esbelteces hasta 13.
% Se obtienen los esfuerzos crı́ticos de pandeo con un factor de longitud efectiva k de 0.5,
0.75 y 1.
for esb=0.1:0.1:13
Y(l)=esb;
for k=1:(size(Gc,2)-1)
Et(k)=((Gc(k+1)-Gc(k))/(ac(k+1)-ac(k)));
Er(k) = (4 ∗ Es ∗ Et(k))/((((Es)( 1/2)) + ((Et(k))( 1/2)))2 );
Gcr(k + 1) = −((((pi)2 ) ∗ (Er(k)))/(16 ∗ ((0.5 ∗ esb)2 )));
Gcr1(k + 1) = −((((pi)2 ) ∗ (Er(k)))/(16 ∗ ((0.75 ∗ esb)2 )));
Gcr2(k + 1) = −((((pi)2 ) ∗ (Er(k)))/(16 ∗ ((1 ∗ esb)2 )));
H2(k+1)=abs(Gcr2(k+1)-Gc(k+1));
H1(k+1)=abs(Gcr1(k+1)-Gc(k+1));
H(k+1)=abs(Gcr(k+1)-Gc(k+1));
end
posicion2=find(H2==min(H2));
posicion1=find(H1==min(H1));
posicion=find(H==min(H));
Z2(l)=-Gc(posicion2);
Z1(l)=-Gc(posicion1);
Z(l)=-Gc(posicion);
def2(l)=df(Z2(l));

def1(l)=df(Z1(l));
def(l)=df(Z(l));
Erp(l)=Er(posicion2-1);
l=l+1;
end
% Se introducen los datos experimentales de esfuerzos crı́ticos de pandeo y deformaciones

crı́ticas de pandeo de los ensayos de Botero.
CSx=[4 4 4 6 6 6 8 8];
CAx=[4 4 6 6 8 8];
CSyc=[8486 8000 7930 6900 6681 6484 5204 4750];
CAyc=[7300 7250 6030 5956 5186 5088];
CSydc=[0.07216 0.05 0.05 0.039 0.031 0.0211 0.0123 0.0097];
CAydc=[0.069 0.069 0.0472 0.0459 0.0236 0.019];
% Por último, se grafica los esfuerzos vs. la relación de esbeltez y las deformaciones vs.
la relación de esbeltez.
subplot(2,1,1),plot(CSx,CSydc,’+r’,CAx,CAydc,’s’,Y,def,Y,def1,Y,def2,’b’);
subplot(2,1,2),plot(CSx,CSyc,’+r’,CAx,CAyc,’s’,Y,Z,Y,Z1,Y,Z2,’b’);

Apéndice C
Codificación del programa Matlab de

la obtención de K mediante
optimización
% Programa para la optimización utilizando loa herramienta fminsearch.

% Las matrices x,x1,...,xm son las esbelteces de los datos a optimizar.
% Las matrices f,f1,...,fm son las deformaciones de pande de los datos a optimizar.
% Las matrices D,Dl1,...,Dlm son los diámetros del refuerzo longitudinal.
% Las matrices De,De1,...,Dem son los diámetro del refuerzo transversal.
x=[];
f=[];
Dl=[];
De=[];
x1=[];
f1=[];
Dl1=[];
De1=[];
.
.
.
xm=[];
fm=[];
Dlm=[];
Dem=[];
% Uso de la herramienta fminsearch para la optimización de m datos en donde el valor a

optimizar es el factor n, minimizando el error.
options = optimset(’TolFun’,1e-15,’Display’,’iter’,’PlotFcns’,@optimplotfval);
[n,fval,exitflag,output] = fminsearch(@(n) Fe(n,x,f,Dl,De,x1,f1,Dl1,De1,...,xm,fm,Dlm,Dem),[1],options)

% Función Fe el cual calcula el error de cada iteración comparando las deformaciones cal-
culadas con las deformaciones de entrada.
function z = Fe(n,x,f,Dl,De,x1,f1,Dl1,De1,...,xm,fm,Dlm,Dem)
yca =(Fs(n,x,f,Dl,De));
yca1 =(Fs(n,x1,f1,Dl1,De1));
.
.
.
ycam =(Fs(n,xm,fm,Dlm,Dem));
dca=Df(yca);
dca1=Df(yca1);
.
.
.
dcam=Df(ycam);
z1=0;
z2=0;
.
.
.
zm=0;
for k=1:size(yca,2)
z1=z1+s(k);
end
for k=1:size(yca1,2)
z1=z1+s1(k);
end
.
.
.
for k=1:size(ycam,2)
zm=zm+sm(k);
end
z=z4+z5+...+zm;
% Función Fs que calcula el esfuerzo de pandeo para cada dato con la variación de la ecua-
ción de Euler para una barra con diferentes restricciones en los apoyos.

function z = Fs(n,x,f,Dl,De)
Er=Erm(f);
for k=1:size(f,2)
f le(n, x(k), Dl(k), De(k), Er(k))) ∗ x(k)).2 ));
end
% Función Erm que calcula el doble módulo reducido con la deformación de cada dato, con
la curva propuesta con las caracterı́sticas de la barra de esbeltez 2.5 ensayada por Botero.
function Erm = Erm(y)
Es=2000000;
Fy=4680;
fsu=10800;
ey=(Fy/Es);
esh=ey;
esu=0.1649;
P=3.474;
Fyc=−F y ∗ ((1 + (ey))2 );
fsuc=-11333;
eyc=-((ey)/(1+ey));
esuc=-0.1963;
et0=0;
ec1=-0.1963;
Pc=2.092;
Esec= (fsuc)/(esuc-et0);
R=0.58;
Eo=Es*(1-dE);
c=(1 + (Eo/Esec)R )( − 1/R);
for i=1:size(y,2)
(f suc)))R )( 1/R)))));
f 2 = −(Eo ∗ ((−y(i) + 0.00001) − et0) ∗ (Q + ((1 − Q)/
((1 + (Eo ∗ (((−y(i) + 0.00001) − et0)/(f suc)))R )( 1/R)))));
Et = (f 1 − f 2)/(y(i) − (y(i) − 0.00001));
Er1(i) = (4 ∗ Es ∗ Et)/((((Es)( 1/2)) + ((Et)( 1/2)))2 );
end
Erm=Er1

% Función K que calcula el factor de longitud efectiva para cada dato con la variación de
n, el cual es el factor que reduce las restricciones de giro .
function k=K(n,x,Dl,De,Er)
L=Dl*x;
I=(pi ∗ Dl4 )/64;
k1=(2000000 ∗ pi ∗ De2 )/(4 ∗ L);
k2=k1;
kr1=((4*Er*I)/L)*n;
kr2=kr1;
syms u
U = vpasolve((((((−k1 + k2) ∗ (u6 /(L6 ))) + ((kr1 ∗ kr2) ∗ (k1 + k2)
+(k1 ∗ k2 ∗ L)) ∗ (u4 /(L4 )) + (k1 ∗ k2 ∗ (kr1 + kr2 − kr1 ∗ kr2 ∗ L)) ∗ (u2 /(L2 )))
∗sin(u) + ((k1 + k2) ∗ (kr1 + kr2) ∗ (u5 /(L5 )) − (k1 ∗ k2 ∗ L) ∗ (kr1 + kr2)
∗(u3 /(L3 )) − (2 ∗ k1 ∗ k2 ∗ kr1 ∗ kr2) ∗ (u/L)) ∗ cos(u) + (2 ∗ k1 ∗ k2 ∗ kr1 ∗ kr2)
∗(u/L))) == 0, u, [17]);
k=pi/U;
% Función Df que devuelve el valor de las deformaciones correspondientes a cada uno de los
esfuerzos calculados con la variación de la ecuación de Euler y las caracterı́sticas de la barra
de esbeltez 2.5 ensayada por Botero.
function Df = Df(y)
Es=2000000;
Fy=4680;
fsu=10800;
ey=(Fy/Es);
esh=ey;
esu=0.1649;
P=3.474;
Fyc=−F y ∗ ((1 + (ey))2 );
fsuc=-11333;
eyc=-((ey)/(1+ey));
esuc=-0.1963;
et0=0;
ec1=-0.1963;
Pc=2.092;
Esec= (fsuc)/(esuc-et0);
R=0.58;
Eo=Es*(1-dE);
c=(1 + (Eo/Esec)R )( − 1/R);

for i=1:size(y,2)
A(i) = vpasolve((y(i) + Eo ∗ (a − et0) ∗ (Q + ((1 − Q)/
((1 + (Eo ∗ ((a − et0)/(f suc)))R )( 1/R))))) == 0, a, [−10]);
df(i)=-A(i);
end


Trabajo de Titulacion Cargas Ciclicas

Cargado por

Copyright:

Formatos disponibles

Trabajo de Titulacion Cargas Ciclicas

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Trabajo de Titulacion Cargas Ciclicas

Cargado por

Copyright:

Formatos disponibles

UNIVERSIDAD DE CUENCA

ESCUELA DE INGENIERÍA CIVIL

ESTUDIO DEL COMPORTAMIENTO E INESTABILIDAD DEL

TESIS PREVIA A LA OBTENCIÓN DEL

ANDREA FERNANDA MONCAYO TELLO

MAURO ANDRÉS VALLEJO BOJORQUE

INGENIERO HOMERO FERNANDO CARRIÓN CABRERA MsC.

Durante sismos, el desprendimiento del recubrimiento del concreto en estructuras de hor-

Andrea Moncayo Tello 2

Keywords: Cyclic loads, reinforced concrete, longitudinal reinforcement, hysteretic cycles,

Andrea Moncayo Tello 3

2. Estado del arte 18

Andrea Moncayo Tello 4

3. Pandeo en barras de acero bajo cargas cı́clicas 38

4. Pandeo en barras de acero del refuerzo longitudinal de estructuras de

A. Codificación del programa Matlab del modelo propuesto del comporta-

B. Codificación del programa Matlab de las curvas esfuerzo crı́tico de pandeo

C. Codificación del programa Matlab de la obtención de K mediante optimi-

Andrea Moncayo Tello 5

2.1. Esfuerzo vs. deformación del acero bajo carga de tracción. . . . . . . . . . . 18

3.1. Curva esfuerzo vs. deformación bajo carga a tracción . . . . . . . . . . . . . 39

Andrea Moncayo Tello 6

4.6. Gráfica relación de excentricidad vs. relación de esbeltez, con k=0.632 . . . . 71

Andrea Moncayo Tello 7

3.1. Parámetros de las probetas ensayadas por Botero [3]. . . . . . . . . . . . . . 39

4.1. Resultados de ensayos experimentales en muros y columnas bajo cargas cı́cli-

Andrea Moncayo Tello 8

En la siguiente lista se encuentra la simbologı́a utilizada en el presente estudio, en el orden

Estado del Arte

εy Deformación del esfuerzo de fluencia.

Andrea Moncayo Tello 9

Mmax Momento máximo resistente a flexión.

Pandeo en barras de acero bajo cargas cı́clicas

α0 Restricción lateral al inicio de la barra.

Pandeo en barras de acero del refuerzo longitudinal de estructuras

K Factor de longitud efectiva optimizado.

Andrea Moncayo Tello 10

Andrea Moncayo Tello 11

Andrea Moncayo Tello 12

Andrea Moncayo Tello 13

Andrea Moncayo Tello 14

La falla por pandeo se ha presentado en varios terremotos como el de México en 1985,

Andrea Moncayo Tello 15

Proponer una relación entre la deformación crı́tica de pandeo y la esbeltez de barras de

Andrea Moncayo Tello 16

Posteriormente, se realiza un análisis de la existencia de excentricidades y su efecto en ensa-

Andrea Moncayo Tello 17

ESTADO DEL ARTE

2.1. Comportamiento del acero bajo carga de tensión

Andrea Moncayo Tello 18

2.1.1.1. Zona elástica lineal

Donde Et , fs y εs son el módulo elástico, el esfuerzo y la deformación en cualquier pun-

2.1.1.2. Zona de fluencia

El valor del esfuerzo en esta zona se obtiene utilizando la siguiente expresión:

Andrea Moncayo Tello 19

2.1.1.3. Zona de endurecimiento por deformación

A continuación, se presentan los distintos modelos revisados:

2.2. Modelos del Acero

Para la obtención de p mediante la ecuación 2.5, se puede utilizar la definición de loga-