Mecánica Avanzada de Estructuras
Mecánica Avanzada de Estructuras
Publicado por:
Mecánica Avanzada de Estructuras
2020
ISBN: 978-84-121101-3-5
A nuestros estudiantes
I
II
Prólogo
“Mecánica Avanzada de Estructuras” tiene el objetivo de proporcionar a los estudiantes
de grado de ingenierı́a una visión y una comprensión extendida del análisis estructural
y de sus aplicaciones, más allá del comportamiento estático y lineal de estructuras de
barras. El planteamiento del texto coincide con el de asignaturas de Mecánica de Es-
tructuras que, con diferentes nombres, desarrollan la materia de Teorı́a de Estructuras
en las diversas titulaciones de grado y master en ingenierı́a. La obra recoge la experien-
cia docente de los autores en asignaturas relacionadas con el Análisis de Estructuras en
la Escuela Técnica Superior de Ingenieros de Caminos, Canales y Puertos de Barcelona.
El libro se presenta como una extensión del cálculo estático de estructuras de barras
desarrollado en el libro “Mecánica de Estructuras”, por lo que se suponen conocimien-
tos previos de elasticidad lineal y teorı́a de estructuras, ası́ como de algebra matricial y
ecuaciones diferenciales. Los desarrollos teóricos se acompañan con abundantes ejem-
plos y figuras explicativas.
En el Capı́tulo 4 se afronta análisis estructural más allá del lı́mite elástico del material,
hasta el colapso de la estructura. El objetivo de este Capı́tulo es determinar tanto la
III
carga última de la estructura como el mecanismo de colapso. Esto implica introducir
conceptos como el momento plástico de la sección y el de rótula plástica. Se presentan
los teoremas fundamentales del análisis lı́mite ası́ como las metodologı́as derivadas de
los mismos.
En el Capı́tulo 5 se desarrolla la formulación de placas delgadas. Las placas son ele-
mentos estructurales en los que existen dos dimensiones predominantes, por lo que el
problema asociado es bidimensional. La formulación desarrollada se aplica a los dos
casos de mayor interés: placas rectangulares y placas circulares.
El Capı́tulo 6 está dedicado al cálculo dinámico. Se estudia un sistema dinámico de un
grado de libertad para introducir y estudiar las caracterı́sticas diferenciadoras del prob-
lema dinámico en relación al problema estático. Los métodos se extienden a sistemas
de múltiples grados de libertad, lo que lleva a la formulación matricial del problema
dinámico.
Los autores agradecen la contribución de los profesores Prof. Benjamı́n Suárez, cuyos
años de docencia son el origen de esta obra, Prof. Ramon Codina, por su ayuda,
comentarios y sugerencias en varios de los capı́tulos, y al Prof. Gabriel Bugeda, por
acometer la revisión del manuscrito. Asimismo, agradecemos a Raúl Giménez por la
delineación de las figuras que complementan el texto y a Irene Cervera por la realización
de la portada.
Agradecemos igualmente a la Escuela Técnica Superior de Ingenieros de Caminos,
Canales y Puertos de Barcelona y al Centro Internacional de Métodos Numéricos en
Ingenierı́a (CIMNE) por el apoyo recibido a través de sus programas de innovación
docente que han permitido la elaboración de este material.
Finalmente, agradecemos a nuestras familias el apoyo y comprensión recibidos, sin los
cuales no hubiera sido posible este trabajo.
IV
Índice
V
1.6.4.1 Sistemas de cargas puntuales . . . . . . . . . . . . . . . 30
1.6.4.2 Cargas distribuidas . . . . . . . . . . . . . . . . . . . . 32
1.7 Envolventes de esfuerzos . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
2 Cables y arcos 49
2.1 Introducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
2.2 Cables . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50
2.2.1 Antecedentes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50
2.2.2 Hipótesis de cálculo . . . . . . . . . . . . . . . . . . . . . . . . . 51
2.2.3 Cable sometido a cargas concentradas . . . . . . . . . . . . . . . 52
2.2.4 Cable sometido a cargas distribuidas . . . . . . . . . . . . . . . . 54
2.3 Arcos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65
2.3.1 Antecedentes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65
2.3.2 Arco indeformable . . . . . . . . . . . . . . . . . . . . . . . . . . 66
2.3.3 Arco deformable . . . . . . . . . . . . . . . . . . . . . . . . . . . 70
2.3.4 Arco triarticulado . . . . . . . . . . . . . . . . . . . . . . . . . . 71
3 Inestabilidad 83
3.1 Introducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 83
3.2 Niveles de aproximación geométrica. Análisis de segundo orden . . . . . 84
3.3 Teorı́a de Euler para columnas esbeltas. Piezas ideales . . . . . . . . . . 85
3.3.1 Carga crı́tica en piezas simples . . . . . . . . . . . . . . . . . . . 86
3.3.1.1 Viga articulada . . . . . . . . . . . . . . . . . . . . . . . 86
3.3.1.2 Viga empotrado-articulada . . . . . . . . . . . . . . . . 88
3.3.1.3 Viga biempotrada . . . . . . . . . . . . . . . . . . . . . 89
3.3.1.4 Viga en ménsula . . . . . . . . . . . . . . . . . . . . . . 91
3.3.1.5 Viga biarticulada con apoyo elástico . . . . . . . . . . . 92
3.3.2 Longitud efectiva de pandeo . . . . . . . . . . . . . . . . . . . . . 94
3.3.3 Verificación estructural a pandeo. Concepto de esbeltez . . . . . 96
3.4 Piezas no ideales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 102
3.4.1 Axil excéntrico . . . . . . . . . . . . . . . . . . . . . . . . . . . . 103
3.4.2 Deformación inicial . . . . . . . . . . . . . . . . . . . . . . . . . . 104
3.4.3 Pieza con carga transversal . . . . . . . . . . . . . . . . . . . . . 106
3.5 Energı́a potencial total estacionaria. Método de Rayleigh-Ritz . . . . . . 110
3.5.1 Pieza biarticulada con carga axial de compresión . . . . . . . . . 110
3.5.2 Pieza con carga transversal . . . . . . . . . . . . . . . . . . . . . 113
3.5.3 Método de Rayleigh-Ritz . . . . . . . . . . . . . . . . . . . . . . 114
3.5.4 Comportamiento real de las estructuras . . . . . . . . . . . . . . 116
VI
3.6 Formulación matricial del problema de inestabilidad . . . . . . . . . . . 117
3.6.1 Estructuras reticuladas de plano medio . . . . . . . . . . . . . . 118
3.6.1.1 Matriz elemental de rigidez . . . . . . . . . . . . . . . . 118
3.6.1.2 Matriz de rigidez aproximada. Matriz geométrica . . . 121
3.6.2 Estructuras articuladas planas . . . . . . . . . . . . . . . . . . . 123
3.6.3 Transformación de coordenadas y ensamblaje de la matriz global 124
3.6.4 Resolución del problema no lineal . . . . . . . . . . . . . . . . . . 125
3.6.4.1 Método de Picard . . . . . . . . . . . . . . . . . . . . . 125
3.6.4.2 Método de Newton-Raphson . . . . . . . . . . . . . . . 125
3.6.4.3 Criterios de Convergencia . . . . . . . . . . . . . . . . . 126
3.6.4.4 Esquema de resolución del problema no lineal . . . . . . 126
3.6.5 Análisis lineal de bifurcación . . . . . . . . . . . . . . . . . . . . 127
VII
4.6.1.1 Teorema de unicidad . . . . . . . . . . . . . . . . . . . 184
4.6.1.2 Teorema estático o del lı́mite inferior . . . . . . . . . . 185
4.6.1.3 Teorema cinemático o del lı́mite superior . . . . . . . . 185
4.6.2 Método cinemático . . . . . . . . . . . . . . . . . . . . . . . . . . 186
4.6.2.1 Viga biarticulada bajo carga puntual . . . . . . . . . . 188
4.6.2.2 Viga biarticulada bajo carga distribuida . . . . . . . . . 188
4.6.2.3 Viga biarticulada con voladizo con carga distribuida . . 189
4.6.2.4 Viga biempotrada bajo cargas puntuales asimétricas . . 191
4.6.2.5 Viga continua de dos vanos . . . . . . . . . . . . . . . . 194
4.6.3 Combinación de mecanismos . . . . . . . . . . . . . . . . . . . . 196
4.6.3.1 Pórtico pilar-dintel con carga de viento . . . . . . . . . 197
4.6.3.2 Pórtico hiperestático asimétrico . . . . . . . . . . . . . 200
VIII
6.2.1 Vibración libre no amortiguada . . . . . . . . . . . . . . . . . . . 247
6.2.2 Vibración libre amortiguada . . . . . . . . . . . . . . . . . . . . . 249
6.2.2.1 Sobreamortiguamiento . . . . . . . . . . . . . . . . . . . 251
6.2.2.2 Amortiguamiento crı́tico . . . . . . . . . . . . . . . . . 251
6.2.2.3 Inframortiguamiento . . . . . . . . . . . . . . . . . . . . 252
6.2.3 Vibración forzada armónica . . . . . . . . . . . . . . . . . . . . . 260
6.2.4 Vibración forzada no armónica . . . . . . . . . . . . . . . . . . . 267
6.2.4.1 Respuesta dinámica frente a funciones de excitación sim-
ples . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 267
6.2.4.2 Respuesta dinámica para una excitación general.
Integral de Duhamel . . . . . . . . . . . . . . . . . . . . 271
6.2.5 Espectros de respuesta para cargas impulsivas . . . . . . . . . . . 273
6.2.6 Espectros de respuesta para cargas rampa-escalón con tiempo
finito de incremento . . . . . . . . . . . . . . . . . . . . . . . . . 277
6.2.7 Espectros sı́smicos de respuesta . . . . . . . . . . . . . . . . . . . 279
6.2.7.1 Espectro de respuesta para una aceleración en la base . 279
6.2.7.2 Pseudo-espectros de respuesta sı́smica . . . . . . . . . . 281
6.2.7.3 Pseudo-espectros sintéticos y espectros de diseño . . . . 282
6.3 Sistemas con múltiples grados de libertad . . . . . . . . . . . . . . . . . 284
6.3.1 Vibración libre. Frecuencias y modos propios de vibración . . . . 284
6.3.2 Vibración libre amortiguada . . . . . . . . . . . . . . . . . . . . . 286
6.3.3 Vibración forzada . . . . . . . . . . . . . . . . . . . . . . . . . . 286
6.4 Trabajos virtuales y minimización de la energı́a potencial . . . . . . . . 287
6.4.1 Principio de Hamilton o de Mı́nima Acción . . . . . . . . . . . . 287
6.4.2 Conservación de la energı́a . . . . . . . . . . . . . . . . . . . . . . 288
6.5 Solución numérica de la ecuación de equilibrio dinámico.
Método de Newmark . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 291
6.5.1 Sistema dinámico con un grado de libertad . . . . . . . . . . . . 291
6.5.2 Sistema dinámico con múltiples grados de libertad . . . . . . . . 297
6.6 Solución modal de la ecuación de equilibrio . . . . . . . . . . . . . . . . 297
6.6.1 Frecuencias y modos propios de vibración . . . . . . . . . . . . . 298
6.6.2 Descomposición y superposición modal . . . . . . . . . . . . . . . 299
6.7 Análisis sı́smico de edificios . . . . . . . . . . . . . . . . . . . . . . . . . 301
6.7.1 Introducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 301
6.7.2 Idealización estructural . . . . . . . . . . . . . . . . . . . . . . . 302
6.7.3 Análisis modal . . . . . . . . . . . . . . . . . . . . . . . . . . . . 303
IX
6.7.3.1 Componentes modales de fuerzas dinámicas, esfuerzos
cortantes y momentos . . . . . . . . . . . . . . . . . . . 304
6.7.3.2 Valores de diseño . . . . . . . . . . . . . . . . . . . . . 305
B Tablas 331
X
XI
XII
1 Análisis estructural y seguridad
1.1 Introducción
• funcionalidad: toda estructura debe servir para aquello para lo que ha sido con-
cebida.
• seguridad: toda estructura debe soportar las cargas a las que se ve sometida
durante su vida útil.
1
2 CAPÍTULO 1. ANÁLISIS ESTRUCTURAL Y SEGURIDAD
analizadas.
Una estructura debe ser proyectada, construida y mantenida para soportar las acciones
que le afectan tanto durante el periodo de construcción como durante la vida útil
prevista en el proyecto, manteniendo su funcionalidad.
Se define la vida útil de una estructura como el periodo de tiempo a partir del final de su
ejecución durante el cual debe garantizar unos requisitos de seguridad y funcionalidad
de proyecto, además de unas condiciones estéticas aceptables. Esto requiere un plan
de seguimiento y mantenimiento durante ese periodo.
La vida útil nominal es un parámetro que debe ser fijado por la propiedad de la es-
tructura al inicio del proyecto. Este valor nunca será inferior a lo indicado en las
reglamentaciones aplicables en cada caso. En la Tabla 1.1 se detallan las vidas útiles
para algunos tipos de estructura según la Instrucción de Hormigón Estructural (EHE-
08) en su capı́tulo I “Principios generales”, artı́culo 5◦ .
1.2.2 Incertidumbres
Tabla 1.1: Vida útil nominal para diversos tipos de estructura (EHE-08)
En los métodos probabilistas se trata la solicitación (S) y la respuesta (R) como varia-
bles aleatorias, de naturaleza incierta, asignándoles distribuciones de probabilidad. Los
coeficientes de seguridad se introducen de manera objetiva para cuantificar la proba-
bilidad de fallo.
Dichos coeficientes parciales de seguridad no tienen en cuenta los posibles errores hu-
manos; éstos deben ser evitados mediante mecanismos de control de calidad sobre el
proyecto, la ejecución, el uso y el mantenimiento de la estructura.
1.4 Acciones
En España, las acciones de diseño de una estructura están definidas en los siguientes
reglamentos vigentes:
1.4. ACCIONES 5
Todas estas normativas están adaptadas a las normas desarrolladas en el ámbito europeo
(Eurocódigos).
En las normativas vigentes las acciones se clasifican de diversos modos, con especial
relevancia a su variación en el tiempo:
- Permanentes (G): Son las acciones que actúan en todo momento y son constantes
en posición y magnitud. Habitualmente comprenden el peso propio y las cargas
muertas.
- Permanentes de valor no constante (G*): Son las acciones que actúan en todo
momento pero cuya magnitud no es constante. Se consideran en esta categorı́a
fuerzas de pretensado, asientos, acciones reológicas, etc.
- Variables (Q): Son las acciones externas que pueden actuar o no sobre la estruc-
tura. Tı́picamente, se consideran como acciones variables las sobrecargas de uso,
las acciones climáticas como viento, nieve, etc.
A continuación se describen con detalle las diferentes acciones de tipo variable o acci-
dental junto con las consideraciones necesarias para su aplicación.
6 CAPÍTULO 1. ANÁLISIS ESTRUCTURAL Y SEGURIDAD
La sobrecarga de uso es el peso de todo lo que puede gravitar sobre la estructura por
razón de su utilización.
En edificios, la sobrecarga de uso puede simularse con una carga uniformemente repar-
tida. El CTE define dicha carga distribuı́da según el uso del edificio a proyectar. Para
comprobaciones locales de capacidad portante se define una carga concentrada actuando
sobre cualquier punto de la zona. La Figura 1.1 muestra algunos valores caracterı́sticos
de sobrecarga de uso definidos por el CTE.
En puentes de carretera, la IAP define las sobrecargas de uso a tener en cuenta como
un tren de cargas que consta de: a) componentes verticales (fuerzas gravitatorias) y b)
componentes horizontales (fuerzas de frenado, arranque y fuerza centrı́fuga).
La distribución y el valor de las presiones que ejerce el viento sobre una estructura
y las fuerzas resultantes dependen de varios factores: la geometrı́a y las dimensiones
de la estructura, las caracterı́sticas y permeabilidad de sus superficie y la dirección,
intensidad y racheo del viento.
qe = qb · ce · cp (1.1)
donde qb es la presión dinámica del viento, que de forma simplificada puede adop-
tarse igual a 0,5 kN/m2 , ce es el coeficiente de exposición, variable con la altura del
punto considerado, y cp es el coeficiente eólico o de presión, dependiente de la forma y
orientación de la superficie respecto al viento (Figura 1.2).
En puentes, la acción del viento puede asimilarse a una carga estática, salvo en estruc-
turas muy flexibles en las que puede originarse fenómenos vibratorios importantes. El
8 CAPÍTULO 1. ANÁLISIS ESTRUCTURAL Y SEGURIDAD
empuje producido por el viento se calcula por separado para cada elemento del puente,
según la expresión (IAP):
1 2
F = CD · A · ρ VC (1.2)
2
Figura 1.3: Coeficientes de arrastre CD para las secciones más usuales (IAP)
La intensidad y distribución de la carga de nieve depende del clima del lugar, del tipo
de precipitación y del relieve del entorno. En el caso de edificios depende de la forma
del edificio o de la cubierta, de los efectos del viento y de los intercambios térmicos en
los paramentos exteriores.
La construcción debe resistir la acción horizontal del sismo en todas las direcciones, lo
que obliga a analizarlo en más de una dirección. Por lo general, basta hacerlo en dos
direcciones ortogonales en planta.
El Documento Básico de Seguridad en caso de Incendios (DB-SI) del CTE fija los
criterios a seguir en caso de incendio en edificios. El aumento de temperatura que se
produce en el incendio de un edificio afecta a su estructura de dos formas diferentes.
Por un lado, los materiales sufren una modificación importante de sus propiedades
mecánicas. Por otro, las deformaciones que sufren los elementos estructurales debido
al incremento de temperatura dan lugar, en general, a tensiones que se suman a las
debidas a otras acciones.
En caso de incendio, las acciones a considerar son las mismas acciones permanentes y
variables que en el cálculo de situación persistente, si es probable que actúen durante el
incendio. Además, si para el cálculo de la resistencia al fuego se emplean los métodos
indicados en el DB-SI, puede tomarse como efecto de la acción del incendio únicamente
el derivado del efecto de la temperatura en la resistencia del elemento estructural.
1.4. ACCIONES 11
En edificios, según el CTE, los valores de cálculo de las fuerzas estáticas equivalentes
debidas al impacto de vehı́culos de hasta 30 kN de peso total son de 50 kN en la
dirección paralela a la vı́a y de 25 kN en la dirección perpendicular, no actuando
simultáneamente.
En el caso de puentes que crucen cursos de agua navegables, es necesario tener en cuenta
la posibilidad de la colisión de una embarcación contra los elementos de la estructura.
Las normativas vigentes proporcionan el valor caracterı́stico (Sk ) de una acción, que
es su valor de referencia a efectos de proyecto. El valor caracterı́stico puede venir
determinado por un valor medio, un valor nominal o, en los casos en que se fije mediante
criterios estadı́sticos, por un valor correspondiente a una determinada probabilidad de
no ser superado durante un perı́odo de referencia, que tiene en cuenta la vida útil de
la estructura y la duración de la acción.
El valor de cálculo (Sd ) de la acción se obtiene a partir del valor representativo modi-
ficado por un coeficiente de ponderación (γi ) en función de la naturaleza de la acción
y de su condición de favorable o desfavorable para el estado lı́mite en estudio.
1.5. MÉTODO DE LOS ESTADOS LÍMITE 13
Situaciones accidentales
P P
j≥1 γG,j Gk,j + γa Ak + γQ,1 ψ1,1 Qk,1 + i>1 γQ,i ψ2,i Qk,i
Las normativas establecen las combinaciones de cargas que han de considerarse para
1.5. MÉTODO DE LOS ESTADOS LÍMITE 15
los estados lı́mite último y de servicio. En el caso de los estados lı́mite últimos las
combinaciones a considerar se muestran en la Tabla 1.6:
Situaciones frecuentes
P P
j≥1 γG,j Gk,j + γQ,1 ψ1,1 Qk,1 + i>1 γQ,i ψ2,i Qk,i
Situaciones cuasipermanentes
P P
j≥1 γG,j Gk,j + i≥1 γQ,i ψ2,i Qk,i
Sea una viga continua formada por un vano de longitud l y dos voladizos de distinta
longitud li , ld . Supóngase por simplicidad que ld > li .
La viga está sometida a una carga permanente g y dos cargas variables, una sobrecarga
de uso distribuida qsc y una sobrecarga de uso puntual Q, que pueden actuar o no en
cualquier parte de la estructura.
16 CAPÍTULO 1. ANÁLISIS ESTRUCTURAL Y SEGURIDAD
(li + l)2
Mest = γGf av g
2
ld2 l2
Mdes = γGdes g + γQ qsc d + γQ Qld
2 2
Sea una viga contnua formada por dos vanos iguales de longitud l, sometida a una
carga distribuida permanente g y una carga distribuida variable q que puede actuar o
no en cualquier parte de la viga.
Determinar las combinaciones de carga más desfavorables para los estados lı́mite últimos
de rotura a momentos flector positivo y negativo, y para el estado lı́mite de servicio de
deformación.
Datos: g (kN/m), q (kN/m)
ELU rotura: γGf av = 1, 0, γGdes = 1, 35, γQ = 1, 50
Para determinar la situación más desfavorable se notan como q1 y q2 las cargas actuantes
en cada uno de los vanos.
Comprobación a Estado lı́mite Último de rotura
La distribución de las cargas más desfavorable para los máximos momentos flectores
positivo y negativo no es la misma, por lo que es necesario identificar dicha distribución
para cada caso.
La ley de momentos flectores en el primer vano en función de q1 y q2 es
x
M12 (x) = [q1 (7l − x) − q2 l] 0≤x≤l
16
La distribución de momentos flectores en el segundo vano se comporta de forma simétrica
al primero, al ser de la misma longitud, intercambiando las cargas y la coordenada x
18 CAPÍTULO 1. ANÁLISIS ESTRUCTURAL Y SEGURIDAD
q1 = γGdes g + γQ q q2 = γG∗ g
es decir
q1 = q2 = γGdes g + γQ qsc
La distribución de cargas más desfavorable para el momento flector negativo es, por
tanto, la mostrada en la Figura 1.9.
La flecha máxima en función de q1 y q2 se localiza en las cercanı́as del punto medio del
vano 12. El cálculo de esta flecha y su ubicación exacta es complejo. Pos simplicidad
de cálculo, se considera la flecha en el centro luz (CL) del vano
max,CL l4
f12 = (7q1 − 3q2 )
764EI
Muchas estructuras se proyectan para resistir cargas móviles además de cargas fijas
tales como el peso propio. En el caso de puentes, por ejemplo, se considera la carga
de un tren o de un camión desplazándose a lo largo de toda la longitud del tablero,
Figuras 1.11a y b.
En estos casos, los efectos estructurales (reacciones, esfuerzos flector y cortante) de-
penden del punto de aplicación de la carga. El análisis de la variación de estos efectos
estructurales en función de la posición de las cargas móviles, ası́ como la determinación
de sus valores máximos se hace mediante las lı́neas de influencia.
Figura 1.13: Ley del momento flector y lı́nea de influencia del mismo en la sección S
En las estructuras articuladas como la mostrada en la Figura 1.14, las lı́neas de in-
fluencia para los axiles en las barras de interés se obtienen considerando una carga
unitaria que se mueve a lo largo de un cordón de la estructura determinando, para
cada posición, el correspondiente valor del axil en dicha barra. Es habitual considerar
sólo cargas aplicadas en los nudos o e interpolar linealmente los valores de influencia
obtenidos.
Para la cercha de la Figura 1.15, determinar las lı́neas de influencia del axil en las
barras 1, 2 y 3 cuando la carga se supone aplicada en los nudos del cordón inferior.
Para hallar las lı́neas de influencia de los axiles se emplea el método de las secciones
cortante la estructura por un plano vertical que incluye a las barras en estudio. La
Figura 1.17 muestra el esquema de cálculo:
1.6. LÍNEAS DE INFLUENCIA 23
Cuando la carga unitaria está en la parte izquierda del plano de corte, se tiene:
El signo negativo indica que el axil en la barra (1) es de compresión. En la Figura 1.18
puede verse la lı́nea de influencia correspondiente al esfuerzo axil en dicha barra.
En la barra (2) vertical, sólo hay axil cuando la carga está aplicada en el nudo A en
cuyo caso el valor del axil coincide con el de la carga aplicada (Figura 1.19).
24 CAPÍTULO 1. ANÁLISIS ESTRUCTURAL Y SEGURIDAD
El axil en la barra (3) es positivo lo que significa que la barra está traccionada. En la
Figura 1.20 se representa la lı́nea de influencia del axil de dicha barra.
Sea un viga continua en voladizo con las dimensiones indicadas en la Figura 1.21.
Determinar las lı́neas de influencia de las reacciones, del momento flector y del esfuerzo
cortante, en la sección A y en la sección B a la izquierda del apoyo 1.
1.6. LÍNEAS DE INFLUENCIA 25
x − 30
MA = 15V2 = 0, 75 (30 − x) TA = −V2 = 0 < x ≤ 15
20
x − 10
MA = 5V1 = 0, 25 (x − 10) TA = V1 = 15 < x ≤ 30
20
Las lı́neas de influencia del momento flector y del cortante en la sección A se representan
en la Figura 1.23.
MB = 10 − x TB = 1 0 < x ≤ 15
Cuando la carga está aplicada en el vano tanto el momento flector como el esfuerzo
cortante en la sección B son nulos.
MB = 0 TB = 0 15 < x ≤ 30
Las lı́neas de influencia del momento flector y del esfuerzo cortante en la sección B
situada a la izquierda del apoyo 1 se muestran en la 1.24.
1.6. LÍNEAS DE INFLUENCIA 27
Sea un viga continua de dos vanos con las dimensiones indicadas en la Figura 1.25.
Determinar las lı́neas de influencia de las reacciones en los tres apoyos, del del momento
flector en el apoyo 2 y del cortante en la sección O donde se encuentra la rótula.
Para determinar las lı́neas de influencia de las reacciones en los tres apoyos se plantean
las condiciones de equilibrio considerando la fuerza unitaria P aplicada a la izquierda
y a la derecha de la rótula O.
V = 1 (x − 18)
12 < x ≤ 18
1
30
28 CAPÍTULO 1. ANÁLISIS ESTRUCTURAL Y SEGURIDAD
V2 = 1 (18 − x)
12 < x ≤ 18
5
en el primer vano, entre los apoyos 1 y 2 mientras que cuando la carga actúa entre
el apoyo 2 y la rótula dicho momento es igual a la carga unitaria multiplicada por su
distancia a 2.
Si se conoce la lı́nea de influencia de una carga puntual, puede conocerse el efecto que
produce un sistema de varias cargas puntuales por aplicación directa del principio de
superposición (Figura 1.29) mediante la expresión:
X
YS = Pi yiS (xi ) (1.3)
i
donde Pi es el valor de cada una de las cargas puntuales que conforman el sistema de
cargas, y yiS es el valor de la lı́nea de influencia para un determinado efecto sobre la
sección S en el punto de aplicación de la carga Pi .
Es habitual que la posición más desfavorable sea la de colocar las tres cargas de manera
que la carga central se sitúe en la sección de influencia máxima mientras las dos cargas
adyacentes se sitúan a ambos lados de la central. Sin embargo, la influencia de las
cargas adyacentes depende de las pendientes de la lı́nea de influencia a ambos lados
del valor máximo. Cuando una de las dos pendientes es mayor del doble de la otra, la
disposición de cargas que maximiza la influencia del tren de cargas es la consistente en
las tres cargas situadas en el tramo de viga con mayor pendiente en la lı́nea de influencia,
con una de las cargas adyacentes coincidente con el valor de máxima influencia. Esta
relación de pendientes se da cuando la sección está situada entre el apoyo y un tercio
de la luz de la viga.
Si se conoce la lı́nea de influencia y(x) debida a una carga puntual unitaria, es posible
conocer el efecto de una carga repartida q (x) que actúa sobre una cierta longitud de la
viga. Ası́ por ejemplo, para la viga de la Figura 1.32 bajo carga repartida q(x) entre A
y B, el momento flector en la sección S puede calcularse utilizando la lı́nea de influencia
del momento flector en S, mediante la expresión:
Z B
S
Y = q (x) y (x) dx (1.4)
A
La viga continua en voladizo de la Figura 1.33 está sometida a una carga uniformemente
distribuida q = 4 kN/m y dos cargas puntuales de 20 kN cada una, separadas 2 m.
Las cargas pueden actuar sobre cualquier parte de la viga.
1.6. LÍNEAS DE INFLUENCIA 33
Determinar la disposición de las cargas para los valores extremos del momento flector
sobre la sección A. En cada caso, determinar las reacciones y los momentos flectores
máximos positivo y negativo en A.
A partir de las lı́neas de influencia de las reacciones, se calcula la lı́nea de influencia del
momento flector en A como
3
MA = 6V2 = 8 (20 − x) 0 < x ≤ 10
M = 10V = 5 (x − 4)
10 < x ≤ 20
A 1
8
34 CAPÍTULO 1. ANÁLISIS ESTRUCTURAL Y SEGURIDAD
Las áreas de influencia positiva y negativa en la lı́nea de influencia del momento flector
son
A+ = 0, 5 · 16 · 3, 75 = 30 A− = 0, 5 · 4 · (−2, 5) = −5
Analizando la lı́nea de influencia, se puede inducir la disposición de las cargas para ma-
ximizar los momentos flectores positivos y negativos suponiendo que las cargas pueden
actuar en cualquier parte de la estructura. Este planteamiento lleva a dos disposiciones
para ambos escenarios que se muestran en la Figura 1.36
Figura 1.36: Disposición de las cargas para máximos momentos flectores en la sección A
(a) flectores positivos (b) flectores negativos
forma que una de ellas se aplica encima de la sección A donde se observa la máxima
influencia mientras que la segunda carga se aplica al lado de la anterior con mayor
influencia, en esta caso a su izquierda. Ası́, las reacciones para este escenario se calculan
como
6 8
V1 = 4 · 0, 5 · 16 · 1 + 20 16 + 16 = 49, 5 kN
10 8
V2 = 4 · 0, 5 · 16 · 1 + 20 16 + 16 = 54, 5 kN
max,(+)
MA = 4 · 0, 5 · 16 · 3, 75 + 20 (3, 75 + 3, 0) = 255 kN m
En el escenario de máximo momento negativo (posición II) todas las cargas actúan
en la parte negativa de la lı́nea de influencia. La carga distribuida actúa en todo el
voladizo mientras que las cargas puntuales se colocan en el extremo, maximizando su
influencia. En este caso, las reacciones en los apoyos son
max,(−)
MA = 4 · 0, 5 · 4 · (−2, 5) + 20 (−2, 5 − 1, 25) = −95 kN m
Ası́, la Figura 1.37 muestra la envolvente del momento flector en una viga biarticulada
para cargas puntuales. La envolvente lo es de las lı́neas de influencia del momento
flector en las distintas secciones de la viga.
36 CAPÍTULO 1. ANÁLISIS ESTRUCTURAL Y SEGURIDAD
Las expresiones halladas para las reacciones dan lugar a las correspondientes lı́neas de
influencia cuya representación es trivial.
En primer lugar, se determinan las lı́neas de influencia del momento flector en las seccio-
nes situadas a intervalos de 1, 5 m siguiendo el procedimiento explicado anteriormente.
1.7. ENVOLVENTES DE ESFUERZOS 37
Estas lı́neas de influencia son bilineales, según la Figura 1.39, con un valor máximo que
corresponde a la carga unitaria aplicada en x = a.
El valor máximo de la influencia para cada sección dada por s = a se calcula como
a(l − a)
y(x)max = y(x)|x=a = aV1 =
l
El valor del momento flector en la sección s = a para el tren de cargas dado puede
calcularse utilizando la lı́nea de influencia del momento flector en dicha sección situando
el tren de cargas en la posición más desfavorable según muestra la Figura 1.41.
Los valores del máximo momento flector para las diferentes posiciones del tren de cargas
(Figura 1.41) se indican en la Tabla 1.8.
Sea una viga continua en voladizo cuyo vano apoyado tiene una longitud l = 12 m con
un voladizo de longitud lv = 4 m. La viga está sometida a la acción de su peso propio,
g = 8 kN/m, y una sobrecarga de uso uniformemente distriuda qsc = 12 kN/m.
l−x x
V1 = V2 =
l l
Se determinan las lı́neas de influencia del momento flector en una sección genérica
situada a una distancia s = a al apoyo izquierdo, para lo cual se distinguen dos casos
en función de la situación de la sección de referencia.
Los valores relevantes de esta lı́nea de influencia con 0 ≤ a ≤ 12, se dan en la propia
sección de referencia y en el extremo del voladizo
a(12 − a) max,(−) a
y(x)max,(+)
x=a = y(x)x=l+lv = −
12 3
La lı́nea de influencia del momento flector para una sección genérica situada en el vano
de la viga se representa en la Figura 1.45.
Figura 1.45: Lı́nea de influencia del momento flector para una sección en el vano
Figura 1.46: Lı́nea de influencia del momento flector para una sección en el voladizo
1
A− = − (16 − a)2
2
Para dibujar las envolventes de los momentos flectores positivos y negativos se calcula
la integral de la lı́nea de influencia del momento flector para cada sección. Las áreas
de influencia de la Tabla 1.9 se obtienen integrando las partes, positiva y negativa, de
la lı́nea de influencia del momento flector para cada sección.
Se calcula el valor del momento flector debido al peso propio, g = 8 kN/m, multipli-
cando dicho peso propio por el área (positiva+negativa) de la lı́nea de influencia del
momento flector en cada sección.
Para cada sección, se obtiene el valor del máximo momento positivo y negativo debido a
la sobrecarga qsc = 12 kN/m, multiplicando dicha carga por el área positiva y negativa,
1.7. ENVOLVENTES DE ESFUERZOS 43
Tabla 1.9: Determinación de los momentos flectores máximos en las secciones de cálculo
Finalmente, por combinación de los valores debidos al peso propio más la sobrecarga
(+)
(positivo y negativo), se obtienen los valores extremos, M max = Mpp + Msc y el
(−)
valor M min = Mpp + Msc que permiten dibujar las envolventes del máximo y mı́nimo
momento flector.
En la Tabla 1.9 se reúnen todos los pasos realizados. La Figura 1.47 representa las
envolventes de flectores positivos y negativos. La zona sombreada de la figura indica el
rango de los momentos flectores posibles.
Sea una viga continua con un tramo central de longitud l y dos voladizos de longitudes
li y ld con las dimensiones que se indican en la Figura 1.48. Las cargas que actúan
sobre la estructura son:
- Una sobrecarga de uso uniforme q (γQ ) que puede actuar sobre cualquier parte de la
estructura
(c) Calcular y dibujar la lı́nea de influencia del momento flector en las secciones A,B
y en las secciones de los apoyos.
(d) Calcular y dibujar la envolvente de momentos en el vano central, acotando los
valores máximos en las secciones 1,A,B,2.
Datos: g = 3 kN/m, q = 5 kN/m, Q = 15 kN , li = 3 m, l = 14m, ld = 4 m a = 4 m
ELU equilibrio: γGf av = 0, 95, γGdes = 1, 10, γQ = 1, 50
ELU rotura: γGf av = 1, , γGdes = 1, 35, γQ = 1, 50
Para analizar el estado lı́mite de equilibrio respecto al vuelco alrededor del apoyo 2 se
toman momentos respecto del apoyo 2 con las cargas mayoradas en función de si el
momento es estabilizador o volcador según se muestra en la Figura 1.49.
172
Mest =0, 95 · 3 · = 411, 8 kN m
2
42 42
Mdes =1, 1 · 3 · + 1, 5 · 5 · + 1, 5 · 15 · 4 = 176, 4 kN m
2 2
Puesto que Mest > Mdes no es necesario anclar la viga en el apoyo 1. El vuelco respecto
del apoyo 1 es menos crı́tico.
1.7. ENVOLVENTES DE ESFUERZOS 45
Las lı́neas de influencia de las reacciones en los apoyos se muestran en la Figura 1.50.
La lı́nea de influencia del momento flector en la sección justo a la izquierda del apoyo
1 es independiente de las reacciones en los apoyos y vale
En este caso únicamente existe área negativa A− = −4.5. Los valores extremos de los
momentos flectores en 1− se calculan como
max,(+)
M1 − = 1, 0 · 3 · A− = −13, 5 kN m
max,(−)
M1 − = 1, 35 · 3 · A− + 1, 5 · 5 · A− + 1, 5 · 15 · (−3) = −119, 5 kN m
Análisis de la sección A
Las áreas positiva y negativa en la lı́nea de influencia del momento flector en la sección
A son A+ = 20 y A− = −5.5. Los valores extremos de los momentos flectores en A se
46 CAPÍTULO 1. ANÁLISIS ESTRUCTURAL Y SEGURIDAD
calculan como
max,(+)
MA = 1, 35 · 3 · (A+ + A− ) + 1, 5 · 5 · A+ + 1, 5 · 15 · (2, 86) = 273, 0 kN m
max,(−)
MA = 1, 0 · 3 · (A+ + A− ) + 1, 5 · 5 · A− + 1, 5 · 15 · (−2, 14) = −46, 0 kN m
Análisis de la sección B
Las áreas positiva y negativa en las lı́neas de influencia del momento en la sección A
son A+ = 24 y A− = −6.5. Los valores extremos de los momentos flectores en B se
calculan como
max,(+)
MB = 1, 35 · 3 · (A+ + A− ) + 1, 5 · 5 · A+ + 1, 5 · 15 · (3, 43) = 328, 0 kN m
max,(−)
MB = 1, 0 · 3 · (A+ + A− ) + 1, 5 · 5 · A− + 1, 5 · 15 · (−2, 29) = −47, 7 kN m
M2+ (x) = 21 − x 17 ≤ x ≤ 21
En este caso únicamente existe área negativa A− = −8. Los valores extremos de los
momentos flectores en 2+ se calculan como
max,(+)
M2+ = 1, 0 · 3 · A− = −24, 0 kN m
max,(−)
M2+ = 1, 35 · 3 · A− + 1, 5 · 5 · A− + 1, 5 · 15 · (−4) = −182, 4 kN m
Figura 1.51: Lı́neas de influencia para el flector en el apoyo 1, secciones A,B y apoyo 2
2.1 Introducción
Una estructura funicular (o antifunicular) es aquella que, bajo una determinada carga,
se ve sometida únicamente a esfuerzo axil de tracción (o de compresión). Dentro de
las estructuras funiculares lineales se distinguen básicamente dos tipologı́as: cables y
arcos. Los cables y las estructuras suspendidas adoptan la forma adecuada en función
de la carga, de tal forma que se alcance el equilibrio en su configuración deformada,con
todos los elementos en tracción.
49
50 CAPÍTULO 2. CABLES Y ARCOS
y además altamente eficiente. Por ello las estructuras funiculares se han usado en
ingenierı́a y arquitectura desde la más remota antigüedad hasta en los más innovadores
diseños.
Un cable del que cuelga una fuerza puntual se deforma en dos tramos rectos. Si soporta
n cargas puntuales la deformada define una poligonal de n + 1 tramos. Si el cable se ve
sometido a una carga repartida el cable se deforma según una parábola o una catenaria.
2.2 Cables
2.2.1 Antecedentes
Los cables flexibles se han utilizado para formar sistemas estructurales desde tiempos
prehistóricos hasta la actualidad. Algunas de las primeras estructuras hechas por el
hombre fueron puentes colgantes construidos con materiales naturales para atravesar
rı́os y barrancos. Los cables se han utilizado en puentes militares desde el siglo XV.
En tiempos modernos, los cables de acero se emplean en la construcción de puentes de
gran luz (Figuras 2.2 y 2.3). Hoy en dı́a, los records de luz y de altura se baten con
puentes y torres atirantados (Figura 2.4).
Figura 2.2: Puentes atirantados (a) Puente de Brooklyn (b) Golden Gate
2.2. CABLES 51
Figura 2.3: Estructuras atirantadas (a) Puente de Russki Island (b) Torre de Collserola
Un cable de peso despreciable sometido a cargas puntuales adopta una forma poligonal
de tramos rectos cada uno de los cuales está sometido a esfuerzo axil constante. Un
cable fijado en sus extremos sometido a n cargas puntuales da lugar a n + 1 tramos
rectos que configuran la deformada del cable y a n + 2 nudos incluyendo los apoyos, en
los que se puede plantear el equilibrio de fuerzas.
Por un lado, las incógnitas son las 4 reacciones, los axiles en los n+1 tramos y las flechas
en los n nudos, lo que supone un total de 2n + 5 incógnitas; por otro lado, en cada
nudo puede establecerse el equilbrio de fuerzas horizontales y verticales, lo que supone
un total de 2n + 4 ecuaciones. Por tanto, en general, el sistema no está determinado y
para cada sistema de cargas existe una familia de formas funiculares.
Para poder resolver la estructura es necesario añadir una condición adicional que de-
termina una única solución. Esta condición puede ser la longitud del cable, el axil en
alguno de los tramos del cable o la flecha en algún nudo.
Un cable está fijado en dos puntos y soporta dos cargas de valor P1 y P2 . La deformada
del cable se muestra la Figura 2.6.
Datos: a = 2 m, b = 2 m, P1 = 3 kN , P2 = 8 kN
2.2. CABLES 53
El problema tiene 8 incógnitas en fuerzas {H1 , V1 , H4 , V4 , T12 , T23 , T34 , h}. El número
de nodos permite establecer 8 ecuaciones de equilibrio por lo que el problema tiene
solución única. Los ángulos de inclinación de los distintos tramos son función de h:
h 4−h 2
α12 = arctan α23 = arctan α34 = arctan = 53, 1◦
2 2 1, 5
El equilibrio del nudo 4 relaciona las reacciones en dicho nudo y el axil T34 :
T34 = 6, 79 kN H1 = −4, 07 kN V1 = 5, 57 kN H4 = 4, 07 kN V4 = 5, 43 kN
54 CAPÍTULO 2. CABLES Y ARCOS
Conocido el valor de h puede obtenerse el ángulo α12 = 53.87◦ . Por equilibrio del nudo
2 se tiene
X
FH = T12 cos α12 − T23 cos α23 = 0
X
FV = T12 sin α12 − T23 sin α23 − 3 = 0
Se considera un cable sometido a una carga distribuida por unidad de proyección ho-
rizontal. La distancia entre apoyos es 2l y la flecha en el punto medio es f según se
muestra en la Figura 2.7. Se sitúa el origen de coordenadas en el punto del cable cuya
pendiente es nula, con el eje de ordenadas hacia arriba y el eje de abscisas hacia la
derecha.
Se considera un elemento diferencial del cable, en cuyos extremos actúa un cierto axil
T y T + dT con inclinaciones θ y θ + dθ.
y despreciando los términos de segundo orden, se obtienen las siguientes relaciones entre
el axil T y la carga q(x)
d
[T cos θ] = 0 (2.4a)
dx
d
[T sin θ] = q(x) (2.4b)
dx
que indican que la componente horizontal del axil es constante mientras que la compo-
nente vertical varı́a linealmente si la carga q(x) es constante. En tal caso, las ecuaciones
que relacionan la tensión y el ángulo θ son:
T cos θ = H (2.5a)
T sin θ = qx (2.5b)
El ángulo θ es la pendiente del cable, igual al giro de cada sección con respecto a la
horizontal. De las expresiones (2.5a) y (2.5b) se tiene
dy qx qx2
tan θ = = ⇒ y(x) = +C (2.6)
dx H 2H
y(x)|x=0 = C = 0 (2.7a)
ql2 ql2
y(x)|x=l = =f ⇒H= (2.7b)
2H 2f
56 CAPÍTULO 2. CABLES Y ARCOS
d
[T cos θ] = 0 (2.12a)
dx
d
[T sin θ] = q(s) (2.12b)
ds
2.2. CABLES 57
En este caso la componente horizontal del axil permanece constante mientras que la
componente vertical varı́a linealmente con la directriz del cable a través de la longitud
s. Integrando las ecuaciones anteriores e imponiendo las correspondientes condiciones
de apoyo q(s) = q, con s(x)|x=0 = 0, se tiene
T cos θ = H (2.13a)
T sin θ = qs (2.13b)
de donde se obtiene
dy qs
tan θ = = (2.14)
dx H
El parámetro arco s puede relacionarse con las coordenadas cartesianas x, y a nivel
diferencial a través de la siguiente expresión:
s 2
Z Z xp Z x
2 2
dy
s = ds = (dx) + (dy) = 1+ dx (2.15)
0 0 dx
Derivando (2.14) se obtiene la ecuación diferencial del cable actuando bajo su propio
peso: s 2
d2 y q ds q dy
2
= = 1+ (2.16)
dx H dx H dx
El axil del cable y sus valores extremos pueden calcularse en función de cos θ.
s 2
H p dy qx
T = = H 1 + tan2 θ = H 1+ = Hcosh = qy + H (2.21)
cos θ dx H
ql2
f=
8H
2.2. CABLES 59
Por su parte, si se considera la carga distribuida por unidad de longitud del cable, la
deformada es una catenaria. La expresión para la flecha en este caso es
H ql
f= cosh −1
q 2H
Anclaje (kN) fpar (m) f /l Anclaje (kN) fcat (m) f /l fcat /fpar
7,5 16,7 1/6 7,5 17,3 1/6 1,038
10 12,5 1/8 10 12,8 1/8 1,021
15 8,3 1/12 15 8,4 1/11 1,009
25 5,0 1/20 25 5,0 1/19 1,003
50 2,5 1/40 50 2,5 1/39 1,001
Se observa que para relaciones f /l < 1/12 la diferencia de flechas es inferior al 1%,
mientras que para valores 1/8 < f /l < 1/6 dicha diferencia alcanza valores de hasta el
5%. A nivel ingenieril es habitual suponer que ambas curvas coinciden.
Debido a la carga distribuida, la deformada del cable será una parábola de ecuación:
qx2
y(x) =
2H
H
T max = = 447, 2 kN
cos θmax
T max
σ max = = 178, 9 MPa
A
Un cable está soportado en dos extremos situados a una distancia l con un cierto
desnivel, según muestra la Figura 2.11.
Determinar los axiles en los apoyos y en el punto más bajo del cable. Dimensionar la
sección del mismo si la tensión admisible es σ adm .
Figura 2.11: Cable sustentado en apoyos a distinto nivel bajo carga uniformemente distribuida
Puesto que los apoyos no están a la misma altura, la parábola no es simétrica. El punto
más bajo se sitúa a una distancia a del apoyo izquierdo. En dicho punto de pendiente
nula se sitúa el origen de coordenadas con lo que las flechas son conocidas en los apoyos.
Esto permite plantear dos ecuaciones en las que las incógnitas son la distancia a y la
reacción horizontal H:
q(l − a)2
f1 = = h1 = 40 m
2H
qa2
f2 = = h1 − h2 = 20 m
2H
H
T =
cos θ
62 CAPÍTULO 2. CABLES Y ARCOS
Un cable se encuentra sostenido por dos mástiles de altura h entre dos puntos 2 y 3
cubriendo una luz total l. La flecha máxima en el centro del vano central es f . El
esquema se completa con dos cables 12 y 34 que hacen las veces de tirantes de retenida
anclándose en unos muertos de cimentación colocados a una distancia d de los mástiles,
según muestra la Figura 2.12. El peso propio de los cables se considera despreciable. La
estructura debe resistir una carga q uniformemente distribuida en proyección horizontal
entre ambos mástiles.
Determinar el trazado del cable y el axil máximo en el tramo 23.
Calcular la fuerza necesaria en los tirantes de retenida de tal manera que los mástiles
reciban únicamente cargas verticales en los puntos 1 y 4.
Datos: l = 40 m, h = 15 m, d = 15 m, f = 3 m, q = 10 kN/m
La sección con pendiente nula se sitúa en el punto medio del tramo 23 por la simetrı́a
del problema.
La reacción horizontal H es
ql2
H= = 666, 7 kN
2f
El axil máximo se dará en los puntos más alejados del origen de coordenadas, es decir,
en los soportes, y se calcula como
H
T max = = 696 kN
cos θmax
Finalmente, para que los soportes reciban únicamente una carga vertical, es necesario
que los cables de retenida, con inclinación α = 45◦ , proporcionen una componente
horizontal que equilibre la componente horizontal H en los anclajes (Figura 2.13).
H
T0 cos α = H ⇒ T0 = = 942, 8 kN
cos α
64 CAPÍTULO 2. CABLES Y ARCOS
El peso propio de los cables se considera negligible. La estructura debe resistir una
carga q uniformemente distribuida entre ambos mástiles.
Calcular la fuerza necesaria en los tirantes de retenida de tal manera que los mástiles
reciban únicamente cargas verticales en los puntos 1 y 4.
Datos: l = 30 m, h1 = 15 m, h2 = 10 m, d1 = 15 m, d2 = 10 m, f = 6 m, q = 10 kN/m
Puesto que los apoyos no están a la misma altura, la parábola no es simétrica. El punto
más bajo del cable en el que la pendiente es nula se sitúa a una distancia a del apoyo
izquierdo desconocida. En él se sitúa el origen para el cálculo de flechas, con lo que
éstas se conocen en los apoyos. Esto permite plantear dos ecuaciones en las que las
incógnitas son la distancia a y la reacción horizontal H.
en cuyo vértice se sitúa el origen de coordenadas, quedando éste situado a 21,3 m del
apoyo izquierdo y a 8,7 m del apoyo derecho. En dicho vértice se da una tensión igual
a H en la dirección tangente a dicha parábola, es decir, horizontal.
La tensión en cualquier punto del cable se calcula a partir del valor de H por lo que
para conocerla en los apoyos basta con obtener el valor del ángulo θ en dichos puntos.
Esto se consigue con la ley de giros derivando la ecuación de la parábola.
(
dy qx y20 = tan θ2 = 0, 56 ⇒ θ2 = 29, 4◦ ⇒ T2 = 434, 1 kN
= = 0, 02644x →
dx H y20 = tan θ3 = 0, 23 ⇒ θ3 = 13, 0◦ ⇒ T3 = 388, 1 kN
H
T0 cos α = H ⇒ T0 = = 534, 8 kN
cos α
2.3 Arcos
2.3.1 Antecedentes
Los antiguos romanos fueron los primeros en usar los arcos en la construcción de puentes
y acueductos, empleando la forma semicircular debido a la facilidad de su construcción.
Los arcos romanos y románicos se construı́an con bloques de piedra o dovelas, y con
juntas secas o con mortero. Durante la Edad Media, los arcos góticos se generalizan en la
construcción de las grandes catedrales europeas, combinados con el uso de contrafuertes
que proporcionan el necesario empuje horizontal en los apoyos. En el siglo XVIII
se utilizaron arcos de mamposterı́a para construir numerosos puentes, muchos de los
cuales sobreviven hoy en dı́a sobrellevando cargas muy superiores a las de diseño. En
la actualidad, los arcos son de acero o de hormigón armado o pretensado y puede
soportar cargas tanto de tracción como de compresión, salvando luces de varios cientos
de metros.
66 CAPÍTULO 2. CABLES Y ARCOS
Se puede considerar las estructuras en arco como sistemas simétricos respecto a los ca-
bles. En los arcos, la directriz coincide con la directriz funicular de un cable sometido
a cargas de signo opuesto, por lo que cuando nos referimos a arcos hablamos de estruc-
turas antifuniculares de las cargas.
- Los arcos son rı́gidos, por lo que cada configuración es antifunicular de un único es-
tado de carga. En consecuencia, cualquier variación de carga implica la aparición
de esfuerzos flectores en el arco.
Se considera una arco biarticulado sometido a una carga repartida de valor q(x) cualquiera,
como se muestra en la Figura 2.16. La estructura es hiperestática de grado 1 por tener
4 reacciones incógnita en los apoyos frente a las tres ecuaciones de la estática en una
estructura plana.
2.3. ARCOS 67
En el estado I las reacciones horizontales son nulas. En dicho arco curvo sometido a
una carga distribuida arbitraria la resultante de fuerzas verticales y el momento flector
{QV , MV } en cada sección son equivalentes a los correspondientes a una viga recta
de longitud igual a la luz del arco. Esto significa que los esfuerzos en cada sección del
arco {Nq , Qq , Mq } pueden calcularse proyectando las resultantes de fuerzas y momentos
según el ángulo θ que determina la orientación en dicha sección.
Figura 2.17: Convenio de signos para la definición de esfuerzos sobre una sección de arco
(a) Definición frontal y dorsal (b) Variación del ángulo que define la orientación de la sección
2.3. ARCOS 69
H
N = QV sin θ + H cos θ = (2.28a)
cos θ
Q = QV cos θ − H sin θ = 0 (2.28b)
4f 4f
y(x) = x (l − x) y 0 (x) = tan θ(x) = (l − 2x) (2.29)
l2 l2
70 CAPÍTULO 2. CABLES Y ARCOS
La traza de la generatriz del arco viene dada por una función y(x) y la orientación de
cada sección θ(x) se determina a través de su derivada tan θ = y 0 (x): La determinación
de la ley de momentos flectores involucra la geometrı́a del arco y en particular la función
y(x) que permite situar una sección con respecto a su origen. Por su parte, las leyes
de axil y cortante se refieren a la resultante de fuerzas en dirección normal y paralela
a cada sección usando las expresiones
1 1
sin θ = √ cos θ = √ (2.30)
1+ tan−2 θ 1 + tan2 θ
Si la directriz del arco coincide con la lı́nea de empujes el arco es antifunicular, pero
la variación de la directriz con respecto a la lı́nea de empujes provoca la aparición de
flexiones.
Por ejemplo, considerando la deformación por axil, disminuye la reacción H con res-
pecto al valor de la reacción antifunicular. En este caso los esfuerzos flectores en la
2.3. ARCOS 71
estructura son:
M = MV − Hy 6= 0 (2.31)
La introducción de una rótula en una de las secciones intermedias del arco implica
que la estructura queda estáticamente determinada (Figura 2.20), es decir la reacción
H se puede determinar exclusivamente a partir de las ecuaciones de equilibrio y de la
condición de rótula:
X
FH = H1 − H3 = 0 (2.33a)
X X q
FV = V1 + V3 − FV = 0 (2.33b)
X X q
M1 = V 3 l − M1 = 0 (2.33c)
X l X q
M2− = V1 − Hf + M2 − = 0 (2.33d)
2
La primera ecuación determina que ambas reacciones horizontales son iguales a un valor
H que puede obtenerse de la ecuación de momento flector nulo en la rótula. Por su
parte, las reacciones verticales V1 y V3 se obtienen de las ecuaciones de equilibrio global
de fuerzas verticales y de momentos.
72 CAPÍTULO 2. CABLES Y ARCOS
M = MV − Hy = 0 (2.35)
ya que MV es una ley parabólica con valor máximo ql2 /8. Por tanto, un arco parabólico
con apoyos al mismo nivel sometido a carga uniformemente distribuida es siempre un
arco antifunicular y, por tanto, trabaja exclusivamente a axil de compresión.
Datos: l = 90 m, f = 13 m, q = 30 kN/m
El arco triarticulado es una estructura isostática, por lo que las reacciones pueden
2.3. ARCOS 73
H1 = 2337 kN H3 = −2337 kN
V1 = 1350 kN V3 = 1350 kN
Puesto que se trata de un arco antifunicular, la única ley que es necesario definir es la
de axiles. Para ello se calcula el valor del axil en las secciones relevantes teniendo en
cuanta la simetrı́a: sección de apoyo (1), sección a cuartos de luz (A) y sección de la
rótula (2).
Determinar el trazado de la lı́nea directriz del arco y las leyes de esfuerzos de la estruc-
tura.
Datos: l = 20 m, f = 3, 5 m, q = 10 kN/m
La carga que soporta el arco aplicada a una viga recta darı́a lugar a una ley momentos
lineal en la mitad izquierda y parabólica en su mitad derecha. La geometrı́a parabólica
incumple la condición de antifunicularidad por lo que la respuesta del arco incluirá los
tres esfuerzos: axil, cortante y momento flector.
por equilibrio:
H1 = 71, 4 kN H3 = −71, 4 kN
V1 = 25 kN V3 = 75 kN
Para dibujar las leyes de esfuerzos se calculan éstos en secciones relevantes del arco:
apoyos, secciones a cuartos de luz y rótula central.
Sección apoyo izquierdo (1):
Las coordenadas (x, y) de la sección son (0;0).
La pendiente de la directriz en 1 es tan θ1 = 0, 7 y su inclinación es θ1 = 35, 0◦ .
Los esfuerzos en la sección 1 son
Figura 2.24: Leyes de esfuerzos: (a) axil (b) cortante (c) momento flector
2.3. ARCOS 77
La carga que soporta el arco aplicada a una viga recta darı́a lugar a una ley momentos
lineal-constante-lineal en la luz del arco. La geometrı́a parabólica incumple la condición
de antifunicularidad por lo que la respuesta del arco incluirá los tres esfuerzos: axil,
cortante y momento flector.
Las reacciones pueden obtenerse mediante las ecuaciones de equilibrio:
H1 = 14, 3 kN H3 = −14, 3 kN
V1 = 10 kN V3 = 10 kN
Para dibujar las leyes de esfuerzos se calculan éstos en secciones relevantes del arco:
apoyos, secciones a cuartos de luz y rótula central.
Sección apoyo izquierdo (1):
Las coordenadas (x, y) de la sección son (0;0).
La pendiente de la directriz en 1 es tan θ1 = 0, 7 y su inclinación es θ1 = 35, 0◦ .
78 CAPÍTULO 2. CABLES Y ARCOS
Figura 2.26: Leyes de esfuerzos: (a) axil (b) cortante (c) momento flector
Determinar el trazado de la lı́nea directriz del arco y las leyes de esfuerzos de la estruc-
tura.
Datos: l = 20 m, f = 3, 5 m, P = 100 kN
La carga que soporta el arco aplicada a una viga recta darı́a lugar a una ley momentos
bilineal en la luz del arco. La geometrı́a parabólica incumple la condición de antifu-
nicularidad por lo que la respuesta del arco incluirá los tres esfuerzos: axil, cortante y
momento flector.
80 CAPÍTULO 2. CABLES Y ARCOS
La estructura es isostática gracias a la rótula situada en el centro del arco, por lo que
las reacciones se calculan por equilibrio:
H1 = 71, 4 kN H3 = −71, 4 kN
V1 = 25 kN V3 = 75 kN
Para dibujar las leyes de esfuerzos se calculan éstos en secciones relevantes del arco:
apoyos, sección a cuartos de luz y rótula central.
Sección apoyo izquierdo (1):
Las coordenadas (x, y) de la sección son (0;0).
La pendiente de la directriz en 1 es tan θ1 = 0, 7 y su inclinación es θ1 = 35, 0◦ .
Los esfuerzos en la sección 1 son
N1 = 25 sin θ1 + 71, 4 cos θ1 = 72, 9 kN (→←)
Q1 = 25 cos θ1 − 71, 4 sin θ1 = −20, 5 kN (↓↑)
M1 = 0
La ley de momentos flectores en la mitad izquierda del arco puede expresarse como:
(
M1B (x) = 2, 5x(x − 10) 0 ≤ x ≤ 15
M (x) = 2
MB2 (x) = 2, 5x[(x − 25) − 325] 15 ≤ x ≤ 20
Figura 2.28: Leyes de esfuerzos: (a) axil (b) cortante (c) momento flector
3 Inestabilidad
3.1 Introducción
Figura 3.1: Pandeo de una viga esbelta sometida a carga axial de compresión
83
84 CAPÍTULO 3. INESTABILIDAD
o pandeo.
Dicho fenómeno se desarrolla dentro del rango elástico, por lo que la pieza recupera su
forma recta inicial cuando se retira la carga axial. Esta constatación pone de manifiesto
las limitaciones del análisis de primer orden de la resistencia de materiales, en particular,
el principio de rigidez que desprecia la deformación de la pieza en el planteamiento de
las condiciones de equilibrio.
Cuanto más corto (menos esbelto) es una pieza en relación a sus dimensiones transver-
sales, más probable es que falle por resistencia y menos probable es su fallo por pandeo.
Recı́procamente, cuando más esbelta es la pieza, más posible es que pandee antes de
que se produzca el agotamiento del material. Existe un valor lı́mite de esbeltez que
separa el fallo por resistencia del fallo por inestabilidad.
M (x)
y 00 (x) = (3.1)
EI
Se considera una viga biarticulada de longitud l sometida a una carga axial N , según
muestra la Figura 3.2, dando lugar a una ley de flechas y(x).
La ley de esfuerzos en este caso es M (x) = −N y(x) por lo que la ecuación diferencial
de equilibrio de la viga es
y 00 (x) + k 2 y(x) = 0 (3.3)
y(x)|x=0 = 0 → B = 0 (3.5a)
y(x)|x=l = 0 → A sin(kl) = 0 (3.5b)
EI
NE := π 2 (3.7)
l2
Figura 3.3: Situaciones de equilibrio (a) Estable (b) Neutro (c) Inestable
88 CAPÍTULO 3. INESTABILIDAD
La carga crı́tica para la pieza biarticulada se toma como referencia para las cargas
crı́ticas de vigas simples con diferentes condiciones de apoyo. A este efecto se define
el factor τ = Ncr /NE como la relación entre la carga Ncr en un determinado caso y la
carga crı́tica de Euler.
V1
y 00 (x) + k 2 y(x) − k 2 x=0 (3.9)
N
V1
y(x) = A sin(kx) + B cos(kx) + x (3.10)
N
V1
y 0 (x) = Ak cos(kx) − Bk sin(kx) + (3.11)
N
y(x)|x=0 = 0 ⇒ B = 0 (3.12a)
N
y(x)|x=l = 0 ⇒ V1 = −A sin(kl) (3.12b)
l
y 0 (x)x=l = 0 ⇒ A [kl cos(kl) − sin(kl)] = 0
(3.12c)
EI EI
kl = z ∗ ⇒ Ncr = (z ∗ )2 = 20, 19 2 (3.13)
l2 l
Se considera ahora el caso de una pieza biempotrada sometida a esfuerzo axil N , con
reacciones M0 en ambos extremos empotrados y cuya flecha se nota por y(x).
90 CAPÍTULO 3. INESTABILIDAD
Una viga biempotrada con carga simétrica como es el caso, es una estructura hiper-
estática de grado 1. En consecuencia, los esfuerzos pueden expresarse en función de
una incógnita hiperestática; por ejemplo, el momento reacción M0 en los apoyos. La
ley de momentos expresada en función de M0 es
M0
y 00 (x) + k 2 y(x) − k 2 =0 (3.17)
N
M0
y(x) = A sin(kx) + B cos(kx) + (3.18)
N
y 0 (x) = Ak cos(kx) − Bk sin(kx) (3.19)
M0
y(x)|x=0 = 0 ⇒ B = − (3.20a)
N
0
y (x) x=0 = 0 ⇒ Ak = 0
(3.20b)
EI
Ncr = 4m2 π 2 , m = 1, 2, . . . (3.24)
l2
M0 l 2
2mπ
y(x) = 1 − cos x , m = 1, 2, . . . (3.25)
4EIm2 π 2 l
La segunda ecuación tiene dos posibles soluciones, la primera de las cuales lleva a la
solución trivial con axil nulo, mientras que la segunda solución implica que A = 0. La
función de flechas y(x) puede por tanto escribirse como
Existen por tanto infinitas soluciones para la ecuación de equilibrio y sus condiciones
de apoyo. La carga crı́tica mı́nima se obtiene para m = 1 como múltiplo de la carga
crı́tica de Euler con un factor τ = 1/4:
1 EI 1
Ncr = π 2 2 = NE (3.34)
4 l 4
Se considera una Viga biarticulada de longitud l, con un apoyo elástico que restringe
el giro en el apoyo derecho, cuya constante es km . Dicha viga se encuentra sometida a
un esfuerzo axil N , según muestra la Figura 3.7.
3.3. TEORÍA DE EULER PARA COLUMNAS ESBELTAS. PIEZAS IDEALES 93
Figura 3.7: Pieza biarticulada con apoyo elástico. Equilibrio en la configuración deformada
Las reacciones pueden calcularse en función del momento reacción M que es desconocido
ya que depende del giro en el extremo 2 φ2 .
M
V1 = = −V2 (3.35)
l
y(x)|x=0 = 0 → B = 0 (3.39a)
V1 l
y(x)|x=l = 0 → A sin(kl) + =0 (3.39b)
N
Esta ecuación tiene dos posibles soluciones. La solución trivial es A = 0 con lo que
la flecha es nula. La segunda solución es que se anule el término que multiplica a la
constante A, lo que se transforma en un ecuación implı́cita:
Nl
kl = tan(kl) 1 − (3.43)
km
La menor de estas cargas crı́ticas es Ncr = 2, 046NE , la carga crı́tica de una viga
empotrado-articulada, hallada en 3.3.1.2.
Dependiendo de las condiciones de apoyo, la carga crı́tica en piezas simples puede ex-
presarse en función de la carga crı́tica de Euler para la pieza biarticulada. Se introduce
ası́ el concepto de longitud efectiva de pandeo aplicable a toda pieza simple en función
de sus condiciones de apoyo.
Dada una pieza de longitud l sujeta a unas determinadas condiciones de apoyo, la
longitud efectiva de pandeo le se define como la longitud de una pieza biapoyada cuya
3.3. TEORÍA DE EULER PARA COLUMNAS ESBELTAS. PIEZAS IDEALES 95
En la Figura 3.8 se muestran las longitudes efectivas de pandeo para distintas piezas
simples.
La Tabla 3.1 resume las cargas crı́ticas para las distintas vigas simples ası́ como su
longitud de pandeo asociada.
Figura 3.10: Zona admisible de tensión de una pieza a compresión en análisis de segundo orden
Sea una viga biarticulada de lontigud l. Determinar la carga crı́tica de la viga para dos
secciones transversales, según muestra la Figura 3.11:
a) Sección rectangular de hormigón de dimensiones b y h.
b) Sección tubular de acero formada por la unión de dos perfiles UPN200, de los cuales
se conocen sus propiedades mecánicas y la posición de su centro de gravedad (G).
Datos: l = 8 m, b = 30 cm, h = 20 cm, Eh = 20 GP a, btub = 30 cm, htub = 20 cm,
Ea = 200 GP a
98 CAPÍTULO 3. INESTABILIDAD
Figura 3.11: Viga biarticulada y secciones transversales (a) rectangular (b) tubular
La carga crı́tica de una viga biarticulada es la carga crı́tica de Euler, que depende de
la elasticidad del material, de la longitud de la viga y de la inercia de la sección.
EI
Ncr = π 2
l2
Conocidas las dos primeras variables, cada sección transversal presenta dos inercias
principales. La carga crı́tica se calcula con la menor de ambas.
En la sección rectangular los momentos de inercia son
1
Iy = 20 · 303 = 45 000 cm4
12
1
Iz = 30 · 203 = 20 000 cm4
12
La inercia mı́nima es Iz = 20 000 cm4 , con lo que el pandeo se produce alrededor del
eje horizontal z y la carga crı́tica de pandeo la viga (a) es
EI
Ncr = π 2 = 616, 9 kN
l2
La inercia mı́nima es I = 2237 cm4 , con lo que el pandeo se produce alrededor del eje
vertical y y la carga crı́tica de pandeo de la viga (b) es
EI
Ncr = π 2 2 = 441, 6 kN
l
La carga crı́tica de la viga con sección tubular (b) es menor que la carga crı́tica de la
viga con sección rectangular (a) por lo que el pandeo se producirá antes en el caso (b)
Para que la viga pandee dicho axil ha de alcanzar el valor crı́tico, que se expresa como un
múltiplo de la carga crı́tica de Euler. En el caso de una viga biempotrada, Ncr = 4NE .
En consecuencia, el incremento de temperatura necesario para provocar el pandeo de
la viga será
4NE I 4π 2
N = Ncr = 4NE ⇒ ∆T = = 4π 2 =
EAα Aαl2 αλ2
En el caso particular de que la sección transversal sea circular de radio R la esbeltez y
la carga crı́tica son
1 πR 2
2l
λ= ⇒ ∆T =
R α l
100 CAPÍTULO 3. INESTABILIDAD
Sea una estructura articulada formada por dos barras de sección transversal A, cons-
truidas en un material cuya tensión lı́mite es σ adm . Sobre el nudo 2 de dicha estructura,
actúa una carga P vertical, según muestra la Figura 3.13.
Determinar carga de colapso de la estructura considerando el dimensionamiento por
resistencia y por pandeo.
Datos: l = 2 m, α = 20◦ , E = 200 GP a, A = 16 cm2 , I = 21, 33 cm4 , σ adm =
300 M P a
La barra 1 está traccionada por lo que únicamente puede fallar por tensión cuando ésta
alcance el valor lı́mite σ adm :
N1 P
σ1 = = = σ adm ⇒ Pres,1
max
= Aσ adm tan α = 174, 7 kN
A Atan α
La barra 2 está comprimida por lo que puede fallar por tensión cuando ésta alcance
el valor lı́mite σ adm o bien por pandeo si el axil alcanza el valor crı́tico. Al ser una
estructura articulada, la carga crı́tica por pandeo es la carga de Euler NE :
N2 P
σ2 = = = σ adm ⇒ Pres,2
max
= Aσ adm sin α = 164, 2 kN
A A sin α
P max EI
N2 = = NE ⇒ Ppan = NE sin α = π 2 2 cos2 α sin α = 31, 8 kN
sin α l2
3.3. TEORÍA DE EULER PARA COLUMNAS ESBELTAS. PIEZAS IDEALES 101
Sea una estructura articulada formada por tres barras de sección circular con radios
R1 , R2 y R3 construida con un material cuya tensión lı́mite es σ adm , sobre la que actúa
una carga P vertical según indica la Figura 3.14.
Determinar la carga crı́tica de pandeo de la estructura.
Datos: h = 2 m, α = 60◦ , E = 200 GP a, σ adm = 300 M P a, R1 = R3 = 4 cm,
√
R2 = 2 2 cm
P 2
N1 = N3 = − N2 = − P
3 3
Todas las barras están comprimidas por lo que pueden fallar por tensión o por pandeo.
Los valores de la carga P máxima para la barra 1 (y por simetrı́a la barra 3) para
ambos supuestos son las siguientes:
N1 P
σ1 = = = σ adm ⇒ Pres,1
max
= 3A1 σ adm = 4524 kN
A1 3A1
P max EI1
N1 = = NE ⇒ Ppan,1 = 3NE = π 2 2 cos2 α = 744 kN
3 l1
En el caso de la barra central 2 los valores para la carga P máximos son los siguientes:
N2 2P 3
σ2 = = = σ adm ⇒ Pres,2
max
= A2 σ adm = 1131 kN
A2 3A2 2
2 max 3 EI2
N2 = P = NE ⇒ Ppan,2 = NE = 3π 2 2 = 372 kN
3 2 2l2
Una vez analizado el caso de la pieza ideal, la teorı́a de Euler se extiende a condiciones
distintas de las establecidas para piezas ideales. Se consideran los casos en los que la
carga axial sobre la viga no se aplica de forma centrada, la viga presenta una curvatura
inicial o bien la pieza está sometida a carga transversal.
Figura 3.15: Equilibrio en la configuración deformada de una pieza con axil excéntrico
La solución de la ecuación es
y(x) = e [1 − tan(β)sin(kx) − cos(kx)] (3.54)
- La carga crı́tica Ncr es la misma que en el caso de la pieza ideal. Una vez más,
la carga crı́tica depende únicamente de las condiciones de apoyo.
- Para un axil diferente al crı́tico, la flecha no es nula sino que evoluciona de forma
continua con la relación ρ = N/NE . La Figura 3.16 muestra la evolución del
término tan β en función de ρ.
Figura 3.17: Equilibrio en la configuración deformada de una pieza con curvatura inicial
y(x)|x=0 = 0 ⇒ B = 0 (3.61a)
y(x)|x=l = 0 ⇒ Asin(kl) = 0 (3.61b)
En consecuencia, para axiles diferentes del valor crı́tico la solución para la posición final
de cada sección de la pieza ŷ(x) es
1 x
ŷ(x) = y0 sin π (3.63)
1−ρ l
De esta solución se extrae un valor de la carga crı́tica de pandeo para el axil aplicado
N igual a la carga crı́tica de Euler (NE ), coincidente con la carga crı́tica de una pieza
biarticulada ideal. La flecha en este caso es de amplitud A indeterminada.
- Cuando el axil tiende al valor crı́tico Ncr , la flecha tiende a una ası́ntota vertical.
Esta respuesta estructural se conoce como divergencia de equilibrio. La Figura
3.18 muestra la evolución del factor f (ρ) que gobierna la evolución de la flecha
para una sección de la viga dada.
Al problema de pandeo de una viga con deformación inicial también se le conoce como
pandeo con curvatura inicial χ0 (x). La relación entre deformación y curvatura es
χ0 (x) = y 00 (x).
Mq (x)
y 00 (x) + k 2 y(x) − k 2 =0 (3.66)
N
En el caso particular de que la carga sea uniforme (Figura 3.20) la resolución del
problema de segundo orden es la siguiente:
108 CAPÍTULO 3. INESTABILIDAD
qx(l − x) q
y(x) = Asin(kx) + Bcos(kx) + + (3.70)
2N N k2
Aplicando las condiciones de apoyo en los extremos se obtienen las constantes de inte-
gración:
q
y(x)|x=0 = 0 ⇒ B = − (3.71a)
N k2
q kl
y(x)|x=l =0⇒A=− 2
tan (3.71b)
Nk 2
ql4 h x x
2 x
x i
y(x) = − tan(β)sin 2β + cos 2β − 2β 1 − − 1 (3.73)
16EIβ 4 l l l l
3.4. PIEZAS NO IDEALES 109
El término tan β merece especial atención. El paraḿetro β depende del axil N pero su
tangente no está acotada, por lo que existe un valor de β para el cual dicha función, y
en consecuencia la flecha y(x), tiende a infinito.
π
tan β → ∞ ⇒ β = m , m = 1, 3, . . . ⇒ Ncr = NE (3.74)
2
Por tanto, existe un valor crı́tico del axil de compresión, para el cual la flecha tiende
a infinito. Este axil crı́tico es igual al hallado en la misma viga biarticulada sin carga
repartida, lo que permite extraer algunas conclusiones relevantes:
- La solución en flechas del problema es una función no nula que tiende a una
ası́ntota vertical cuando el axil aplicado tiende a la carga crı́tica en la cual la
flecha tiende a un valor infinito.
5ql4 12 2cos−1 β − β 2 − 2
y(x)|x=l/2 = − (3.75)
384EI 5β 4
La flecha en la sección central puede expresarse como el producto del resultado obtenido
en un análisis lineal de primer orden y un factor amplificador dependiente de β, es
decir, de la relación entre el axil aplicado N y la carga crı́tica de Euler NE en una viga
biarticulada, lo que pone de manifiesto la influencia del análisis de segundo orden.
y(x)|x=l/2 ≡ f (β)y(x)|lineal
x=l/2 (3.76)
1 5ql4
y(x)|x=l/2 ≈ − (3.77)
1 − ρ 384EI
110 CAPÍTULO 3. INESTABILIDAD
Es notoria la similitud entre los casos de pieza con axil excéntrico, con curvatura inicial
o con carga transversal. Todos estos casos coaparición de un término de flexión.
Es asimismo notorio que la carga crı́tica es idéntica en todos los casos, lo que permite
afirmar que la expresión de dicha carga no depende de las imperfecciones en cuanto
a la forma de la directriz o el punto de aplicación del axil, sino únicamente de las
condiciones de apoyo.
Los estados de equilibrio neutro de una pieza elástica que pandea se corresponden
con estados de energı́a potencial total estacionaria. Este hecho se puede utilizar para
determinar la carga crı́tica de la pieza, lo que resulta particularmente útil cuando se
desconoce la deformada de la pieza pandeada.
Se considera la pieza biarticulada recta de la Figura 3.21, sometida a una fuerza ho-
rizontal de valor N , que produce un esfuerzo axil centrado. La longitud inicial de la
viga es l. Si se desprecia la deformación por axil, una vez deformada la viga debido al
pandeo, la proyección horizontal de la viga es l0 = l − δ, siendo δ el desplazamiento del
apoyo móvil.
Π := U + H (3.78)
M (x)
y 00 (x) = (3.80)
EI
H = −N δ (3.81)
1 dy 2
Z l
H=− N dx (3.85)
0 2 dx
Los estados de equilibrio son los puntos estacionarios de este potencial. Sin embargo,
la ecuación (3.86) es integral y no puede resolverse para hallar la deformada y(x). Para
imponer la condición de estacionariedad, es necesario suponer una función y(x) como
aproximación a la solución. Supóngase que la función y(x) se desarrolla en serie de
Fourier en la forma
∞
X x
y(x) = λm sin mπ (3.87)
l
m=1
Nótese que la función y(x) queda parametrizada en función de sus coeficientes de Fourier
λm , por lo que determinar la función de flecha y(x) se reduce a obtener tales coeficientes.
Sustituyendo la expresión de la flecha (ec. 3.87), en la expresión de la energı́a potencial
(ec. 3.86) e integrando se tiene
∞ ∞
π4 X
4 2 π2 X 2 2
Π=U +H = EI m λ m − N m λm (3.89)
4l3 4l
m=1 m=1
π2 π 2 EI
∂Π
= m2 2 − N λm = 0 (3.90)
∂λm 2l l
3.5. ENERGÍA POTENCIAL TOTAL ESTACIONARIA. MÉTODO DE RAYLEIGH-RITZ 113
Se observa que cada término del desarrollo en serie de y(x) representa una función
armónica, y viene asociado a un valor crı́tico del axil. La Figura 3.22 representa los
tres primeros términos de dicho desarrollo en serie. Nótese que las cargas crı́ticas son
múltiplos de la carga crı́tica de Euler. En particular, para el primier término, m = 1,
la carga crı́tica es Ncr = NE . Obsérvese que el resultado obtenido es el mismo por el
método diferencial y por el método integral.
transversal. Z l
H = −N δ − q(x)y(x)dx (3.92)
0
Nótese que el término debido a la carga transversal no depende del axil N por lo que
no modifica la parte izquierda de la ecuación de equilibrio, sino la parte derecha. Esta
observación lleva a la conclusión de que la carga crı́tica de la viga no se ve modificada
por la carga transversal q(x), sino que depende únicamente de las condiciones de apoyo.
En este caso también se concluye que la amplitud λ1 tiende a infinito.
Si el axil no alcanza el valor crı́tico, la amplitud λ1 de los modos de deformación sı́ está
definida y depende de las cargas transvesales.
Las conclusiones halladas mediante el planteamiento energético (integral), coinciden
con las obtenidas en 3.4.3.
y(x)|x=0 = 0 y(x)|x=l = yl
0
(3.95)
y (x)|x=0 = 0 M (x)|x=l = EIy 00 (x)|x=l = 0
Las cuatro condiciones de apoyo anteriores permiten definir la función cúbica de aprox-
imación de la solución en función del parámetro yl , ası́ como sus derivadas:
yl 2 3yl 3yl
y(x) = x (3l − x) y 0 (x) = x (2l − x) y 00 (x) = (l − x) (3.96)
2l3 2l3 l3
Este resultado es una aproximación al valor real, Ncr = 0, 25NE , hallado en 3.3.1.4,
con un error de sólo el 0, 16 %.
116 CAPÍTULO 3. INESTABILIDAD
Las cargas crı́ticas obtenidas por este método son siempre superiores a los valores
exactos. Esto se debe a que el error de aproximación de la deformada se traduce en
una aproximación por exceso de la rigidez de la pieza, con el consiguiente incremento
en la estimación de la carga crı́tica.
Las estructuras reales pueden colapsar para valores de cargas muy diferentes a las cargas
crı́ticas predichas por el análisis lineal de bifurcación de pandeo. En este apartado se
describen algunos tipos de colapso por pandeo suponiendo comportamiento elástico
lineal del material y aplicación gradual de la carga.
El problema más allá de la carga de colapso se vuelve dinámico; la curva adquiere una
pendiente descendente, situación en la cual el sistema libera energı́a elástica en forma
de energı́a cinética.
primaria por debajo de la carga crı́tica es no lineal. Una vez alcanzado el punto de
bifurcación la pieza no ideal sigue la solución bifurcada estable en la que disminuye la
carga aplicada, por lo que el la estructura no resiste mayor carga más allá del pandeo.
Por el contrario, en la rama inestable la carga aumenta hasta alcanzar una carga lı́mite
que provoca el colapso de la estructura.
Se considera una pieza recta de una estructura reticulada de plano medio, de sección
transversal A y rı́gidez a flexión EI, ambas uniformes a lo largo de la pieza. La longitud
de la pieza es l y sus extremos se notan como a y b.
donde v(x) es la ley de flechas en la pieza. Nótese que en esta sección se considera
positivo el axil de compresión.
En el caso de la ley de axiles, si el axil N (x) = N es uniforme, el alargamiento de la
pieza se expresa a partir de la relación elástica
Fxb l
δ = ub − ua = (3.103)
EA
y, por tanto, las ecuaciones elásticas de los axiles en los extremos de la pieza son las
conocidas expresiones lineales:
EA EA
Fxa = (ua − ub ) Fxb = − (ua − ub ) (3.104)
l l
En el caso de la deformación por flexión, obtener la ley de flechas v(x) exige integrar
la ecuación de equilibrio
k2
v 00 + k 2 v + (Ma − Fya x − N va ) = 0 (3.105)
N
con las condiciones de apoyo en flechas y giros en los extremos que, de forma general,
se expresan como
v(x)|x=0 = va v(x)|x=l = vb
(3.106)
v 0 (x)|x=0 = θa v 0 (x)|x=l = θb
120 CAPÍTULO 3. INESTABILIDAD
Se define el parámetro ρ = N/NE como la relación entre el axil actuante sobre la pieza
N la carga crı́tica de Euler NE .
La solución de la ecuación (3.105) depende del signo del axil, por lo que se distinguen
dos casos:
+ +
2α 2α
N > 0 (ρ > 0) → v(x) = A+ sin x + B + cos x + vp (x) (3.107a)
l l
− −
− 2α − 2α
N < 0 (ρ < 0) → v(x) = A sinh x + B cosh x + vp (x) (3.107b)
l l
√ √
donde α+ = π2 ρ y α− = π2 −ρ. Las constantes A+ , A− , B + y B − se obtienen
imponiendo las condiciones de apoyo v1 , v2 , θ1 , θ2 en los extremos.
Nótese que para el caso N = 0, es decir, ρ = 0, se recupera el caso lineal.
Por su parte, la solución particular de la ecuación (3.105) es
1
vp (x) = (Fyb x + N va − Ma ) (3.108)
N
Una vez obtenida la solución v(x), empleando la ley de momentos flectores, se pueden
evaluar explı́citamente las relaciones entre los momentos Ma y Mb con los movimientos
en los extremos. Seguidamente, por equilibrio, se determinan también Fya y Fyb . Estos
resultados pueden expresarse matricialmente en la forma
f 0 = K0 d0 (3.109)
2α − sin(2α)
+
C (α) = si ρ > 0
sin(2α) − 2αcos(2α)
C(ρ) = (3.115)
2β − sinh(2β)
C − (β) =
si ρ < 0
sinh(2β) − 2βcosh(2β)
∂K0
K0G |ρ6=0 ≈ K0G |ρ=0 +ρ (3.116)
∂ρ ρ=0
De esta forma la matriz de rigidez para el análisis de segundo orden puede descompo-
nerse en dos matrices:
K0 = K00 + ∆K0 (3.117)
122 CAPÍTULO 3. INESTABILIDAD
donde K00 es la matriz de rigidez del análisis lineal y ∆K0 es una matriz dependiente
del esfuerzo axil N , que explı́citamente se calcula como
∂K0
K0G =ρ (3.118)
∂ρ ρ=0
donde K0G es la denominada matriz geométrica de la pieza. Nótese que las filas y
columnas primera y cuarta de la matriz K0G son nulas, puesto que la componente de
fuerza-desplazamiento en la dirección de la directriz de la pieza no añade términos a la
inestabilidad con respecto al problema lineal.
Idéntico resultado aproximado se puede obtener siguiendo método de Rayleigh-Ritz tal
como se expuso en 3.5.3. El desarrollo completo de este punto puede consultarse en el
apartado A.3.3 del apéndice.
3.6. FORMULACIÓN MATRICIAL DEL PROBLEMA DE INESTABILIDAD 123
Las ecuaciones de equilibrio permiten encontrar las relaciones entre esfuerzos y movimien-
tos en los extremos:
K0 d0 = K00 + K0G d0 = f 0
(3.124)
donde K00 es la matriz de rigidez del problema lineal y K0G es la matriz geométrica de
una pieza articulada, ambas explicitadas a continuación:
1 0 1 0 0 0 0 0
EA 0 0 0 0 N 0 −1 0 1
K00 = K0G = (3.125)
l −1 0 −1 0 l
0 0 0 0
0 0 0 0 0 1 0 −1
f = Tf 0 → f 0 = TT f 0 d = Td0 → d0 = TT d0 (3.126)
K = TK0 TT (3.127)
La matriz de rigidez global de la estructura depende del valor del axil N de cada pieza,
que a su vez depende de los movimientos (d) en los extremos de las piezas. Nótese
que esta dependencia se produce a través de los desplazamientos locales en la dirección
de la pieza indeformada (ua , ub ). Esto significa que la matriz de rigidez depende de la
solución del problema, por lo que el sistema de ecuaciones puede expresarse
K(d)d = f (3.128)
De esta expresión puede obtenerse ∆dk , lo que permite actualizar el vector de despla-
zamientos d de forma iterativa hasta alcanzar una nivel de tolerancia dada.
Es habitual emplear como aproximación inicial d0 = 0, lo que implica que la primera
iteración d1 equivale a la solución del problema lineal. El método de Picard converge
linealmente a la solución con el número de iteraciones, es decir, la función log k dk −
dk−1 k decrece linealmente con k.
Finalmente, la resolución de un análisis de segundo orden sobre una pieza o una estruc-
tura requiere un esquema iterativo en desplazamientos d. Seguidamente se muestra
3.6. FORMULACIÓN MATRICIAL DEL PROBLEMA DE INESTABILIDAD 127
1) Inicializar el iterador k = 0.
3.2) Hallar ∆dk para actualizar la solución en desplazamientos. Esto supone re-
solver un sistema de ecuaciones como el presentado en (3.129) o bien en (3.132).
3.4) Calcular los esfuerzos axiles en las barras a partir de los desplazamientos actu-
alizados dk .
3.5) A partir de los axiles hallados en las barras de la estructura, calcular las ma-
trices de rigidez de cada barra K dk en particular la matriz geométrica, que
Sea una estructura de la que se conoce su geometrı́a y las propieades mecánicas de sus
componentes y sea un sistema de cargas expresado en el sistema global, definido por el
vector f , que no provoca el pandeo de la estructura.
128 CAPÍTULO 3. INESTABILIDAD
(K0 + KG ) d = f (3.133)
Existe un valor del factor de carga λc que da lugar a una matriz de rigidez de la
estructura singular, es decir, para λ = λc se cumple
siendo µc = 1/λc los autovalores del problema redefinido, donde ahora la matriz K0
sı́ es invertible. Los autovalores µc pueden hallarse mediante el método de iteración
vectorial directa o de las potencias. Los factores de carga del problema de inestabilidad
λc son los inversos de los hallados numéricamente µc , mientras que los autovectores
no varı́an. La mı́nima carga crı́tica dada por el factor λc viene dada por el máximo
autovalor µc y su autovector asociado.
Es preciso remarcar que la carga crı́tica en una estructura se obtiene como un factor
multiplicador de carga respecto de un sistema de cargas de referencia, y se obtiene
suponiendo dos aproximaciones: la linealización de la matriz rigidez de segundo orden
y la relación lineal de la propia matriz geométrica con respecto a los desplazamientos.
Puede demostrarse que, en general, el resultado obtenido mediante este procedimiento
supone una sobrestimación de la carga crı́tica, quedando por tanto del lado de la inse-
guridad.
2. Fijar k = 0.
3. Calcular d0 = −KG0 δ 0 .
4. Actualizar k = k + 1.
130 CAPÍTULO 3. INESTABILIDAD
5. Resolver K0 δ k = dk−1 .
La carga crı́tica de la viga depende del estado de cargas. En este caso la viga trabaja
exclusivamente a esfuerzo axil N , con N = P . El análisis de inestabilidad de la viga
se lleva a cabo a través de las matrices de rigidez y goemétrica de la estructura. Si
3.6. FORMULACIÓN MATRICIAL DEL PROBLEMA DE INESTABILIDAD 131
se considera la viga formada por una única barra la matriz de rigidez global de la
estructura es
EA
− EA
l 0 0 l 0 0
12EI
0
l3
− N 5l6 6EI
l2
1
− N 10 0 − l3 + N 5l6
12EI 6EI
l2
1
− N 10
6EI 1 4EI 2l
0
l2
− N 10 l − N 15 0 − 6EI
l2
+ N 10 1 2EI
l + N 30
l
K = EA
− l EA
0 0 l 0 0
− 12EI + N 5l6 − 6EI 1 12EI
N 5l6 − 6EI 1
0
l3 l2
+ N 10 0 l3
− l2
+ N 10
6EI 1 2EI l
0 l2
− N 10 l + N 30 0 − 6EI
l2
+ N 10 1 4EI
l − N 15
2l
30EI
Ncr = = 3, 040NE
l2
Una vez determinadas las matrices locales y realizadas las operaciones estándar de
ensamblaje y aplicación de las condiciones de apoyo, la matriz reducida del sistema es
de orden 5. El sistema de ecuaciones a resolver es el siguiente:
EA
− EA
l 0 0 l 0 u2 0
192EI
0
l3
− N 24
5l 0 0 24EI
l2
− N 101
ν2
0
16EI 2l 4EI l
0 0 − N 15 0 l + N 60 θ2 = 0
EA l
− l EA
0 0 l 0 u3 −P
24EI 1 4EI l 8EI l
0 l2
− N 10 l + N 60 0 l − N 15 θ3 0
Nótese que a pesar de que el orden de la matriz reducida del sistema es 5, el poli-
nomio caracterı́stico es de orden 3 ya que dos de los grados de libertad corresponden
al movimiento horizontal en los nudos 2 y 3. Existen por tanto tres soluciones para el
problema de autovalores:
EI EI EI
Ncr,1 = 20, 7 Ncr,2 = 75, 1 Ncr,3 = 197, 5
l2 l2 l2
EI
Ncr = 20, 7 = 2, 098NE
l2
Al subdividir la pieza original en dos barras se aumentan los grados de libertad del
problema y en consecuencia la flexibilidad de la función de aproximación lo que deriva
en una aproximación más precisa a la solución exacta y consecuentemente a la carga
crı́tica de la estructura.
Sea un pórtico formado por dos vigas con dimensiones según muestra la Figura 3.29,
con sección transversal A e inercia I respecto al eje de flexión uniformes y módulo
resistente E. Dicho pórtico está sometido a una carga V vertical en su vértice superior
y a una carga H horizontal en el apoyo derecho, con la relación V = 2, 5H.
Realizar un análisis lineal de inestabilidad de la estructura mediante formulación ma-
tricial para determinar su carga crı́tica.
Datos: h = 5 m, α = 30◦ , EA = 106 kN , EI = 104 kN m2
Como paso previo a la solución del problema lineal se obtiene la solución del problema
lineal para las cargas de referencia f0 , lo que permite obtener los axiles en ambas barras
para dicha carga de referencia.
El ensamblaje de las matrices elementales da lugar a la matriz de rigidez lineal de la
estructura. El pórtico está constituido por tres nudos con tres grados de libertad por
nodo, luego ambas matrices de rigidez son cuadradas y simétricas de orden 9.
La matriz de rigidez K0 de la estructura no depende de las cargas, ni por tanto de los
axiles en las barras, y por tanto no varı́a en todo el procedimiento:
960 0 −2400 −960 0 −2400 0 0 0
0 200000 0 0 −200000 0 0 0 0
−2400 0 8000 2400 0 4000 0 0 0
−960 0 2400 75990 −43249 2700 −75030 43249 300
K0 =
0 −200000 0 −43249 225090 520 43249 −25090 520
−2400 0 4000 2700 520 12000 −300 −520 2000
0 0 0 −75030 43249 −300 75030 −43249 −300
−25090 −520 −43249 −520
0 0 0 43249 25090
0 0 0 300 520 2000 −300 −520 4000
Los axiles obtenidos en la solución del problema lineal son los siguientes:
N1 = 2, 10 kN (→←) N2 = 1, 07 kN (→←)
Estos axiles determinan las matrices geométricas elementales de ambas barras. Me-
diante las operaciones de rotación de ejes y ensamblaje apropiadas, se obtiene la matriz
geométrica KG de la estructura, para el estado de cargas de referencia f0 , que incorpora
los axiles N1 y N2 obtenidos en ambas barras:
−0, 5030 0 0, 2096 0, 5030 −0, 0000 0, 2096 0 0 0
0 0 0 0 0 0 0 0 0
0, 2096 0 −1, 3972 −0, 2096 0, 0000 0, 3493 0 0 0
0, 5030 0 −0, 2096 −0, 5350 −0, 0555 −0, 2630 0, 0320 0, 0555 −0, 0534
KG0 =
0 0 0 −0, 0555 −0, 0961 −0, 0925 0, 0555 0, 0961 −0, 0925
0, 2096 0 0, 3493 −0, 2630 −0, 0925 −2, 8214 0, 0534 0, 0925 0, 3560
0 0 0 0, 0320 0, 0555 0, 0534 −0, 0320 −0, 0555 0, 0534
−0, 0555 −0, 0961
0 0 0 0, 0555 0, 0961 0, 0925 0, 0925
0 0 0 −0, 0534 −0, 0925 0, 3560 0, 0534 0, 0925 −1, 4242
Una vez aplicadas las condiciones de apoyo se obtiene el sistema de ecuaciones reducido
que de forma sintética adquiere la forma
h i
K̂0 + λK̂G0 d̂ = λf̂0
Problema de autovalores
Los factores de carga hallados permiten establecer los estados de carga que provocan
la inestabilidad de la estructura. Asimismo, permiten aproximar los axiles en ambas
barras en pandeo. En la Tabla 3.3 se detallan todas estas magnitudes expresdas en
(kN):
El estado de carga crı́tico que provoca el pandeo de la estructura está formado por una
carga horizonal H = 698, 7 kN y una carga vertical V = 1746, 9 kN .
Merece la pena analizar los axiles encontrados en el estado crı́tico con respecto los
valores de la carga crı́tica en los casos simples estudiados en este capı́tulo. La Tabla
3.4 muestra las cargas crı́ticas de Euler de cada barra (NE,i ) y las cargas crı́ticas de
cada barra atendiendo a sus condiciones de apoyo (Ncr,i ). En este caso, la barra (1)
se aproxima a una viga empotrado-articulada mientras que la barra (2) se aproxima a
una viga biarticulada. Esta comparación ofrece una medida de la influencia tanto del
grado de empotramiento en los nudos como de la influencia del resto de la estructura
sobre cada barra.
3.6. FORMULACIÓN MATRICIAL DEL PROBLEMA DE INESTABILIDAD 137
Tabla 3.4: Análisis de bifurcación. Axiles crı́ticos en relación a los casos simples
4.1 Introducción
En las estructuras dimensionadas por métodos elásticos, las tensiones debidas a las
cargas mayoradas no sobrepasan en ningún punto de la estructura el lı́mite elástico del
material. Este proceder tiene dos serios inconvenientes.
Por un lado, el análisis elástico es muy sensible a las imperfecciones estructurales (por
ejemplo, apoyos y enlaces no perfectamente rı́gidos o elásticos) o a acciones de difı́cil
determinación (variaciones térmicas, descensos de apoyos, etc). Por otro lado, la su-
peración local del lı́mite elástico no suele llevar al colapso estructural; bien al contrario,
la ruina de la estructura se da generalmente para cargas muy superiores a las cargas
últimas obtenidas en régimen elástico.
Ambos inconvenientes son tanto más serios cuanto mayor es el grado de hiperestatismo
de la estructura, y cuanto más dúctil es el material empleado. Si el material conserva
la capacidad de transmitir tensiones más allá del lı́mite elástico, la estructura puede
soportar cargas superiores gracias a un proceso de redistribución de tensiones; este
proceso tiende a ser tanto más extensivo cuanto mayor es el grado de indeterminación
estática de la estructura.
Por estas consideraciones surgen los métodos de análisis elasto-plástico y análisis lı́mite
plástico objeto de este Capı́tulo en el que se trata de determinar la carga última de
la estructura y el mecanismo de colapso asociado considerando un comportamiento
elasto-plástico del material.
139
140 CAPÍTULO 4. ANÁLISIS PLÁSTICO
• Tramo elástico no lineal (AB): En este tramo la respuesta del material conserva
su reversibilidad - elasticidad - por lo que ante una descarga, la trayectoria tenso-
deformacional es la misma que en carga, llegando en último término hasta la
4.2. RELACIÓN TENSO-DEFORMACIONAL. COMPORTAMIENTO PLÁSTICO 141
La diversidad de materiales existentes hace que este esquema sea general, por lo que
cada material puede presentar algunos de los comportamientos antes descritos. En
ingenierı́a es relevante la clasificación de los materiales en función de su comportamiento
en plasticidad. De esta forma se habla de materiales dúctiles y frágiles.
Los materiales dúctiles se caracterizan por un amplio rango plástico, alcanzando la ro-
142 CAPÍTULO 4. ANÁLISIS PLÁSTICO
tura después de un nivel de deformación muy elevado. Esto suele deberse a la existencia
de un tramo de fluencia en el cual la deformación crece a tensión constante. Ası́, en este
tipo de materiales εD ≈ 20εA y εF ≈ 200εA . En estos materiales en general pueden
observarse todos los tramos especificados. Son materiales dúctiles el acero o el cobre,
y en general los metales.
Por el contrario, en los materiales frágiles los tramos elástico no lineal (AB) y de
fluencia (CD) son mucho más pequeños o no se encuentran bien definidos, por lo
que habitualmente no se consideran. La rotura del material por tanto, se produce a
deformaciones notablemente menores que en el caso de materiales dúctiles εF ≈ 5εA
por lo que la estructura “no avisa” antes de colapsar. Por convenio, en estos materiales
se toma como tensión de fluencia aquella que da lugar a una deformación permanente
del 0, 2%. Son materiales frágiles el hormigón, el vidrio o la fundición.
En la Figura 4.2 se muestran las curvas tensión-deformación caracterı́sticas para mate-
riales dúctiles y frágiles.
con lı́mites resistentes más altos y más frágiles a tracción con resistencias del orden de
un 10% más bajas.
Las estructuras articuladas son estructuras de barras enlazadas entre sı́ mediante ar-
ticulaciones y que trabajan fundamentalmente sometidas a esfuerzo axil.
En las estructuras articuladas con cargas en los nudos, el axil es uniforme en cada una
de las barras. Esto da lugar a una respuesta estructural en la cual la plastificación
se da en cada barra, en todas las secciones y en todos los puntos de la sección en el
mismo instante, por lo que el mecanismo de colapso de la estructura se traduce en una
secuencia de plastificación referida al orden en que las barras alcanzan la tensión de
fluencia σF o equivalentemente el axil de plastificación.
Cuando una de las barras alcanza el esfuerzo axil de plastificación, el axil en dicha barra
es fijo y conocido, es decir, deja de ser una incógnita en las ecuaciones de equilibrio,
reduciendo ası́ el grado de indeterminación estática. El colapso de la estructura se
produce cuando la estructura ha reducido su indeterminación estática lo suficiente
como para convertirse en un mecanismo, es decir, que la estructura se encuentra en
un estadio en el que no puede garantizar el equilibrio estable, por lo que se produce el
colapso.
Notación: En este capı́tulo se notan las cargas últimas de colapso de una estructura
como Pu o bien qu según la naturaleza puntual P o distribuida q de la carga. El valor de
carga que da lugar a una plastificación parcial de la estructura sin alcanzar el colapso
de la misma se nota como Pp(i) siendo i el paso asociado a dicha carga.
La estructura articulada de la Figura 4.4 está formada por tres barras, una barra
vertical y dos barras inclinadas un ángulo α = 45◦ . La estructura está sometida a una
carga vertical de valor P aplicada en el vértice en el que confluyen las tres barras.
En el rango elástico, todos los axiles son proporcionales a la carga P , al igual que las
tensiones.
La plastificación de una barra se da cuando la tensión en un punto de la sección alcanza
el valor de la tensión de fluencia σF o equivalentemente el axil alcanza el valor del axil
de plastificación NP = AσF . Al ser el axil uniforme en toda la barra, el proceso de
plastificación se da en toda la barras según la magnitud del axil.
La barra (2) plastifica cuando el axil alcanza el valor de plastificación NP lo que sucede
para un valor de la carga P calculado como
√ !
2P 2
N2 = √ = NP ⇒ Pp1 = 1 + NP
2+ 2 2
Para un valor de la carga Pp1 los axiles en las barras de la estructura son
1
N1,1 = N3,1 = NP N2,1 = NP
2
146 CAPÍTULO 4. ANÁLISIS PLÁSTICO
Los esfuerzos totales se calculan sumando los obtenidos en el paso 1 a los incrementos
calculados
1 ∆P
N1 = N3 = NP + √ N2 = NP
2 2
La plastificación de las barras (1) y (3) se produce cuando el axil en estas barras alcanza
el valor de plastificación.
√
1 ∆P 2
N1 = N3 = NP + √ ⇒ ∆P = NP
2 2 2
4.3. PLASTIFICACIÓN EN ESTRUCTURAS ARTICULADAS 147
Alcanzada la carga Pp2 , la estructura no puede admitir más carga ya que todas las
barras han plastificado, lo que supone el colapso de la misma. Por tanto, la carga de
colapso de la estructura es Pu = Pp2 .
Las barras (1) y (3) evolucionan en régimen elástico durante todo el proceso, por lo
que su alargamiento se calcula a través de la relación elástica en cada paso.
√ √
N1,1 2h 2NP h
δ1,1 = δ3,1 = =
EA√ √2EA
N1,2 2h 2NP h
δ1,2 = δ3,2 = =
EA EA
La barra (2), plastificada en primer lugar, se mueve por el rango elástico hasta el
instante 1 y en el rango plástico entre los instantes 1 y 2. En el primer intervalo, el
alargamiento de la barra (2) se calcula mediante la relación elástica.
N2,1 h NP h
δ2,1 = =
EA EA
2NP h
δ2,2 =
EA
En la Figura 4.6 se muestran las trayectorias del axil normalizado η = N/NP en todas
las barras, todas ellas traccionadas (a), en función de la carga normalizada al axil de
plastificación ρ = P/NP y (b) en función del alargamiento normalizado ν = δ/δ0 hasta
el colapso de la estructura, con δ0 = NP h/EA.
148 CAPÍTULO 4. ANÁLISIS PLÁSTICO
En la figura se aprecia cómo la barra (2), que plastifica en primer lugar, pierde la
proporcionalidad entre axil y alargamiento, mientras que las barras (1) y (3) mantienen
la linealidad durante todo el proceso de carga.
El marco articulado de la Figura 4.7 está formada por cinco barras de sección constante
A, dos laterales verticales, una superior horizontal y dos diagonales que se cruzan sin
cortarse, con una inclinación α = 45◦ . La estructura soporta una carga horizontal P
aplicada en el el nudo superior izquierdo.
Siendo todas las barras de igual sección transversal A y del mismo material, la plasti-
ficación de cualquiera de ellas se alcanza para un axil de plastificación
NP = AσF = 150 kN
En el rango lineal del problema todos los axiles conservan la proporcionalidad con
respecto a la carga. El máximo axil se da en la barra (5), que se encuentra comprimida.
El valor de la carga que ocasiona la plastificación de la barra (5) es
23
N5 = −NP ⇒ Pp1 = √ NP = 1, 268NP
4 + 10 2
150 CAPÍTULO 4. ANÁLISIS PLÁSTICO
Para la carga de plastificación Pp1 los axiles en todas las barras son:
√
1 2−3 2
N1,1 = √ NP = 0, 707NP N2,1 = N3,1 = NP = −0, 561NP
2 4
√
3− 2
N4,1 = NP = 0, 793NP N5,1 = −NP
2
Una vez que la barra (5) ha plastificado, el valor del axil en dicha barra permanece
constante con valor NP comprimiendo la barra. Al conocerse el axil N5 el grado de
hiperestatismo se reduce, por lo que el problema resultante es isostático.
Paso 2: Una vez que se produce la primera plastificación el problema pierde su linea-
lidad, es decir, la relación entre axiles y alargamientos no viene dada por la expresión
obtenida de la elasticidad. Sı́ se conservan las relaciones de equilibrio y compatibilidad
de movimientos, en particular las relaciones entre alargamientos (locales) y movimientos
(globales).
√
2−1
∆P = √ NP
2 2
Para este incremento de carga los axiles en las barras en el paso 2 son
NP NP
N1,2 = √ N2,2 = N3,2 = − √ N4,2 = NP N5,2 = −NP
2 2
Una vez han plastificado las barras (5) y (4) la estructura se convierte en un mecanismo
hipoestático y no puede soportar más carga. La carga de colapso será entonces
√
Pu = Pp1 + ∆P = 2NP = 212, 14 kN
La Figura 4.9 muestra las trayectorias del axil normalizado η = N/NP en todas las
barras (a) en función de la carga normalizada al axil de plastificación ρ = P/NP y (b)
en función del alargamiento normalizado ν = δ/δ0 , hasta el colapso de la estructura.
152 CAPÍTULO 4. ANÁLISIS PLÁSTICO
Por su parte, en la Figura 4.9b se puede apreciar la pérdida de linealidad entre axil y
deformación en la barra (5), linealidad que sı́ se conserva en el resto de barras. Cuando
la barra (4) plastifica, se produce el colapso estructural.
Este ejemplo pone de manifiesto cómo la resolución del problema plástico, esencial-
mente no lineal, se afronta mediante la solución de una secuencia de problemas lineales
incrementales con los condicionantes estructurales adecuados.
según la cual las secciones transversales de una viga permanecen planas y normales a
la directriz de la viga deformada.
En flexión simple recta, el eje neutro pasa por el centro de gravedad y es perpendicular
al plano de solicitación a flexión. Para un momento flector Mz la distribución de
tensiones y deformaciones en régimen elástico es lineal con la coordenada y que mide
la distancia al eje neutro,
Mz
σz = − y (4.1)
Iz
ME = W z σ F (4.2)
En un análisis elastoplástico las fibras que alcanzan este nivel tensional no pueden
incrementarlo y entran en la rama plástica de la curva de comportamiento. Sin embargo,
las fibras con un nivel inferior de tensión sı́ pueden aumentarla hasta alcanzar la tensión
lı́mite. Ası́, el momento plástico de la sección respecto al eje z es
siendo mez (Ω̂) el momento estático de la parte de la sección por encima del eje z de
flexión.
La sección se considera agotada cuando todas las fibras alcanzan la tensión de fluencia
σF . El momento flector asociado a esta distribución de tensiones se denomina mo-
mento plástico (MP ). La Figura 4.10 muestra la evolución de la sección hasta su total
plastificción.
Figura 4.10: Evolución de las tensiones en una sección durante la plastificación por flexión
154 CAPÍTULO 4. ANÁLISIS PLÁSTICO
El factor de forma f es una medida del intervalo de momento flector entre el comienzo
de la plastificación y el agotamiento de la sección.
La Figura 4.11 muestra los factores de forma para las secciones más habituales.
Sea una sección rectangular de ancho b y altura h como muestra la Figura 4.13, que
soporta un momento flector en dirección z.
El área de la sección es A = bh, la inercia con respecto al eje horizontal es Iz = bh3 /12,
mientras que la coordenada de su fibra más alejada del eje de flexión es |y max | = h/2.
El momento elástico de la sección es el que produce la tensión de fluencia σF en el
punto de máxima tensión y se calcula como
Mz max Mz bh2
σ= y = = σF ⇒ ME = Wz σF = σF
Iz Wz 6
Sea una sección sección romboidal de la que se conocen sus diagonales, horizontal b y
vertical h sometida a un momento flector en dirección z.
El área de la sección es A = bh/2, la inercia con respecto al eje horizontal es Iz = bh3 /48,
mientras que la coordenada de su fibra más alejada del eje de flexión es |y max | = h/2.
bh2
ME = Wz σF = σF
24
Z h
bh2
2 2σF 2y
ME = 2 y·y·b 1− dy = σF
0 h h 24
Sea una sección en corona con radio exterior R y radio interior r, siendo el espesor
e = R − r, sometida a un momento flector en dirección z.
El área de la sección, la coordenada de su fibra más alejada del eje de flexión es y max
y la inercia de la sección respecto de cualquier eje son
πR4
r 2 r 4
2 max
A = πR 1 − |y |=R I= 1−
R 4 R
Para el caso particular de una sección circular maciza, esto es r = 0, los momentos
elástico, plástico y el factor de forma se particularizan en
πR3 4 16
ME = σF MP = R3 σF f= = 1, 70
4 3 3π
Sea una sección sección en doble T de altura h, ancho de las alas b y espesor e, sometida
a un momento flector en dirección z.
El área de la sección es A = (2b + h)e, la coordenada de la fibra más alejada del eje
neutro es |y max | = h/2 y la inercia de la sección respecto de su eje principal de inercia
es
h3 e be3 beh2 h2 e
Iz = + + ≈ (h + 6b)
12 6 2 12
El momento elástico de la sección puede hallarse mediante la forma general
h2 e
6b
M E = Wz σF = 1+ σF
6 h
Por su parte el momento plástico se calcula integrando el momento que producen las
tensiones normales σ = σF
h2 e
h h h 4b
MP = 2 σF eb + σF e = 1+ σF
2 2 4 4 h
Para el caso particular en que la relación canto/ancho sea (b/h = 1/2) se tiene un factor
de forma de f = 1, 125.
Sea una sección sección en doble T asimétrica con respecto al eje normal a su alma. La
sección tiene una altura h, el ancho de las alas es b y b0 para el ala inferior y superior
respectivamente, mientras que el espesor e es uniforme en toda la sección. La sección
se encuentra sometida a un momento flector en dirección z.
Las caracterı́sticas geométricas de las sección, área, posición del dentro de gravedad
y distancia desde el centro de gravedad a la fibra superior más alejada (y max ) son las
siguientes:
h2 h2
0 b0 h + 2 max bh + 2
A = e(b + b + h) yG = y =
b+ b0 + h2 b + b0 + h
2
160 CAPÍTULO 4. ANÁLISIS PLÁSTICO
La fibra inferior más alejada se sitúa a una distancia h − y max . Si b > b0 , el centro de
gravedad se sitúa por debajo de la mitad de la altura, por lo que la tensión máxima
se da en la fibra superior con un valor igual a la tensión de fluencia σ sup = σF . En la
fibra inferior, la tensión correspondiente es
h
σ inf = σF
y max −1
Cuando una sección está sometida a flexión compuesta, es decir, a la acción simultánea
de un axil N y un momento flector M , la combinación de esfuerzos determina el valor
lı́mite de los esfuerzos que provocan la plastificación total de la sección.
Se plantea en este apartado el estudio de la influencia del axil en la plastificación de
la sección;los casos de axil puro y flexión simple son casos lı́mite de ésta. En la Figura
4.18 se muestran las distribuciones de tensión en plastificación para los tres casos: axil
puro, flexión compuesta y flexión simple.
A modo de ejemplo, consideremos el caso de una sección rectangular de dimensiones
hxb. Los resultados obtenidos son extrapolables para otras secciones.
En una sección rectangular el axil de plastificación NP y el momento plástico MP son
bh2
NP = bhσF MP = σF (4.5)
4
En flexión compuesta el eje neutro no pasa por el centro de gravedad sino que se sitúa
a una distancia y0 según se aprecia en la Figura 4.18b. A partir de las distribuciones
4.4. PLASTIFICACIÓN EN ESTRUCTURAS RETICULADAS 161
y0 h NP,F C
NP,F C = 2b NP ⇒ y0 = (4.7)
h 2 NP
Por tanto, la actuación del axil reduce el momento flector que plastifica la sección en
la medida que el axil actuante se aproxima al axil de plastificación NP . La Figura 4.19
muestra la influencia del axil sobre el momento flector en el estado de plastificación
MP,F C para una sección rectangular.
Esto significa que en una sección en doble T sometida a un axil menor al 15% del axil
de plastificación el momento de plastificación puede ser considerado en su totalidad.
162 CAPÍTULO 4. ANÁLISIS PLÁSTICO
Esto es debido a que el alma de la sección soporta parte de la carga axial mientras que
su contribución al momento plástico es mı́nima.
El esfuerzo cortante puede igualmente reducir el valor del momento plástico de una
sección en algunos casos. En caso de considerar la interacción con el cortante se requiere
una superficie de interacción tridimensional.
El momento plástico de una sección marca el momento flector máximo que puede
soportar dicha sección una vez se ha producido la plastificación de todas sus fibras.
Para momentos mayores al momento plástico, todas las fibras fluyen sin resistencia, lo
que implica el giro libre de dicha sección. Esta situación se traduce en el concepto de
rótula plástica necesario para determinar el mecanismo de colapso y la carga última en
vigas y en pórticos.
En el caso de los pórticos reticulados, los esfuerzos axiles reducen el momento de plas-
tificación de las secciones, aunque su influencia suele despreciarse en la práctica.
Las rótulas plásticas van apareciendo en diferentes secciones hasta que se produce el
colapso plástico de la estructura. Este colapso se produce cuando el equilibrio estático
no puede garantizarse. Se denomina mecanismo de colapso a la configuración deformada
de la estructura en la cual no se garantiza el equilibrio por carecer de las suficientes
coacciones internas y/o externas (estado hipoestático o mecanismo). Se denomina carga
última al valor de la carga que provoca el colapso de la estructura.
El proceso de formación de las n+1 rótulas plásticas es un proceso secuencial que puede
determinarse analı́ticamente hasta llegar a la formación del mecanismo de colapso. En
cada uno de los pasos, el momento flector en las secciones en las que ya se ha alcanzado
el momento plástico MP no puede aumentar, mientras que sı́ puede hacerlo en las
demás secciones de la estructura. A este proceso se le denomina redistribución plástica
de esfuerzos. Este proceso es no lineal, es decir, las leyes de esfuerzos que se producen
en los diferentes estadios que conducen al mecanismo de colapso no son proporcionales
a la carga actuante.
• Dicha sección se convierte en una rótula plástica sobre la cual actúa un momento
plástico constante MP , que pasa a tratarse como un valor conocido de la ley de
flectores. La rótula disminuye en uno el grado de hiperestatismo de la estructura.
rótulas plásticas.
Sea una viga biarticulada de longitud l cuya sección transversal tiene un momento
plástico MP , sometida a una carga puntual P en su punto medio.
Conocido el valor del momento máximo y el valor del momento plástico de la sección
puede calcularse la carga de colapso de la viga:
Pl 4MP
M max = = MP ⇒ P u =
4 l
Figura 4.21: Mecanismo de colapso en una viga biarticulada con carga puntual
Sea una viga biarticulada sometida a una carga distribuida de valor q cuya sección
transversal tiene un momento plástico MP .
ql
M (x) = x(l − x) 0≤x≤l
2
Conocido el valor del momento máximo y el valor del momento plástico de la sección,
puede calcularse la carga de colapso de la viga.
ql2 8MP
M max = = MP ⇒ q u = 2
8 l
Figura 4.23: Mecanismo de colapso de una viga biarticulada con carga distribuida
Sea una viga biarticulada formada por un vano de longitud l y un voladizo de longitud
αl con 0 ≤ α ≤ 1 sometida a una carga repartida q en toda su longitud. Considerar el
caso particular α = 0.25. La viga se encuentra sometida a una carga distribuida q en
toda su longitud.
ql ql
V1 = (1 − α2 ) V2 = (1 + α)2
2 2
ql (1 − α2 )x − qx2 0 ≤ x ≤ l
2 2
M (x) =
q [(1 + α)l − x]2 l ≤ x ≤ (1 + α)l
2
ql2 l
M max,(+) = (1 − α2 )2 en xmax = (1 − α2 )
8 2
ql2 α2
M max,(−) = en x = l
2
Ambos valores dependen de la longitud relativa del voladizo expresada en función del
parámetro α. Existen dos situaciones de colapso posibles en función de la magnitud
relativa de ambos máximos, según se muestra en la Figura 4.25:
168 CAPÍTULO 4. ANÁLISIS PLÁSTICO
√
Los valores de α que igualan ambos momentos máximos son α0 = 2 ± 1. La situación
del máximo positivo M max,(+) está limitada por la condición |α| < 1. El valor α =
√
2 + 1 da lugar a una ley de momentos flectores cuyo máximo flector positivo se sitúa
√
fuera de la viga. De modo que el único valor de α con sentido fı́sico es α = α0 = 2 − 1.
Para valores de α < α0 se tiene que M max,(+) > M max,(−) y por tanto el valor del
momento plástico se alcanza en el vano, siendo la carga de colapso:
ql2 8 MP
M max,(+) = (1 − α2 )2 = MP ⇒ qu =
8 (1 − α ) l2
2 2
Para valores de α > α0 se tiene que M max,(−) > M max,(+) y por tanto, el valor del
momento plástico se alcanza antes en la sección del apoyo, y la carga de colapso es:
ql2 α2 2 MP
M max,(−) = = MP ⇒ q u = 2 2
2 α l
En el caso lı́mite en el cual α = α0 se producen dos rótulas plásticas simultáneamente,
dando lugar a un mecanismo de colapso supercompleto.
El caso particular de α = 0, 25 se encuentra en el primer caso (α < α0 ), por lo que
el colapso se produce por la aparición de una rótula plástica en el vano. La carga de
colapso es
2048 MP 15
qu = 2
en x = l
225 l 32
4.4. PLASTIFICACIÓN EN ESTRUCTURAS RETICULADAS 169
Nótese que la rótula plástica no se produce exactamente en el centro del vano sino
ligeramente a la izquierda de este punto.
ql2 8MP
M max = = MP ⇒ qp1 = 2
8 l
Paso 2: En un segundo paso se resuelve la viga fijando el valor del momento flector
en el empotramiento, ya que éste no puede superar el valor del momento plástico MP .
Nótese que la rótula plástica fija el valor del momento flector pero permite el giro lo que
lleva a un nuevo esquema estructural equivalente formado por una viga biarticulada
con un momento MP aplicado en el extremo que ha plastificado. La ley de momentos
flectores en este nuevo esquema según muestra la Figura 4.28(2) se expresa como:
qx2
ql MP
M (x) = −MP + + x− 0≤x≤l
2 l 2
170 CAPÍTULO 4. ANÁLISIS PLÁSTICO
El nuevo momento flector máximo se sitúa allı́ donde se anula la derivada primera de
la ley de momentos, es decir
dM (x) ql MP l MP
= + − qx = 0 ⇒ xmax = +
dx 2 l 2 ql
Nuevamente se calcula el valor de la carga repartida q para alcanzar el valor del mo-
mento plástico en la sección xmax :
√ MP qp2a = 0, 34 Ml2P
max
M = MP ⇒ qp2 = (6 ± 4 2) 2 →
l qp2b = 11, 65 Ml2P
La primera de las soluciones es menor que la carga obtenida en el primer paso qp1 , por
lo que la carga de colapso de la viga es qu = qp2b .
El mecanismo de colapso de una viga empotrado-articulada se muestra en la Figura
4.29.
Figura 4.29: Mecanismo de colapso en una viga empotrado-articulada bajo carga distribuida
4.4. PLASTIFICACIÓN EN ESTRUCTURAS RETICULADAS 171
Sea una viga biempotrada sometida a dos cargas puntuales iguales de valor P situadas
de forma simétrica, cuya sección transversal presenta un momento plástico MP .
2P l 9MP
M max = = MP ⇒ Pp1 =
9 2l
Este valor de la carga no produce el colapso de la viga, sino la aparición de una rótula
plástica en cada uno de los empotramientos (A,D), en los cuales el valor del momento
flector es MP .
mecanismo de colapso está conformado por cuatro rótulas y la carga asociada a este
escenario es
3MP 6MP
M max = 0, 5MP + ∆P l/3 = MP ⇒ ∆P = ⇒ Pp2 =
2l l
Esta configuración provoca el colapso de la estructura por lo que la carga última para
este problema es Pu = Pp2 .
Figura 4.32: Mecanismo de colapso en una viga biempotrada bajo cargas puntuales simétricas
4.4. PLASTIFICACIÓN EN ESTRUCTURAS RETICULADAS 173
Sea una viga biempotrada sometida a dos cargas puntuales P y 2P situadas según se
indica en la Figura 4.33, cuya sección transversal presenta un momento plástico MP .
La estructura es dos veces hiperestática por lo que requieren nr = 3 rótulas para llegar
al colapso. Existen cuatro posibles rótulas plásticas, dos en los empotramientos (A,D)
y dos más en los puntos de aplicación de las cargas (B,C).
La resolución del problema elástico da lugar a una ley de flectores como se muestra en
la Figura 4.34, en la cual el momento máximo se da en el empotramiento derecho, lo
que permite determinar la carga necesaria para producir la pimera rótula plástica.
30P l 27MP MP
M max = = MP ⇒ Pp1 = = 2, 7
81 10l l
Tras la formación de la primera rótula la viga pasa comportarse como una viga empotrado-
articulada con un valor fijo del momento flector MP en el extremo derecho (D). Al haber
plastificado una sección, la proporcionalidad entre carga y esfuerzo flector se pierde. El
comportamiento de la viga en el rango plástico se determina superponiendo a la dis-
tribución de flectores en el instante de la plastificación de (D), un estado incremental
que tiene en cuenta la variación del momento flector en la viga sin incrementar su valor
en la sección plastificada. En la nueva configuración el valor máximo del momento se
da en el empotramiento izquierdo (A). El valor de la carga que ocasiona la segunda
rótula plástica en dicha sección es Pp2 :
MP
Pu = Pp3 = 3, 6
l
En la Figura 4.34 se muestran las sucesivas leyes de momentos flectores a medida que
se forman las rótulas, en las que se aprecian los valores de los momentos flectores en
los puntos crı́ticos, hasta el colapso final de la viga.
Figura 4.35: Mecanismo de colapso en una viga biempotrada bajo cargas puntuales
asimétricas
Sea una viga continua formada por dos vanos iguales de longitud l, sometida a dos
cargas puntuales P y 2P según se indica en la Figura 4.36. Los momentos plásticos de
las secciones son 2MP y MP en el primer y segundo vano respectivamente, incluyendo
éste la sección en el nudo 2.
Paso 1: La resolución del problema elástico da lugar a una ley de momentos flectores
en tramos lineales, como se muestra en la Figura 4.37(0), en la cual el máximo se da
en el punto de aplicación de la carga mayor (A). Sin embargo, al ser los momentos
plásticos diferentes en ambos vanos no se garantiza que la primera rótula se produzca
en la sección de momento máximo. De hecho:
23P l 128MP
MA = = 2MP ⇒ Pp1,A =
64 23l
18P l 64MP
MB = = MP ⇒ Pp1,B =
64 18l
176 CAPÍTULO 4. ANÁLISIS PLÁSTICO
En consecuencia, la primera rótula se produce en el apoyo central (B) para una carga
64MP
Pp1 =
18l
La ley de flectores para esta carga se muestra en 4.37(1). Esta carga no produce el
colapso de la viga continua.
23 l 13MP
MA = MP + ∆P = 2MP ⇒ ∆PA =
18 2 9l
7 l 22MP
MC = MP + ∆P = MP ⇒ ∆PC =
18 4 9l
5MP
Pu = Pp1 + ∆P =
l
1) Se define la geometrı́a, las barras con sus propiedades mecánicas y las condiciones
de apoyo.
5) Se ejecuta un bucle que se repite tantas veces como sea necesario hasta llegar al
colapso de la estructura.
5.2) Se define el estado incremental (∆k) como el problema estructural que incluye
las k rótulas plásticas en las secciones en las que ya se ha alcanzado el momento
plástico en los pasos anteriores. En este estado incremental las rótulas plásticas
se consideran rótulas de momento nulo ya que la sección correspondiente no
admite incremento de momento flector.
5.3) Se resuelve dicho estado incremental como problema elástico sometido a las
fuerzas de referencia {P0 } obteniendo los esfuerzos incrementales {∆E}k co-
rrespondientes.
5.4) Los esfuerzos totales se calculan superponiendo los esfuerzos {E}k asociados
4.5. IMPLEMENTACIÓN COMPUTACIONAL DEL ANÁLISIS ELASTO-PLÁSTICO 179
La Figura 4.39 muestra un pórtico plano formado por dos pilares de longitudes l y 3l/4
unidos por un dintel de longitud l. La sección transversal es igual para todas las barras
y se caracteriza por un momento plástico de valor MP . La estructura está sometida a
una carga horizontal P aplicada en el nudo 2 y otra carga puntual de valor 4P aplicada
en el centro del dintel.
Determinar mediante el cálculo elasto-plástico, el mecanismo de colapso y la carga
última de la estructura.
La estructura presenta un grado de hiperestatismo de n = 2, por lo que se requieren
nr = 3 rótulas para alcanzar el colapso. Las posibles rótulas plásticas se sitúan en el
empotramiento (A), en los nudos rı́gidos (B,D) y en el punto de aplicación de la carga
4P (C), lo que supone un total de 4 posibles rótulas plásticas para formar el mecanismo
de colapso.
180 CAPÍTULO 4. ANÁLISIS PLÁSTICO
El sistema de cargas de referencia es el formado por las cargas {P0 , 4P0 } = {1, 4} con
P0 = 1 carga de referencia. El incremento de las cargas mantiene la relación entre ellas
por lo que el sistema de cargas aplicado siempre puede expresarse como αP0 {1, 4}.
Conocidos los resultados del nuevo estado incremental puede definirse el estado 3 como
superposición del estado (2) y el estado incremental (∆2) multiplicado por un factor
α3 .
Estado [3] = Estado [2] + α3 · Estado [∆2]
ESTADO MA MB MC MD α [ MP /P0 l ]
Estado de referencia -0,23Pl -0,04P0 l 0,68P0 l -0,61P0 l α1 = 1, 479
Estado 1 -0,34MP -0,06MP MP -0,90MP
Estado incremental ∆1 0,07P0 l -0,69P0 l 0 -1,31P0 l α2 = 0, 077
Estado 2 -0,35MP -0,11MP MP MP
Estado incremental ∆2 -3P0 l -2P0 l 0 0 α3 = 0, 222
Estado 3 MP -0,55MP MP MP
Estado incremental ∆3 Mecanismo
Nótese igualmente que el cálculo se realiza sobre el estado inicial y los estados incre-
mentales (∆1, ∆2, ∆3) y su objetivo es evaluar la variación de la ley de momentos
flectores para el sistema de cargas de referencia. Entre los sucesivos estados incremen-
tales únicamente varı́an las coacciones internas, al añadir una rótula de momento nulo
en cada incremento, lo que hace que el proceso sea fácilmente sistematizable.
Estos teoremas son: el teorema de unicidad, el teorema estático o del lı́mite inferior y
el teorema cinemático o del lı́mite superior.
1) condición de equilibrio, tal que las fuerzas internas están en equilibrio con las cargas
externas.
4.6. ANÁLISIS LÍMITE PLÁSTICO 185
El teorema estático o del lı́mite inferior establece que cualquier distribución de esfuerzos
compatible con las ecuaciones de equilibrio, en la cual no se supera el axil/momento
plástico en ninguna sección, está originada por una carga igual o inferior a la carga
última de colapso es decir la carga que verifica el teorema estático cumple
Pestatico ≤ Pultima
El teorema cinemático o del lı́mite superior establece que cualquier mecanismo de co-
lapso cinemáticamente admisible, obtenido mediante la distribución arbitraria de sec-
ciones/rótulas plásticas en distintas secciones de la estructura, es provocado por una
carga de colapso igual o mayor a la carga última de colapso.
Es decir que, conociendo todos los posibles mecanismos de colapso de la estructura, que
cumplan las condiciones de compatibilidad, y conocidas las cargas de rotura asociadas
a dichos mecanismos, la carga última es la menor de todas ellas. Como consecuencia
de ello, puede afirmarse que cada una de las cargas obtenidas en cada uno de los
mencionados mecanismos cumple
Pcinematico ≥ Pultima
En el caso de una carga distribuida q(x), el trabajo virtual de las fuerzas externas se
obtiene integrando el trabajo realizado por un diferencial de carga en una sección dx
Z l
We = q(x)δ(x)dx (4.14)
0
Es habitual que el valor de la carga repartida sea uniforme, con lo que la expresión
anterior se simplifica calculando el trabajo virtual externo We como el valor de la carga
q por el área delimitada por la función de desplazamientos virtuales. Si se supone una
rótula en el centro, la deformada de la viga es de forma triangular con desplazamiento
virtual máximo δ max , y
Z l
l
We = q δ(x)dx = qδ max (4.15)
0 2
Una vez obtenida la mı́nima carga de colapso de los mecanismos de colapso disponibles,
únicamente puede asegurarse que dicha carga es mayor o igual a la carga de colapso
de la estructura según el enunciado del teorema cinemático. Por tanto es preciso una
verificación adicional para lo cual se emplea el teorema estático. Se trata de encon-
trar para el mecanismo de colapso y para su carga de colapso asociada, una ley de
esfuerzos flectores que cumpla las condiciones de equilibrio, junto con las condiciones
de rótula plástica implı́citas, en las cuales el momento flector es un valor conocido igual
al momento plástico de la sección.
Sea la viga biarticulada bajo carga puntual analizada en 4.4.3.1 cuyo mecanismo de
colapso se muestra en la Figura 4.41.
Puede comprobarse que la carga última obtenida por el méodo incremental y el método
cinemático coinciden para el mismo mecanismo de colapso.
Sea la viga biarticulada bajo carga distribuida analizada en 4.4.3.2 cuyo mecanismo de
colapso se muestra en la Figura 4.42.
We = ql2 δ
8MP
⇒ qu =
Wi = MP 2θ = MP 4δl l2
4.6. ANÁLISIS LÍMITE PLÁSTICO 189
Sea la viga biapoyada con voladizo con carga distribuida analizada en 4.4.3.3 mediante
el método incremental y reproducida en la Figura 4.43.
30 MP 50 MP
V1 = (↑) V2 = (↑)
7 l 7 l
We = q 21 lδ − q 12 αlδ l/2
lα
= qlδ 12 − α2
4 MP 64 MP
⇒ qu,m1 = =
Wi = MP 2 l/2δ
= 4MlP δ 2 2
1/2 − α l α= 1
7 l2
4
We = qlδα
2 2MP MP
δ ⇒ qu,m2 = 2 2 = 32 2
Wi = MP αl α l α= 1 l
4
La ley de flectores verifica que MA− = MA+ = MP , mientras que el momento flector
en el apoyo 2 es M2 = 2MP /7 < MP .
Se cumple por tanto que la ley de momentos flectores no supera el valor del momento
plástico en ninguna sección, por lo que el teorema estático permite concluir que la carga
de colapso qu,m1 es menor o igual a la carga última.
64MP 1
método cinemático qu = αcr = √ ≈ 0, 408
7l2 6
Sea la viga biempotrada sometida a dos cargas puntuales analizada en la sección 4.4.3.6
y reproducida en la Figura 4.46.
Existen cuatro rótulas posibles que combinadas dan lugar a cuatro mecanismos ele-
mentales de colapso, según muestra la Figura 4.47, cuyas cargas de colapso asociadas
se calculan mediante el principio de los trabajos virtuales, según se detalla en la Tabla
4.3.
192 CAPÍTULO 4. ANÁLISIS PLÁSTICO
Figura 4.47: Mecanismos de colapso en viga biempotrada bajo cargas puntuales asimétricas
Mecanismo 1
We = 2P δ 6MP
3δ 6δ 3δ 12MP δ ⇒ Pu,m1 =
W i = MP l + MP l + MP l = l l
Pδ Mecanismo 2
We = 2P δ + 2 = 5P2 δ 18MP MP
9MP δ ⇒ Pu,m2 = = 3, 6
Wi = MP 3δ
2l + MP 3δ
2l + 3δ
l + MP 3δl = l 5l l
2P δ Mecanismo 3
We = P δ + 2 =
2P δ 9MP MP
9MP δ ⇒ Pu,m3 = = 4, 5
Wi = MP 3δ
2l + MP 3δ
2l + 3δ
l + MP 3δl = l 2l l
Mecanismo 4
We = P δ 12MP
3δ 6δ 3δ 12MP δ ⇒ Pu,m4 =
W i = MP l + MP l + MP l = l l
Conocidas las cargas últimas de todos los mecanismos elementales de colapso posibles,
la menor de las cargas últimas es la carga última de colapso de la viga, lo que sucede
para el mecanismo 2.
4.6. ANÁLISIS LÍMITE PLÁSTICO 193
Esta conclusión debe ser corroborada por el teorema estático, para lo cual se establecen
las ecuaciones de equilibrio para el mecanismo 2 cuya carga de colapso es Pu,m2 =
18MP /5l:
X 54MP
FV = V1 + V2 − =0
5l
X 6MP
M1 = V2 l − =0
l
FV = V1 + V2 − 27M
P
2l
P
=0 12MP 15MP
P → V1 = (↑) V2 = (↑)
M1 = V2 l − 6MP = 0 2l 2l
Sea la viga continua formada por dos vanos iguales de longitud l, sometida a dos cargas
puntuales P y 2P según se indica en la Figura 4.49. Los momentos plásticos son 2MP
y MP en el primer y segundo vano respectivamente, incluyendo éste la sección en el
nudo 2.
Las cargas últimas de los tres mecanismos de colapso se determinan mediante el prin-
cipio de trabajos virtuales como se muestra en la Tabla 4.4.
que en aplicación del teorema estático permite concluir que la carga de colapso Pu,m1
es menor o igual a la carga última.
La ley de flectores para el mecanismo 1 de colapso se muestra en la Figura 4.51.
Sea un pórtico formado por dos barras de longitud l y sometido a una carga repartida
en su pilar de valor q según se indica en la Figura 4.52. Las secciones transversales de
la estructura tienen un momento plástico MP .
Las posibles rótulas se dan en el empotramiento (A), el punto medio del pilar (B) y en
el nudo de unión entre el pilar y el dintel (C). Para este conjunto de rótulas los posibles
mecanismos elementales de colapso plástico se muestran en la Figura 4.53.
Para cada mecanismo elemental se calcula la carga última aplicando el principio de los
trabajos virtuales según se muestra en la Tabla 4.5.
We = 21 qlδ ql2
4MP
δ δ → = −MA + MC ⇒ qu,m1 = 2
Wi = −MA l + MC l 2 l
We = 34 qlδ 3ql2
16MP
→ = −2MA + 2MB ⇒ qu,m2 =
Wi = −MA 2δl 2δ
+ MB l 4 3l2
We = 14 qlδ ql2
16MP
→ = −2MB + 2MC ⇒ qu,m3 =
Wi = −MB 2δl 2δ
+ MC l 4 l2
El valor del momento flector en el extremo superior del pilar es conocido de valor
M2 = MP . Derivando la expresión de la ley de momentos, se tiene que en y = l se
encuentra el extremo relativo de la parábola, por lo que el valor del momento plástico
MP no se supera en ninguna sección, cumpliendo ası́ el teorema estático. Por esta
razón puede asegurarse que la carga última que ocasiona esta distribución de flectores
es inferior o igual a la carga última de colapso de la estructura.
Sea un pórtico formado por dos pilares de diferente longitud {l, 0.75l} y un dintel de
longitud l, sometido a dos cargas puntuales P y 4P según indica la Figura 4.55. Las
secciones transversales de la estructura se caracterizan por un momento plástico MP .
Para cada mecanismo elemental se calcula la carga última a partir del principio de los
trabajos virtuales, cuyos resultados se muestran en la Tabla 4.7.
4.6. ANÁLISIS LÍMITE PLÁSTICO 201
Mecanismo 2: Dintel
We = 4P δ 2MP
→ 4P l = −2MB + 4MC − 2MD ⇒ Pu,m2 =
Wi = −MB 2δl 2δ 2δ
+ MC l 2 − MD l l
We = P δ + 4P2 δ = 3P δ
7 16MP
δ δ 4δ δ
→ 3P l = −MA +2MC − MD ⇒ Pu,m3 =
Wi = −MA l + MC l 2 − MD 3l + l 3 9l
We = 4P 2δ 5
3 − Pδ = 3Pδ 5 7 8
→ P l = MA − MB + MC ⇒ Pu,m4 =
18MP
δ 4δ δ 8δ
Wi = MA l − MB 3l + l + MC 3l 3 3 3 5l
Por otro lado, la menor de las cargas obtenidas es Pu,m3 . Por el teorema cinemático se
sabe que esta carga es mayor o igual a la carga última de colapso de la estructura.
Para verificar que efectivamente es la carga última, se recurre al teorema estático. Las
ecuaciones de equilibrio de la estructura, ası́ como las condiciones de rótula plástica
que permiten determinar la ley de momentos flectores del mecanismo de colapso, son
las siguientes:
X 16MP
FH = H1 + H4 + =0
9l
X 64MP
FV = V1 + V4 − =0
9l
X 16MP
M1 = MP + V4 l + H4 l −
3l
X 32MP
M3− = −MP + V1 l − H1 l − = −MP
9l
X 3l
M3+ = −H4 = −MP
4
El valor del momento flector en el punto superior del dintel es M2− = −5MP /9, inferior
al momento plástico. Por tanto, el valor del momento plástico no se supera en valor
absoluto en ninguna sección. Por esta razón puede asegurarse que la carga Pu,m3 que
ocasiona esta distribución de flectores es inferior o igual a la carga de colapso de la
estructura Pu .
Combinando las conclusiones obtenidas de los teoremas estático y cinemático puede
afirmarse que la carga Pu,m3 es la carga de colapso de la estructura.
16MP
Pu =
9l
Este resultado coincide con el calculado por el método incremental en 4.5, con l = 2 m.
204 CAPÍTULO 4. ANÁLISIS PLÁSTICO
5 Flexión de placas delgadas
5.1 Generalidades
Las placas son elementos estructurales planos muy utilizados en ingenierı́a civil y en
edificación. Ası́, en la Figura 5.1 se muestra el tablero de hormigón armado y pretensado
de un puente, y en las Figuras 5.2 y 5.3 edificios con forjados rectangulares y circulares,
también de hormigón armado.
Una placa plana es un elemento estructural que, como las vigas, soporta cargas transver-
sales trabajando principalmente a flexión y torsión. En una viga que soporta cargas
transversales actúan un momento flector y un esfuerzo cortante. En una placa, actúan
los momentos flectores, un momento torsor y dos esfuerzos cortantes según se muestra
en la Figura 5.4. A estos esfuerzos pueden superponerse los debidos a las cargas en el
plano de la placa, llamados esfuerzos de membrana, que no se tratarán aquı́.
205
206 CAPÍTULO 5. FLEXIÓN DE PLACAS DELGADAS
Figura 5.2: Forjados rectangulares (a) Forjado de hormigón (b) Forjado mixto
Figura 5.3: Forjados circulares (a) Marina Towers (Chicago) (b) Torres Blancas (Madrid)
parados con el espesor, y las pendientes del plano medio deformado son pequeñas.
La hipótesis 2 se considera válida para espesores menores a 1/10 de la luz menor. Para
espesores mayores, la deformación debida a los cortantes cobra importancia y se ha
de tener en cuenta en la teorı́a de placas moderadamente gruesas de Reissner-Mindlin
(1945-1951).
Las hipótesis 4,5 y 6 implican que todos los puntos situados en la normal al plano
medio de la placa tienen la misma flecha. Por tanto, la deformación de la placa está
completamente determinada por la deformada (ley de flechas) de la superficie media.
La hipótesis 4 implica que las deformaciones angulares transversales son nulas. Esta es
la hipótesis de Kirchhoff, especı́fica de las placas delgadas.
Las hipótesis 3,4 y 6 implican que las distribuciones de tensiones a través del espesor
son lineales, con valores proporcionales a la distancia al plano medio.
donde z es la coordenada normal al plano medio en el punto (x, y). La Figura 5.6
muestra una interpretación geométrica de las expresiones anteriores.
210 CAPÍTULO 5. FLEXIÓN DE PLACAS DELGADAS
Se observa como la hipótesis de normalidad (4) da lugar a distorsiones γxz y γyz nulas en
todos los puntos de la placa. En consecuencia una placa delgada presenta únicamente
tres componentes de deformaciı́on {εx , εy , γxy } que pueden expresarse en forma de vec-
tor:
∂2w
∂x 2
ε x
∂2w
χ x
ε= εy = −z = −z χ y = −zχ (5.4)
γ ∂y 2
χ
2
2∂ w
xy
xy
∂x∂y
σ = Cε (5.6)
donde C es la matriz constitutiva elástica en tensión plana. Nótese que la placa delgada
está a la vez en su estado de deformación plana y de tensión plana; esta inconsistencia
proviene de despreciar el efecto Poisson en la dirección transversal.
Los esfuerzos actuantes en la placa, cortantes (Qx ,Qy ), momentos flectores (Mx ,My ) y
momento torsor (Mxy ) adoptan el convenio de signos descrito en la Figura 5.7.
Los esfuerzos en una placa se obtienen integrando las tensiones a lo largo del espesor t
en cada punto del plano medio. Ası́:
Mx
Z t
2
Z t
2 t3
M= My =− zCεdz = z 2 Cχdz = Cχ = Dχ (5.7)
M
− 2t − 2t 12
xy
212 CAPÍTULO 5. FLEXIÓN DE PLACAS DELGADAS
1 ν 0
t3 Et3
D= C= ν 1 0 (5.8)
12 12(1 − ν 2 )
1−ν
0 0 2
Et3
D= (5.9)
12(1 − ν 2 )
La rigidez de una placa presenta una forma similar a la rigidez a flexión de una viga
rectangular de ancho unidad y canto igual al espesor de la placa. El término de Poison
(1 − ν 2 ) en el denominador resulta en una rigidez superior a la rigidez de la viga.
Sustituyendo las curvaturas de la ec. 5.4 en la definición de los momentos (ec. 5.7), los
momentos flectores y torsor son:
t
∂2w ∂2w
Z
2
Mx = − σx zdz = D + ν (5.10a)
− 2t ∂x2 ∂y 2
t 2
∂ w ∂2w
Z
2
My = − σx zdz = D ν 2 + (5.10b)
− 2t ∂x ∂y 2
t
∂2w
Z
2
Mxy = − τxy zdz = D(1 − ν) (5.10c)
− 2t ∂x∂y
Nótese que al integrar sobre el espesor [m] los esfuerzos en una placa se expresan como
momento por unidad de anchura [N m/m].
∂Qx ∂Mx
Q0x = Qx + dx Mx0 = Mx + dx (5.11)
∂x ∂x
5.3. PLACAS RECTANGULARES 213
∂Qx ∂Qy
+ +q =0 (5.12)
∂x ∂y
P
La ecuación de equilibrio de momentos Mx = 0 es
∂My ∂Mxy ∂Qy ∂Qx dy dx
dy dx + dx dy + Qy + dy dx + dx dy + qdxdy =0
∂y ∂x ∂y ∂x 2 2
Esta ecuación se usa por ejemplo en el cálculo de una placa apoyada sobre un terreno
elástico, lo que es habitual en losas de cimentación en edificación.
5.3. PLACAS RECTANGULARES 215
Sea una placa rectangular delgada cuyos contornos son paralelos a los ejes cartesianos
y se encuentran situados en x = x0 , x = xf y y = y0 , y = yf . Las condiciones de
contorno de una placa rectangular se aplican en cada uno de sus bordes en función
de los movimientos y/o esfuerzos prescritos en sus ellos. Se consideran tres tipos de
condiciones de contorno: borde empotrado, borde apoyado o borde libre.
Navier propuso en 1820 una solución para la ecuación diferencial de equilibrio de una
placa rectangular simplemente apoyada basada en la teorı́a de desarrollos en serie de
funciones de Fourier.
Sea una placa rectangular apoyada en sus cuatro bordes, de dimensiones a, b en las
direcciones x, y respectivamente según muestra la Figura 5.12. La placa está sometida
a una carga distribuida q(x, y) en cada punto.
Navier propuso una solución para la flecha w(x, y) a partir de un desarrollo en doble
serie de Fourier:
∞ X ∞
X x y
w(x, y) = wmn sin mπ sin nπ (5.26)
a b
m=1 n=1
218 CAPÍTULO 5. FLEXIÓN DE PLACAS DELGADAS
con lo que la solución para la flecha queda totalmente determinada en función de los
coeficientes de la carga qmn
∞ ∞
1 XX qmn x y
w(x, y) = sin mπ sin nπ (5.28)
Dπ 4 2 2 a b
2
m=1 n=1 m n
+
a2 b2
Una vez determinada la flecha w(x, y) el resto de magnitudes puede obtenerse derivando
según las expresiones obtenidas anteriormente. Los momentos flectores pueden expre-
sarse como:
∞ ∞ m2 2
1 XX a2
+ ν nb2 x y
Mx (x, y) = 2 qmn h i2 sin mπ sin nπ (5.29a)
π m2 2 a b
m=1 n=1
a2
+ nb2
∞ ∞ 2 2
1 XX νma2
+ nb2 x y
My (x, y) = q mn h sin mπ sin nπ (5.29b)
π2 2 2 a b
i
m2 n
m=1 n=1
a2
+ b2
∞ ∞ mn
1−ν X X ab
x y
Mxy (x, y) = q mn 2 cos mπ cos nπ (5.29c)
π2 a b
h i
m2 n2
m=1 n=1
a2
+ b2
La Figura 5.14 representa las funciones de carga q(x, y), y de flecha w(x, y) solución
del problema.
En este punto es interesante verificar la hipótesis inicial según la cual, la tensión normal
al plano medio de la placa σz es despreciable frente a las tensiones normales en el plano
de la placa σ y σy . La tensión σz es del orden del la carga repartida q.
Según la expresión 5.5 la tensión σx puede explicitarse en función de la flecha w como
2
∂2w
E ∂ w
σx = −z +ν 2 =
1 − ν 2 ∂x2 ∂y
E X mπ 2 nπ 2 x y
=z w mn + sin mπ sin nπ
1 − ν2 a b a b
se establece para placas delgadas en al menos un orden de magnitud a/t > 10, por
lo que la relación entre las tensiones en el plano de placa σx , σy y la tensión normal
al plano medio σz es de al menos dos órdenes de magnitud, confirmando la hipótesis
inicial.
El coeficiente de Fourier está definido para todos los valores de m naturales, mientras
que los términos en n pares se anulan. Nótese que esto último es debido a la simetrı́a
horizontal de la carga.
8q
(−1)m+1 n = 1, 3, 5 . . .
mnπ 2
La Figura 5.16 representa las funciones de carga q(x, y), y de flecha w(x, y) solución
del problema.
222 CAPÍTULO 5. FLEXIÓN DE PLACAS DELGADAS
b
Para el caso particular en la que la carga en banda actúe en el intervalo 2 ≤ y ≤ b, el
coeficiente de Fourier es
8q h π i
qmn = cos n − cos (nπ) m = 1, 3, 5 . . .
mnπ 2 2
La Figura 5.18 representa las funciones de carga q(x, y), y de flecha w(x, y) solución
del problema.
b0 sin nπ b
8q 8q b0
qmn = lim sin nπ = sin nπ m = 1, 3, 5 . . .
→0 bmπ b nπ b bmπ b
Para el caso particular en el que la carga esté situada en la lı́nea media b0 = b/2, se
obtiene
n−1 8q
qmn = (−1) 2 m, n = 1, 3, 5 . . .
bmπ
La Figura 5.20 representa las funciones de carga q(x, y), y de flecha w(x, y) solución
del problema.
Se considera una carga repartida sobre una superficie rectangular de dimensiones 2ξ,2
según x, y respectivamente, centrada en el punto de aplicación de la carga, cuya resul-
tante es P , es decir se define q(x, y) de la forma
P
q(x, y) = xP − ξ ≤ x ≤ xP + ξ, yP − ≤ y ≤ yP +
2ξ2
5.3. PLACAS RECTANGULARES 225
La carga puntual es el caso lı́mite del anterior para ξ, → 0 por lo que los coeficientes
de Fourier se calculan como
ξ
sin mπ sin nπ
4P xP a yP
b
qmn = lim sin mπ sin nπ
ξ→0 ab a ξ b nπ
→0 mπ
| {za } | {zb }
1 1
4P xP
y
P
= sin mπ sin nπ
ab a b
La Figura 5.22 representa las funciones de carga q(x, y), y flecha w(x, y) solución del
problema para una carga puntual en un punto (x, y) arbitrario de la placa.
226 CAPÍTULO 5. FLEXIÓN DE PLACAS DELGADAS
Existen también soluciones analı́ticas para placas rectangulares por el método del doble
desarrollo en serie de Fourier para otras condiciones de apoyo, empotradas, libres y
combinaciones de éstas. Ası́ por ejemplo, la solución de Hencky (1913) para placas
rectangulares con los cuatro bordes empotrados.
Algo más general es el método de Lèvy (1899), que propone una solución en serie de
Fourier de la forma
∞
X x
w(x, y) = Ym (y) sin mπ (5.34)
a
m=1
que satisface de forma natural las condiciones de contorno simplemente apoyado, flecha
y momento flector Mx nulos en los contornos x = 0 y x = a. De las derivadas de
las funciones Ym (y) se obtienen ecuaciones diferenciales ordinarias de cuarto orden que
satisfagan las condiciones de contorno adecuadas en los bordes y = 0 y y = b.
Ninguna de estas soluciones resulta tan fácil de obtener como el caso simplemente apo-
yado resuelto mediante la solución de Navier, para el cual las condiciones de contorno se
satisfacen a priori, pero la mayorı́a de los casos de interés práctico pueden encontrarse
en la literatura.
5.4. PLACAS CIRCULARES 227
El estudio de placas delgadas circulares se lleva a cabo con las hipótesis de Kirchhoff-
Love enunciadas en la sección 5.2, por lo que las ecuaciones de gobierno son las mismas
obtenidas allı́.
La geometrı́a de revolución del problema hace que sea recomendable utilizar coorde-
nadas cilı́ndricas (r, ϕ, z) tales que las coordenadas (r, z) describen la superficie media
de la placa y el eje z es el eje perpendicular a dicha superficie media.
El cambio de coordenadas entre las coordenadas (x, y) y las coordenadas (r, ϕ) según
la Figura 5.23 es el siguiente:
( ( p
x = rsin ϕ r = x2 +y 2
(x, y) ≡ ⇒ (r, ϕ) ≡ y (5.35)
y = rcos ϕ ϕ = arctan
x
Las derivadas del cambio de variable entre las coordenadas cartesianas y cilı́ndricas se
detalla en el apéndice A. La ecuación biarmónica de flexión de placas se expresa en
coordenadas cilı́ndricas como
q(r, ϕ)
∆2 w(r, ϕ) = (5.36)
D
siendo el operador Laplaciano en coordenadas cilı́ndricas
∂2 1 ∂ 1 ∂2
∆= + + (5.37)
∂r2 r ∂r r2 ∂ϕ2
228 CAPÍTULO 5. FLEXIÓN DE PLACAS DELGADAS
∂2w 1 ∂ 2 w 1 ∂w
Mr = D +ν + (5.38a)
∂r2 r2 ∂ϕ2 r ∂r
2 2
∂ w 1 ∂ w 1 ∂w
Mϕ = D ν 2 + + (5.38b)
∂r r2 ∂ϕ2 r ∂r
2
1 ∂ w 1 ∂w
Mrϕ = (1 − ν)D − (5.38c)
r ∂r∂ϕ r2 ∂ϕ
∂
Qr = D ∆w(r, ϕ) (5.39a)
∂r
∂
Qϕ = D ∆w(r, ϕ) (5.39b)
∂ϕ
5.4. PLACAS CIRCULARES 229
1 ∂2w
∂ 1−ν ∂ 1 ∂w
Vr = D ∆w + − (5.40a)
∂r r ∂ϕ r ∂r∂w r2 ∂ϕ
∂ 1 ∂2w
1 ∂ 1 ∂w
Vϕ = D ∆w + (1 − ν) − (5.40b)
r ∂ϕ ∂r r ∂r∂w r2 ∂ϕ
Si tanto la geometrı́a como la carga son de revolución la solución del problema mantiene
dicha simetrı́a por lo que las expresiones anteriores se simplifican notablemente ya que
la función q = q(r, ϕ) y la función de flecha w = w(r, ϕ) son únicamente función de r.
En estas condiciones, las derivadas respecto a ϕ se anulan y la ecuación de equilibrio
es
2 1 d d 1 d dw(r) q(r)
∆ w(r) = r r = (5.41)
r dr dr r dr dr D
o equivalentemente
d 1 d dw
Qr = D r (5.42)
dr r dr dr
F = πR2 q (5.46)
Sea una placa circular de radio R, espesor t y rigidez D sometida a una carga distribuida
q uniforme en toda su superficie.
1
2πrQr = qπr2 ⇒ Qr = qr 0≤r≤R
2
qr3 r C2
φ(r) = + C1 +
16D 2 r
qr4 r2
w(r) = + C1 + C2 ln r + C3
64D 4
Nótese que la ley de cortantes no depende de las condiciones de contorno sino de las
condiciones de carga por lo que el resultado anterior permite particularizar la solución
para distintas condiciones de contorno de la placa. Se plantea la solución del problema
para dos casos: borde empotrado y borde apoyado.
Sea una placa circular con el borde empotrado según se muestra en la Figura 5.25.
5.4. PLACAS CIRCULARES 231
La condición de empotramieto implica que las flechas y los giros son nulos en el contorno,
es decir,
w(r)|r=R = 0 φ(r)|r=R = 0
−qR2 qR4
C1 = C2 = 0 C3 =
8D 64D
qR2
1 − ρ2
χ(r) = −
16D
qR3
ρ 1 − ρ2
φ(r) = −
16D
qR4 2
w(r) = 1 − ρ2
64D
qR4
ω(r)|r=0 =
64D
qR2
M (r)| = − (1 + ν)
2
ϕ r=0
qR 2
16
Mϕ (r) = − (1 + ν) − (1 + 3ν)ρ → 2
16
Mϕ (r)|r=R = 2qR ν
16
Las flechas y los giros junto con los momentos radiales y circunferenciales se representan
en la Figura 5.26.
Sea una placa circular con el borde apoyado según se muestra en la Figura 5.27
La condición de apoyo implica que las flechas y los momentos radiales son nulos en el
contorno, es decir,
w(r)|r=R = 0 Mr (r)|r=R = 0
qR2 3 + ν
χ(r) = − − 3ρ2
16D 1 + ν
qR3
3+ν 2
φ(r) = − ρ −ρ
16D 1 + ν
qR4 4
3+ν 2 5+ν
w(r) = ρ −2 ρ +
64D 1+ν 1+ν
qR3
φ(r)|r=R = −
8D(1 + ν)
qR4 5 + ν
ω(r)|r=0 =
64D 1 + ν
qR2 qR2
(3 + ν) 1 − ρ2 → Mr (r)|r=0 = −
Mr (r) = − (3 + ν)
16 16
qR2
M (r)| = − (3 + ν)
2
ϕ r=0
qR 1 + 3ν 2
16
Mϕ (r) = − (3 + ν) 1 − ρ → 2
16 3+ν
Mϕ (r)|r=R = − 2qR (1 − ν)
16
Las flechas y los giros junto con los momentos radiales y circunferenciales se representan
en la Figura 5.28.
234 CAPÍTULO 5. FLEXIÓN DE PLACAS DELGADAS
Sea una placa circular de radio R, espesor t y rigidez D sometida a una carga disribuida
q(r) lineal ascendente con la coordenada r de valor máximo q en su contorno.
r
q(r) = q 0≤r≤R
R
r r 2 qr2
2πrQr = q 2π r ⇒ Qr = 0≤r≤R
2 R 3 3R
qr4 r C2
φ(r) = + C1 +
45DR 2 r
qr5 r2
w(r) = + C1 + C2 ln r + C3
225DR 4
Sea una placa circular con el borde empotrado según se muestra en la Figura 5.29.
La condición de empotramieto implica que las flechas y los giros son nulos en el contorno,
es decir,
w(r)|r=R = 0 φ(r)|r=R = 0
2qR2 qR4
C1 = − C2 = 0 C3 =
45D 150D
qR2
1 − 4ρ3
χ(r) = −
45D
qR3
ρ 1 − ρ3
φ(r) = −
45D
qR4 5
2ρ − 5ρ2 + 3
w(r) =
450D
qR4
ω(r)|r=0 =
150D
qR2
M (r)| = − (1 + ν)
2
ϕ r=0
qR 1 + 4ν 3
45
Mϕ (r) = − (1 + ν) 1 − ρ → 2
45 1+ν
Mϕ (r)|r=R = qR ν
15
Las flechas y los giros junto con los momentos radiales y circunferenciales se representan
en la Figura 5.30.
Sea una placa circular con el borde apoyado según se muestra en la Figura 5.31
La condición de apoyo implica que las flechas y los momentos radiales son nulos en el
contorno, es decir,
w(r)|r=R = 0 Mr (r)|r=R = 0
Las flechas y los giros junto con los momentos radiales y circunferenciales se representan
en la Figura 5.32.
Sea una placa circular de radio R, espesor t y rigidez D sometida a una carga puntual
P aplicada en el centro de la placa.
Integrando la ley de cortantes, se obtienen las expresiones para el giro φ(r) y la flecha
w(r):
Pr 1 r
φ(r) = ln r − + C1 + C2
4πD 2 2
P r2 r2
w(r) = [ln r − 1] + C1 + C2 ln r + C3
8πD 4
Sea una placa circular con el borde empotrado según se muestra en la Figura 5.33
La condición de empotramieto implica que las flechas y los giros son nulos en el contorno,
5.4. PLACAS CIRCULARES 239
es decir,
w(r)|r=R = 0 φ(r)|r=R = 0
P R2
P 1
C1 = − ln R − C2 = 0 C3 =
2πD 2 16πD
P
χ(r) = (ln ρ + 1)
4πD
PR
φ(r) = ρln ρ
4πD
P R2
2 1
w(r) = 2ρ ln ρ − +1
16πD 2
P R2
ω(r)|r=0 =
16πD
Las flechas y los giros junto con los momentos radiales y circunferenciales se representan
en la Figura 5.34.
Sea una placa circular con el borde apoyado según se muestra en la Figura 5.35
La condición de apoyo implica que las flechas y los momentos radiales son nulos en el
contorno, es decir,
w(r)|r=R = 0 Mr (r)|r=R = 0
240 CAPÍTULO 5. FLEXIÓN DE PLACAS DELGADAS
P R2 3 + ν
P 1−ν
C1 = − 2lnR + C2 = 0 C3 =
4πD 1+ν 16πD 1 + ν
PR
φ(r)|r=R = −
4πD(1 + ν)
P R2 3 + ν
ω(r)|r=0 =
16πD 1 + ν
P
Mr (r) = (1 + ν)ln ρ → Mr (r)|r→0 → −∞
4π
P Mϕ (r)|r→0 → −∞
Mϕ (r) = [(1 + ν)ln ρ − (1 − ν)] → P
4π Mϕ (r)|r=R = − (1 − ν)
4π
Las flechas y los giros junto con los momentos radiales y circunferenciales se representan
en la Figura 5.36.
Sea una placa circular de radio R, espesor t y rigidez D sometida a un momento radial
distribuido en el borde exterior de valor constante M según muestra la Figura 5.37.
242 CAPÍTULO 5. FLEXIÓN DE PLACAS DELGADAS
2πrQr = 0 ⇒ Qr = 0 0≤r≤R
Integrando la ley de cortantes, se obtienen las expresiones para el giro φ(r) y la flecha
w(r)
C1 C2
χ(r) = − 2
2 r
r C2
φ(r) = C1 +
2 r
r 2
w(r) = C1 + C2 ln r + C3
4
w(r)|r=R = 0 Mr (r)|r=R = M
Aplicando las condiciones de contorno para el borde apoyado, y teniendo en cuenta que
la flecha en el centro es finita, las constantes de integración son
2M M R2
C1 = C2 = 0 C3 = −
D(1 + ν) 2D(1 + ν)
5.4. PLACAS CIRCULARES 243
M R2
ω(r)|r=0 = −
2D(1 + ν)
Mr (r) = Mϕ (r) = M
Las flechas y los giros junto con los momentos radiales y circunferenciales se representan
en la Figura 5.38.
6.1 Introducción
A menudo, las fuerzas que actúan sobre las estructuras son variables y dependientes del
tiempo. De hecho, prácticamente todas las estructuras se ven sometidas a una u otra
forma de carga dinámica a lo largo de su vida útil, ya sean éstas debidas a fenómenos
naturales o a la actividad humana. Ası́, en la Figura 6.1 se muestran acciones de natu-
raleza periódica regular debidas a equipos mecánicos en funcionamiento o a la acción
del oleaje. Las funciones regulares pueden ser armónicas simples o más complejas. En
la Figura 6.2 se muestran acciones de naturaleza no periódica de larga duración, tales
como las acciones sı́smicas o el viento, mientras que en la Figura 6.3 se indican cargas
impulsivas tales como las debidas a impactos y explosiones.
245
246 CAPÍTULO 6. DINÁMICA DE ESTRUCTURAS
Por otro lado, como consecuencia de la consideración de estas derivadas temporales, las
ecuaciones de equilibrio de fuerzas se convierten en ecuaciones diferenciales respecto al
tiempo, cuya solución requiere una interpretación en términos fı́sicos y estructurales.
Todo esto hace que el cálculo dinámico de estructuras precise de técnicas de resolución
distintas del cálculo estático y que la interpretación de los resultados sea también
diferente.
o bien
mẍ(t) + kx(t) = 0 (6.2)
Escrito de esta manera, el término mẍ(t) se puede interpretar com una fuerza de inercia
que se opone al movimiento y tiene, por tanto, el mismo sentido positivo que la fuerza
elástica (principio de D’Alembert). Ası́ el sistema dinámico se puede estudiar como un
sistema estático con fuerzas de inercia.
248 CAPÍTULO 6. DINÁMICA DE ESTRUCTURAS
Dividiendo por la masa, la ecuación (6.2) puede reescribirse obteniendo la ecuación del
movimiento dinámico:
ẍ(t) + ω 2 x(t) = 0 (6.3)
p
siendo ω := k/m la denominada frecuencia natural del sistema expresada en rad/s.
La solución de la ecuación diferencial de segundo orden es de la forma x(t) = Aert ,
donde r es un número complejo. Derivando, sustituyendo en la ecuación del movimiento
y operando, la ecuación de equilibrio queda de la siguiente forma:
Aert r2 + ω 2 = 0 → r = ±iω
(6.4)
Las relación entre las constantes de integración A, ϕ y las condiciones iniciales del
movimiento son:
s 2
2 ẋ0 ωx0
A = x0 + tan ϕ = (6.7)
ω ẋ0
En la Figura 6.5 se representan las tres funciones halladas para el desplazamiento, ve-
locidad y aceleración en un sistema dinámico que vibra libremente sin amortiguamiento.
Aert r2 + 2βr + ω 2 = 0
(6.11)
6.2.2.1 Sobreamortiguamiento
El coeficiente de amortiguamiento es ν > 1 y por tanto las dos raı́ces son reales
p
r1,2 = −ω ν ± ν 2 − 1 (6.14)
Considerando las condiciones iniciales x(0) = x0 y ẋ(0) = ẋ0 , se obtienen valores para
las constantes A,B
A = x0 B = ẋ0 + ωx0 (6.18)
6.2.2.3 Inframortiguamiento
El coeficiente de amortiguamiento es ν < 1 y, por tanto, las dos raı́ces obtenidas son
complejas p
r1,2 = −ων ± iω 1 − ν 2 = −ων ± iων (6.20)
6.2. SISTEMAS CON UN GRADO DE LIBERTAD 253
√
siendo ων := ω 1 − ν 2 la denominada frecuencia amortiguada del sistema.
A1 e−νωt √2πν
= −νω(t+T ) = eνωTν = e 1−ν 2 (6.23)
A2 e ν
Sea un peso W que oscila sobre un soporte fabricado en acero de módulo E, con una
longitud l y una sección transversal cuadrada de lado a. El desplazamiento inicial del
sistema es x0 y la velocidad inicial ẋ0 . El coeficiente de amortiguamiento del sistema
es ν.
a4
I= = 1, 33 · 10−8 m4
12
Se deja caer una masa m desde una altura H sobre una estructura asimilable a una
viga biarticulada de hormigón de longitud l en su punto medio. La sección transversal
de la viga es rectangular de ancho b y canto h. El sistema tiene un amortiguamiento
de valor c.
Al caer libremente la masa impacta con la viga a una velocidad v0 que supone la
velocidad inicial del movimiento vibratorio
p
v0 = 2gh = 4, 43 m/s
P l3 P 48EI
δ= = ⇒ k = 3 = 2500 kN/m
48EI k l
La amplitud A calculada para el caso amortiguado resulta ser A = 88, 6 · 10−3 m. Las
funciones de desplazamiento y velocidad quedan completamente determinadas:
Sea un pórtico formado por dos pilares iguales de altura h y separados una distancia
b y un forjado rı́gido de masa m, que se desplaza horizontalmente. Los pilares son de
hormigón y de sección cuadrada de lado a. Las condiciones iniciales son x0 y ẋ0 . El
amortiguamiento del sistema es ν.
En primer lugar se estudian las caracterı́sticas dinámicas del pórtico para oscilar la-
teralmente. Para calcular la rigidez del sistema es preciso conocer la inercia de los
pilares:
a4
I= = 6, 75 · 10−4 m4
12
La rigidez del sistema es la relación entre el desplazamiento estático producido por una
carga en el esquema estructural abordado. Para evaluarla se calcula la rigidez de cada
pilar suponiendo que está empotrado en ambos extremos y suponiendo que no existe
desplazamiento horizontal relativo entre los nudos del pórtico.
F h3 F 12EI
u= = ⇒ kpilar =
12EI kpilar h3
La solución particular del problema tiene en cuenta la función forzada y es del tipo
El oscilador actúa por tanto, como un filtro de frecuencias, porque responde dinámicamente
a frecuencias forzadas similares a la frecuencia natural del sistema, anula la respuesta
6.2. SISTEMAS CON UN GRADO DE LIBERTAD 263
Volviendo a la ecuación (6.27), la respuesta dinámica del sistema a una excitación ex-
terna es la combinación de las dos soluciones, homogénea y particular. La solución
homogénea o transitoria depende de las condiciones iniciales y su influencia es decre-
ciente de forma exponencial, de modo que a tiempo infinito la solución es únicamente
la solución particular o permanente.
En la Figura 6.16 se muestran las soluciones homogénea, particular y total para una
vibración forzada, en la cual se puede apreciar cómo la solución homogénea (transitoria)
se amortigua disminuyendo su amplitud tendiendo a anularse en t → ∞, mientras que
la solución particular conserva su amplitud. En consecuencia, a medida que transcurre
el tiempo la solución coincide con la solución particular.
Un depósito de sección transversal S contiene agua hasta una altura h. El agua se eleva
hasta el depósito mediante una bomba cuya masa es mb . Dicho depósito está situado
sobre un soporte de acero, de longitud l y de sección transversal circular tubular de
diámetro exterior D y espesor e.
La bomba ejerce una fuerza F0 girando a una frecuencia θ. El sistema tiene un coefi-
ciente de amortiguamiento ν.
Determinar la altura de agua para la cual se igualan las frecuencias natural y forzada.
264 CAPÍTULO 6. DINÁMICA DE ESTRUCTURAS
Explicar las condiciones necesarias para que el sistema entre en resonancia y calcular
la amplitud del movimiento dinámico.
Datos: mb = 500 kg, S = 1 m2 , l = 8 m, D = 25 cm, e = 1 cm, Eacero = 200 GP a,
F0 = 100 N , θ = 5 rad/s, ν = 2, 5%
En primer lugar se calculan las caracterı́sticas dinámicas del sistema en el que la masa de
agua se supone concentrada y el soporte aporta la rigidez que determina el movimiento
vibratorio.
La inercia de la sección transversal del pilar es
π 4
R − (R − e)4 = 5, 44 · 10−5 m4
I=
4
A partir de la inercia a flexión del pilar, la rigidez del sistema estructural se calcula
suponiendo que éste se comporta como una ménsula
3EI
k= = 63750 N/m
l3
La masa del sistema depende de la altura de agua e incluye la masa de la bomba y por
tanto puede calcularse como m = 500 + 1000h. La frecuencia natural de vibración del
sistema se calcula como r
63750
ω=
500 + 1000h
6.2. SISTEMAS CON UN GRADO DE LIBERTAD 265
F0 1 F0 1
Adin = q = · = 0, 031 m = 31 mm
k 2 2 2 2 k 2ν
(1 − γ ) + 4ν γ
bh3
I= = 2, 083 · 10−6 m4
12
266 CAPÍTULO 6. DINÁMICA DE ESTRUCTURAS
La rigidez del sistema se calcula para un esquema estructural de viga biarticulada con
carga centrada:
P l3 P 48EI
f= = ⇒ k = 3 = 7, 78 · 105 N/m
48EI k l
El motor, incluyendo la masa excéntrica, tiene una masa M = 8154, 9 kg, que es la
masa oscilante. Por su parte, la excéntrica tiene una masa me = 51 kg y su movimiento
es diferente al del motor excluyendo la excéntrica. Por tanto, la aportación inercial de
la masa que vibra pero no rota, M − me , y la masa que vibra y además rota, me , se
han de considerar por separado.
La función de movimiento de la masa excéntrica ye (t) está relacionada con la vibración
del fleje y(t). La función de movimiento de la masa excéntrica y sus derivadas se
expresan como:
La ecuación del movimiento oscilatorio del sistema se obtiene sumando las fuerzas de
inercia de las masas M −me y me , las fuerzas viscosas y las fuerzas elásticas que actúan
sobre la masa oscilante M , según se muestra en la Figura 6.19
A partir del esquema anterior, se establece la ecuación de equilibrio dinámico del sis-
tema:
Se consideran las funciones elementales impulso, escalón y rampa como funciones tem-
porales para la excitación forzada de un sistema dinámico con frecuencia natural ω,
representadas en la Figura 6.20 y definidas de la siguiente forma:
Figura 6.20: Funciones elementales de excitación. (a) Impulso (b) Escalón (c) Rampa
siendo δ(t) la función delta-Dirac y H(t) la función de Heaviside standard cuya amplitud
es unitaria. Ası́ definidas, la tres funciones se relacionan entre sı́:
Según se vio en 6.2.3, el equilibrio de fuerzas para una sistema dinámico de masa m
y rigidez k, con k = mω 2 , sobre el que actúa una fuerza exterior F (t), da lugar a la
ecuación del movimiento
F (t)
ẍ(t) + 2νω ẋ(t) + ω 2 x(t) = (6.35)
m
esto es, xesc,p = xest la solución estática para una fuerza unidad.
Dadas las condiciones iniciales, x0 = ẋ0 = 0, las constantes de integración son:
√
1 1 − ν2
A= √ tan ϕ = (6.39)
k 1 − ν2 ν
Nótese que la solución dinámica obtenida se expresa como la solución estática bajo
una fuerza constante, xest , multiplicada por un factor de amplificación dependiente del
tiempo y de los parámetros dinámicos.
Siendo la función impulso la derivada de la función de Heaviside, la respuesta dinámica
ante una excitación forzada dada por un impulso es la derivada de la solución xesc (t):
ω
ximp (t) = ẋesc (t) = xest √ e−ωνt sin (ων t) (6.42)
1 − ν2
e−ωνt 2ν 2 − 1
2ν
xram (t) = xest t − + 2νcos (ων t) + √ sin (ων t) (6.44)
ω ω 1 − ν2
es decir, a tiempo infinito, la respuesta dinámica replica la función rampa que excita el
sistema con un retraso en el tiempo tr = 2ν/ω.
En la Figura 6.21 se muestran las respuestas obtenidas para las tres funciones de ex-
citación simple obtenidas, normalizadas a la respuesta estática.
Se muestran las respuestas para dos coeficientes de amortiguamiento distintos, ν = 1% y
ν = 5%, durante 10 periodos de oscilación, donde se puede observar tanto la naturaleza
oscilatoria de la respuesta como su tendencia a tiempo infinito.
En los desarrollos anteriores se han empleado las siguientes relaciones trigonométricas
√
consecuentes con el valor de tan ϕ del que se obtiene que sin ϕ = 1 − ν 2 y cos ϕ = ν:
p
sin (ων t + ϕ) = ν sin (ων t) + 1 − ν 2 cos (ων t)
p (6.46)
cos (ων t + ϕ) = ν cos (ων t) − 1 − ν 2 sin (ων t)
Sea una función de excitación F (t) general que se discretiza en una serie de impulsos
sucesivos de duración ∆τ , tal como se muestra en la Figura 6.22.
que se expresa en función de la solución estática xest = F0 /k, donde k es la rigidez del
sistema, con k = mω 2
siendo las constantes de integración para la solución ası́ definida A, ϕ, las correspondien-
tes a la vibración libre obtenidas en 6.2.1, expresiones (6.7), aplicadas a las condiciones
274 CAPÍTULO 6. DINÁMICA DE ESTRUCTURAS
Se considera ahora la función de impulso triangular decreciente definida como Ftri (t) =
F0 (1 − t/tf ) , 0 ≤ t ≤ tf , actuando sobre un oscilador simple no amortiguado, según
se muestra en la Figura 6.26.
Los espectros de respuesta permiten evaluar directamente los máximos efectos espe-
rados para una carga dinámica que actúa sobre un sistema dinámico de un grado de
libertad.
El espectro de respuesta dinámica está formado por el valor máximo de cada una de
las curvas anteriores, dependiente de la relación tr /T , tal como se muestra en la Figura
6.31.
Nótese que la respuesta dinámica, y el correspondiente espectro, correspondiente a una
función de excitación formada por tramos lineales puede obtenerse por superposición
de las respuestas de sucesivas funciones rampa-escalón.
6.2. SISTEMAS CON UN GRADO DE LIBERTAD 279
La solución para una función de aceleración a(t) general se obtiene mediante la integral
de Duhamel: Z t
1
x(t) = − a(τ )e−νω(t−τ ) sin(ων (t − τ ))dτ (6.64)
ων 0
Durante muchos años el mayor terremoto del cual se disponı́an registros fue el terremoto
de El Centro (California, 1940), por lo que este evento se emplea a menudo como
referencia para la respuesta sı́smica de un sistema estructural. La Figura 6.32 muestra
la componente Norte-Sur de la aceleración registrada en el terremoto del El Centro.
es mayor en amplitud para las frecuencias más cercanas a la frecuencia natural del
sistema. De este modo, dicha respuesta x(t) y en particular su valor máximo max|x(t)|
están fuertemente condicionadas por la frecuencia natural del sistema ω.
Dado un acelerograma asociado a un evento sı́smico y un sistema dinámico caracte-
rizado por su frecuencia natural ω y su amortiguamiento ν, se obtiene una respuesta
en desplazamientos x(t), de la cual se considera el desplazamiento máximo. La curva
de desplazamientos máximos obtenida en la respuesta de un sistema dinámico Sd (ω, ν)
dado un valor de ν se denomina desplazamiento espectral :
Figura 6.33: Espectro de respuesta en aceleración para el sismo de El Centro (California, 1940)
con lo que basta conocer uno de los tres para obtener los restantes.
A partir de los años 70 del pasado siglo el volumen de información registrada acerca
de los eventos sı́smicos en todo el mundo empieza a ser los suficientemente importante
como para plasmar la influencia del tipo de suelo en la respuesta de la estructura. La
Figura 6.35 muestra cuatro espectros de aceleración normalizados correspondientes a
un cierto número de terremotos registrados en distintos tipos de suelo, para un total de
104 registros. Los tipos de suelo considerados son: rocas y suelos duros (tipo I), suelos
profundos sin cohesión (tipo II) y arcillas/arenas blandas o medias (tipo III). La curva
A-B-C-D define un espectro de respuesta conservador para suelos duros de tipo I.
6.2. SISTEMAS CON UN GRADO DE LIBERTAD 283
Figura 6.35: Espectro sintético de respuesta en función de las caracterı́sticas del suelo
Sea una estructura cuya respuesta dinámica se modela mediante un único grado de
libertad con periodo natural de oscilación T = 1 s. Determinar el máximo absoluto
de la aceleración, el máximo desplazamiento relativo y la máxima pseudo-velocidad
relativa en dos casos:
(b) un espectro de diseño sı́smico con una aceleración máxima del suelo igual a 0, 32g.
Suponer un amortiguamiento de la estructura del 10% del amortiguamiento crı́tico.
En el caso del espectro sintético de respuesta, los valores espectrales se obtienen di-
rectamente de la gráfica para las tres magnitudes de respuesta. Estos valores han de
escalarse al factor dado por la máxima aceleración del suelo con respecto a la gravedad,
es decir, f = 0, 32. La velocidad de diseño se obtiene a partir del valor espectral de
la Figura 6.36b como Sv = 0, 32 · 1, 530 = 0, 490 m/s. El desplazamiento máximo
se obtiene como Sd = 0, 32 · 24, 9 = 7, 98 cm mientras que la aceleración máxima es
Sa = 0, 32 · 0, 899 = 0, 288g.
En la Tabla 6.2 se resumen los valores de diseño obtenidos por ambos métodos.
donde F(t) es un vector de componentes Fi (t) que dependen del tiempo y que adoptan
la forma
donde los tres primeros sumandos son las fuerzas internas (de inercia, amortiguamiento
y resistente, respectivamente) y el cuarto sumando son a las fuerzas externas.
Sea δx una variación arbitraria del campo de desplazamientos, compatible con las
condiciones de apoyo. Sea un intervalo de tiempo [t1 , t2 ], tal que δx(t = t1 ) = δx(t =
t2 ) = 0.
Multiplicando la ecuación (6.79) por δxT , integrando en el tiempo en el intervalo [t1 , t2 ]
y operando, se tiene que
Z t2 Z t2
1 1 T T
δ xT Cẋ dt = 0
δ ẋMẋ − δ x Kx − x F dt + (6.80)
t1 2 2 t1
δV = δU + δH = 0 (6.83)
o sea
d 1 T 1 T
ẋ Mẋ + x Kx − x F + ẋT Cx = 0
T
(6.85)
dt 2 2
K̇ + Ḣ = 0 (6.87)
t/T0 Ref x ẋ U K H E
k A,A’ 0 0 0 0 0 0
k 1 1 2 1 2
2 + 4 B,B’ xest ωxest 2 kxest 2 kxest −kx2est 0
k + 12 C 2xest 0 2kx2est 0 −2kx2est 0
La solución del problema dinámico implica resolver una ecuación diferencial de se-
gundo orden con condiciones iniciales en desplazamientos y velocidades. A menudo, la
obtención de la solución analı́tica del problema dinámico es imposible en la práctica
y resulta necesaria la solución numérica del mismo. Para ello es preciso discretizar el
dominio temporal en intervalos finitos de tiempo [tn , tn+1 ] de amplitud ∆t aproximando
la solución en el instante tn+1 a partir de la solución conocida en el instante tn .
La integración numérica puede plantearse para un único grado de libertad o bien para
n grados de libertad. Se presenta a continuación el método de Newmark para la inte-
gración numérica de la ecuación de equilibrio dinámico en ambos casos.
En el caso de una aceleración conocida de la base del sistema a(t), la fuerza actuante
es F (t) = −ma(t) y el movimiento x(t) es relativo a la base del sistema dinámico.
donde xPn+1 y ẋPn+1 son predicciones de los valores de xn+1 y ẋn+1 calculados a partir
de los valores conocidos en tn .
Sustituyendo las aproximaciones (6.95b) en la ecuación (6.94) escrita para tiempo tn+1 ,
se tiene
1 γ P P
m + c + k ∆x = F n+1 − cẋ n+1 + kx n+1 (6.96)
β∆t2 β∆t
o bien
k ∗ ∆x = Fn+1
∗
(6.97)
siendo
1 γ
k∗ = m+ c+k (6.98)
β∆t2 β∆t
∗
=Fn+1 − cẋPn+1 + kxPn+1
Fn+1 (6.99)
6.5. SOLUCIÓN NUMÉRICA DE LA ECUACIÓN DE EQUILIBRIO DINÁMICO.
MÉTODO DE NEWMARK 293
∗ /k ∗ y ẋ
De la ecuación (6.97) se obtiene ∆x = Fn+1 n+1 y ẍn+1 sustituyendo ∆x en la
ecuación (6.95b).
siendo xtri,0 = F0 /kωtp una amplitud de referencia. La función de velocidad del sistema
294 CAPÍTULO 6. DINÁMICA DE ESTRUCTURAS
dinámico es
1 − cos ωt 0 ≤ t ≤ tp
ẋtri (t)
= −1 + 2cos ω(t − tp ) − cos ωt tp < t ≤ 2tp
ẋtri,0
−cos ω(t − 2t ) + 2cos ω(t − t ) − cos ωt 2t < t
p p p
siendo ẋtri,0 = F0 /ktp un factor de amplitud para la velocidad. Por su parte, la función
de aceleración de sistema se expresa como
sin ωt 0 ≤ t ≤ tp
ẍtri (t)
= −2sin ω(t − tp ) + sin ωt tp < t ≤ 2tp
ẍtri,0
sin ω(t − 2t ) − 2sin ω(t − t ) + sin ωt 2t < t
p p p
La solución puede ser aproximada mediante integración numérica tal como se explica en
la sección 6.5. Se utilizan los parámetros de integración γ = 1/2 y β = 1/4. Nótese que
para los parámetros escogidos el método de Newmark es incondicionalmente estable.
La integración numérica se realiza con diversos incrementos de tiempo ∆t en relación
al periodo natural T del sistema, ∆t/T = {1/5, 1/10, 1/20} para evaluar su grado de
aproximación.
La Figura 6.40 muestra la solución exacta (rojo) y las tres aproximaciones numéricas
para los tres pasos de tiempo. Se observa en ella que las soluciones aproximadas mues-
tran un cierto error tanto en amplitud como en fase, respecto de la solucióon analı́tica.El
error se reduce de forma cuadrática a medida que se reduce el incremento de tiempo
utilizado en el cálculo.
(·)tp (·)2tp
Desplazamientos
Exacta 3,142 0,000
Newmark ∆t/T = 1/5 2,684 1,183
Newmark ∆t/T = 1/10 3,044 0,389
Newmark ∆t/T = 1/20 3,116 0,102
Velocidades
Exacta 2,000 -4,000
Newmark ∆t/T = 1/5 1,516 -3,380
Newmark ∆t/T = 1/10 1,995 -3,971
Newmark ∆t/T = 1/20 2,000 -3,998
Aceleraciones
Exacta 0,000 0,000
Newmark ∆t/T = 1/5 -0,171 -1,183
Newmark ∆t/T = 1/10 0,097 -0,389
Newmark ∆t/T = 1/20 0,025 -0,102
Esta ecuación se puede resolver por el método de Newmark escribiendo las ecuaciones
(6.94) y (6.97) en forma vectorial. Esto es
con
1
xPn+1 = xn + ∆tẋn + − β ∆t2 ẍn (6.102a)
2
ẋPn+1 = ẋn + (1 − β) ∆tẍn (6.102b)
K∗ ∆x = F∗n+1 (6.103)
siendo
1 γ
K∗ = 2
M+ C+K (6.104)
β∆t β∆t
F∗n+1 = Fn − CẋPn+1 − KxPn+1 (6.105)
Las ecuaciones matriciales (6.73), (6.74) y (6.76) son ecuaciones diferenciales de segundo
orden en el tiempo que se resuelven conocidas las correspondientes condiciones iniciales
en desplazamientos x(t)|t=0 = x0 y ẋ(t)|t=0 = ẋ0 .
Dado que la matriz de rigidez K y en su caso la matriz de amortiguamiento C no
son diagonales, la respuesta de los distintos grados de libertad está acoplada. Para la
resolución de la ecuación diferencial de equilibrio existe una alternativa basada en la
obtención de las frecuencias propias y los modos propios de vibración del sistema. Estos
298 CAPÍTULO 6. DINÁMICA DE ESTRUCTURAS
ser nulo:
K − ω 2 M = 0 (6.108)
Un problema dinámico con n grados de libertad tiene, por tanto n frecuencias propias
de vibración ωi con sus correspondientes periodos de vibración propios Ti = 2π/ωi .
Cada una de las frecuencias propias tiene asociado un vector propio ai que define el
modo de vibración asociado a dicha frecuencia ωi .
Se define un vector propio normalizado ϕi como el vector propio dividido por su primera
componente, es decir
1
ϕi = ai (6.110)
ai1
Los vectores propios presentan algunas propiedades destacables con respecto a las ma-
trices de rigidez K y de masa M:
• Los vectores propios son ortogonales entre sı́, es decir ϕ̂Ti · ϕ̂j = 0, i 6= j.
• Los vectores propios son ortogonales con respecto a la matriz de masa, es decir
ϕ̂Ti Mϕ̂j = 0, i 6= j.
Gracias a las propiedades de los autovectores o modos propios de vibración del sistema
ϕ̂i , i = 1, . . . , n, linealmente independientes, la solución x = x(t) puede escribirse como
x = Ωz (6.115)
donde
1 α
ξi = + βωi (6.117)
2 ωi
Calcular la matriz modal Ω implica conocer los n modos del sistema. A menudo
esto es innecesario. En la mayorı́a de los problemas, los modos correspondientes a
las frecuencias más altas participan poco en la respuesta estructural y su contribución
puede despreciarse. En la práctica, sólo se necesita conocer la respuesta modal de
un reducido número de modos de baja frecuencia. Por tanto, únicamente se requiere
plantear y resolver las m primeras de las n ecuaciones modales (6.116).
6.7.1 Introducción
• Se considera que los forjados son mucho más rı́gidos que los pilares, por lo que se
suponen nulos los giros en los nudos que los unen.
Las anteriores hipótesis definen un sistema de masas concentradas y rigideces para los
grados de libertad del edificio, que de forma esquemática pueden observarse en la Figura
6.41.
Según la idealización estructural adoptada (Figura 6.41), cada piso i del edificio está
asociado a un ńico grado de libertad de movimiento horizontal.
modo m, S(ωm ):
xmax
i,m =βm φi,m Sd (ωm ) (6.120a)
ẋmax
i,m =βm φi,m Sv (ωm ) (6.120b)
ẍmax
i,m =βm φi,m Sa (ωm ) (6.120c)
La fuerza sı́smica Fi,m que actúa sobre el piso i debido al modo m se calcula mediante
la segunda ley de Newton como el producto de la masa en cada grado de libertad mi y
la aceleración máxima asociada al modo m:
n
X
Fi,m = mi ẍi,m (6.122)
m=1
6.7. ANÁLISIS SÍSMICO DE EDIFICIOS 305
El cálculo de la componente modal del momento flector Mi,m en cada piso i del edificio
es la resultante de los momentos producidos por las fuerzas sı́smicas Fi,m en los niveles
superiores. Si hi es la altura de cada piso respecto del suelo, la componente del momento
Mi,m se calcula como
Xn
Mi,m = Fk,m (hk − hi ) (6.125)
k=i+1
siendo ei la excentricidad del esfuerzo cortante Qi,m con respecto al centro de torsión
de la masa mi . Nótese que la excentricidad ei tendrá en general dos componentes en el
plano del forjado. La normativa ASCE 7-16 requiere desplazar el centro de la masa mi
un 5% de la dimensión del edificio en dirección normal a cada dirección considerada.
Existen diversos métodos o criterios para determinar los valores de diseño en un edificio
a partir de las contribuciones modales de las magnitudes anteriores. Habitualmente se
emplea el método SRSS (Square Root of Sum of Squares) que define tales valores de
diseño a partir de la norma ecuclı́dea aplicada a las componentes modales. En el caso de
306 CAPÍTULO 6. DINÁMICA DE ESTRUCTURAS
Se emplean igualmente como valores de diseño las fuerzas dinámicas y los cortantes en
cada piso, además del cortante en la base del edificio, todos ellos calculados como
v v v
u n u n u n
X 2
X 2
uX
(Q0,m )2
u u
Fi = t (Fi,m ) Qi = t (Qi,m ) Q0 = t (6.129)
m=1 m=1 m=1
Finalmente, también se define como valor de diseño el momento en cada piso, el mo-
mento en la base del edificio y el momento torsional, calculados como
v v v
u n u n u n
uX 2
X 2
uX
(M ti,m )2
u
Mi = t (Mi,m ) M0 = t (M0,m ) M ti = t (6.130)
m=1 m=1 m=1
Plantear igualmente el análisis de la respuesta dinámica del pórtico frente a una ex-
citación armónica de amplitud F0 y frecuencia θ. Analizar igualmente si se dan las
condiciones para que el sistema entre en resonancia.
Posteriormente se añade el segundo piso al pórtico inicial, con las caracterı́sticas men-
cionadas. Calcular las frecuencias propias y los modos naturales de vibración del sis-
tema y dibujar esquemáticamente éstos últimos.
La inercia de los pilares, la rigidez de cada pilar individual y la rigidez del primer piso
son:
I1 = 1, 25 · 10−3 m4
12EI1
k1 = 2kpilar = 2 = 27, 8 · 106 N/m
h31
308 CAPÍTULO 6. DINÁMICA DE ESTRUCTURAS
Para el estudio del edificio con dos pisos, se adopta un sistema con 2 grados de libertad
horizontales x1 y x2 en los forjados.
Las masas son conocidas, mientras que las rigideces dependen de la inercia de la sección
transversal y de la altura del piso. En este caso las inercias son iguales para ambos
grados de libertad, pero no lo son las alturas de ambos pisos. Las rigideces para cada
grado de libertad del sistema dinámico son
12EI 12EI
k1 = 2 = 27, 8 · 106 N/m k2 = 2 = 11, 7 · 106 N/m
h31 h32
6.7. ANÁLISIS SÍSMICO DE EDIFICIOS 309
Los autovectores asociados a cada autovalor que definen los modos de vibración de la
estructura se calculan como sigue y se representan en la Figura 6.44.
Mϕ1 Mϕ2
β1 = = 0, 389 β2 = = 0, 611
ϕT1 Mϕ1 ϕT2 Mϕ2
Determinar las frecuencias propias y los modos naturales de vibración del edificio y
dibujar estos últimos. Hallar además los factores de participación de cada modo de
vibración.
La masa de los pilares se calcula a partir del peso especı́fico del hormigón γh y el
volumen de material en cada forjado, que para una superficie A constante y un espesor
variable ei da lugar a una masa de cada forjado calculada como
mi = γh Ati
12EIi
ki =
Hi3
Finalmente, la rigidez de cada piso es la suma de las rigideces de los pilares que lo
sustentan (vs A.5.4) ki = 2kpilar,i . La Tabla 6.9 muestra los resultados obtenidos.
312 CAPÍTULO 6. DINÁMICA DE ESTRUCTURAS
Para el estudio del edificio se adopta un sistema con 3 grados de libertad horizontales
x1 , x2 y x3 en orden creciente con los pisos. El sistema de ecuaciones diferenciales a
resolver es:
0, 940 0 0 ẍ1 148, 2 −47, 4 0 x1 0
106 0 822, 5 0 ẍ2 +106 −47, 4 73, 3 −25, 9 x2 = 0
0 0 705, 0 ẍ3 0 −25, 9 25, 9 x3 0
Los modos de vibración del edificio se representan en la Figura 6.46 a partir de los
autovectores normalizados a su primera componente.
6.7. ANÁLISIS SÍSMICO DE EDIFICIOS 313
Conocidas las frecuencias y los modos propios de vibración, junto con los factores de
participación de cada modo, se realiza un análisis modal de la respuesta del edificio.
Para ello, en primer lugar se determinan los valores espectrales de desplazamiento (Sd ),
velocidad (Sv ) y aceleración (Sa ) a partir del espectro de respuesta en velocidades dado.
Los valores de Sv dependen del periodo natural en cada modo de vibración (Tm ). Los
valores espectrales Sd y Sa se calculan a partir del anterior empleando la frecuencia
natural (ωm ) de cada modo de vibración (o equivalentemente, su periodo Tm ).
En la tercera parte de la Tabla 6.11 se presentan los resultados del análisis modal en
esfuerzos. En primer lugar se determina la contribución en cada piso i de cada modo m
a las fuerzas sı́smicas Fi,m . A partir de estos valores se determinan los cortantes Qi y los
momentos flectores Mi,m , ası́ como los correspondientes valores Qi,m y Mi,m , obtenidos
componiendo las contribuciones modales, según las expresiones de 6.7.3. Nótese la
contribución determinante del primer modo en los esfuerzos máximos obtenidos.
314 CAPÍTULO 6. DINÁMICA DE ESTRUCTURAS
Los valores máximos en cada magnitud se emplean para el diseño y comprobación del
edificio frente a acción sı́smica, tanto en rotura como en deformación
A Fórmulas y desarrollos matemáticos
Sea una viga continua de dos vanos simplemente apoyados, de longitudes l1 y l2 que
soportan cargas distribuidas de valor q1 y q2 respectivamente, como muestra la Figura
A.1. Se supone por simplicidad que l1 ≥ l2 sin pérdida de generalidad del problema,
asumiendo que el estudio de los efectos más desfavorables se da en el vano de mayor
longitud.
Figura A.1: Distribución de cargas distribuidas sobre viga continua de dos vanos
Se calculan las reacciones y los valores máximos de la ley de momentos flectores para
el caso general. Se definen para ello dos constantes para la relación de longitudes y de
cargas como
l1 q2
α := β := (A.1)
l2 q1
En cada caso se particulariza el resultado para el caso habitual de vanos iguales (l1 = l2
o α = 1).
315
316 APÉNDICE A. FÓRMULAS Y DESARROLLOS MATEMÁTICOS
La ley de momentos flectores presenta dos tramos parabólicos con dos máximos flectores
positivos. En el primer vano la ley de momentos es M12 (x) con un valor positivo máximo
max,(+)
M12 situado en xmax
12 :
2
q1 l12 4 + 3α − βα−2 q1 l 2
max,(+)
M12 = = (7 − β)2 (A.3)
128 1+α 512
l1 =l2
−2
max 1 α + βα 1 1+β
x12 = l1 − = l1 − (A.4)
2 8 (1 + α) 2 16
l1 =l2
En el segundo vano la ley de momentos es M23 (x) con un valor positivo máximo
max(+)
M23 situado en xmax
23 :
2
q1 l12 −4 − 3α−1 + β −1 α2 q1 l 2
max,(+) 2
M23 = = 7 − β −1 (A.5)
128 1+α 512
l1 =l2
−1
3 1 1 + β −1
max 1 1+β α
x23 = l1 + l2 + = l1 + l2 + (A.6)
2 8 (1 + α) 2 16
l1 =l2
El resultado muestra que el valor del máximo flector positivo en el vano mayor (12)
es decreciente con el valor de β con lo cual el valor máximo se da cuando β = 0 es
decir cuando q2 = 0 por tanto sin carga en el segundo vano. Nótese igualmente que el
máximo valor del flector positivo siempre se da en el vano más cargado, es decir cuando
β > 1.
El máximo momento flector negativo se da sobre el apoyo central con un valor
q1 l2 α + βα−2 q1 l 2
M max,(−) = − 1 =− (1 + β) (A.7)
8 1 + α l1 =l2 16
ex − e−x ex + e−x
sinh(x) = cosh(x) = (A.9)
2 2
d d
[sinh(x)] = cosh(x) [cosh(x)] = sinh(x) (A.10)
dx dx
A.3 Inestabilidad
1 l 1 l
Z Z Z l
00
2 0
2
Π= EI v (x) dx − N v (x) dx − q(x)v(x)dx = 0 (A.17)
2 0 2 0 0
v̂(x) = C0 + C1 x + C2 x2 + C3 x3 (A.18)
La función v̂(x) viene caracterizada por los cuatro coeficientes del polinomio de tercer
grado que quedan determinadas en función de los movimientos perpendiculares a la
directriz y los giros, en ambos extremos de la viga d0 = {v1 , θ1 , v2 , θ2 } entendidos
como condiciones de apoyo y que constituyen los parámetros incógnita de la función
de aproximación. La expresión de la deformada en términos de los movimientos y los
giros d0 en los extremos a, b es la siguiente:
3(va − vb ) 2θa + θb 2 θa + θb 2(va − vb )
v̂(x) = va + θa x − + x + + x3 (A.19)
l2 l l2 l3
La solución para la deformada v̂(x) ası́ como su primera derivada – ley de giros θ(x) –
y su segunda derivada – ley de curvaturas χ(x) – pueden expresarse matricialmente en
función del vector de movimientos en los extremos (d0 ).
v10
h
3x2 2x3 2x2 x3 3x2 2x3 x3 x2
i θ1
v̂(x) = C(x)δ 0 = 1− + x− + − −
l2 l3 l l2 l2 l3 l2 l v20
θ2
(A.21)
v10
dv̂(x) h
6x2 3x2 6x2 3x2
i θ
1
θ̂(x) = = C0 (x)δ 0 = − 6x
l2
+ l3
1− 4x
l + l2
6x
l2
− l3
− 2x
l + l2
0
dx v2
θ2
(A.22)
v10
d2 v̂(x) 00 0
h i θ1
χ̂(x) = = C (x)δ = − l62 + 12x
− 4l + 6x 6
− 12x
− 2l + 6x
dx2 l3 l2 l2 l3 l2
0
v2
θ2
(A.23)
Z l Z l Z l
1 1
Π = d0T C00 (x)T EIC00 (x)dx − N C0 (x)T C0 (x)dx − d0T q(x)C(x)dx d0
2 0 2 0 0
(A.24)
K00 − N K0G d0 = f 0
(A.25)
en ejes locales
Z l
K00 = C00 (x)T EIC00 (x)dx (A.26a)
0
Z l
K0G = C0 (x)T C0 (x)dx (A.26b)
0
Z l
f0 = q(x)C(x)dx (A.26c)
0
EA
− EA
l 0 0 l 0 0
12EI 6EI
0
l3 l2
0 − 12EI
l3
6EI
l2
6EI 4EI
0
0
l2 l 0 − 6EI
l2
2EI
l
K0 = EA (A.27)
− l EA
0 0 l 0 0
− 12EI − 6EI 12EI
− 6EI
0 0
l3 l2 l3 l2
6EI 2EI
0 l2 l 0 − 6EI
l2
4EI
l
0 0 0 0 0 0
0 − 5l6 1
− 10 0 5l6 − 101
1 2l 1 l
0
0 − 10 − 15 0 10 30
KG = N (A.28)
0 0 0 0 0 0
0 5l6 1
0 − 5l6 1
10 10
1 l 1 2l
0 − 10 30 0 10 − 15
Sea f (t) una función de variable real integrable en el intervalo [t0 − T /2, t0 + T /2]. En
este caso puede obtenerse el desarrollo en serie de Fourier de dicha función f (t) como
∞
a0 X 2nπ 2nπ
f (t) ≈ + an cos t + bn sin t (A.29)
2 T T
n=1
Z T /2 Z T /2 Z T /2
2 2 2
a0 := f (t)dt an := f (t)cos (wn t) dt bn := f (t)sin (wn t) dt
T −T /2 T −T /2 T −T /2
(A.30)
2nπ
donde wn := T es la frecuencia asociada a cada armónico.
Aplicaciones especı́ficas:
a
a cos mπ x
Z mπ
a a 2a
sin x dx = − = (1 − cos(mπ)) = , m impar (A.37)
mπ
0 a a
mπ mπ
0
a2
Z mπ
x sin x dx = (−1)m+1 (A.38)
a mπ
Z
e −at
e−at sin(bt)dt = − 2 [asin(at) + bcos(bt)] (A.39)
a + b2
rk+1
Z
k 1
r ln rdr = ln r − (A.40)
k+1 k+1
Z
Para el caso k=0 → ln rdr = r (ln r − 1) (A.41a)
r2
Z
1
Para el caso k=1 → rln rdr = ln r − (A.41b)
2 2
ln r r−1 1 rk
lim rk ln r = lim = lim = lim = lim − =0 (A.42)
r→0 r→0 r −k r→0 −kr −k−1 r→0 −kr −k r→0 k
A.4. CÁLCULO DE PLACAS 325
Las transformaciones de variables entre coordenadas cartesianas y polares son las sigu-
ientes:
( ( p
x = rsinϕ r = x2 + y 2
(x, y) ≡ ⇒ (r, ϕ) ≡ (A.43)
ϕ = arctan xy
y = rcosϕ
Las derivadas de las coordenadas cartesianas con respecto a las coordenadas polares
son las siguientes:
∂r ∂r ∂ϕ ∂ϕ
∂x = cosϕ ∂y = sinϕ ∂x = − 1r sinϕ ∂y = 1r cosϕ (A.44)
En el segundo paso se han separado los términos con exponente par (2n) e impar (2n+1)
sin variar el contador lo que no supone pérdida de generalidad en el desarrollo.
Ambos sumatorios se identifican con el desarrollo en serie de Taylor de las funciones
cos(x) y sin(x) respectivamente definidas al principio.
Asin (ωt + ϕ) = (Acos ϕ) sin (ωt) + (Asin ϕ) cos (ωt) = C1 sin (ωt) + C2 cos (ωt) (A.50)
A.5. CÁLCULO DINÁMICO 327
Sea x1 (t) una solución de la ecuación diferencial de segundo orden homogénea del tipo
Se define x2 = u(t)x1 (t) siendo u(t) una función continua y derivable en el intervalo
de definición. Derivando sucesivamente la función x2 (t) y sustituyendo en la ecuación
diferencial se obtiene la expresión para u(t) que permite definir x2 (t):
R
e− p(t)dt
Z
u(t) = dt (A.53)
x21 (t)
Se verifica por tanto que el Wronskiano de las funciones x1 (t) y x2 (t) es no nulo en
todo el dominio de definición de las funciones, por lo que ambas funciones conforman
un sistema fundamental de soluciones para la ecuación diferencial.
Sea un pórtico formado por dos pilares de altura h con inercias I1 , I2 separados por
un forjado de ancho l de inercia If , sometido a una carga P horizontal sobre el forjado
asumiendo éste como inextensible, según se muestra en la Figura A.2.
En estas condiciones, y asumiendo que ambos pilares tienen la misma sección transversal
y por tanto I1 = I2 = I, la rigidez del sistema es
24EI
K= (A.60)
h3
330 APÉNDICE A. FÓRMULAS Y DESARROLLOS MATEMÁTICOS
B Tablas
VIGAS EN VOLADIZO
CARGAS yB (+ ↑) φB (+ ) ELÁSTICA y(x) ( + ↑)
P l3 P l2 P
− − (x3 − 3lx2 )
3EI 2EI 6EI
pl4 pl3 p
− − − (x4 − 4lx3 + 6l2 x2 )
8EI 6EI 24EI
M x2
0≤x≤a
Ma Ma 2EI
(2l − a)
2EI EI
Ma
(2x − a) a≤x≤l
2EI
Tabla B.1: Flechas y giros en vigas en voladizo
331
332 APÉNDICE B. TABLAS
VIGAS BIARTICULADAS
CARGAS GIROS (+ ) ELÁSTICA y(x) ( + ↑ )
P ab
φA = − (l + b)
6EIl P bx 2
− (l − b2 − x2 ) 0 ≤ x ≤ a
P ab 6EIl
φB = (l + a)
6EIl
pl3
φA = −
24EI px
− (x3 − 2lx2 + l3 )
24EI
pl3
φB =
24EI
3pl3 px
16x3 −24lx2 +9l3
φA = − − 384EI
128EI 0 ≤ x ≤ l/2
pl
8x3 −24lx2 + + 17l2 x − l3
− 384EI
7pl3
φB = l/2 ≤ x ≤ l
384EI
pl3
φA = −
45EI pl3 x
x2 x4
− 7 − 10 2 + 3 4
360EI l l
7pl3
φB =
360EI
b2 x2
M l b2 −M lx
6EI 1 − 3 l2 − l2
φA = (3 − 1)
6EI l2
h0 ≤ x ≤ a i
M l(l−x) a2 l−x 2
Ml a2 6EI 1 − 3 l2 − l
φB = (3 − 1)
6EI l2 a≤x≤l
Ml
φA = −
6EI Mx 2
− (l − x2 )
Ml 6EIl
φB =
3EI
VIGAS BIEMPOTRADAS
CARGAS M0AB (+ ) M0BA (+ )
P ab2 P a2 b
Pl Pl
2
= − 2 =−
l
a=b 8 l a=b 8
pl2 pl2
−
12 12
pl2 pl2
−
20 30
5pl2 5pl2
−
96 96
Mb M Ma M
2
(2a − b) = 2
(2b − a) =
l a=b 4 l a=b 4
6EI 6EI
∆ ∆
l2 l2
VIGAS EMPOTRADO-ARTICULADAS
CARGAS (M0AB )∗ (+ ) φB (+ )
P a2 b P l2
Pb 2 2
3P l
2
(l − b ) = =
2l
a=b 16 4EIl a=b 32EI
pl2 pl3
8 48EI
pl2 pl3
15 120EI
5pl2 5pl3
64 384EI
M 2 2
M 2 2
l − 3b = l − 4lb + 3b =
2l2 a=b 4EIl a=b
M Ml
= =−
8 16EI
3EI 3∆
∆ −
l2 l
Múltiplos y submúltiplos en el SI
Prefijo Sı́mbolo Factor multiplicativo
tera T 1012 = 1 000 000 000 000
giga G 109 = 1 000 000 000
mega M 106 = 1 000 000
kilo k 103 = 1 000
hecto h 102 = 100
deca da 101 = 10
deci d 10−1 = 0,1
centi c 10−2 = 0,01
mili m 10−3 = 0,001
micro µ 10−6 = 0,000001
nano n 10−9 = 0,000 000 001
pico p 10−12 = 0,000 000000 001
ALFABETO GRIEGO
A α Alpha I ι Iota P ρ Rho
B β Beta K κ Kappa Σ σ Sigma
Γ γ Gamma Λ λ Lambda T τ Tau
∆ δ Delta M µ Mu Υ υ Upsilon
E Epsilon N ν Nu Φ φ Phi
Z ζ Zeta Ξ ξ Xi X χ Chi
H η Eta O o Omicron Ψ ψ Psi
Θ θ Theta Π π Pi Ω ω Omega
[1] K. J. Bathe. Finite Element procedures. Prentice Hall, Pearson Education, Inc.,
2nd edition, 2014.
[2] Neal B.G. The plastic methods of structural analysis. Chapman & Hall, 3rd edition,
1977.
[3] J. Blaauwendraad. Plates and FEM. Surprises and Pitfalls. Springer, 2010.
[11] J. Heyman. Beams and Framed Structures. Pergamon Press, 2nd edition, 1974.
337
338 BIBLIOGRAFÍA
La respuesta dinámica de un pórtico de hormigón ante vibraciones libres y amortiguadas se determina considerando varios parámetros clave: la frecuencia natural del sistema, el amortiguamiento, y la respuesta espectral. La frecuencia natural se obtiene a partir del periodo de oscilación y la rigidez de la estructura . El grado de amortiguamiento, que afecta la reducción de las vibraciones, generalmente se expresa como un porcentaje del amortiguamiento crítico y depende del tipo de material y las conexiones estructurales . Para evaluar la respuesta bajo excitaciones dinámicas, se utilizan espectros de respuesta, que indican desplazamientos máximos, velocidades y aceleraciones relativas al movimiento del suelo debido a un evento sísmico . La solución de la ecuación diferencial del movimiento, que incluye fuerzas de inercia y amortiguamiento, permite predecir la respuesta del sistema bajo vibraciones libres . En la práctica, se adoptan modelos simplificados como espectros de diseño para evaluar las acciones sísmicas .
La frecuencia natural es crucial en la planificación de estructuras dinámicas porque influye directamente en su respuesta ante cargas variables, como las sismos. Las estructuras exhiben una respuesta máxima cuando la frecuencia de la carga es cercana a su frecuencia natural, lo que puede llevar a fenómenos de resonancia, aumentando significativamente las amplitudes de vibración . Además, en la práctica, solo un número limitado de modos de baja frecuencia tienen una contribución significativa en la respuesta estructural, permitiendo simplificar el análisis modal . La correcta identificación y consideración de la frecuencia natural es fundamental para diseñar estructuras seguras y funcionales al minimizar los riesgos de resonancia y garantizar que soporten las cargas dinámicas durante su vida útil .
El momento flector y el esfuerzo cortante en las secciones A y B de una viga continua de dos vanos se determinan utilizando las líneas de influencia y las condiciones de carga aplicadas. Para la sección A, la línea de influencia del momento flector es MA = 15V2 = 0,75 (30 − x) y el esfuerzo cortante es TA = −V2 = x − 30 / 20 para 0 < x ≤ 15. Para la sección B, el esfuerzo cortante cuando la carga está en el voladizo es igual a la unidad, y el momento flector es MB = 10 − x, mientras que en el vano central ambos son nulos . Para obtener estos valores, se consideran las condiciones de equilibrio de la viga, donde las reacciones en los apoyos se equilibran con las cargas aplicadas, y estos esfuerzos extremos se maximizan según las posiciones de carga más desfavorables utilizando las envolventes de esfuerzos . Además, en vigas hiperestáticas es importante considerar la plasticidad para determinar si se alcanzan las condiciones de colapso, influenciando así la distribución de los momentos flectores a lo largo de la estructura . El cálculo del momento flector y del esfuerzo cortante también se beneficia del uso de matrices de rigidez y métodos incrementarles donde las secciones críticas pueden plastificar y convertirse en rótulas plásticas, permitiendo redistribuir el momento a otras secciones hasta alcanzar el estado de colapso . Este enfoque se fundamenta en la superposición de patrones de carga incrementales y en el análisis de puntos de rótula para definir la ley de momentos flectores a través de los vanos .
El criterio para seleccionar el modelo de análisis límite plástico en estructuras reticuladas con respecto a los mecanismos de colapso se basa en la identificación de los mecanismos que presentan menores cargas últimas de colapso. Se deben determinar las secciones susceptibles de formar rótulas plásticas y elegir aquellos mecanismos de colapso elementales o combinados que activen un número suficiente de éstas para inducir el colapso de la estructura . La carga de colapso se calcula mediante el principio de trabajos virtuales, eligiendo la carga mínima entre los mecanismos evaluados . El modelo de análisis asegura que la carga última de colapso es mayor o igual a la carga última de la estructura, cumpliendo el teorema cinemático y verificándose con el teorema estático .
El proceso de formación de un mecanismo de colapso en el análisis plástico de estructuras reticuladas implica varias etapas. Primero, se identifican las secciones donde pueden desarrollarse rótulas plásticas, que generalmente incluyen apoyos, nudos y puntos de aplicación de cargas . A medida que las cargas aumentan, las rótulas plásticas empiezan a formarse en estas secciones críticas hasta que se alcanzan al menos n+1 rótulas para provocar el colapso en estructuras con un grado de hiperestatismo n, lo que se conoce como un mecanismo completo . Este proceso de formación de rótulas se realiza mediante la redistribución de momentos flexores, manteniendo constante el momento plástico MP en las secciones donde ya se ha alcanzado . Para cada posible mecanismo de colapso identificado, se calcula la carga última utilizando el principio de los trabajos virtuales y se mantiene como carga última de colapso la menor de todas las calculadas . Finalmente, se verifica si la carga última calculada cumple el teorema estático, asegurando que la sumatoria de momentos flexores en todas las secciones no supera el valor de MP .
La estabilidad de una viga sometida a un esfuerzo axial, según la teoría de Euler, se evalúa mediante el cálculo de la carga crítica de pandeo, conocida como carga crítica de Euler (NE). Esta carga se determina considerando la longitud de la viga, la rigidez a flexión (EI) y las condiciones de apoyo. La ecuación diferencial de equilibrio y las condiciones de apoyos determinan la forma de la curva de pandeo y el valor de la carga crítica en diferentes configuraciones de la viga (articulada, empotrada, biempotrada, etc.). Para una viga biarticulada, la carga crítica se expresa como NE = π²EI/l² . Cuando la carga axial aplicada (N) alcanza o supera NE, la viga experimenta inestabilidad, manifestada como pandeo, es decir, una deformación lateral significativa sin un incremento adicional en la carga aplicada . La carga crítica no depende de cargas transversales sino del axil aplicado y las condiciones de apoyo .
La carga distribuida uniformemente afecta el cálculo de momentos flectores en una viga continua con voladizo al generar un momento negativo máximo en la sección del apoyo y un momento positivo máximo en el vano, dependiendo de la longitud del voladizo relativo al vano . Para una viga compuesta por un vano de longitud l y un voladizo de longitud αl bajo carga q, la ubicación y magnitud de los momentos flectores máximos varían con α. El colapso puede ocurrir por la aparición de una rótula plástica en el vano o en el apoyo, con momentos dependientes de α. El momento máximo positivo Mmax,(+) se sitúa en el vano mientras que el negativo Mmax,(−) se sitúa en el apoyo, con cargas de colapso dadas por la relación entre estos momentos y el momento plástico de la sección MP .
El análisis dinámico de estructuras difiere entre sistemas no amortiguados y amortiguados principalmente en la consideración de las fuerzas de amortiguamiento. En un sistema no amortiguado, las ecuaciones de movimiento incluyen solo la fuerza elástica y la de inercia, resultando en oscilaciones sostenidas a la frecuencia natural del sistema, sin pérdida de energía . En contraste, un sistema amortiguado incorpora una fuerza adicional proporcional a la velocidad del sistema, que disipa la energía vibratoria del sistema, haciendo que las oscilaciones disminuyan gradualmente con el tiempo . El efecto del amortiguamiento es crucial para prever una disminución de la amplitud de las oscilaciones, afectando así el análisis y diseño estructural, particularmente en la respuesta sísmica . Un sistema amortiguado es más representativo de las condiciones reales donde siempre existe alguna forma de pérdida de energía .
La teoría de Euler se aplica al diseño de columnas esbeltas bajo cargas de compresión calculando la carga crítica de pandeo, la cual es aquella a la que una columna larga se vuelve inestable y falla lateralmente antes de que el material alcance su resistencia máxima. La carga crítica de pandeo está relacionada con la esbeltez de la columna, su longitud y sus condiciones de apoyo. Columnas más esbeltas son más propensas al pandeo antes de agotar el material . El fenómeno ocurre dentro del rango elástico, lo que significa que, al retirar la carga, la columna puede volver a su forma original .