C.1.5MDOC Análisis de Estructuras
C.1.5MDOC Análisis de Estructuras
C.1.5MDOC Análisis de Estructuras
ELABORACIÓN
REVISIÓN
COORDINACIÓN
En la realización de este capítulo del MDOC-CFE colaboraron personal técnico1 y asesor externo4
de la Comisión Federal de Electricidad, personal académico del Instituto de Ingeniería – UNAM2 y
el Subdirector de Proyectos y Construcción de CFE cuando se elaboró este capítulo3.
AGRADECIMIENTOS
CONTENIDO
Presentación..................................................................................................XV
1. INTRODUCCIÓN ......................................................................................... 1
1.1 Objetivos del análisis estructural ............................................................ 1
1.2 Elementos, formas y sistemas estructurales ............................................. 2
1.3 Grados de libertad ................................................................................ 3
1.4 Linealidad y no linealidad de estructuras .................................................. 4
1.5 Principio de superposición ...................................................................... 5
1.6 Modelado de la estructura ...................................................................... 6
viii
MANUAL DE DISEÑO DE OBRAS CIVILES-CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
ix
C.1.5 ANÁLISIS DE ESTRUCTURAS
CONTENIDO
x
MANUAL DE DISEÑO DE OBRAS CIVILES-CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
xi
PRÓLOGO
En esta coyuntura de importantes retos para el Sector Eléctrico Nacional, pero también
de grandes oportunidades, se inscribe la presente obra. El Manual de Diseño de Obras
Civiles contribuye al cumplimiento de las ambiciosas metas planteadas por la reforma y
por el Programa Nacional de Infraestructura 2014-2018, que buscan impulsar no sólo el
desarrollo de más y mejor infraestructura, sino también la formación de cuadros técnicos
altamente capacitados en el diseño y construcción de estos proyectos.
Así, esta nueva edición del Manual elaborado por la CFE y la Universidad Nacional
Autónoma de México provee los lineamientos de diseño de todas las obras de ingeniería
civil, y lo hace incorporando los avances tecnológicos y la experiencia técnica acumulada
por la ingeniería mexicana desde 1983, año en que se editó el último manual.
Sin duda, esta obra, que conjunta elementos teóricos y prácticos de la ingeniería en los
campos de la Hidrotecnia, Geotecnia y Estructuras, constituye una referencia técnica
obligada no sólo para los ingenieros de la CFE, sino para todos aquellos encargados de
construir obras de ingeniería civil en PEMEX y dependencias del sector público federal y
local. Al igual que su predecesor, este manual único en su tipo será de gran ayuda técnica
para los ingenieros de toda América Latina.
Quiero destacar que con el Manual de Diseño de Obras Civiles, la CFE refrenda su
compromiso con la formación de nuevos cuadros de ingeniería. En efecto, esta
publicación está pensada para ser consultada y utilizada en las aulas por nuestros
estudiantes de ingeniería, no sólo en el ramo civil, sino en diversas disciplinas.
La presente edición actualiza las versiones previas e incorpora capítulos nuevos, con una
visión en la que se incluyen los múltiples avances tecnológicos alcanzados en las últimas
décadas. El signo de cambio desde 1983 está asociado al uso de las computadoras en la vida
cotidiana. Así, sin menoscabo de la exposición de los conocimientos invariantes y los métodos
de diseño convencionales, en la presente versión del Manual se incluyen los procedimientos
basados en los métodos numéricos computacionales, discutiendo sus capacidades y
limitaciones.
El Manual actual constituye una herramienta de utilidad para quienes se dedican al diseño de
las estructuras civiles que forman parte de una central de generación eléctrica, o de un
sistema de distribución eléctrica. Si bien este Manual sirve fundamentalmente para este fin
como una ayuda práctica para los ingenieros de CFE, se considera que, en general, también
será útil para los ingenieros civiles que ejercen su profesión en otros campos. Adicionalmente,
en la elaboración del Manual se ha procurado que cumpla con una misión didáctica, al
constituirse en un documento auxiliar para que el estudiante de ingeniería se familiarice con
los temas que desarrollará en su vida profesional.
El Manual fue elaborado por académicos del Instituto de Ingeniería de la Universidad Nacional
Autónoma de México, en colaboración con ingenieros de la Subdirección de Proyectos y
Construcción a la que pertenecen especialistas de las Coordinaciones de Proyectos
Hidroeléctricos, Proyectos Termoeléctricos, Líneas de Transmisión y Transformación, y la
Gerencia de Estudios de Ingeniería Civil de la Comisión Federal de Electricidad. Es de interés
para la CFE difundir esta obra que se ha consolidado como el manual de consulta en ingeniería
civil más requerido no sólo en México, sino también por los ingenieros de otros países de
habla hispana.
Con la estrecha y fructífera cooperación lograda entre una de las más prestigiadas
instituciones dedicadas a la investigación y la docencia del país, y la Comisión Federal de
Electricidad, cuya misión es el suministro de electricidad con altos estándares de calidad al
menor costo posible, para beneficio de México, se ha logrado actualizar y complementar este
documento de alto valor técnico y académico, que será de mucha utilidad para el desarrollo
de nuestra infraestructura eléctrica.
Ing. César Fernando Fuentes Estrada
Dr. Luis Agustín Álvarez - Icaza Longoria
SIMBOLOGÍA
𝑎𝑒𝑓 altura del elemento finito
𝑎𝑙 longitud medida del origen a la posición de la fuerza P
𝑎𝑚𝑒𝑓 base menor de elemento finito trapezoidal
𝐚 término que hace que el campo de desplazamientos sea
continuo
𝐴 área de la sección transversal
𝐴𝑐 área efectiva de cortante
𝐴𝑛 , 𝐵𝑛 constantes de integración
𝐴𝑚 área promedio de la sección transversal
𝐴𝑣a área de la sección transversal de una viga de sección variable
en el nodo a
𝐴(𝑥) área de la sección transversal de una viga a una distancia x
𝐴1 área de la sección transversal en el nodo 1
𝐴2 área de la sección transversal en el nodo 2
𝐀 matriz de autovectores
𝑏 fuerza de cuerpo
𝑏𝑒𝑓 ancho del elemento finito
𝑏𝑓 ancho del elemento en la dirección normal a éste
𝑏𝑙 longitud medida de la posición de la carga P al nodo final del
elemento
𝑏𝑚𝑒𝑓 base mayor de elemento finito trapezoidal
𝑏𝑝𝑥 longitud de un elemento plano en dirección x
𝑏𝑝𝑦 longitud de un elemento plano en dirección y
𝑏𝑣 ancho de la sección transversal de una viga
𝑏𝑥 , 𝑏𝑦 , 𝑏𝑧 fuerzas de cuerpo en dirección 𝑥, 𝑦, 𝑧
𝐛 vector de fuerzas de cuerpo
𝐵1 𝑦 𝐵2 factores de amplificación
𝐁 matriz de transformación
𝑐 mecanismo de amortiguamiento
𝑐11 , 𝑐12 , 𝑐13 , 𝑐14 … 𝑐66 coeficientes de la matriz constitutiva simétrica del material
CG carga gravitacional
𝐶𝑚 coeficiente que se determina considerando que no hay
traslación lateral del marco
𝐂 matriz constitutiva simétrica del material
𝑪 matriz de amortiguamiento
C.1.5 ANÁLISIS DE ESTRUCTURAS
SIMBOLOGÍA
d𝑥1 , d𝑦1 , d𝑧1 , d𝑥2 , d𝑦2 , d𝑧2 coeficientes del vector global de desplazamientos de una
armadura 3D
1
𝑑𝑎𝑧 desplazamiento de azotea, debido al primer incremento de
carga lateral
2
𝑑𝑎𝑧 desplazamiento de azotea, debido al segundo incremento de
carga lateral
𝑑𝑛𝑥 , 𝑑𝑛𝑦 , 𝑑𝑛𝑧 … 𝑑𝑛𝑓 coeficientes del vector nodal de desplazamientos
′
𝑑𝑥1 ′
, 𝑑𝑦1 ′
, 𝑑𝑧1 ′
, 𝜃𝑥1 ′
, 𝜃𝑦1 ′
, 𝜃𝑧1 ′
, 𝑑𝑥2 ′
, 𝑑𝑦2 ′
, 𝑑𝑧2 ′
, 𝜃𝑥2 ′
, 𝜃𝑦2 ′
, 𝜃𝑧2 coeficientes del vector local de
desplazamientos de un elemento cable,
viga, marco y marco 3D
′
𝑑1′ , 𝑑2′ … 𝑑𝑛𝑑 coeficientes del vector local de desplazamientos de un
elemento de 𝑛𝑑 nodos
𝐝 vector global de desplazamientos
𝐝′ vector local de desplazamientos
𝐝𝑔 vector espectral global de grados de libertad
𝒅𝒏 vector nodal de desplazamientos
𝐷 módulo de rigidez a flexión equivalente de la placa
𝐃𝑓 vector de desplazamientos de nodos no restringidos
𝐃𝑠 vector de desplazamientos de nodos restringidos
𝐸 módulo de elasticidad o de Young
𝐸𝐼𝑍 rigidez elástica a flexión
𝐸𝑖,max elementos mecánicos máximos del modo 𝑖
𝐸max elementos mecánicos máximos combinados de todos los
modos
𝐸𝑄1 , 𝐸𝑄2 , 𝐸𝑄3 vector de cargas laterales necesario para provocar la primera
articulación plástica
f𝑐 matriz de coeficientes de flexibilidad de un elemento de la
estructura completa
𝑓ℎ flujo de calor por unidad de área
𝑓ℎ𝑦 , 𝑓ℎ𝑥 flujo de calor por unidad de área paralelo al gradiente de
temperatura
𝑓𝑘𝑖 distribución de rigidez del nodo i
𝑓𝑠 fuerza de superficie
𝑓𝑣1 fuerza vertical generada debido a la tensión del cable en el
nodo 1
𝑓𝑣2 fuerza vertical generada debido a la tensión del cable en el
nodo 2
𝑓𝑝𝑥1 , 𝑓𝑝𝑦1 , 𝑓𝑝𝑥2 , 𝑓𝑝𝑦2 , 𝑓𝑝𝑥3 , 𝑓𝑝𝑦4 , 𝑓𝑝𝑥4 , 𝑓𝑝𝑦4 coeficientes de un vector de fuerzas de un elemento
plano, en N
𝑓𝑥1 , 𝑓𝑦1 , 𝑓𝑧1 , 𝑀𝑥1 , 𝑀𝑦1 , 𝑀𝑧1 , 𝑓𝑥2 , 𝑓𝑦2 , 𝑓𝑧2 , 𝑀𝑥2 , 𝑀𝑦2 , 𝑀𝑧2 coeficientes del vector local de fuerzas de
un elemento cable, viga, marco y marco
3D
𝑓𝐷 (𝑡) fuerza de amortiguamiento en función del tiempo
𝑓𝑆 (𝑡) fuerza del resorte en función del tiempo
𝑓11 , 𝑓12 , 𝑓13 , 𝑓14 , 𝑓15 … 𝑓66 coeficientes de la matriz de flexibilidades
𝑓𝐼 (𝑡) fuerza de inercia en función del tiempo
xviii
MANUAL DE DISEÑO DE OBRAS CIVILES-CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
xix
C.1.5 ANÁLISIS DE ESTRUCTURAS
SIMBOLOGÍA
ℎ𝑣 peralte de la viga
ℎ𝑣 (𝑥) altura de la sección transversal de una viga de sección
variable a una distancia 𝑥
𝐻𝑢̆ suma de un desplazamiento en un campo discontinuo más un
desplazamiento en un campo continuo
ℋ𝛤𝑑 término que hace que el campo de desplazamientos sea
discontinuo
𝐼𝑥 , 𝐼𝑦 y 𝐼𝑧 momentos de inercia de la sección transversal del elemento
𝐼𝑧 (𝑥) momento de inercia de la sección transversal de una viga a una
distancia x
𝐼𝐼 , 𝐼𝐼𝐼 , 𝐼𝐼𝐼𝐼 invariantes del tensor de esfuerzos
𝚰 vector de invariantes
𝑰 matriz identidad
𝐽 constante de torsión de la sección
𝑘𝑛 , 𝑘𝑡 constantes que representan la rigidez de la interfaz para el
modo I y II, respectivamente
𝑘𝑝 coeficiente teórico crítico de pandeo
k𝑟 rigidez de un elemento resorte
𝑘𝑟 coeficiente de la matriz local de un elemento resorte
𝑘𝑟𝑛 rigidez rotacional del elemento n
𝑘𝑟𝑜 rigidez de un resorte elástico de un oscilador de un grado de
libertad
k matriz local de rigidez
𝐤𝐺 matriz local de rigidez geométrica
𝐤𝑇 matriz local de rigidez total en coordenadas globales
𝒌𝑡𝑒 conductividad térmica del material
KL longitud efectiva de la columna
𝑘11 , 𝑘12 , 𝑘13 , 𝑘14 … 𝑘𝑛𝑛 coeficientes de la matriz global de rigidez
K11 , K12 , K 21 , K 22 submatrices de la matriz global de rigidez
𝐴
K11 𝐵
, K11 𝐶
, K11 𝐷
, K11 submatrices de la matriz global de rigidez de los elementos
estructurales A, B, C, D.
𝐾∗ módulo volumétrico
̅
𝐾 matriz de rigidez efectiva
̅𝑡
𝐾 matriz de rigidez efectiva triangular
K matriz global de rigidez
𝑲 matriz de rigidez de la ecuación de equilibrio dinámico
𝐊 (𝑛) matriz global de rigidez para para una enésima iteración
𝐊 (𝑡) matriz global de rigidez para un tiempo determinado
𝐊 (𝑛−1) matriz global de rigidez para una iteración 𝑛 − 1
𝐊 𝑐𝑎 matriz de rigidez de un elemento cascarón
𝐊𝑓𝑓 sub-matriz de rigidez con condiciones de apoyo no
restringidos
𝐊𝑓𝑠 , 𝐊 𝑠𝑓 sub-matriz de rigideces con condiciones de apoyos restringido
y no restringidos
𝐊𝒍 matriz de rigidez lateral
𝐊 𝑠𝑠 sub-matriz de rigideces con condiciones de apoyo restringido
xx
MANUAL DE DISEÑO DE OBRAS CIVILES-CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
xxi
C.1.5 ANÁLISIS DE ESTRUCTURAS
SIMBOLOGÍA
xxii
MANUAL DE DISEÑO DE OBRAS CIVILES-CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
xxiii
C.1.5 ANÁLISIS DE ESTRUCTURAS
SIMBOLOGÍA
2
𝑃𝑖,𝑒𝑞 fuerza equivalente del nivel 2
3
𝑃𝑖,𝑒𝑞 fuerza equivalente del nivel 3
4
𝑃𝑖,𝑒𝑞 fuerza equivalente del nivel 4
𝑃𝑙𝑡 fuerza axial de primer orden, causada únicamente por
traslación lateral del marco
𝑃𝑛𝑡 fuerza axial de primer orden, considerando que no hay
traslación lateral del marco
𝑃𝑟 resistencia axial de primer orden requerida
𝑃𝑢 resistencia axial de segundo orden requerida
𝑃𝑥1 , 𝑃𝑦1 , 𝑃𝑧1 , 𝑀𝑥1 , 𝑀𝑦1 , 𝑀𝑧1 fuerzas en el nodo 1 en un elemento marco 3D
𝑃𝑥2 , 𝑃𝑦2 , 𝑃𝑧2 , 𝑀𝑥2 , 𝑀𝑦2 , 𝑀𝑧2 fuerzas en el nodo 2 en un elemento marco 3D
𝑃𝐶 carga crítica
𝑃𝐸 carga de Euler
𝑃1 , 𝑃2 coeficiente del vector local de fuerzas de un elemento resorte
P𝐼 , P𝐼𝐼 fuerza 𝐼 y 𝐼𝐼
𝑃∗ carga unitaria aplicada a una distancia x
𝐏 vector de fuerzas
𝐏𝑓 vector de fuerza aplicada en nodo no restringido
𝐏𝑠 vector de fuerza aplicada en nodo restringido
𝑷𝑒𝑞 vector de cargas equivalente
𝑞 carga uniformemente distribuida
𝑞(𝑥) carga uniformemente distribuida a una distancia x
𝑞𝑛 (𝑡) desplazamiento con respecto al tiempo de un oscilador de un
grado de libertad
𝑄 flujo de calor absorbido del medio
𝐐 matriz de anulación
𝐐𝑻 vector de cargas térmicas
𝑟𝑥 𝑋, 𝑟𝑥 𝑌, 𝑟𝑥 𝑍 … 𝑟𝑧 𝑍 coeficientes de la submatriz 𝐫𝒐
𝐫, 𝐫′ radio vector con componentes
𝐫𝒐 submatriz del tensor de transformación
R radio del segmento del arco
𝑅 residuo
𝑅𝑟𝑖 fuerza total acumulada del edificio
𝑅(𝜹𝑎𝑣 ) cociente de Rayleigh
𝑅𝑒𝑖 suma de las fuerzas equivalentes de los entrepisos
𝑅𝐴 reacción en el apoyo A
𝑅𝑀 factor de arriostramiento
𝐑 𝑑𝑒 vector de fuerza desbalanceada
𝐑 matriz de rotación
𝐑o vector de fuerza desbalanceada inicial
𝑆𝑎 seudoaceleración
𝑆𝑑 desplazamiento espectral
𝐒𝑔 matriz espectral global
𝑡∗ tracciones prescritas
𝑡 tiempo
𝑡𝑐 tracciones
xxiv
MANUAL DE DISEÑO DE OBRAS CIVILES-CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
xxv
C.1.5 ANÁLISIS DE ESTRUCTURAS
SIMBOLOGÍA
u1 desplazamiento debido a P1
u2 desplazamiento debido a P 2
u̇ (𝑡) velocidad en función del tiempo de un sistema de un grado de
libertad
u̇ 𝑖 (𝑡) velocidad del modo 𝑖
u̇ 0 velocidad inicial
u̇ 1 velocidad inicial en un tiempo 1
ü aceleración de la gravedad
ü (𝑡) aceleración en función del tiempo de un sistema de un grado
de libertad
ü os aceleración máxima de la coordenada modal
ü 𝑖 (𝑡) aceleración del modo 𝑖
ü𝑥 ,ü𝑦 ,ü𝑧 aceleración de la gravedad en dirección 𝑥, 𝑦, 𝑧
ü 0 aceleración inicial de la gravedad
ü 1 aceleración inicial en un tiempo 1
𝑢, 𝑣, 𝑤, 𝑢𝑥 , 𝑢𝑦 , 𝑢𝑧 coeficiente del vector de desplazamiento en coordenadas 𝑥, 𝑦,
𝑧,
𝑢(𝑟, 𝑧) coeficientes del vector de desplazamientos radiales de un
elemento axisimétrico
𝑢ei desplazamiento en el entrepiso i
𝑢𝑖,𝑚𝑎𝑥 desplazamientos máximos asociado al modo 𝑖, en m
𝑢𝑥 , 𝑢𝑦 , 𝑢𝑧 coeficiente del vector de desplazamiento
𝑢1 , 𝑢2 coeficiente del vector local de desplazamientos de un
elemento resorte y armadura
𝑢1 , 𝑣1 , 𝑤1 , 𝑢2 , 𝑣2 , 𝑤2 , 𝑢3 , 𝑣3 , 𝑤3 , 𝑢4 , 𝑣4 , 𝑤4 desplazamientos del nodo 1 a 4 de un elemento
cascarón
𝑢∗ desplazamiento prescrito
𝑢 ′
desplazamiento campo discontinuo
𝑢̆ desplazamiento campo continuo
𝑢̈ 𝑆 histograma de las aceleraciones en el tiempo
𝐮 vector de desplazamientos
𝐮ℎ desplazamientos discontinuos
𝐮𝑡 vector de desplazamientos con respecto al tiempo
𝐮𝑡−Δ𝑡 vector de desplazamientos con respecto a un incremento de
tiempo anterior
𝐮 (𝑛)
vector desplazamiento para una enésima iteración
𝐮(𝑡) vector de desplazamiento para un tiempo determinado
𝐮 (𝑡+Δ𝑡)
vector de desplazamiento para un n-iesimo tiempo
𝐮(𝑛−1) vector desplazamiento para una iteración 𝑛 − 1
𝐮̇ 𝑡 vector de velocidades con respecto al tiempo
𝐮̇ 𝑡−Δ𝑡 vector de velocidades con respecto a un incremento de
tiempo anterior
𝐮̈ 𝒕 vector de aceleraciones con respecto al tiempo
𝐮̈ 𝑡−Δ𝑡 vector de aceleraciones con respecto a un incremento de
tiempo anterior
1 2 3 𝑛
𝐮 ,𝐮 ,𝐮 ,𝐮 vector de desplazamiento para diferentes incrementos de
fuerza
xxvi
MANUAL DE DISEÑO DE OBRAS CIVILES-CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
xxvii
C.1.5 ANÁLISIS DE ESTRUCTURAS
SIMBOLOGÍA
xxix
C.1.5 ANÁLISIS DE ESTRUCTURAS
SIMBOLOGÍA
xxxi
MANUAL DE DISEÑO DE OBRAS CIVILES-CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
1. INTRODUCCIÓN
1
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.1 INTRODUCCIÓN
2
MANUAL DE DISEÑO DE OBRAS CIVILES-CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
3
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.1 INTRODUCCIÓN
x2
2
Posición
ux1 deformada
1
uy2
uy1
Nodo 1 Nodo 2
Posición original
donde
u𝑥1 desplazamiento del nodo 1 en dirección x, en m
u𝑦1 desplazamiento del nodo 1 en dirección y, en m
u𝑥2 desplazamiento del nodo 2 en dirección x, en m
u𝑦2 desplazamiento del nodo 2 en dirección y, en m
θ1 rotación del nodo 1, en rad y/o grados
θ2 rotación del nodo 2, en rad y/o grados
No obstante de que los resultados obtenidos de análisis elásticos lineales pueden ofrecer
resultados satisfactorios, particularmente cuando la estructura está bajo una acción de
4
MANUAL DE DISEÑO DE OBRAS CIVILES-CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
u𝑡 = u1 + u2 (1.1)
1
𝑝𝑡 = 𝑝𝑖𝑛 2
+ 𝑝𝑖𝑛 (1.2)
donde
P1 , P 2 sistema de fuerzas aplicadas a un elemento y/o estructura, en N
1
𝑝𝑖𝑛 fuerza interna debido a la fuerza P1 , en N
2
𝑝𝑖𝑛 fuerza interna debido a la fuerza P 2 , en N
𝑝𝑡 fuerza interna total, en N
u1 desplazamiento debido a P1 , en m
u2 desplazamiento debido a P 2 , en m
u𝑡 desplazamiento total, en m
5
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.1 INTRODUCCIÓN
6
MANUAL DE DISEÑO DE OBRAS CIVILES-CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
2.1 Introducción
7
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.2 MÉTODOS MODERNOS DE ANÁLISIS ESTRUCTURAL
ecuaciones diferenciales que gobiernan el fenómeno físico asociado. Por todo lo anterior
en el presente capítulo se describen brevemente los fundamentos, suposiciones y bases
teóricas de la mecánica estructural que el analista debe de conocer para llevar a cabo el
modelado de una estructura.
fs
ƒs
b
y
x
u*
donde
𝑏 fuerza de cuerpo, en N
𝑓𝑠 fuerza de superficie, en N
𝑡∗ tracciones prescritas, en MPa
𝑢∗ desplazamiento prescrito, en m
Γ frontera del sólido, en m2
Γ𝜎 frontera correspondiente a la región donde aplican cargas, en m2
𝛺 volumen de un sólido, en m³
2.2.1 Desplazamientos
Bajo la acción de las fuerzas aplicadas, un cuerpo sólido cambia, en general, de forma y
volumen, es decir, exhibe deformación hasta cierto punto. La deformación de este
cuerpo se describe matemáticamente de la siguiente manera. La posición de cualquier
punto en el cuerpo se define mediante su radio vector 𝐫 (con componentes 𝑥, 𝑦, 𝑧).
8
MANUAL DE DISEÑO DE OBRAS CIVILES-CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
Cuando el cuerpo se deforma, cada punto en él se desplaza a una nueva posición definida
por su radio vector 𝐫’ (con componentes 𝑥′ , 𝑦′ , 𝑧′ ). Debido a la deformación, el
desplazamiento está dado por el vector 𝐫′ − 𝐫, que en este manual se denota por 𝐮. en
este caso definido por sus tres componentes, 𝑢𝑥 , 𝑢𝑦 e 𝑢𝑧 , ec 2.1.
𝑢𝑥 = 𝑢𝑥 (𝑥, 𝑦, 𝑧)
𝑢𝑦 = 𝑢𝑦 (𝑥, 𝑦, 𝑧) (2.1)
𝑢𝑧 = 𝑢𝑧 (𝑥, 𝑦, 𝑧)
𝑢𝑥
𝐮 = [ 𝑦 ] = [𝑢 𝑥
𝑢 𝑢𝑦 𝑢𝑧 ]𝑇 = [𝑢 𝑣 𝑤]𝑇 (2.2)
𝑢𝑧
donde
𝐫, 𝐫′ radio vector con componentes 𝑥, 𝑦, 𝑧 y 𝑥′, 𝑦′, 𝑧′
u vector de desplazamientos
𝑢, 𝑣, 𝑤, 𝑢𝑥 , 𝑢𝑦 , 𝑢𝑧 coeficiente del vector de desplazamiento en coordenadas 𝑥,
𝑦, 𝑧, en m
Cuando un cuerpo cambia su forma y volumen debido a acciones externas, las distancias
entre dos puntos en él, cambian antes y después de la deformación. Este cambio se
mide mediante el estado de deformaciones representado por un tensor simétrico de
segundo orden, 𝛆 , que dentro del contexto de análisis estructural mediante el método
de los elementos finitos, es común representarlo en forma compacta por la matriz
simétrica de orden 3x3 de la ec 2.3,
𝜀𝑥 𝜀𝑥𝑦 𝜀𝑥𝑧
𝜀
𝛆 = [ 𝑦𝑥 𝜀𝑦 𝜀𝑦𝑧 ] (2.3)
𝜀𝑧𝑥 𝜀𝑧𝑦 𝜀𝑧
donde
𝛆 tensor simétrico de deformaciones de segundo orden
𝜺 vector de deformaciones
𝜀x , 𝜀y , 𝜀z coeficientes del tensor simétrico de deformaciones en su
diagonal principal, adimensional
𝜀𝑥𝑦 , 𝜀𝑥𝑧 , 𝜀𝑦𝑧 coeficientes del tensor de deformaciones en dirección angular,
fuera de la diagonal principal, adimensional
9
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.2 MÉTODOS MODERNOS DE ANÁLISIS ESTRUCTURAL
𝜕𝑢 𝜕𝑣 𝜕𝑤
𝜀𝑥 = ; 𝜀𝑦 = ; 𝜀𝑧 = (2.5)
𝜕𝑥 𝜕𝑦 𝜕𝑧
1 𝜕𝑢 𝜕𝑣 1 𝜕𝑢 𝜕𝑤 1 𝜕𝑣 𝜕𝑤
𝜀𝑥𝑦 = ( + ) ; 𝜀𝑥𝑧 = ( + ); 𝜀𝑦𝑧 = ( + ) (2.6)
2 𝜕𝑦 𝜕𝑥 2 𝜕𝑧 𝜕𝑥 2 𝜕𝑧 𝜕𝑦
Utilizando una notación matricial, las seis componentes del “vector de deformaciones”
se definen mediante la ec 2.7.
𝛆 = 𝐋𝐮 (2.7)
donde
𝐋 matriz que contiene operadores diferenciales, adimensional
𝜕
0 0
𝜕𝑥
𝜕
0 0
𝜕𝑦
𝜕
0 0
𝜕𝑧
𝐋= 1 𝜕 1 𝜕 (2.8)
0
2 𝜕𝑧 2 𝜕𝑦
1 𝜕 1 𝜕
0
2 𝜕𝑧 2 𝜕𝑥
1 𝜕 1 𝜕
0
[ 2 𝜕𝑦 2 𝜕𝑥 ]
lugar a fuerzas sobre el plano de corte que divide al cuerpo, mismas que tienden a
devolver su equilibrio. Estas fuerzas internas que se producen cuando un cuerpo se
deforma, expresadas como un vector fuerzas por unidad de área, se denominan
esfuerzos o vector de esfuerzos asociados a ese plano de corte. Es importante enfatizar
que en condiciones de equilibrio térmico, si no existe deformación, no hay esfuerzos.
𝜎𝑥 𝜎𝑥𝑦 𝜎𝑥𝑧
𝛔 = [𝜎𝑦𝑥 𝜎𝑦 𝜎𝑦𝑥 ] (2.9)
𝜎𝑧𝑥 𝜎𝑧𝑦 𝜎𝑧
donde
𝛔 tensor de esfuerzos de orden 2
𝜎𝑥 , 𝜎𝑦 , 𝜎𝑧 componentes del tensor de esfuerzos en la diagonal
principal, en MPa
𝜎𝑦𝑧 , 𝜎𝑧𝑦 , 𝜎𝑧𝑥 , 𝜎𝑥𝑧 , 𝜎𝑥𝑦 , 𝜎𝑦𝑥 componentes del tensor esfuerzos fuera de la
diagonal principal, en MPa
z
zx zy
x
xy yz
yx
y xz y
xz
yx
yz xy
x
zy zx y
z
11
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.2 MÉTODOS MODERNOS DE ANÁLISIS ESTRUCTURAL
𝜎𝑦𝑧 = 𝜎𝑧𝑦
𝜎𝑧𝑥 = 𝜎𝑥𝑧 (2.10)
𝜎𝑥𝑦 = 𝜎𝑦𝑥
donde
𝝈 vector de esfuerzos
donde
𝐼𝐼 invariante del tensor de esfuerzo, en MPa
𝐼𝐼𝐼 invariante del tensor de esfuerzo, en (MPa)²
𝐼𝐼𝐼𝐼 invariante del tensor esfuerzo, en (MPa)³
λ constante de Lamé, en MPa
𝚻 tensor de deformaciones o de esfuerzos principales
𝚰 vector de invariantes
donde
𝐼𝐼 = 𝜎𝑥 + 𝜎𝑦 + 𝜎𝑧
2
𝐼𝐼𝐼 = 𝜎𝑥 𝜎𝑦 − 𝜎𝑥𝑦 2
+ 𝜎𝑦 𝜎𝑧 − 𝜎𝑦𝑧 2
+ 𝜎𝑥 𝜎𝑧 − 𝜎𝑥𝑧 (2.13)
2 2 2
𝐼𝐼𝐼𝐼 = 𝜎𝑥 𝜎𝑦 𝜎𝑧 − 𝜎𝑥 𝜎𝑦𝑧 − 𝜎𝑦 𝜎𝑥𝑧 − 𝜎𝑧 𝜎𝑥𝑦 + 2𝜎𝑥𝑦 𝜎𝑥𝑧 𝜎𝑦𝑧
donde
𝜎𝐼 , 𝜎𝐼𝐼 , 𝜎𝐼𝐼𝐼 esfuerzos principales en dirección 𝐼, 𝐼𝐼 𝑦 𝐼𝐼𝐼, en MPa
𝜀𝐼 , 𝜀𝐼𝐼 , 𝜀𝐼𝐼𝐼 deformaciones principales en dirección 𝐼, 𝐼𝐼 𝑦 𝐼𝐼𝐼, adimensional
P
P
o y
donde
P fuerza en N
donde
𝜌 densidad de masa, N/m3
𝑏𝑥 , 𝑏𝑦 , 𝑏𝑧 fuerzas de cuerpo en dirección 𝑥, 𝑦, 𝑧 , en N
ü𝑥 ü𝑦 , ü𝑧 aceleración de la gravedad en dirección 𝑥, 𝑦, 𝑧, en m/s2
13
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.2 MÉTODOS MODERNOS DE ANÁLISIS ESTRUCTURAL
𝜎𝑥 l + 𝜎𝑥𝑦 m + 𝜎𝑥𝑧 n = 𝑡𝑥
𝜎𝑧𝑥 l + 𝜎𝑧𝑦 m + 𝜎𝑧 n = 𝑡𝑧
donde
l, m, n cosenos directores, adimensional, rad y/o grados
𝑡𝑐 tracciones, en MPa
𝑡𝑥 , 𝑡𝑦 , 𝑡𝑧 tracciones en dirección 𝑥, 𝑦, 𝑧, en MPa
Se considera que los sólidos elásticos lineales están hechos de una clase de material
muy específico con un comportamiento basado en las siguientes suposiciones:
14
MANUAL DE DISEÑO DE OBRAS CIVILES-CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
o matricialmente
𝝈 = 𝐂𝜺 (2.18)
donde
𝐂 matriz constitutiva simétrica del material
𝑐11 , 𝑐12 , 𝑐13 , 𝑐14 … 𝑐66 coeficientes de la matriz constitutiva simétrica del
material, en MPa
Los elementos de esta matriz son diferentes para cada tipo de comportamiento mecánico
del material. Para un material elástico lineal e isótropo, la matriz constitutiva se obtiene
a partir de la ley de Hooke generalizada, ec. 2.19.
(1 − 𝜇) 𝜇 𝜇 0 0 0
𝜇 (1 − 𝜇) 𝜇 0 0 0
𝜇 𝜇 (1 − 𝜇) 0 0 0
𝐸 (1 − 2𝜇)
𝐂= 0 0 0 0 0 (2.19)
(1 + 𝜇)(1 − 2𝜇) 2
(1 − 2𝜇)
0 0 0 0 0
2
(1 − 2𝜇)
[ 0 0 0 0 0 ]
2
donde
𝐸 módulo de elasticidad o de Young, en MPa
𝜇 relación de Poisson, adimensional
𝐸
𝐺= (2.20)
2(1 + 𝜇)
donde
𝐺 módulo de cortante, en MPa
15
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.2 MÉTODOS MODERNOS DE ANÁLISIS ESTRUCTURAL
𝜇𝐸
λ= (2.21)
(1 + 𝜇)(1 − 2𝜇)
𝐸
𝐾∗ = (2.22)
3(1 − 2𝜇)
donde
𝐾 ∗ módulo volumétrico, en MPa
Para problemas en los que existen deformaciones iniciales, la relación constitutiva para
un sólido elástico está dada por la ec 2.23,
𝝈 = 𝐂(𝜺 − 𝜺𝟎 ) (2.23)
donde
𝜺𝟎 vector de deformaciones iniciales
1
1
1
𝜺𝟎 = 𝛼Δ𝛵 (2.24)
0
0
{ 0 }
donde
𝛼 coeficiente de expansión térmica del material, adimensional
Δ𝛵 incremento de temperatura, en °C y/o K
16
MANUAL DE DISEÑO DE OBRAS CIVILES-CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
Matemáticamente esta ley establece que el trabajo realizado por un sistema mecánico
por las fuerzas externas actuando en sólido, más el flujo de calor absorbido o liberado
por el medio es igual al incremento de energía cinética más el incremento de energía de
deformación.
𝑊 + 𝑄 = ∆𝑇 ∗ + ∆𝑈𝑜∗ (2.25)
donde
𝑊 trabajo realizado en un sistema mecánico por una fuerza, en N-m
𝑄 flujo de calor absorbido del medio, en J/s
∆𝑇 ∗
incremento de energía cinética, N-m/m
∆𝑈𝑜∗ incremento de la energía de deformación interna, N-m/m
1 ∗ 2
∆𝑇 ∗ = (𝑈 ) (2.26)
2 𝑜
donde
𝑈𝑜∗ energía de deformación, en N-m/m
La consecuencia directa de esta ley es que en un cuerpo aislado la energía no puede ser
creada ni destruida, solo puede ser transformada de un estado a otro.
𝑊 ∗ = 𝑈∗ (2.27)
17
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.2 MÉTODOS MODERNOS DE ANÁLISIS ESTRUCTURAL
donde
𝑈∗ trabajo interno, en N-m
𝑊 ∗
trabajo externo, en N-m
En las ecs 2.25 y 2.27 “Trabajo externo = Trabajo interno” lo que implica que el
trabajo realizado por las fuerzas aplicadas sobre un cuerpo más la energía asociada al
flujo de calor, se transforma en la suma de energías cinética y de deformación.
donde
Δ𝑊 ∗ incremento de trabajo externo, N-m
𝐏 vector de fuerzas
𝐏 = [𝑝𝑥 𝑝𝑦 𝑝𝑧 ]𝑇 (2.29)
𝐮 = [𝑢𝑥 𝑢𝑦 𝑢𝑧 ]𝑇 (2.30)
donde
𝑝𝑥 , 𝑝𝑦 , 𝑝𝑧 coeficientes del vector de fuerzas, en N
du
dP
Carga P, en N
Pdu
dPu
W c = ud P
W =Pdu
(trabajo complementario)
(trabajo)
Desplazamiento u, en m Desplazamiento u, en m
(a) (b)
donde
𝑊𝑐 trabajo complementario, en N-m
donde
Δ𝑊𝑐 incremento de trabajo complementario, en N-m
1
𝑊 = 𝑊𝑐 = 𝐏𝐮 (2.32)
2
Wc
W = Wc
Sistema lineal
Carga P, en N
Desplazamiento u, en m
19
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.2 MÉTODOS MODERNOS DE ANÁLISIS ESTRUCTURAL
𝜀
𝑈0 = ∫ 𝛔 𝑑𝜀 (2.33)
0
donde
𝑈0 densidad de energía de deformación, N-m/m3
𝜀 deformación unitaria, adimensional
d
Esfuerzo , en MPa
Esfuerzo , en MPa
du d
P
Uc = d
Uo = d
Densidad de la energía de
Densidad de la energía
deformación complementaria
de deformación
𝑈 = ∫ 𝑈0 𝑑𝛺 (2.34)
donde
𝑈 energía de deformación total, en N-m/m
𝜎
𝑈𝑐 = ∫ 𝜺 𝑑𝝈 (2.35)
0
20
MANUAL DE DISEÑO DE OBRAS CIVILES-CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
𝑈𝑐∗ = ∫ 𝑈𝑐 𝑑𝛺 (2.36)
donde
𝑈𝑐∗ densidad de energía de deformación complementaria total, en N-m/m³
Para un cuerpo con un comportamiento elástico lineal, Fig 2.7, la densidad de energía
de deformación y la densidad de energía complementaria de deformación son
equivalentes, ec 2.37.
1 T
𝑈𝑐 = 𝑈𝑜 = 𝝈 𝜺 (2.37)
2
Uc
Esfuerzo , en MPa
Uc = Uo
Ley constitutiva
Uo lineal del material
Deformación , adimensional
1 𝟏
𝐏𝐮 = ∫ 𝝈T 𝜺𝑑𝛺 (2.38)
2 𝟐 𝛺
Bajo este contexto, este principio definido en términos del trabajo real efectuado por un
sistema de fuerzas actuando en un cuerpo podría utilizarse para determinar el
desplazamiento en dirección de la fuerza.
21
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.2 MÉTODOS MODERNOS DE ANÁLISIS ESTRUCTURAL
el sistema estructural, llevando en la generalidad de los casos a tener una ecuación con
más de una incógnita.
La definición de cada uno de estos principios, se basa en conceptos que los hacen ser
principios duales. El primero se relaciona con el trabajo debido a una fuerza real
moviéndose en la dirección de un desplazamiento virtual y el segundo se relaciona con
el trabajo debido a una fuerza virtual moviéndose en la dirección de un desplazamiento
real. Gráficamente ambos principios están correspondientemente representados por el
área bajo o sobre la curva en las Figs 2.8a y 2.8b.
Estos dos principios duales pueden ser planteados para una partícula, para un cuerpo
rígido, cuerpo no deformable y para un cuerpo deformable. Por la importancia que tienen
estos principios en los métodos de análisis estructural en este manual solo se presentan
los principios aplicables a cuerpos deformables con un comportamiento constitutivo
general.
22
MANUAL DE DISEÑO DE OBRAS CIVILES-CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
u Pu = W c
P
Carga P, en N
Carga P, en N
P u = W
Desplazamiento u, en m Desplazamiento u, en m
(a) (b)
donde
𝛿𝐮 vector de desplazamiento virtual
𝛿𝐏 vector de carga virtual
𝛿𝑊 trabajo debido a un desplazamiento virtual, en N-m
𝛿𝑊𝑐 trabajo complementario debido a una carga virtual, en N-m
t*
b
u
u*
donde
Γ𝑢 frontera correspondiente a la región con desplazamientos prescritos
23
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.2 MÉTODOS MODERNOS DE ANÁLISIS ESTRUCTURAL
Para este cuerpo sólido, el principio de los desplazamientos virtuales se establece a partir
de sus ecuaciones de equilibrio interno y en su frontera. Los desplazamientos virtuales
impuestos, aun cuando no están relacionados con las fuerzas reales impuestas, deben
de satisfacer los requisitos de continuidad dentro del cuerpo, además de ser igual a cero
sobre la frontera donde se imponen condiciones cinemáticas.
𝛿𝑢 = 0, 𝛿𝑣 = 0, 𝛿𝑤 = 0 𝑒𝑛 Γ𝑢
𝛿𝑈 = 𝛿𝑊 (2.40)
donde
𝛿𝑈 trabajo virtual interno, en N-m
𝛿𝑊 trabajo debido a un desplazamiento virtual, en N-m
𝐭𝒓 vector de tracciones
𝐛 vector de fuerzas de cuerpo
𝛿𝑢, 𝛿𝑣, 𝛿𝑤 desplazamientos virtuales
𝛿𝐮 vector de desplazamientos virtuales
Este principio establece que un sólido continuo bajo la acción de un sistema de fuerzas
está en equilibrio, si para, desplazamientos virtuales compatibles con las condiciones de
frontera impuestas, el trabajo virtual externo es igual al trabajo virtual interno.
24
MANUAL DE DISEÑO DE OBRAS CIVILES-CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
Para ilustrar este principio básico en el análisis estructural considérese un cuerpo sólido
en equilibrio como se ilustra en la Figs 2.10a y b, al que se aplican dos sistemas de
fuerzas P𝐼 y P𝐼𝐼 , que producen correspondientemente dos configuraciones deformadas, u𝐼
y u𝐼𝐼 . El principio del trabajo recíproco o de Maxwell-Betti establece que el trabajo
producido por el sistema P𝐼 sobre los desplazamientos u𝐼𝐼 .es igual al trabajo ejercido por
el sistema P𝐼𝐼 sobre los desplazamientos u𝐼 , lo anterior se expresa como,
donde
𝑊𝐼/𝐼𝐼 , 𝑊𝐼𝐼/𝐼 trabajo recíproco de Maxwell-Betty, en N-m
P𝐼 , P𝐼𝐼 fuerza 𝐼 y 𝐼𝐼, en N
u𝐼 desplazamiento debido a la fuerza P𝐼𝐼 , en m
u𝐼𝐼 desplazamiento debido a la fuerza P𝐼 , en m
25
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.2 MÉTODOS MODERNOS DE ANÁLISIS ESTRUCTURAL
PII
PI
PII
PI
(a) (b)
El caso particular del principio de los desplazamientos virtuales aplicado a sólidos con un
comportamiento elástico, lineal o no—lineal, éste da lugar al principio de la energía
potencial mínima. Para establecer este principio variacional se acepta, para un cuerpo
elástico deformado en condiciones adiabáticas, la existencia de una funcional positiva
definida, denominada densidad de energía de deformación 𝑈0 , función de los elementos
del tensor deformación, tal que,
𝜕𝑈0
𝝈= (2.43)
𝜕𝜺
𝜕𝑈0
∫ 𝐛𝛿𝐮𝑑𝛺 + ∫ 𝐭 𝒓 𝛿𝐮𝑑Γ = ∫ 𝛿𝛆𝑑𝛺 (2.44)
𝛺 Γ1 𝛺 𝜕𝛆
En la que
𝜕𝑈0
∫ 𝛿𝛆𝑑𝛺 = ∫ 𝛔𝛿𝛆𝑑𝛺 = 𝛿 ∫ 𝑈0 𝑑𝛺 = 𝛿𝑈 (2.46)
𝛺 𝜕𝛆 𝛺 𝛺
26
MANUAL DE DISEÑO DE OBRAS CIVILES-CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
donde
Π energía potencial del cuerpo, N-m/m
𝜕𝑈𝑜∗
𝜺= (2.48)
𝜕𝝈
Con base en la ec 2.48 y en las ecuaciones con las que se formula el principio de las
fuerzas virtuales se encuentra que
𝛿(𝑈 ∗ + 𝑊 ∗ ) = 𝛿Π ∗ = 0 (2.49)
donde
Π∗ energía potencial complementaria, en N-m/m
La ec 2.49 establece que: De todos los estados de esfuerzo admisibles que satisfacen
las ecuaciones de equilibrio y las tracciones prescritas en la frontera, aquel que hace
27
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.2 MÉTODOS MODERNOS DE ANÁLISIS ESTRUCTURAL
1. Las fuerzas de cuerpo y las tracciones prescritas deben cumplir las ecuaciones de
equilibrio.
2. Los esfuerzos relacionados con el campo de deformaciones a través de una ley
constitutiva deben minimizar la energía potencial complementaria con respecto a
sus correspondientes campos.
donde
Π𝐻𝑅 funcional de energía, N-m/m
𝐭𝒓 vector de tracciones
En este principio se establece que los desplazamientos y esfuerzos que producen el valor
estacionario de la funcional dada por la ec 2.49, satisfacen las ecuaciones de
compatibilidad y de equilibrio. Es posible demostrar que el valor extremo del funcional
se obtiene de,
𝛿Π𝐻𝑅 = 0 (2.51)
𝑡2 𝑡2 𝑡2
− ∫ ∫ 𝛿𝑢𝑇 𝜌ü 𝑑𝛺 𝑑𝑡 = ∫ 𝛿𝑈 𝑑𝑡 + ∫ 𝛿𝑊𝑐 𝑑𝑡 (2.52)
𝑡1 𝛺 𝑡1 𝑡1
o en su forma compacta
28
MANUAL DE DISEÑO DE OBRAS CIVILES-CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
𝑡2
𝛿 (∫ (Δ𝑇 ∗ − 𝑈 ∗ − 𝑊 ∗ )𝑑𝑡 ) = 0 (2.53)
𝑡1
𝐿∗ = Δ𝑇 ∗ − 𝑈 ∗ − 𝑊 ∗ (2.54)
donde
ü aceleración de la gravedad, cm/s²
𝐿 ∗
función Lagrangiana, N-m/m
𝑡1 , 𝑡2 tiempo 1 y tiempo 2, en s
La ec 2.5.2 define el principio de Hamilton que establece que las posibles configuraciones
que una estructura puede adoptar a lo largo de un intervalo de tiempo [𝑡1 , 𝑡2 ], aquella
que satisface el equilibrio dinámico es la que hace estacionaria la integral de 𝐿∗ durante
dicho intervalo de tiempo.
Se debe entender por elemento finito como una subregión de un medio continuo dividido
en varias partes, elementos. Este elemento es de tamaño finito y normalmente su
geometría es simple, líneas, triángulos, rectángulos, tetraedros, cubos; en comparación
con la geometría del problema completo. Desde el punto de vista algebraico, el método
transforma un problema continuo, en uno discreto. Esto se logra debido a que el MEF
transforma los sistemas de ecuaciones diferenciales en sistemas de ecuaciones
algebraicas.
29
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.2 MÉTODOS MODERNOS DE ANÁLISIS ESTRUCTURAL
(a) (b)
(c) (d)
(c) (d)
Fig 2.11 Estructuras reticulares
30
MANUAL DE DISEÑO DE OBRAS CIVILES-CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
(a) (b)
Para un problema general de la Mecánica de Sólidos, que está caracterizado por tener
un modelo en tres dimensiones, se divide el medio continuo en elementos de volumen
simple, Fig 2.13.
(a) (b)
31
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.2 MÉTODOS MODERNOS DE ANÁLISIS ESTRUCTURAL
bajo la acción de cargas externas aplicadas de forma normal a su plano, Fig 2.14, teoría
de flexión en placas.
(a) (b)
Una de las grandes ventajas de utilizar el MEF, es que en un mismo problema es posible
utilizar una malla que considere la interacción entre diferentes tipos de elementos finitos;
aplicables al comportamiento estructural particular de una región del problema. Esto es
posible considerando la compatibilidad de los grados de libertad que convergen a un
mismo nodo, Fig 2.16.
32
MANUAL DE DISEÑO DE OBRAS CIVILES-CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
Pre—proceso
33
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.2 MÉTODOS MODERNOS DE ANÁLISIS ESTRUCTURAL
Solución
Post—proceso
En este proceso de mallado del problema, se crean nodos elementos y nodos en los
cuales convergen uno o más elementos y que definen las conectividades para cada
elemento. El número de elementos de malla está en función de la aproximación deseada;
a mayor número de elementos, mayor exactitud en la solución.
34
MANUAL DE DISEÑO DE OBRAS CIVILES-CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
Dependiente del comportamiento del elemento finito, es posible tener diferentes tipos
de grados de libertad para cada problema. Para elementos sólidos, las incógnitas son los
desplazamientos en las direcciones 𝑥, 𝑦, 𝑧. Para problemas en los que los elementos
finitos involucran flexión, se puede tener hasta tres grados de libertad más asociados al
giro alrededor los ejes cartesianos.
Para obtener los desplazamientos, dentro de cada elemento finito, se utilizan polinomios
simples que satisfagan las condiciones de frontera del problema,
𝑛𝑑
donde
𝐝′ vector local de desplazamientos
𝒅𝒏 vector nodal de desplazamientos
𝐍 matriz de funciones de forma
𝑛𝑑 numero de nodos de un elemento
𝑇
𝒅𝒏 = [𝑑𝑛𝑥 , 𝑑𝑛𝑦 , 𝑑𝑛𝑧 … 𝑑𝑛𝑓 ] (2.56)
donde
𝑑𝑛𝑥 , 𝑑𝑛𝑦 , 𝑑𝑛𝑧 … 𝑑𝑛𝑓 coeficientes del vector nodal de desplazamientos, en m
𝑢
𝐝′ = 𝐮 = [ 𝑣 ] (2.57)
𝑤
𝑑1′
𝑑′
𝐝′ = [ 2 ] (2.58)
⋮
′
𝑑𝑛𝑑
35
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.2 MÉTODOS MODERNOS DE ANÁLISIS ESTRUCTURAL
donde
′
𝑑1′ , 𝑑2′ … 𝑑𝑛𝑑 coeficientes del vector local de desplazamientos de un elemento
𝑛𝑓 número de grados de libertad de un nodo
𝑛𝑑 número de nodos de un elemento
𝑁1 0 0
𝐍=[0 𝑁2 0] (2.60)
0 0 𝑁𝑛
donde
𝑁1 , 𝑁2 . . . . 𝑁𝑛 coeficientes de la matriz que contiene las funciones de forma,
adimensional
Una vez definidas las funciones de forma, se procede a determinar la matriz de rigidez
elemental por medio del principio de la energía potencial, ec 2.61.
1 1 1
Π = ∫ 𝜺𝑇 𝐂𝜺 𝑑𝛺 = ∫ 𝐝𝑇 𝐁 T 𝐂𝐁𝑑𝛺 = 𝐝𝑇 ∫ 𝐁 T 𝐂𝐁𝑑Ω (2.61)
2 𝛺 2 𝛺 2 𝛺
donde
𝐁 matriz de transformación
𝐁 = 𝐋𝐍 (2.62)
donde
𝐋 matriz que contiene los operadores diferenciales
36
MANUAL DE DISEÑO DE OBRAS CIVILES-CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
𝐤 = ∫ 𝐁 𝑇 𝐂𝐁 𝑑𝛺 (2.63)
𝛺
donde
k matriz local de rigidez
𝑇 𝑇 𝑇
𝑊𝑃 = 𝐝′ 𝐟𝑏 + 𝐝′ 𝐟𝑠 = 𝐝′ 𝐟𝑒 (2.64)
𝐟𝑒 = 𝐛 + 𝑭𝑠 (2.65)
donde
𝑊𝑃 trabajo realizado por las cargas externas, en N-m
𝐟𝑏 vector local de fuerzas nodales debidos a las fuerzas de cuerpo
𝐟𝑠 vector local de fuerzas nodales debido a las fuerzas de superficie
𝐟𝑒 suma de vectores de fuerzas de cuerpo y fuerzas de superficie
Todas las ecuaciones definidas anteriormente, están formuladas para cada elemento en
coordenadas locales; por lo que es necesario realizar las transformaciones de
coordenadas necesarias de acuerdo al tipo de elemento finito a utilizar, ver Fig 2.18.
37
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.2 MÉTODOS MODERNOS DE ANÁLISIS ESTRUCTURAL
Después de obtener las ecuaciones para todos los elementos que conforman la malla, y
su transformación a coordenadas globales, se ensamblan en el sistema global para
formar las ecuaciones del problema completo. Al realizar el ensamble se está
reestableciendo la continuidad del problema,
𝐊𝐝 = 𝐅 (2.66)
K K12
𝐅 = [ 11 ]𝐝 (2.67)
K 21 K 22
donde
𝐅 vector global de fuerzas
K matriz global de rigidez
K11 , K12 , K 21 , K 22 submatrices de la matriz global de rigidez
Existen varias formas para imponer los apoyos; sin embargo, la más común es modificar
la matriz de rigidez y el vector de cargas, removiendo los renglones y las columnas
correspondientes a los grados de libertad en los que se tiene un apoyo. Después de este
procedimiento, es posible obtener la inversa de la matriz de rigidez para determinar el
vector de desplazamientos nodales totales del problema,
𝐝 = 𝐊 −𝟏 𝐅 (2.68)
En este apartado se establecen los fundamentos básicos del método directo de las
rigideces, para construir las ecuaciones de equilibrio global de una estructura; en el que
los desplazamientos en los nodos
38
MANUAL DE DISEÑO DE OBRAS CIVILES-CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
desplazamientos en todos y cada uno de los puntos nodales del modelo; satisfaciendo la
compatibilidad cinemática de desplazamiento en los apoyos de la estructura.
un6
un5
F
un8 un2
B
un7 C
un1
E
G A
D
un3
un9 un4
un10
donde
un1 , un2 , un3 … . un10 desplazamientos en los nodos, en m
Se establece que las matrices de rigidez de cada elemento (A – G, ver Fig 2.19) y sus
fuerzas de sistema algebraico global de la estructura. Bajo este punto, partamos del
hecho de que se dispone de tal sistema de ecuaciones algebraicas dado por,
F3 u3
𝑘 𝑘32 𝑘33 𝑘34 … 𝑘3𝑛
= 31 u4 (2.69)
F4 𝑘41 𝑘42 𝑘43 𝑘44 … 𝑘4𝑛
⋮
⋮ ⋮ ⋮ ⋮ ⋮ ⋱ ⋮ u
{F𝑛 } [𝑘𝑛1 𝑘𝑛2 𝑘𝑛3 𝑘𝑛3 … 𝑘𝑛𝑛 ] { 𝑛 }
donde
F1 , F2 , F3 , F4 … F𝑛 coeficientes del vector global de fuerzas, en N
u1 , u2 , u3 , u4 … u𝑛 coeficientes del vector global de desplazamientos, m
𝑘11 , 𝑘12 , 𝑘13 , 𝑘14 … 𝑘𝑛𝑛 coeficientes de la matriz global de rigidez, en N/m
Los grados de libertad de la estructura están identificados por los subíndices 𝑛𝑛, que
corresponden a una numeración global. Se entiende que el sistema de numeración global
es independiente de la numeración local adoptada para cada elemento en particular. La
39
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.2 MÉTODOS MODERNOS DE ANÁLISIS ESTRUCTURAL
Para mostrar la aplicación de este método, adicionalmente se considera que los sistemas
de ecuaciones algebraicas han sido resueltas numéricamente y que se dispone de los
desplazamientos, fuerzas internas, y cargas asociadas a cada grado de libertad; por lo
que es posible evaluar estas mismas cantidades para todos y cada uno de los elementos
que conforman la estructura.
El proceso de aplicación del método, se aplica a todos y cada uno de los grados de
libertad que comprende la estructura; sin embargo, con el fin de ejemplificar,
únicamente consideremos como ejemplo el nodo 𝑖, para el grado de libertad horizontal,
en el que las cargas aplicadas y las fuerzas internas, de los elementos que concurren en
el nodo, deben de satisfacer el equilibrio, y donde la compatibilidad de desplazamientos
se ha satisfecho a priori en el proceso de ensamble del sistema algebraico global de la
estructura, Fig 2.20. Las cantidades de interés asociadas al nodo y grado de libertad
considerado, se identifican mediante el subíndice 𝑖.
FC FB
P
F A , F B, F C , F D
i
FD FA
donde
𝐹 𝐴 , 𝐹𝐵 , 𝐹𝐶 , 𝐹𝐷 fuerzas internas de cada uno de los elementos que concurren en el
nodo 𝑖, en N
P fuerza, en N
40
MANUAL DE DISEÑO DE OBRAS CIVILES-CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
Para que la estructura esté en equilibrio, se debe de cumplir que parte de ella también
lo esté. A nivel nodal esto se puede mostrar al hacer la sumatoria de fuerzas en el nodo
y grado de libertad sea igual a cero; esto conduce a establecer la siguiente ecuación de
equilibrio,
𝑃 = 𝐹 𝐴 + 𝐹𝐵 + 𝐹𝐶 + 𝐹𝐷 (2.70)
Por la condición de compatibilidad los desplazamientos que concurren al nodo 𝑖 son los
mismos para todos los elementos, por lo que la podemos reescribir la ec 2.71 como
u𝑖 = u𝑖𝐴 = u𝐵𝑖 = u𝐶𝑖 = u𝐷
𝑖 (2.71)
donde
u𝑖𝐴 , u𝐵𝑖 , u𝐶𝑖 , u𝐷
𝑖 desplazamiento nodal de los elementos estructurales A, B, C y
D, en m
La ecuación que relaciona las fuerzas con los desplazamientos, se deriva a partir del
sistema algebraico dado por la ec 2.67. Por lo que el equilibrio en el grado de libertad
en consideración se reescribe como,
𝐴
𝐅 = K11 𝐵
u𝑖 + K11 𝐶
u𝑖 + K11 𝐷
u𝑖 + K11 u𝑖 (2.72)
donde
𝐴
K11 𝐵
, K11 𝐶
, K11 𝐷
, K11 submatrices de la matriz global de rigidez de los
elementos estructurales A, B, C, D.
u𝑖 desplazamiento en el nodo 𝑖, en m
41
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.2 MÉTODOS MODERNOS DE ANÁLISIS ESTRUCTURAL
j-ésimo
un grado de
u1 libertad un
F1 k11
i-ésima kij
Fuerza
knn
Fn
donde
u1 , u2 , u3 , u4 … u𝑛 coeficientes del vector global de desplazamientos, en m
F1 , F2 , F3 , F4 … F𝑛 coeficientes del vector global de fuerzas, en N
𝑘11 , 𝑘12 , 𝑘13 , 𝑘14 … 𝑘𝑛𝑛 coeficiente de la matriz global de rigidez, en N/m
42
MANUAL DE DISEÑO DE OBRAS CIVILES-CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
a los grados de libertad que están restringidos. Para esto, se puede dividir el sistema
como.
𝐏𝑓 𝐊𝑓𝑓 𝐊𝑓𝑠 𝐃𝑓
[ ]=[ ][ ] (2.74)
𝐏𝑠 𝐊 𝑠𝑓 𝐊 𝑠𝑠 𝐃𝑠
donde
𝐏𝑓 vector de fuerza aplicada en nodo no restringido
𝐏𝑠 vector de fuerza aplicada en nodo restringido
𝐃𝑓 vector de desplazamientos de nodos no restringidos
𝐊𝑓𝑓 sub-matriz de rigidez con condiciones de apoyo no restringidos
𝐊 𝑠𝑠 sub-matriz de rigideces con condiciones de apoyo restringido
𝐊𝑓𝑠 , 𝐊 𝑠𝑓 sub-matriz de rigideces con condiciones de apoyos restringido y no
restringidos
𝐃𝑠 vector de desplazamientos de nodos restringidos
−1
𝐃𝑓 = [𝐊𝑓𝑓 ] 𝐏𝑓 (2.75)
𝐏𝑠 = 𝐊 𝑠𝑓 𝐃𝑓 (2.76)
6. Finalmente, para obtener las fuerzas internas para cada uno de los elementos, se
obtiene el vector de desplazamientos asociado a cada elemento y se multiplica por
su respectiva matriz de rigidez. Las fuerzas así obtenidas, están dadas en
coordenadas globales; las cuales se utilizan para establecer el equilibrio en cada
uno de los nodos de la estructura. Las fuerzas internas para el elemento A
(Fig 2.20) se estiman como
𝐅𝐴 = 𝐊 𝐴 𝐝 𝐴 (2.77)
o en coordenadas locales
𝐟 𝐴 = 𝐑𝐅𝐴 (2.78)
donde
𝐑 matriz de rotación
𝐟 𝐴
vector local de fuerzas del elemento A
𝐅𝐴 vector global de fuerzas del elemento A
43
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.2 MÉTODOS MODERNOS DE ANÁLISIS ESTRUCTURAL
1. La matriz es simétrica
2. No existe un término sobre la diagonal principal igual o menor a cero, todos son
positivos. Esto hace que la matriz sea definida positiva.
Existen muchos términos fuera de la diagonal principal que son iguales a cero. Este
hecho se deriva en que, estrictamente hablando, el problema se está resolviendo de
manera discreta. A este tipo de matriz se les conoce como porosa.
El método de las flexibilidades deriva su nombre del hecho de que la relación fuerza—
desplazamiento tanto de la estructura como de los elementos están en términos de
flexibilidades. Empezando con la relación de flexibilidad entre las fuerzas y los
desplazamientos de los elementos estructurales, se utilizan las relaciones de equilibrio y
de compatibilidad del sistema para generar un sistema de 𝑛 ecuaciones en términos de
𝑛 fuerzas desconocidas. Una vez planteadas estas ecuaciones, es posible resolver para
las fuerzas desconocidas. Después, es posible utilizar esas fuerzas para encontrar los
desplazamientos restantes.
Aunque los pasos más importantes del método de las flexibilidades son análogos a los
del método de las rigideces, los detalles y las sucesiones de las operaciones son, de
alguna manera, más complejos. Como consecuencia de la complejidad del método, no
es fácil extender los pasos del método para generar un método general de flexibilidades,
similar al método de las rigideces.
44
MANUAL DE DISEÑO DE OBRAS CIVILES-CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
Paso 2 Utilizando las ecuaciones de equilibrio, determinar las fuerzas en todos los
miembros, debidas a las fuerzas primarias P𝑃 y redundantes P𝑅 .
M, V, N = f(P𝑃 , P𝑅 ) (2.79)
u = f(P, P𝑅 ) = u𝑃 , u𝑅 (2.80)
uI = u𝑃 + u𝑅 (2.81)
Paso 6 Una vez conocidas las fuerzas redundantes P𝑅 , se puede regresar al paso
2 para encontrar las fuerzas totales en todos los miembros de la
estructura.
donde
P𝑃 fuerzas primarias, en N
P𝑅 fuerzas redundantes, en N
M momento flexionante, en N-m
V cortante, en N
N fuerza normal, en N
P fuerza, en N
u desplazamiento, en m
u𝑃 desplazamiento debido a las fuerzas primarias, en m
u𝑅 desplazamiento debido a las fuerzas redundantes, en m
uI desplazamiento impuesto, en m
45
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.2 MÉTODOS MODERNOS DE ANÁLISIS ESTRUCTURAL
A pesar de ser el método de las flexibilidades una herramienta de análisis muy poderosa
para la solución de problemas en los que es fundamental cumplir con el equilibrio en
todos y cada uno de los puntos del continuo, en este manual se aplica el método de las
rigideces para formular el método de los elementos finitos. Sin embargo, es posible
obtener la matriz de flexibilidades a través de primero calcular la matriz de rigidez.
Py1 Py1
Nodo 1 Nodo 1
Px1 Px1
Px3 Nodo 3
Nodo 3
y, v Py3 Py2 y, v
Nodo 2
x, u x, u
(a) (b)
donde
P𝑥1 , P𝑥2 , P𝑥3 fuerza en el nodo 1, 2 y 3, dirección 𝑥, en N
P𝑦1 , P𝑦2 , P𝑦3 fuerza en el nodo 1, 2 y 3, dirección 𝑦, en N
u𝑥1 , u𝑦1 , u𝑥2 desplazamientos en los nodos 1 y 2, dirección 𝑥, 𝑦, en m
Fig 2.22 Estructura con tres elementos, a) sin apoyos y b) con apoyos
Las cantidades relacionadas con los apoyos se referencian con un subíndice 𝑠 y aquellas
debidas a los restantes grados de libertad no restringidos, se utiliza el subíndice 𝑓.
Entonces, la matriz de rigidez de un elemento se particiona como la ec 2.74.
𝑻
𝐏𝑓 = [𝑃𝑥1 𝑃𝑦1 𝑃𝑥2 ] (2.82)
46
MANUAL DE DISEÑO DE OBRAS CIVILES-CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
𝑻
𝐃𝑓 = [𝑢𝑥1 𝑢𝑦1 𝑢𝑥2 ] (2.83)
𝐏𝑓 𝐊𝑓𝑓
[ ]=[ ] [𝐃𝑓 ] (2.84)
𝐏𝑠 𝐊 𝑠𝑓
𝐃𝑓 = 𝐟𝐟 𝐏𝑓 (2.85)
−1
f𝑐 = [𝐊𝑓𝑓 ] (2.86)
donde
𝐟𝐟 matriz de flexibilidades
f𝑐 matriz de coeficientes de flexibilidad de un elemento o de la estructura
completa
donde
𝐏𝑠 = 𝐁𝐏𝑓 (2.89)
𝐏𝑠 = 𝐁𝐊𝑓𝑓 𝐃𝑓 (2.90)
47
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.2 MÉTODOS MODERNOS DE ANÁLISIS ESTRUCTURAL
o, de la ec 2.76,
𝐏𝑠 = 𝐊 𝑠𝑓 𝐃𝑓 (2.91)
en donde
Para establecer la sub—matriz 𝐊𝑓𝑠 , se utiliza el teorema del trabajo recíproco, que en
forma matricial se puede escribir como
𝑇
𝐊𝑓𝑠 = [𝐊 𝑠𝑓 ] = 𝐁[𝐟𝐟 ]−1 (2.93)
𝐊 𝑠𝑠 = 𝐁𝐊𝑓𝑠 (2.94)
La ec 2.95 es una fórmula general para transformar la matriz de flexibilidad a una matriz
de rigidez; que incluye movimientos de cuerpo rígido, y se puede aplicar a cualquier tipo
de estructura o elemento.
El MEE puede ser considerado como una combinación de diferentes métodos numéricos,
i.e., método del elemento finito, método de la matriz de rigidez dinámica y el método de
análisis espectral.
48
MANUAL DE DISEÑO DE OBRAS CIVILES-CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
donde
refinamiento h es una malla más fina que la original
refinamiento P es una malla elaborada utilizando con elementos cuadráticos
donde
𝐒𝑔 matriz espectral global
𝜔𝑛𝑎𝑡 frecuecia natural, en hertz
𝐝𝑔 vector espectral global de grados de libertad
𝐟𝑔 vector espectral global de fuerzas
Cuarto. Resolver la ec 2.97 para los grados de libertad global 𝐝𝑔 y las frecuencias
discretas. A partir de ello, determinar las componentes espectrales de los
desplazamientos.
49
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.2 MÉTODOS MODERNOS DE ANÁLISIS ESTRUCTURAL
donde
TRF transformada rápida de Fourier
ITRF inversa de la transformada rápida de Fourier
MAE método de análisis espectral
MDR método directo de rigideces
50
MANUAL DE DISEÑO DE OBRAS CIVILES-CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
𝐮ℎ = 𝐍𝐚 + ℋΓd 𝐍𝐛 (2.98)
donde
𝐮ℎ desplazamientos discontinuos, en m
𝐚 término que hace que el campo de desplazamientos sea continuo,
en m
ℋ𝛤𝑑 término que hace que el campo de desplazamientos sea discontinuo
en m
51
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.2 MÉTODOS MODERNOS DE ANÁLISIS ESTRUCTURAL
52
MANUAL DE DISEÑO DE OBRAS CIVILES-CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
u'
Líneas
de
+ desplazamiento
u
=
Hu
donde
𝑢′ desplazamiento campo discontinuo, en m
𝑢̆ desplazamiento campo continuo, en m
𝐻𝑢̆ suma de un desplazamiento en un campo discontinuo más un desplazamiento
en un campo continuo, en m
53
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.2 MÉTODOS MODERNOS DE ANÁLISIS ESTRUCTURAL
En los métodos libres de malla, MLM, la solución se obtiene con colocar únicamente
nodos, o partículas, dentro del dominio del problema, Fig 2.29 a) y 2.29 b); eliminando
las restricciones derivadas de tener una malla de elementos finitos.
54
MANUAL DE DISEÑO DE OBRAS CIVILES-CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
(a) (b)
Fig 2.29 a) Malla de elementos finitos y b) Discretización del dominio con MLM
55
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.2 MÉTODOS MODERNOS DE ANÁLISIS ESTRUCTURAL
a una dimensión, como se observa en la Fig 2.30a, debido a que solo es necesario
realizar un mallado de las fronteras del problema. Al resolver el mismo problema con el
método de los elementos finitos, Fig 2.30b, el mallado tiene que realizarse en todo el
dominio del problema, no solo en las fronteras; utilizando elementos finitos en dos
dimensiones.
Sólo la frontera es
dicretizada
( a )a (b) b
56
MANUAL DE DISEÑO DE OBRAS CIVILES-CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
Una comparativa del método estándar de elementos finitos se presenta en la Tabla 2.1.
57
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.2 MÉTODOS MODERNOS DE ANÁLISIS ESTRUCTURAL
58
MANUAL DE DISEÑO DE OBRAS CIVILES-CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
59
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
3.1 Introducción
61
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.3 ELEMENTOS ESTRUCTURALES PARA MODELAJE
Un resorte, es uno de los elementos estructurales más simples, solo resiste cargas
axiales de compresión o tensión; dando lugar a un alargamiento o acortamiento de su
longitud y que es directamente proporcional a la fuerza axial aplicada, en un rango lineal,
i.e., comportamiento lineal.
u1 u2
P Nodo 1 kr Nodo 2 P x
donde
k𝑟 rigidez de un elemento resorte, en N/m
P1 , P2 fuerza aplicada en el nodo 1 y 2, en N
u1 , u2 desplazamiento del nodo 1 y 2, en m
El resorte mostrado en la Fig 3.1 consiste de dos nodos localizados en los extremos de
su longitud. Los desplazamientos, u1 y u2 , se obtienen en estos puntos nodales. Es posible
que en algunos puntos nodales de una estructura completa existan cargas aplicadas
referidas como P
62
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
donde
δ alargamiento, en m
N = k 𝑟 δ = k 𝑟 (u2 − u1 ) (3.2)
donde
N fuerza normal, en N
En la Fig 3.2 se muestra la relación elástica lineal entre la deformación axial del resorte
y la fuerza actuante.
Fuerza P1,2, en N
kr
Alargamiento , en m
donde
P fuerza, en N
Del equilibrio del elemento, se obtiene la siguiente relación matricial entre los
desplazamientos y las fuerzas nodales,
𝑘𝑟 −𝑘𝑟 𝑢1 𝑃
[ ] { } = { 1} (3.3)
−𝑘𝑟 𝑘𝑟 𝑢 2 𝑃2
donde
𝑢1 , 𝑢2 coeficiente del vector local de desplazamientos de un elemento resorte
y armadura, en m
𝑃1 , 𝑃2 coeficiente del vector local de fuerzas de un elemento resorte, en N
𝑘𝑟 coeficiente de la matriz local de un elemento resorte, en N/m
63
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.3 ELEMENTOS ESTRUCTURALES PARA MODELAJE
O de forma matricial
𝐤𝐝′ = 𝐟 (3.4)
𝑘 −𝑘𝑟
𝐤=[ 𝑟 ] (3.5)
−𝑘𝑟 𝑘𝑟
donde
𝐤 matriz local de rigidez
𝐟 vector local de fuerzas
𝐝′ vector local de desplazamientos
El elemento armadura, geométricamente, se caracteriza por ser una barra de eje recto
en la que solo se aplican cargas axiales. De acuerdo a esto, el elemento solo tiene
deformación axial, acortamiento o alargamiento.
A
Nodo 1
P
l
l
Nodo 2
donde
δ alargamiento, en m
P fuerza, en N
𝐴 área de la sección transversal del elemento estructural, en m2
𝑙 longitud de un elemento estructural, barra, viga, marco y cable, en m
Para dar la formulación del elemento finito tipo armadura, considere un elemento
perteneciente a un sistema estructural, como el que se muestra en la Fig 3.3, con sección
transversal constante 𝐴. El elemento se define por los dos nodos extremos 1 y 2. La
longitud del elemento es 𝑙. El sistema local de referencias coincide con el eje axial del
elemento, en el origen en el nodo 1. En este sistema de referencia, únicamente se tiene
un grado de libertad por nodo, desplazamiento axial, y el número total de grados de
libertad del elemento es 2, (𝑛𝑓 = 2).
𝐮 = 𝐍𝐝´ (3.6)
donde
𝐮 vector de desplazamientos
𝐍 matriz que contiene las funciones de forma
donde
𝐍 = [𝑁1 𝑁2 ] (3.7)
𝑢1
𝐝´ = {𝑢 } (3.8)
2
𝑥 𝑥
𝑁1 = 1 − , 𝑁2 = (3.9)
𝑙 𝑙
donde
𝑁1 , 𝑁2 coeficientes de la matriz que contiene las funciones de forma,
adimensional
𝑥 coordenada cartesiana en dirección x, en m
𝜕𝑢 u2 − u1
ε𝑥 = = (3.10)
𝜕𝑥 𝑙
65
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.3 ELEMENTOS ESTRUCTURALES PARA MODELAJE
donde
ε𝑥 deformación unitaria en dirección x
𝜕𝑢
𝛆= = 𝐋𝐍𝐝´ = 𝐁𝐝´ (3.11)
𝜕𝑥
donde
𝛆 tensor simétrico de deformaciones de segundo orden
𝐋 matriz que contiene operadores diferenciales
𝐁 matriz de transformación
𝜕 1 1
𝐁 = 𝐋𝐍 = [𝑁1 𝑁2 ] = [− ] (3.12)
𝜕𝑥 𝑙 𝑙
donde
σ𝑥 esfuerzo normal a la dirección x, en MPa
𝐸 módulo de elasticidad, en MPa
d) Matriz de rigidez. Para obtener la matriz local de rigidez, se sustituyen las ecs
3.12 y 3.13 en la ec 2.63, entonces
1
2 −
1 1
𝐤 = ∫ 𝐁 𝑇 𝐂𝐁𝑑𝛺 = ∫ [ 𝑙 ] 𝐸𝐴 [− ] 𝑑𝑥 (3.14)
𝛺 1
1 𝑙 𝑙
𝑙
𝐸𝐴 1 −1
𝐤= [ ] (3.15)
𝑙 2 −1 1
donde
𝐤 matriz local de rigidez
66
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
Para el caso de elementos con sección variable se definirá un área inicial 𝐴1 , una área
final 𝐴2 y una área media 𝐴𝑚 como se indica en ec 3.15 y la Fig 3.4. La matriz de rigidez
es igual a la que se obtiene con la ec 3.16, con 𝐴 = 𝐴𝑚 .
𝐴1 + 𝐴2
𝐴𝑚 = (3.16)
2
Nodo 1 Nodo 2
donde
𝐴𝑚 área promedio de la sección transversal, en m²
𝐴1 área de la sección transversal en el nodo 1, en m²
𝐴2 área de la sección transversal en el nodo 2, en m²
e) Vector de fuerzas nodales. El vector local de fuerzas tiene la misma forma que
el vector local de desplazamientos, ec 3.8, es decir
𝑃
𝐟 = 𝐟𝑏 + 𝐟𝑠 = { 1 } (3.17)
𝑃2
donde
𝐟𝑏 vector local de fuerzas nodales debidos a las fuerzas de cuerpo
𝐟𝑠 vector local de fuerzas nodales debido a las fuerzas de superficie
g) Elemento armadura tres dimensiones. Se asume que los nodos locales 1 y 2 del
elemento coinciden con los nodos globales, como se muestra en la Fig 3.3. El
desplazamiento axial de cada nodo tiene tres componentes en el sistema global, 𝑥, 𝑦, 𝑧.
67
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.3 ELEMENTOS ESTRUCTURALES PARA MODELAJE
𝐝´ = 𝐓𝐨 𝐝 (3.18)
d𝑥1
d𝑦1
d
𝐝 = 𝑧1 (3.19)
d𝑥2
d𝑦2
[ d𝑧2 ]
l m n 0 0 0
𝐓𝐨 = [ ] (3.20)
0 0 0 l m n
donde
𝐝 vector global de desplazamientos
d𝑥1 , d𝑦1 , d𝑧1 , d𝑥2 , d𝑦2 , d𝑧2 coeficientes del vector global de desplazamientos de
una armadura 3D, en m
𝐓𝐨 tensor de transformación
l, m, n cosenos directores, adimensional
𝑥1 − 𝑥2
l = cos(𝑥´, 𝑥) =
𝑙
𝑦1 − 𝑦2
m = cos(𝑦´, 𝑦) = (3.21)
𝑙
𝑧1 − 𝑧2
n = cos(𝑧´, 𝑧) =
𝑙
donde
𝑥, 𝑦, 𝑧 ejes coordenados globales
𝑥′, 𝑦′, 𝑧′ ejes coordenados locales
Esta misma transformación también aplica entre los vectores de fuerzas nodales, locales
𝐟 y globales, 𝐅.
𝐟 = 𝐓𝐨 𝐅 (3.22)
donde
𝐅 vector global de fuerzas
68
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
F𝑥1
F𝑦1
F
𝐅 = 𝑧1 (3.23)
F𝑥2
F𝑦2
[ F𝑧2 ]
donde
F𝑥1 , F𝑦1 , F𝑧1 , F𝑥2 , F𝑦2 , F𝑧2 coeficientes del vector global de fuerzas de una
armadura 3D
Para transformar la matriz de rigidez elemental del sistema local al sistema global, se
sustituye la ec 3.18 en la ec 2.61 y se aplica el principio de la energía potencial total.
donde
𝐊 matriz global de rigidez
l m 0 0
𝐓𝐨 = [ ] (3.25)
0 0 l m
l2 lm −l2 −lm
𝐴𝐸 lm m2 −lm −m2 ]
𝐊 = 𝐓𝐨𝐓 𝐤𝐓𝐨 = [ (3.26)
𝑙 −l2 −lm l2 lm
−lm −m2 lm m2
69
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.3 ELEMENTOS ESTRUCTURALES PARA MODELAJE
Esfuerzo de
compresión Eje Neutro
Esfuerzo de tensión
1. Eje neutro. Se considera que el eje longitudinal es recto y que las secciones
transversales a lo largo de la viga son simétricas con respecto a este eje.
70
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
q
y, v
y, v
x, u z, w
Eje
l
neutro
71
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.3 ELEMENTOS ESTRUCTURALES PARA MODELAJE
z, w
x, u
z, w
x, u
donde
𝜃 giro de la elástica, en rad
𝑑𝑧
𝜃= (3.27)
𝑑𝑥
y
𝑢(𝑥, 𝑦) = −𝑦ℎ 𝜃
yh
donde
𝑦ℎ distancia medida a partir del eje neutro, en m
72
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
𝜕𝑢 𝜕2𝑣 𝑑2 𝑣
ε𝑥 = = −𝑦ℎ 2 = −𝑦ℎ 2 = −𝑦ℎ 𝜅 (3.28)
𝜕𝑥 𝜕𝑥 𝑑𝑥
donde
𝜅 curvatura del eje deformado de un elemento estructural, en m
𝑑2𝑣
𝜅= (3.29)
𝑑𝑥 2
𝑑2𝑣
σ𝑥 = 𝐸ε𝑥 = −𝐸𝑦ℎ = −𝐸𝑦ℎ 𝜅 (3.30)
𝑑𝑥 2
Mz = - x yh dA
A
x
z Vy
x
Mz
donde
M𝑧 momento flexionante alrededor del eje “z”, en N-m
V𝑦 cortante en dirección “y”, en N
𝑑2𝑣
M𝑧 = ∫ −𝑦ℎ σ𝑥 𝑑𝐴 = ∫ 𝑦ℎ2 𝐸 𝑑𝐴 = 𝐸𝐼𝑧 𝜅 (3.31)
𝐴 𝐴 𝑑𝑥 2
donde
𝐼𝑧 momento de inercia de la sección transversal, en m4
2
y, v
V2
1 EIz
V1
Nodo 1 Nodo 2
x, u
donde
θ1 , θ2 giro en el nodo 1 y 2, en rad
V1 , V2 cortante del nodo 1 y 2 en N
𝑣 = 𝐍𝐮 (3.32)
74
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
3𝑥 2 2𝑥 3
𝑁𝑣1 = 1 − + 3
𝑙2 𝑙
2𝑥 2 𝑥 3
𝑁𝑤1 = 𝑥 − + 2 (3.34)
𝑙 𝑙
3𝑥 2 2𝑥 3
𝑁𝑣2 = − 3
𝑙2 𝑙
𝑥2 𝑥3
𝑁𝑤2 = − + 2
𝑙 𝑙
donde
𝑁𝑣1 , 𝑁𝑤1 , 𝑁𝑣2 , 𝑁𝑤2 coeficientes de la matriz que contiene las funciones de forma
de un elemento viga, adimensional
ε = 𝐁𝐮 (3.36)
1
= 𝐁𝐮 (3.38)
𝜌𝑐
donde
𝜌𝑐 radio de curvatura, en rad
σ𝑥 = 𝐸𝐁𝐮 (3.39)
Al igual que para la deformación, se relacionan los esfuerzos con el momento flexionante,
que en forma discreta se da mediante
M𝑧 = 𝐸𝐼𝑧 κ𝐮 (3.40)
75
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.3 ELEMENTOS ESTRUCTURALES PARA MODELAJE
𝑙
1 1 𝑙 𝑙
Π = ∫ σε𝑑𝛺 − ∫ 𝑞𝑣𝑑𝑥 = ∫ M𝑧 𝜅𝑑𝑥 − ∫ 𝑞𝑣𝑑𝑥 (3.41)
2 𝛺 0 2 0 0
donde
Π energía potencial del cuerpo, N-m/m
12 6𝑙 −12 6𝑙
𝐸 𝐼𝑧 4𝑙 2 2𝑙 2 ]
𝐤 = 3 [ 6𝑙 −6𝑙 (3.42)
𝑙 −12 −6𝑙 12 −6𝑙
6𝑙 2𝑙 2 −6𝑙 2𝑙 2
La matriz de rigidez dada por la ec 3.42, corresponde al elemento finito viga de acuerdo
a la teoría de Euler—Bernoulli; cuya característica fundamental, consiste en tomar en
cuenta la deformación debido al momento flexionante.
𝑓𝑦1
𝑀𝑧1
𝐟= (3.43)
𝑓𝑦2
[𝑀𝑧2 ]
𝑑′𝑦1
𝜃′𝑧1
𝐝′ = (3.44)
𝑑′𝑦2
[ 𝜃′𝑧2 ]
donde
𝑓𝑦1 , 𝑀𝑧1 , 𝑓𝑦2 , 𝑀𝑧2 coeficientes del vector local de fuerzas para un elemento
viga, N y N-m
′
𝑑𝑦1 ′
, 𝜃𝑧1 ′
, 𝑑𝑦2 ′
, 𝜃𝑧2 coeficientes del vector local de desplazamientos de un
elemento viga, en m
76
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
Los marcos son estructuras reticulares formadas por elementos conectados entre sí de
forma rígida, y que transmite las fuerzas aplicadas a la estructura, externamente, por
medio de esfuerzos axiales, cortante y de flexión o torsión.
Un elemento marco es una barra recta en la que las dimensiones de la sección transversal
son mucho menores que la longitud del elemento. Pueden ser marcos en dos o tres
dimensiones. Para que un marco sea plano, Fig 3.12a, se deben cumplir las siguientes
condiciones:
1 Los ejes geométricos de todas las barras están contenidos en un mismo plano, que
es el plano de la estructura. Se asume que este plano es el 𝑥 − 𝑦.
2 Las fuerzas exteriores aplicadas a la estructura, están contenidas en el mismo plano
de la estructura, y los momentos son perpendiculares a él.
posible ya que ambos efectos están desacoplados, i.e., uno es independiente del otro. De
acuerdo a esto, las ecuaciones de equilibrio interno de un elemento marco corresponden a
las ecuaciones de un elemento armadura más las del elemento viga.
2 5
3 6
x
1 4
Nodo 1 Nodo 2
𝐤𝐝′ = 𝐟 (3.45)
donde
𝐸𝐴 𝐸𝐴
0 0 − 0 0
𝑙 𝑙
12𝐸𝐼𝑧 6𝐸𝐼𝑧 12𝐸𝐼𝑧 6𝐸𝐼𝑧
0 0 −
𝑙3 𝑙2 𝑙3 𝑙2
6𝐸𝐼𝑧 4𝐸𝐼𝑧 6𝐸𝐼𝑧 2𝐸𝐼𝑧
0 0 − 2
𝐤= 𝑙2 𝑙 𝑙 𝑙 (3.46)
𝐸𝐴 𝐸𝐴
− 0 0 0 0
𝑙 𝑙
12𝐸𝐼𝑧 12𝐸𝐼 12𝐸𝐼𝑧 6𝐸𝐼𝑧
0 − − 0 −
𝑙3 𝑙3 𝑙3 𝑙2
6𝐸𝐼𝑧 2𝐸𝐼𝑧 6𝐸𝐼𝑧 4𝐸𝐼𝑧
[ 0 𝑙2 𝑙
0 − 2
𝑙 𝑙 ]
′ ′ ′ ′ ′ ′ 𝑇
𝐝′ = [𝑑𝑥1 𝑑𝑦1 𝜃𝑧1 𝑑𝑥2 𝑑𝑦2 𝜃𝑧2 ] (3.47)
donde
′
𝑑𝑥1 ′
, 𝑑𝑦1 ′
, 𝜃𝑧1 ′
, 𝑑𝑥2 ′
, 𝑑𝑦2 ′
, 𝜃𝑧2 coeficientes del vector local de desplazamientos de
un elemento marco, en m
𝑓x1 , 𝑓y1 , 𝑀z1 , 𝑓x2 , 𝑓y2 , 𝑀z2 coeficientes del vector local de fuerzas de un
elemento marco, en N
Sistemas de referencia
están alineados uno después del otro, misma orientación, formando una línea recta que no
necesariamente es una viga continua.
Para un caso general en el que los elementos pueden tener diferentes orientaciones los
sistemas de referencia locales de los elementos no coinciden con el sistema de referencia
global de la estructura, Fig 3.14.
y'
x'
de es
as cal 2 3
m lo x'
ste a
y Si enci
er
ref
x'
y'
x 1 4
Sistema de
referencia global
y'
x'
Una vez obtenidas las ecuaciones para los elementos, es necesario transformarlas del
sistema local al global mediante un Matriz de rotación, Fig 3.15, que se calcula para cada
uno de los elementos de acuerdo a su orientación.
y
y'
x'
'
x
79
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.3 ELEMENTOS ESTRUCTURALES PARA MODELAJE
𝑐𝑜𝑠𝛼′ −𝑠𝑒𝑛𝛼′ 0
𝐑 = [ 𝑠𝑒𝑛𝛼′ 𝑐𝑜𝑠𝛼′ 0] (3.49)
0 0 1
donde
𝛼’ ángulo entre los ejes coordenados globales y locales, en rad
𝐑 matriz de rotación
𝐝´ = 𝐑𝐝 (3.50)
𝐟 = 𝐑𝐅 (3.51)
Sustituyendo las ecs 3.50 y 3.51 en la ecuación de energía potencial total del sistema se
tiene
1
Π = (𝐝′ )𝑻 𝐤𝐝´ − (𝐝′ )𝑻 𝐟 (3.52)
2
1
Π = (𝐑𝐝)𝑇 𝐤(𝐑𝐝) − (𝐑𝐝)𝑇 𝐟 (3.53)
2
𝐊𝐝 = 𝐅 (3.54)
𝐊 = 𝐑𝑇 𝐤𝐑 (3.55)
80
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
y 8
7
11
10
Nodo 2
9
12
z x
2
Eje neutro
5
Nodo 1
1 3
6
y' 4
x'
z'
donde
1, 2, 3…..12 grados de libertad de un elemento marco 3D, adimensional
𝐴 𝐴
0 0 0 0 0 0 0 0 0 0
𝑙 𝑙
12𝐼𝑧 6𝐼𝑧 12𝐼𝑧 6𝐼𝑧
0 0 0 0 0 − 0 0 0
𝑙3 𝑙2 𝑙3 𝑙2
12𝐼𝑦 6𝐼𝑦 12𝐼𝑦 6𝐼𝑦
0 0 0 − 0 0 0 − 0 0
𝑙3 𝑙2 𝑙3 𝑙2
𝐽 𝐽
0 0 0 0 0 0 0 0 − 0 0
2(1 + 𝜇)𝑙 2(1 + 𝜇)𝑙
6𝐼𝑦 4𝐼𝑦 6𝐼𝑦 2𝐼𝑦
0 0 − 0 0 0 0 0 0
𝑙2 𝑙 𝑙2 𝑙
6𝐼𝑧 4𝐼𝑧 6𝐼𝑧 2𝐼𝑧
0 0 0 0 0 0 0 0
𝐤=𝐸 𝑙2 𝑙 𝑙2 𝑙
𝐴 𝐴
− 0 0 0 0 0 0 0 0 0 0
𝑙 𝑙
12𝐼𝑧 6𝐼𝑧 12𝐼𝑧 6𝐼𝑧
0 − 3 0 0 0 0 0 0 0 −
𝑙 𝑙2 𝑙3 𝑙2
12𝐼𝑦 6𝐼𝑦 12𝐼𝑦 6𝐼𝑦
0 0 − 0 0 0 0 0 0
𝑙3 𝑙2 𝑙3 𝑙2
𝐽 𝐽
0 0 0 − 0 0 0 0 0 0 0
2(1 + 𝜇)𝑙 2(1 + 𝜇)𝑙
6𝐼𝑦 2𝐼𝑦 6𝐼𝑦 4𝐼𝑦
0 0 0 0 0 0 0 0
𝑙2 𝑙 𝑙2 𝑙
6𝐼𝑧 2𝐼𝑧 6𝐼𝑧 4𝐼𝑧
[ 0 𝑙2
0 0 0
𝑙
0 −
𝑙2
0 0 0
𝑙 ]
(3.56)
𝐝´ = [𝑑𝑥1 𝑑𝑦1 𝑑𝑧1 𝜃𝑥1 𝜃𝑦1 𝜃𝑧1 𝑑𝑥2 𝑑𝑦2 𝑑𝑧2 𝜃𝑥2 𝜃𝑦2 𝜃𝑧2 ]𝑻 (3.57)
81
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.3 ELEMENTOS ESTRUCTURALES PARA MODELAJE
𝐟 = [𝑓𝑥1 𝑓𝑦1 𝑓𝑧1 𝑀𝑥1 𝑀𝑦1 𝑀𝑧1 𝑓𝑥2 𝑓𝑦𝑥 𝑓𝑧2 𝑀𝑥2 𝑀𝑦2 𝑀𝑧2 ]𝑻 (3.58)
donde
′
𝑑𝑥1 ′
, 𝑑𝑦1 ′
, 𝑑𝑧1 ′
, 𝜃𝑥1 ′
, 𝜃𝑦1 ′
, 𝜃𝑧1 ′
, 𝑑𝑥2 ′
, 𝑑𝑦2 ′
, 𝑑𝑧2 ′
, 𝜃𝑥2 ′
, 𝜃𝑦2 ′
, 𝜃𝑧2 coeficientes del vector de
desplazamientos de un
elemento marco 3D, en m y rad
𝑓𝑥1 , 𝑓𝑦1 , 𝑓𝑧1 , 𝑀𝑥1 , 𝑀𝑦1 , 𝑀𝑧1 , 𝑓𝑥2 , 𝑓𝑦2 , 𝑓𝑧2 , 𝑀𝑥2 , 𝑀𝑦2 , 𝑀𝑧2 coeficientes del vector local de
fuerzas de un elemento marco
3D, en N y N-m
𝐼𝑥 , 𝐼𝑦 y 𝐼𝑧 momentos de inercia de la
sección transversal del
elemento, en m4
𝐽 constante de torsión de la
sección, en m4
La rotación de las ecs 3.56 a 3.58 al sistema de referencia global, Fig 3.17, se realiza de la
misma manera que para el caso de un elemento marco, utilizando la siguiente matriz de
rotación,
𝐫𝒐 0 0 0
0 𝐫𝒐 0 0
𝐓𝐨 = [ ] (3.59)
0 0 𝐫𝒐 0
0 0 0 𝐫𝒐
donde
𝐫𝒐 submatriz del tensor de transformación
M y2
P y2 x'
y'
M y1 Nodo 2
Px2 M x2
l Pz2
P y1
M x1 Px1 Mz2
Nodo 1 x
Pz1
Mz1
z'
donde
𝑃𝑥1 , 𝑃𝑦1 , 𝑃𝑧1 , 𝑀𝑥1 , 𝑀𝑦1 , 𝑀𝑧1 fuerzas en el nodo 1 en un elemento marco 3D, en N y N-m
𝑃𝑥2 , 𝑃𝑦2 , 𝑃𝑧2 , 𝑀𝑥2 , 𝑀𝑦2 , 𝑀𝑧2 fuerzas en el nodo 2 en un elemento marco 3D, en N y N-m
𝑟𝑥 𝑋 𝑟𝑥 𝑌 𝑟𝑥 𝑍
𝐫𝒐 = [ 𝑦 𝑋
𝑟 𝑟𝑦 𝑌 𝑟𝑦 𝑍] (3.60)
𝑟𝑧 𝑋 𝑟𝑧 𝑌 𝑟𝑧 𝑍
donde
𝑟𝑥 𝑋, 𝑟𝑥 𝑌, 𝑟𝑥 𝑍 … 𝑟𝑧 𝑍 coeficientes de la submatriz 𝐫𝒐
El primer renglón de la matriz dada por la ec 3.60 corresponde a los cosenos directores del
eje local 𝑥. Los renglones dos y tres son, respectivamente, los cosenos directores de los
ejes locales 𝑦 y 𝑧.
Un arco es una estructura plana constituida por un elemento curvo de sección transversal
despreciable frente a su longitud, y cuya curvatura es pequeña comparada con su sección
transversal. Los dos puntos extremos pueden estar apoyados de distintas formas y las
cargas exteriores están aplicadas por lo general de manera vertical.
Los arcos son una de las estructuras más utilizadas desde la antigüedad. Ello es debido a
que, si su geometría es adecuada, soportan grandes cargas transversales y las transmiten
a los apoyos extremos trabajando básicamente a compresión, con muy poco esfuerzo de
flexión. Esto permite utilizar en su construcción, materiales que no soportan grandes
niveles de tensión, como el concreto o simplemente ladrillos o bloques de roca, Fig 3.18.
Los arcos están normalmente sometidos a grandes niveles de carga vertical, aplicadas bien
desde la parte superior del arco o desde la parte inferior, como se muestra en la Fig 3.19,
83
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.3 ELEMENTOS ESTRUCTURALES PARA MODELAJE
así como a cargas horizontales debidas a cargas de viento, frenado de vehículos, etc. Son
también frecuentes las cargas térmicas o las debidas a los asentamientos de los apoyos,
que pueden ser importantes en arcos de gran tamaño.
V + V
s M + M
M
N N + N
V
donde
∆V incremento de fuerza cortante, en N
∆M incremento de momento flexionante, en N-m
∆N incremento de fuerza normal, en N
Δ𝑠 elemento arco aislado, en m
Δ𝜑 ángulo correspondiente a un segmento de arco aislado, en rad
M momento flexionante, en N-m
V cortante, en N
𝑑V N
=𝑞+ (3.61)
𝑑𝑥 R
y el equilibrio de momentos
𝑑M
= −V (3.62)
𝑑𝑥
donde
R radio del segmento del arco, en m
P
P P
Fuerza (P)
Cable suspendido
Desplazamiento (u)
1. Analíticos
2. Simplificados
3. Numéricos
85
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.3 ELEMENTOS ESTRUCTURALES PARA MODELAJE
Dentro del tercer grupo de métodos utilizados para analizar las estructuras a base de
cables se encuentran, principalmente, el método de las diferencias finitas y el método
de los elementos finitos. El enfoque adoptado en este manual, es el correspondiente a
este último método.
donde
𝑥1 , 𝑦1 coordenada del nodo 1 de un elemento cable, en m
𝑥2 , 𝑦2 coordenada del nodo 2 de un elemento cable, en m
86
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
y2
Nodo 2
y1
Nodo 1
x x x
𝑥1 − 𝑦1 𝑥2 − 𝑦2
l= , m= (3.64)
𝑙 𝑙
a) Cargas. Las cargas aplicadas en los nodos 1,2 del elemento, referidas al sistema
global de coordenadas, son:
donde
𝑓𝑥1 , 𝑓𝑦1 , 𝑓𝑥2 , 𝑓𝑦2 coeficientes de del vector local de fuerzas de un elemento
cable, en N
87
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.3 ELEMENTOS ESTRUCTURALES PARA MODELAJE
x'
Nodo 2
y'
Nodo 1
y
f y2
Nodo 2
f x2
f y1
Nodo 1 f x1
𝑇
𝐝′ = [𝑑𝑥1
′ ′
, 𝑑𝑦1 ′
, 𝑑𝑥2 ′
, 𝑑𝑦2 ] (3.66)
donde
′
𝑑𝑥1 ′
, 𝑑𝑦1 ′
, 𝑑𝑥2 ′
, 𝑑𝑦2 coeficientes de un vector local de desplazamientos de un
elemento cable, en m
88
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
Nodo 2'
u y2
ux2
Nodo 2
Nodo 1'
u y1
Nodo 1
ux1
donde
u𝑥1 , u𝑥2 , u𝑦1 , u𝑦2 desplazamiento del nodo 1 y 2 en dirección 𝑥, 𝑦, en m
La matriz global de rigidez de un elemento cable se obtiene entre la relación del vector
global de fuerzas 𝐅 y desplazamientos 𝐝, ec 3.67
𝐊𝐝 = 𝐅 (3.67)
89
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.3 ELEMENTOS ESTRUCTURALES PARA MODELAJE
y, v y, v
z, w x, u
x, u
z, w
Fig 3.26 Elementos planos para esfuerzo plano (a) y deformación plana (b)
P
bpy
P bpx
x
t pl
donde
𝑏𝑝𝑥 longitud de un elemento plano en dirección x, en m
𝑏𝑝𝑦 longitud de un elemento plano en dirección y, en m
𝑡𝑝𝑙 espesor de un elemento plano, en m
90
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
Nodo 3
Nodo 2 ux3
ux2 u y3
u y2
Nodo 1
Nodo 4
ux1
u y1 ux4
u y4
donde
u𝑥1 , u𝑦1 , u𝑥2 , u𝑦2 , u𝑥3 , u𝑦3 , u𝑥4 , u𝑦4 desplazamiento del nodo 1, 2, 3, 4 en dirección 𝑥
y 𝑦, en m
El vector de desplazamientos nodales para este tipo de elemento, está dado de manera
general como,
𝑢𝑥
𝐮 = [𝑢 ] (3.68)
𝑦
El tamaño del vector de la ecuación anterior, dependerá del tipo de elemento de esfuerzo
plano a utilizar. Para un elemento de tres nodos, el tamaño es de seis por uno y para el
elemento de cuatro nodos es de ocho por uno. El tamaño de este vector se obtiene como
el producto del número de nodos del elemento por número de grados de libertad por
nodo.
donde
𝜇 relación de Poisson, adimensional
91
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.3 ELEMENTOS ESTRUCTURALES PARA MODELAJE
𝜎𝑥
𝜎𝑦
𝝈 = {𝜎 = 0 } (3.71)
𝑧
𝜎𝑥𝑦 = 𝜎𝑦𝑥
𝜎𝑥 1 𝜇 0 𝜀𝑥
𝜎
[ 𝑦] 𝐸 𝜇 1 0 [ 𝜀𝑦 ]
𝜎𝑥𝑦 = 1 − 𝜇 2
[ ]
1 − 𝜇 𝛾𝑥𝑦 (3.72)
0 0
⏟ ⏟ 2 ⏟
σ ε
.
donde
𝜎𝑥 , 𝜎𝑦 , 𝜎𝑥𝑦 coeficientes de un tensor de esfuerzos de un elemento plano,
en MPa
𝜀𝑥 , 𝜀𝑦 , 𝛾𝑥𝑦 coeficientes de tensor simétrico deformaciones, adimensional
𝐟 = ∫ 𝐍 𝑇 𝐟Γ 𝑑Γ (3.74)
Γ
donde
𝐟Γ vector de fuerzas distribuidas en la frontera
92
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
Derivado del hecho de que en el MEF todas las fuerzas deben de transformarse para ser
aplicadas en los nodos de la estructura; por lo general las fuerzas externas se evalúan
directamente de las cargas aplicadas en el cuerpo y no es necesario aplicar la ec 3.74
cuando las cargas son puntuales. En caso de tener cargas distribuidas en los lados del
elemento, es necesario resolver la integral indicada, numéricamente, para obtener el
vector de cargas en coordenadas locales del elemento plano correspondientes a cada
uno de los nodos, ec 3.75.
𝑓𝑝𝑥1
𝑓𝑝𝑦1
𝑓𝑝𝑥2
𝑓𝑝𝑦2
𝐟= (3.75)
𝑓𝑝𝑥3
𝑓𝑝𝑦3
𝑓𝑝𝑥4
𝑓
[ 𝑝𝑦4 ]
donde
𝑓𝑝𝑥1 , 𝑓𝑝𝑦1 , 𝑓𝑝𝑥2 , 𝑓𝑝𝑦2 , 𝑓𝑝𝑥3 , 𝑓𝑝𝑦4 , 𝑓𝑝𝑥4 , 𝑓𝑝𝑦4 coeficientes de un vector de fuerzas de un
elemento plano, en N
t pl=1
z x
93
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.3 ELEMENTOS ESTRUCTURALES PARA MODELAJE
Este tipo de elemento tiene como característica principal que su espesor se considera
unitaria; a diferencia del elemento de esfuerzo plano donde el espesor puede variar. La
cinemática de deformaciones establece que la componente 𝜀𝑧 = 0. Dentro de este tipo de
modelo se pueden considerar problemas en los que la dimensión perpendicular al plano
del elemento es muy grande en comparación con sus dimensiones en planta.
(1 − 𝜇) 𝜇 𝜇 0
𝜎𝑥 𝜀𝑥
𝜇 (1 − 𝜇) 𝜇 0
𝜎𝑦 𝐸 𝜀𝑦
{𝜎 }= (1 − 𝜇) 0 { 𝜀 } (3.76)
𝑧 (1 + 𝜇)(1 − 2𝜇) 𝜇 𝜇 𝑧
𝜎𝑥𝑦 (1 − 2𝜇) 𝛾
0 0 0 𝑥𝑦
[ 2 ]
Para este caso, deformación plana, se puede observar que a pesar de que la componente
de deformación 𝜀𝑧 = 0, si existe el esfuerzo normal en esta dirección; a menos de que la
relación de poisson se igual a cero.
𝐟 = ∫ 𝐍 𝑇 𝐟Γ 𝑑Γ (3.77)
Γ
Las placas son estructuras continuas cuya geometría se caracteriza por tener un espesor
menor a las otras dos dimensiones en planta del elemento. Se considera que, de manera
normal al plano actúan fuerzas o momentos, ver Fig 3.30. Las cargas externas provocan
un comportamiento estructural del elemento dominado por flexión y deformación
transversal. En este tipo de elemento no se consideran cargas contenidas en el plano.
94
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
Pz
y
w
y tp
x
x
My
donde
𝑡𝑝 espesor del elemento placa, en m
M𝑥 , M𝑦 momentos flexionantes alrededor de x, y
θ𝑥 giro alrededor del eje x, en rad
θ𝑦 giro alrededor del eje y, en rad
P𝑧 fuerza en dirección z, en N
La formulación de un elemento finito placa que cumpla con estas dos consideraciones
básicas, conduce a un elemento cuya aplicación da soluciones exactas; de acuerdo a la
teoría de flexión de placas.
95
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.3 ELEMENTOS ESTRUCTURALES PARA MODELAJE
consideran despreciables. Asimismo, se admite que los esfuerzos cortantes 𝜎𝑥𝑦 varían
linealmente con 𝑡𝑝 y que los esfuerzos cortantes 𝜎𝑦𝑧 y 𝜎𝑧𝑥 varían de forma cuadrática.
Vy
y y
M yx
My
P Vx
z z
Mxy Mx
xz
Mxy
Mx
x
Vx
xy
tp x x
yz yx My
y M yx Vy
(a) (b)
donde
P𝑧 fuerza en dirección z, en N
V𝑥 cortante en dirección x, en N
V𝑦 cortante en dirección y, en N
M𝑦𝑥 cortante en dirección x, normal a y, en N
M𝑥𝑦 cortante en dirección y, normal a x, en N
M𝑥 momento flexionante en x, en N-m
M𝑦 momento flexionante en y, en N-m
σ𝑥 esfuerzo normal dirección x, en MPa
σ𝑦 esfuerzo normal dirección y, en MPa
σ𝑦𝑧 esfuerzo cortante en dirección z normal a y, en MPa
σ𝑦𝑥 esfuerzo cortante dirección x normal a y, en MPa
σ𝑥𝑦 esfuerzo cortante dirección y normal a x, en MPa
σ𝑥𝑧 esfuerzo cortante dirección z normal a x, en MPa
Fig 3.31 (a) configuración de esfuerzos en los ejes del elemento diferencial
(b) momentos y fuerzas internas en el mismo elemento
En los problemas de flexión es preferible utilizar, en lugar de los esfuerzos, los momentos
flexionantes y las fuerzas correspondientes a los esfuerzos; que se obtienen al integrar
los esfuerzos en todo el espesor 𝑡𝑝 , mediante:
1 1 1
𝑡 𝑡 𝑡
2 𝑝 2 𝑝 2 𝑝
M𝑥 = ∫ 𝑧σ𝑥 𝑑𝑧, M𝑦 = ∫ 𝑧σ𝑦 𝑑𝑧 M𝑥𝑦 = ∫ 𝑧σ𝑥𝑦 𝑑𝑧 (3.78)
1 1 1
− 𝑡𝑝 − 𝑡𝑝 − 𝑡𝑝
2 2 2
𝑀𝑥 1
𝑡 𝜎𝑥 1
𝑡
2 𝑝 2 𝑝
𝐌𝒑 = [ 𝑀𝑦 ] = ∫ [ 𝜎𝑦 ] 𝑧𝑑𝑧 ; 𝐌𝒑 = ∫ 𝝈𝑧𝑑𝑧 (3.79)
1 1
𝑀𝑥𝑦 − 𝑡𝑝 𝜎𝑥𝑦 − 𝑡𝑝
2 2
96
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
1 1
𝑡 𝑡
𝑉𝑧𝑥 2 𝑝 𝜎𝑧𝑥 2 𝑝
[𝑉 ] = ∫ [𝜎 ] 𝑑𝓏 ; 𝐕𝒑 = ∫ 𝝉𝑑𝑧 (3.80)
𝑧𝑦 1 𝑧𝑦 1
− 𝑡 − 𝑡𝑝
2 𝑝 2
donde
𝝉 vector de esfuerzo cortante
𝐌𝒑 vector de momentos flexionantes de un elemento placa
𝐕𝒑 vector de fuerzas cortantes de un elemento placa
𝑀𝑥 , 𝑀𝑦 , 𝑀𝑥𝑦 , 𝜎𝑥 , 𝜎𝑦 , 𝜎𝑥𝑦 coeficientes del vector de momentos flexionantes de un
elemento placa, en N-m y Mpa
𝑉𝑧𝑥 , 𝑉𝑧𝑦 , 𝜎𝑧𝑥 , 𝜎𝑧𝑦 coeficientes del vector de fuerzas cortantes de un
elemento placa, en N-m y Mpa
𝜕 2 𝑀𝑥 𝜕 2 𝑀𝑥𝑦 𝜕 2 𝑀𝑦
+ 2 + +𝑞 =0 (3.81)
𝜕𝑥 2 𝜕𝑥𝜕𝑦 𝜕𝑦 2
donde
𝑞 carga uniformemente distribuida, en N/m
Son placas en las que los esfuerzos debido al efecto membrana son muy pequeños en
comparación con los esfuerzos de flexión. Una placa delgada es cuando la dimensión
menor excede 20 veces su espesor.
97
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.3 ELEMENTOS ESTRUCTURALES PARA MODELAJE
y u
a
Superficie media
v
p= w
x
b
donde
𝛼𝑝 ángulo de la superficie rotada de una placa, en grados y/o radianes
3. Se desprecian los efectos que pudieran tener los esfuerzos normales 𝜎𝑧 , por lo que
𝜎𝑧 = 0.
𝐸
σ𝑥 = (𝜀 + 𝜇𝜀𝑦 )
1 − 𝜇2 𝑥
𝐸
σ𝑦 = (𝜀 + 𝜇𝜀𝑥 )
1 − 𝜇2 𝑦
σ𝑥𝑦 = 𝐺𝛾𝑥𝑦 (3.82)
donde
𝛾𝑥𝑦 , 𝛾𝑥𝑧 , 𝛾𝑦𝑧 deformaciones angulares, adimensional
98
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
Para tener una formulación consistente con el método de los elementos finitos, se
sustituyen las ecuaciones de esfuerzo y deformación en la ec 3.79
𝐷(1 − 𝜇)
M𝑥 = 𝐷(𝜅𝑥 + 𝜇𝜅𝑦 ), M𝑦 = 𝐷(𝜅𝑦 + 𝜇𝜅𝑥 ), M𝑥𝑦 = 𝜅𝑥𝑦 (3.83)
2
𝐸𝑡𝑝3
𝐷= (3.84)
[12(1 − 𝜇 2 )]
donde
𝐷 módulo de rigidez a flexión equivalente de la placa, en N/m.
𝜅, 𝜅𝑥 , 𝜅𝑥𝑦 curvatura del eje deformado de un elemento viga, placa, entre
otros, en rad
Matricialmente se tiene
𝜕2𝑤 𝜕2𝑤
−𝐷 ( + 𝜇 )
𝜕𝑥 2 𝜕𝑦 2
𝑀𝑥
𝜕2𝑤 𝜕2𝑤
𝐌𝒑 = [ 𝑀𝑦 ] = −𝐷 ( 2 + 𝜇 2 ) (3.85)
𝑀𝑥𝑦 𝜕𝑦 𝜕𝑥
𝜕2𝑤
−𝐷(1 − 𝜇)
[ 𝜕𝑥𝜕𝑦 ]
𝜕4𝑤 2𝜕 4 𝑤 𝜕4𝑤
𝐷( + + )=𝑞 (3.86)
𝜕𝑥 4 𝜕𝑥 2 𝜕𝑦 2 𝜕𝑦 4
99
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.3 ELEMENTOS ESTRUCTURALES PARA MODELAJE
siendo pequeño para que las deformaciones en la dirección normal a la placa puedan ser
despreciadas.
Tangente al
eje medio
Eje medio
y,v x,u
(a) (b)
donde
θ𝑥 giro alrededor del eje x, en rad
θ𝑦 giro alrededor del eje y, en rad
En esta teoría deja de ser válida la hipótesis de Kirchhoff, por lo que las secciones
rectas y normales al plano medio de la placa permanecen rectas, pero no siguen siendo
normales después de la deformación. Es decir que el ángulo que gira una sección plana
cualquiera no tiene que coincidir con la tangente a la deformada. Las tres
consideraciones fundamentales en la teoría de Mindlin, para placas, son:
𝜎𝑥𝑧 𝐺𝛾𝑥𝑧
𝝉𝒑 = [𝜎 ] = [𝐺𝛾 ] = G𝜸 (3.87)
𝑦𝑥 𝑦𝑥
100
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
donde
𝝉𝒑 vector de esfuerzos cortantes de un elemento placa
𝐺 módulo de cortante, en MPa
𝜸 vector de deformación angular
𝜎𝑥𝑧 , 𝜎𝑦𝑥 , 𝐺𝛾𝑥𝑧 , 𝐺𝛾𝑦𝑥 coeficientes del vector de esfuerzos cortantes de un
elemento placa
𝜕θ𝑦 𝜕θ𝑥
𝐷(
−𝜇 )
𝜕𝑥 𝜕𝑦
𝑀𝑥 M0𝑥
𝜕θ𝑦 𝜕θ𝑥
𝐌𝑝 = { 𝑀𝑦 } = 𝐷 (𝜇 − ) + { 𝑀0𝑦 } (3.88)
𝑀𝑥𝑦 𝜕𝑥 𝜕𝑦 𝑀0𝑥𝑦
𝐷(1 − 𝜇) 𝜕θ𝑦 𝜕θ𝑥
( − )
{ 2 𝜕𝑦 𝜕𝑥 }
donde
𝑀0𝑥 , 𝑀0𝑦 , 𝑀0𝑥𝑦 coeficientes del vector de momentos flexionantes
iniciales de un elemento placa, en N-m
las fuerzas cortantes son las resultantes de los esfuerzos cortantes están dados como
1
𝑡
2 𝑝
V=∫ 𝝉𝑝 𝑑𝑧 = 𝑦ℎ 𝐺𝜸 (3.89)
1
− 𝑡𝑝
2
Una membrana es una estructura con una superficie delgada y flexible que soporta cargas,
principalmente, por medio del desarrollo de esfuerzos de tensión. Las superficies de una
membrana tienden a adaptar sus formas de acuerdo a como se aplican las cargas. Este
tipo de estructuras son muy sensibles a los efectos aerodinámicos del viento. En
consecuencia, muchas de las membranas que se utilizan en los edificios son estabilizadas
de alguna manera para que sus formas se mantengan ante la acción de las cargas. En la
Fig 3.34 se muestra una estructura típica construida a base de membranas.
Desde el punto de vista mecánico, la forma de soportar cargas por una combinación de
efectos axiales y de flexión, los cascarones se asemejan a la forma de trabajar de los
marcos y los arcos. En los cascarones, los arcos y los marcos combinan el efecto de
flexión y cortante con el axial para soportar mayores niveles de carga. Por ello, los
cascarones tienen una mayor eficiencia de trabajo; considerando la relación peso de la
estructura contra carga soportada.
102
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
N M
𝜎𝑥 = 𝜎𝑚𝑥 + 𝜎𝑏𝑥 = + 3 (3.90)
𝜎𝑥 𝜎𝑥
12
donde
N fuerza normal, en N
M momento flexionante, N-m
𝜎𝑥 esfuerzo normal, dirección “x”, en MPa
𝜎𝑚𝑥 esfuerzos de membrana, en MPa
𝜎𝑏𝑥 esfuerzos debido a la flexión, en MPa
Se considera que los esfuerzos varían linealmente a lo largo del espesor del cascarón.
Los esfuerzos en la superficie media 𝑧 = 0, son cero. Esto se cumple si la fuerza de
membrana es cero.
1 Cilíndricos
2 Esféricos
3 Elípticos
4 Cónicos
5 Paraboloides
6 Paraboloide hiperbólico
103
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.3 ELEMENTOS ESTRUCTURALES PARA MODELAJE
z, w
Nodo 4 (-1 , +1) Nodo 3 (+1 , +1)
(u4 , v4 , w4 , (u3 , v3 , w3 ,
x4 , y4 , z4) x3 , y3 , z3)
donde
𝑢1 , 𝑣1 , 𝑤1 , 𝑢2 , 𝑣2 , 𝑤2 , 𝑢3 , 𝑣3 , 𝑤3 , 𝑢4 , 𝑣4 , 𝑤4 desplazamientos del nodo 1 a 4 de un
elemento cascarón, en m
θ𝑥1 , θ𝑦1 , θ𝑧1 , θ𝑥2 , θ𝑦2 , θ𝑧2 , θ𝑥3 , θ𝑦3 , θ𝑧3 , θ𝑥4 , θ𝑦4 , θ𝑧4 giros en el nodo 1, 2,3 y 4, alrededor del
eje 𝑥 𝑦 y 𝑧, de un elemento cascarón, en
rad
𝜉, 𝜂 coordenadas locales de un elemento
discretizado mediante el método de los
elementos finito, adimensionales
Fig 3.36 Elemento cascarón con seis grados de libertad por nodo
104
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
x
Mx M N
Nx
y z
dM x
M Mx + dxx dx
dx N x + dNx dx
dx
N
dV
Vx + dxx dx
donde
M𝑥 momento flexionante alrededor del eje x, en N-m
M𝜑 componente del vector de momentos en un elemento cascarón, en N-m
𝜑 ángulo de la curvatura de un elemento cascarón, en rad y/ grados
N𝑥 fuerza normal en dirección x de un elemento cascarón, en N
N𝜑 fuerza normal en dirección φ de un elemento cascarón, en N
Un elemento cualquiera del cascarón se aproxima por un elemento plano de tres o cuatro
nodos, como se muestra en la Fig 3.38, al que se asocia un sistema de coordenadas,
que se denomina sistema local del elemento, y en el que se pueden obtener sus
propiedades geométricas y de rigidez, considerándolo por separado como una membrana
y como una placa.
105
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.3 ELEMENTOS ESTRUCTURALES PARA MODELAJE
z' w
y
v y'
x
u
Y
Elemento plano
X x'
donde
θ𝑥 giro alrededor del eje x, en grados o rad
θ𝑦 giro alrededor del eje y, en grados o rad
Cada uno de los elementos cascarón, así modelado, tiene cinco grados de libertad: dos
de ellos como membrana o elemento de esfuerzo plano, 𝑢, 𝑣, y tres como placa, 𝜔, 𝜃𝑥 , 𝜃𝑦 ;
todos ellos referidos al sistema local del elemento. El sexto grado de libertad 𝜃𝑧 no
aparece en la formulación. Sin embargo, este grado de libertad se incluye al asignarle
una cierta rigidez; con lo que se logra tener un elemento con seis grados de libertad por
nodo, ver Fig 3.39.
y'
w y
y x'
x
v
u z
x
donde
θ𝑧 giro alrededor del eje z, en grados o rad
106
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
𝐊𝑀 0
𝐊 𝑐𝑎 = [ ] (3.91)
0 𝐊𝐹
donde
𝐊 𝑐𝑎 matriz de rigidez de un elemento cascarón
𝐊𝑀 matriz de rigidez debida al efecto membrana
𝐊𝐹 matriz de rigidez debida al efecto de flexión
𝐊 = 𝐓𝒐𝑻 𝐊 𝑐𝑎 𝐓𝐨 (3.92)
𝐅𝑒 = 𝐓 𝑇 𝐟𝑒 (3.93)
donde
𝐓𝐨 tensor de transformación, está dado como,
𝐓𝟑 0 0 0 0 0 0 0
0 𝐓𝟑 0 0 0 0 0 0
0 0 𝐓𝟑 0 0 0 0 0
𝐓𝟑 0 0 0 0
𝐓𝐨 = 0 0 0
0 (3.94)
0 0 0 0 𝐓𝟑 0 0
0 0 0 0 0 𝐓𝟑 0 0
0 0 0 0 0 0 𝐓𝟑 0
[0 0 0 0 0 0 0 𝐓𝟑 ]
l𝑥 m𝑥 n𝑥
𝐓3 = [l𝑦 m𝑦 n𝑦 ] (3.95)
l𝑧 m𝑧 n𝑧
donde
𝐓3 submatriz de cosenos directores
l𝑥 , m𝑥 , n𝑥 , l𝑦 , m𝑦 , n𝑦 , l𝑧 , m𝑧 , n𝑧 componentes de la submatriz de cosenos
directores entre ambos sistemas de referencia,
mostrados en la Fig 3.40.
107
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.3 ELEMENTOS ESTRUCTURALES PARA MODELAJE
z, w
z', w'
x'
Eje medio
x', u' y, v
x, u
Esta técnica es simple de formular, requiere pocos datos y cumple con los criterios de
convergencia, por lo que es una buena aproximación al problema de cascarones. Sin
embargo, tiene ciertos inconvenientes: es un método poco preciso, y requiere el empleo
de gran cantidad de elementos pequeños, tantos más cuanto más curva sea la geometría
del cascarón. Además, el acople entre los efectos de membrana y de flexión se da
únicamente en los nodos y no a nivel elemental, donde estos efectos están desacoplados.
Esto hace que aparezcan momentos flexionantes donde no los hay, por ejemplo, en un
cilindro infinito sometido a presión interna aparecen momentos de flexión en los
elementos, cuando el material está sometido solo a esfuerzo axial.
Otro problema ocurre cuando todos los elementos que llegan a un nodo son coplanares,
ya que en este caso no aparece ninguna rigidez en el giro normal a los elementos 𝜃𝑧 ,
Fig 3.39, y la matriz de rigidez global de la estructura es singular. Para evitar este
problema, se utiliza una rigidez no nula, asociada a este grado de libertad.
108
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
Un primer inconveniente de esta técnica es que involucra una gran cantidad de nodos
en el modelo. Otro problema a tener en cuenta es que, si un elemento es poco grueso,
comparado con las otras dimensiones, su rigidez en dicha dimensión es demasiado
grande, mayor que la real, lo que produce deformaciones incorrectas. Por lo tanto, este
método solo debe utilizarse cuando el cascarón sea grueso y se aproxime más a un
sólido que a un cascarón.
109
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.3 ELEMENTOS ESTRUCTURALES PARA MODELAJE
El término sólido 3D, se utiliza para referirse a problemas tridimensionales en los que no
se hace restricción alguna en lo referente a la geometría, cargas, propiedades del
material, y condiciones de frontera. Debido a esto, se tiene un problema en el que
existen seis componentes de esfuerzo que deben de considerarse en el análisis. De igual
forma, se tienen tres componentes de desplazamiento en cada punto del elemento. Los
elementos más comunes son los tetraedros y los exaedros con tres grados de libertad
por nodo, Fig 3.42.
y
X
Cada de uno de los elementos tiene: 4 nodos y 4 caras, 6 nodos y 5 caras, y 8 nodos y
6 caras, respectivamente, Fig 3.44. Las configuraciones geométricas de estos elementos
corresponden a elementos lineales o de primer orden; derivado de la aproximación
numérica utilizada para definir los desplazamientos dentro del elemento.
111
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.3 ELEMENTOS ESTRUCTURALES PARA MODELAJE
Nodo 4
Nodo 5 Nodo 6
Nodo 1
Nodo 3 Nodo 1
(a) (b)
Nodo 8
Nodo 7
Nodo 5
Nodo 6
Nodo 4
Nodo 3
Nodo 1
Nodo 2
(c)
Para el caso en el que se desea realizar análisis de esfuerzos más refinados y/o
representar geometrías curvas, se puede utilizar estos mismos elementos pero
enriquecidos con nodos en los ejes medios del elemento, Fig 3.45.
112
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
En el contexto del método del elemento finito, los elementos con las características
mostradas en la Fig 3.45, se les denomina elementos sólidos tridimensionales
cuadráticos o de segundo orden; que al igual que para el caso de los elementos de la
Fig 3.44, la aproximación numérica utilizada para definir los desplazamientos da origen
al término “segundo orden”.
Coordenadas
Región Acrónimo* C+E+F† Restricciones
Globales
Tetraedro Tetr 4+6+4 x, y, z Ninguna
Cuña Wedg 6+9+5 x, y, z Ninguna
Pirámide Pyra 5+8+5 x, y, z Ninguna
Wrick Wric 7+11+6 x, y, z Ninguna
Hexaedro Hexa 8+12+6 x, y, z Ninguna
*
Acrónimo puede ser seguido de un número de nodos e.g. Tetr 10, significa un tetraedro de 10 nodos.
†
C: esquinas, E: aristas. F: caras
113
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.3 ELEMENTOS ESTRUCTURALES PARA MODELAJE
3.9.2 Desplazamientos
uz
z
ux uy
y
X
… 𝑁1
𝑢(𝑥, 𝑦, 𝑧) 𝑢1 𝑢2 𝑢𝑛
… 𝑁
𝐮 = [ 𝑣(𝑥, 𝑦, 𝑧) ] = [ 𝑣1 𝑣2 𝑣𝑛 ] [ 2 ] (3.96)
𝑤(𝑥, 𝑦, 𝑧) 𝑤1 𝑤2 … 𝑤𝑛 ⋮
𝑁𝑛
𝐮 = [𝑢1 𝑣1 𝑤1 … … … 𝑢𝑛 𝑣𝑛 𝑤𝑛 ]𝑇 (3.97)
y su tamaño es de tres por el número de nodos del elemento, i.e. número de grados de
libertad.
𝑁1 0 0 … 𝑁𝑛 0 0
𝐍=[0 𝑁1 0 … 0 𝑁𝑛 0] (3.98)
0 0 𝑁1 … 0 0 𝑁𝑛
114
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
3.9.3 Deformaciones
donde
𝛆 tensor simétrico de segundo orden
𝜀x , 𝜀y , 𝜀z coeficientes del tensor simétrico de deformaciones en su
diagonal principal, adimensional
𝜀𝑥𝑦 , 𝜀𝑥𝑧 , 𝜀𝑦𝑧 coeficientes del tensor de deformaciones en dirección angular,
fuera de la diagonal principal, adimensional
𝜕𝑢
𝜕𝑥
𝜕𝑣
𝜕𝑦
𝜕𝑤
𝜕𝑧
𝜺 = 𝜕𝑢 𝜕𝑣 = 𝐁𝐮 (3.100)
+
𝜕𝑦 𝜕𝑥
𝜕𝑢 𝜕𝑤
+
𝜕𝑧 𝜕𝑥
𝜕𝑣 𝜕𝑤
+
[ 𝜕𝑧 𝜕𝑦 ]
𝜕𝑵
0 0
𝜕𝑥
𝜕𝑵
0 0
𝜕𝑦
𝜕𝑵
0 0
𝜕𝑧
𝐁 = 𝜕𝑵 𝜕𝑵 (3.101)
0
𝜕𝑦 𝜕𝑥
𝜕𝑵 𝜕𝑵
0
𝜕𝑧 𝜕𝑥
𝜕𝑵 𝜕𝑵
0
[ 𝜕𝑧 𝜕𝑦 ]
115
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.3 ELEMENTOS ESTRUCTURALES PARA MODELAJE
3.9.4 Esfuerzos
donde
𝝈 tensor de esfuerzos de orden 2
3
donde
33
𝜎𝑥 = 𝜎11
𝜎𝑦 = 𝜎22
𝜎𝑧 = 𝜎33
32 𝜎𝑥𝑧 = 𝜎13
31
23 𝜎𝑦𝑧 = 𝜎23
13
22
12 21
11
x2
x1
Los esfuerzos se relacionan con las deformaciones mediante una ley constitutiva que
define el comportamiento del material. De manera general, esta relación se da como,
o en forma matricial
𝝈 = 𝐂𝛆 (3.104)
116
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
𝜇 𝜇
1 0 0 0
1−𝜇 1−𝜇
𝜈
0 1 0 0 0
1−𝜇
0 0 1 0 0 0
𝐸(1 − 𝜇) 1 − 2𝜇
𝐂= (3.105)
(1 + 𝜇)(1 − 2𝜇) 0 0 0
2(1 − 𝜇)
0 0
1 − 2𝜇
0 0 0 0 0
2(1 − 𝜇)
1 − 2𝜇
0 0 0 0 0
[ 2(1 − 𝜇)]
𝐤 = ∫ 𝐁 𝑇 𝐂𝐁𝑑𝛺 (3.106)
Ω
donde
𝐤 matriz local de rigidez
La matriz de rigidez total se obtiene al integrar en todo el volumen del elemento. Para
fines de cómputo de la matriz de rigidez, generalmente se utiliza un método numérico
para evaluar la integral.
117
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.3 ELEMENTOS ESTRUCTURALES PARA MODELAJE
Coordenada
circunferencial
r
Sección
transversal
meridional
r, u
donde
θ𝑟 coordenada circunferencial, en rad
118
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
3.10.2 Desplazamientos
Los grados de libertad para este elemento son los desplazamientos. Sin embargo, debido
a la característica de simetría axial del problema, los desplazamientos son:
desplazamiento radial y desplazamiento axial. El correspondiente vector de
desplazamientos es,
𝑢(𝑟, 𝑧)
𝐮={ } (3.107)
𝑣(𝑟, 𝑧)
donde
𝑢(𝑟, 𝑧) coeficientes del vector de desplazamientos radiales de un elemento
axisimétrico, en m
𝑣(𝑟, 𝑧) coeficientes del vector de desplzamientos tangenciqales de un
elemento axisimétrico, en m
3.10.3 Deformaciones
𝜕
0
𝜕𝑟
𝜀𝑟 𝜕
𝜀 0
𝜺𝒂 = { 𝜀 𝑧 } = 𝜕𝑧 {𝑢} = 𝐁𝐮 (3.108)
𝜃 1 𝑣
𝜀𝑟𝑧 0
𝑟
𝜕 𝜕
[𝜕𝑧 𝜕𝑟 ]
donde
𝜀𝑟 , 𝜀𝑧 , 𝜀𝜃 , 𝜀𝑟𝑧 coeficientes del vector axisimétrico de deformaciones
𝜺𝒂 vector de deformaciones de un elemento axisimétrico
𝜕
0
𝜕𝑟
1
0
𝐁= 𝑟 (3.109)
𝜕
0
𝜕𝑧
𝜕 𝜕
[𝜕𝑧 𝜕𝑟 ]
119
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.3 ELEMENTOS ESTRUCTURALES PARA MODELAJE
3.10.4 Esfuerzos
donde
𝜎𝑟 , 𝜎𝑧 , 𝜎𝜃 , 𝜎𝑟𝑧 coeficientes del vector de esfuerzos de un elemento
axisimétrico
𝝈𝒂 vector de esfuerzos de un elemento axisimétrico
z rz
r
rz
r
donde
σ𝑟 esfuerzo radial MPa
σ𝑧 esfuerzo normal, dirección z, en MPa
σ𝜃 esfuerzo circunferencial, en MPa
σ𝑟𝑧 esfuerzo tangencial, en MPa
lf bf lf
Entre los problemas comunes que involucran este tipo de comportamiento discontinuo,
se tienen: interacción de diferentes materiales, interacción suelo estructura, interacción
fluido estructura, modelado del concreto reforzado, modelado de las diferentes capaz de
roca o concreto, cimentaciones; entre otros. De manera general, este problema se
presenta en todos los casos en los que están en contacto dos cuerpos.
𝜉=+1
𝐤 = 𝑏𝑓 ∫ 𝐁 𝑇 𝐂𝐁 𝑑𝜉 (3.111)
𝜉=−1
donde
𝑏𝑓 ancho del elemento en la dirección normal a éste, en m
𝑘𝑛 0
𝐂=[ ] (3.112)
0 𝑘𝑡
donde
𝑘𝑛 , 𝑘𝑡 constantes que representan la rigidez de la interfaz para el modo I y
II, respectivamente, adimensional
121
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.3 ELEMENTOS ESTRUCTURALES PARA MODELAJE
1 1
𝑁 = [ (1 − 𝜉) (1 + 𝜉)] (3.113)
2 2
donde
𝑁 función de forma
donde
𝑙𝑓 longitud del elemento interfaz, en m
122
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
4. TIPOS DE ANÁLISIS
123
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.5.1.4 TIPOS DE ANÁLISIS
Tiende a
infinito
Fuerza P, en N
Comportamiento
fundamental lineal
124
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
Para poder aplicar un análisis estático a una estructura se debe cumplir los siguientes
incisos:
Estas consideraciones no solo son incongruentes con la física del problema, sino también
se contradicen unas con otras. Por ejemplo, si la estructura es capaz de soportar
cualquier nivel de carga, es de esperarse que las deformaciones aumenten hasta niveles
en los que ya no son pequeñas. Esto es posible siempre y cuando la estructura se
comporte como cuerpo rígido y no elástico; lo que contradice la primera consideración.
Esta observación hace notar la necesidad de un límite en el rango de aplicación de un
análisis estático lineal.
A pesar de esta limitante, los modelos lineales pueden brindar información muy útil en
el análisis de estructuras que pueden presentar comportamientos no lineales después
de ciertos niveles de carga, como se observa en la Fig 4.2. El comportamiento
inicialmente elástico y lineal de una estructura, representa el nivel de servicio en el que
se requiere que la estructura trabaje, sin sufrir efectos inelásticos.
Punto
crítico
Comportamiento
secundario
Fuerza P, en N
Comportamiento
primario
125
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.5.1.4 TIPOS DE ANÁLISIS
4
P
q
Nodo i
Nodo A Nodo C
1 3
126
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
Nodo D
ME4D
4
MEi4
ME1A MEC3
Nodo i
1 3
Nodo A Nodo C
ME1i MEi3
MEi2
2
MEB2
Nodo B
donde
𝑀𝐸𝐴1 momento de empotramiento perfecto en el nodo A del elemento
1, en N-m
𝑀𝐸𝐵2 momento de empotramiento perfecto en el nodo B del elemento
2, en N-m
𝑀𝐸𝐶3 momento de empotramiento perfecto en el nodo C del elemento
3, en N-m
𝑀𝐸𝐷4 momento de empotramiento perfecto en el nodo D del elemento
4, en N-m
𝑀𝐸𝑖1 , 𝑀𝐸𝑖2 , 𝑀𝐸𝑖3 , 𝑀𝐸𝑖4 momentos de empotramiento perfecto en el nodo i de los
elementos 1, 2, 3 y 4, en N-m
𝑛𝑏
donde
𝑀𝐸𝑖𝑛 momento de empotramiento del elemento 𝑛 en el nodo 𝑖, en N-m
𝑀𝑑𝑖 momento total de desequilibrio en el nodo 𝑖, en N-m
𝑛𝑏 número de barras
127
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.5.1.4 TIPOS DE ANÁLISIS
MD4D
MD1i Elemento
deformado
2
MDB2
donde
𝑀𝐷𝐴1 momento de distribución en el nodo 𝐴 debido a la rigidez
rotacional del elemento 1, en N-m
𝑀𝐷𝐵2 momento de distribución en el nodo 𝐵 debido a la rigidez
rotacional del elemento 2, en N-m
𝑀𝐷𝐶3 momento de distribución en el nodo 𝐶 debido a la rigidez
rotacional del elemento 3, en N-m
𝑀𝐷𝐷4 momento de distribución en el nodo 𝐷 debido a la rigidez
rotacional del elemento 4, en N-m
𝑀𝐷𝑖1 , 𝑀𝐷𝑖2 , 𝑀𝐷𝑖3 , 𝑀𝐷𝑖4 momento de distribución en el nodo 𝑖 debido a la rigidez
rotacional de los elementos 1, 2, 3, 4, en N-m
4𝐸𝐼𝑧
𝑘𝑟𝑛 =
𝑙
(4.3)
3𝐸𝐼𝑧
𝑘𝑟𝑛 =
𝑙
𝑀𝐷1𝑖
𝜃𝑖 = (4.4)
𝑘1𝑟
4𝐸𝐼𝑧
Cuando el extremo opuesto al nodo es empotrado se utiliza ( ), si el extremo es
𝑙
3𝐸𝐼𝑧
articulado se utiliza ( ).
𝑙
128
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
donde
𝑀𝐷𝑖 momento de distribución en el nodo 𝑖, en N-m
𝑀𝐷𝑖𝑛 momento de distribución en el nodo 𝑖 para un elemento n, en N-m
𝑘𝑟𝑛 rigidez rotacional del elemento n, en N/m.
𝑘𝑟1 , 𝑘𝑟2 , 𝑘𝑟3 , 𝑘𝑟4 rigidez rotacional de los elementos 1, 2, 3, 4 en N/m.
𝜃𝑖 rotación del nodo 𝑖, en grados y/o rad
𝐸 módulo de elasticidad, en MPa
𝐼𝑧 momento de inercia de la sección transversal del elemento, en m 4
𝑙 longitud del elemento, en m
Generalizando se llega a:
𝑘𝑟𝑛
𝑀𝐷𝑖𝑛 = 𝑀𝐸𝑖𝑛 (4.5)
∑ 𝑘𝑟𝑛
donde
𝑓𝑘𝑖 distribución de rigidez del nodo i, adimensional
b) Método de Kani. El método de Kani es un proceso iterativo que sirve para calcular
los momentos en los extremos de las barras que conforman estructuras hiperestáticas.
Esta técnica, al igual que el método de Cross, se basa en el método de pendiente—
deflexión, teniendo las siguientes hipótesis:
Considérese la barra 1 de la estructura mostrada en la Fig 4.6. Todas las barras tienen
una rigidez a flexion constante, 𝐸𝐼. Los momentos en los extremos de los elementos son:
4𝐸𝐼𝑧 2𝐸𝐼𝑧
𝑀𝐴1 = 𝑀𝐸𝐴1 + + (4.7)
𝑙 𝑙
2𝐸𝐼𝑧 4𝐸𝐼𝑧
𝑀𝑖1 = 𝑀𝐸𝑖1 + + (4.8)
𝑙 𝑙
donde
𝑀𝐴1 momento flexionante en el nodo A del elemento 1, en N-m
𝑀𝑖1 momento flexionante en el nodo 𝑖 del elemento 1, en N-m
𝑙 longitud de un elemento estructural, barra, viga, marco, cable, en m
130
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
M4D
Mi3
M1A MC3
1 Mi4
Nodo i
Nodo A Nodo C
Mi2 3
M1i Elemento
deformado
2
MB2
donde
𝑀𝐴1 momento flexionante en el nodo A del elemento 1, en N-m
𝑀𝐵2 momento flexionante en el nodo B del elemento 2, en N-m
𝑀𝐶3 momento flexionante en el nodo C del elemento 3, en N-m
𝑀𝐷4 momento flexionante en el nodo D del elemento 4, en N-m
𝑀𝑖1 , 𝑀𝑖2 , 𝑀𝑖3 , 𝑀𝑖4 momentos flexionantes en el nodo 𝑖 de los elementos 1, 2, 3, 4,
en N-m
Fig 4.6 Momentos en los extremos de las barras que conforman la estructura
4𝐸𝐼 3𝐸𝐼
𝑘𝑟1 = ; 𝑘𝑟1 = (4.11)
𝑙 𝑙
donde
𝜃𝐴 rotación nodo A, en grados y/o rad
𝑛𝑏 𝑛𝑏 𝑛𝑏
∑ 𝑀𝐸𝑖 + ∑ 𝑀𝑗 + 2 ∑ 𝑀𝑖 = 0 (4.12)
𝑛=1 𝑛=1 𝑛=1
131
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.5.1.4 TIPOS DE ANÁLISIS
𝑛𝑏 𝑛𝑏
1 𝑘𝑟𝑛
𝑀𝑖 = − ( 𝑛 ) [∑ 𝑀𝐸𝑖 + ∑ 𝑀𝐸𝑗 ] (4.13)
2 ∑ 𝑘𝑟
𝑛=1 𝑛=1
1 𝑘𝑟𝑛
𝑓𝑘𝑖 = − ( 𝑛 ) (4.14)
2 ∑ 𝑘𝑟
donde
𝑁𝐵 𝑁𝐵
𝑗
𝑀𝑖 = 𝑓𝑘𝑖 [∑ 𝑀𝐸𝑛𝑖 + ∑ 𝑀𝐸𝑛 ] (4.15)
𝑛=1 𝑛=1
Los momentos finales en los extremos se calculan utilizando las siguientes expresiones:
donde
𝑛
𝑀𝑓𝑖 momento flexionante final en el nodo 𝑖 del elemento n, en N-m
𝑛
𝑀𝑓𝑗 momento flexionante final en el nodo 𝑗 del elemento n, en N-m
En los siguientes incisos se enumeran y describen los pasos del método de Kani.
132
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
𝐊 = 𝐝𝐟𝑒𝑥𝑡 (4.18)
donde
𝐊 matriz global de rigidez
𝐝 vector global de desplazamientos
𝐟𝑒𝑥𝑡 vector de fuerzas exteriores debido a un incremento de carga
133
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.5.1.4 TIPOS DE ANÁLISIS
(𝑡+Δ𝑡) (𝑡)
𝐊 (𝑡) 𝚫u = 𝐟𝑒𝑥𝑡 − 𝐟𝑖𝑛𝑡 (4.19)
donde
𝐊 (𝑡) matriz global de rigidez para un tiempo determinado
(𝑡+Δ𝑡)
𝐟𝑒𝑥𝑡 vector de fuerza exterior para un n-iesimo tiempo
(𝑡)
𝐟𝑖𝑛𝑡 vector de fuerza interior para un tiempo determinado
𝚫u incremento del vector de desplazamiento
Δ𝑡 incremento de tiempo, en s
𝑡 tiempo, en s
Las fuerzas internas, obtenidas a partir de los desplazamientos del paso de carga
anterior, se dan para un tiempo anterior que las fuerzas externas, lo que da origen a
una fuerza diferente de cero, conocida como residuo, ec 4.20.
(𝑡+Δ𝑡) (𝑡)
𝑅 = 𝐟𝑒𝑥𝑡 − 𝐟𝑖𝑛𝑡 (4.20)
donde
𝑅 residuo
donde
𝐮(𝑡+Δ𝑡) vector de desplazamiento para un n-iesimo tiempo
𝐮(𝑡) vector de desplazamiento para un tiempo determinado
Como se observa en la Fig 4.8, mediante un análisis incremental, como el aquí descrito,
siempre existirá un residuo diferente de cero; dando lugar a una curva denominada
solución numérica, diferente a la solución de equilibrio; no necesariamente la solución
exacta del problema.
134
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
R0
Solución
Incremento de fuerza, en N
numérica
ƒ Solución de equilibrio
u
Incremento de desplazamiento, en m
donde
𝚫u incremento de del vector desplazamiento, en m
𝚫f incremento de fuerza, en N
𝑅 residuo, en N
(𝑛) (𝑛−1)
𝐊 (𝑛−1) 𝛿𝐮(𝑛) = 𝐟𝑒𝑥𝑡 − 𝐟𝑖𝑛𝑡
} 𝑛 = 1,2,3, … (4.22)
𝐮(𝑛) = 𝐮(𝑛−1) + 𝛿𝐮(𝑛)
donde
𝐊 (𝑛−1) matriz global de rigidez para una iteración 𝑛 − 1
𝐮(𝑛−1) vector desplazamiento para una iteración 𝑛 − 1
𝐮(𝑛) vector desplazamiento para una enésima iteración
(𝑛−1)
𝐟𝑒𝑥𝑡 vector de fuerza exterior para una iteración 𝑛 − 1
(𝑛)
𝐟𝑖𝑛𝑡 vector de fuerza interior para una enésima iteración
𝛿𝐮 (𝑛)
corrección del vector desplazamientos para una enésima iteración
𝚫𝐮𝑡 incremento total del vector de desplazamientos para un tiempo
determinado
135
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.5.1.4 TIPOS DE ANÁLISIS
donde
𝚫𝐮𝑡 incremento total del vector de desplazamientos para un tiempo
determinado
Inicio de
incremento
u = 0
Determinación de la nueva
fuerza interna
No
¿Detener
iteración?
Si
Fin de
incremento
136
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
donde
𝐊 (𝑛) matriz global de rigidez para para una enésima iteración
𝐑 𝑑𝑒 vector de fuerza desbalanceada
𝜕𝑅
𝐊𝑇 = (4.25)
𝜕𝚫u
donde
𝐊𝑇 matriz de rigidez tangente
Las variables que se examinan, durante todo el proceso iterativo, para decidir si la
solución converge o no, son: la fuerza, el desplazamiento y la energía. Estos criterios
son los que se encuentran comúnmente en los programas de computadora.
137
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.5.1.4 TIPOS DE ANÁLISIS
el valor de tolerancia por la norma del residuo para el incremento inicial, es decir,
ec 4.26.
donde
𝐑 𝑑𝑒 vector de fuerza desbalanceada
𝐑o vector de fuerza desbalanceada inicial
𝜖 tolerancia
1 Control de carga
2 Control de desplazamiento
3 Control de la longitud de arco
138
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
t
f ext
Incremento de fuerza, en N
3
f ext
ƒ
2
f ext
1
f ext
u1 u2 u3 un
Desplazamiento, en m
donde
1
𝐟𝑒𝑥𝑡 2
, 𝐟𝑒𝑥𝑡 3
, 𝐟𝑒𝑥𝑡 𝑛
, 𝐟𝑒𝑥𝑡 vector de fuerzas exterior para diferentes incrementos de fuerza o
desplazamiento
𝐮1 , 𝐮2 , 𝐮3 , 𝐮𝑛 vector de desplazamiento para diferentes incrementos de fuerza
139
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.5.1.4 TIPOS DE ANÁLISIS
3
f ext
Incremento de fuerza, en N
n
f ext
2
f ext
1
f ext
u
1 2 3 n
u u u u
Incremento de desplazamientoesplazamiento, en m
donde
Δ1𝑢 , Δ2𝑢 , Δ3𝑢 , Δ𝑛𝑢 incremento de del vector desplazamiento 1, 2, 3, n
1
𝐟𝑒𝑥𝑡 2
, 𝐟𝑒𝑥𝑡 3
, 𝐟𝑒𝑥𝑡 𝑛
, 𝐟𝑒𝑥𝑡 vector de fuerzas exterior para diferentes incrementos de fuerza o
desplazamiento
Los métodos descritos anteriormente, para controlar la respuesta, son eficaces para
casos específicos de análisis; sin embargo, existen problemas en los que la no linealidad
es muy fuerte y es necesario tener un control más estable y robusto.
140
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
La carga se aplica en un solo punto de la estructura, lo cual da origen a una curva como
la que se muestra en la Fig 4.13. Este comportamiento incluye un retroceso en la
respuesta de la estructura; conocido como “snap—back”.
Fuerza, en N
Desplazamiento, en m
El control de la longitud de arco puede eliminar los problemas numéricos asociados a los
métodos de control descritos anteriormente. Este método se basa en restringir la norma
de los desplazamientos incrementales a un valor prestablecido. Esto es adaptando en el
proceso iterativo, el tamaño del incremento.
Los métodos y procedimientos para los análisis no lineales, descritos en las secciones
anteriores, se aplican a cualquier tipo de no linealidad presente en los problemas. Sin
embargo, su uso no significa que la respuesta coincida con la solución exacta del
problema.
Las relaciones no lineales entre las deformaciones y los esfuerzos son la causa más
común de un comportamiento no lineal en las estructuras. Existen muchos factores
pueden afectar las propiedades de los materiales (módulo de elasticidad, relación de
poisson, viscosidad, cohesión, fricción, energía de fractura, esfuerzo de fluencia) que
relacionan los esfuerzos con las deformaciones: historia de carga (elastoplasticidad),
condiciones del ambiente (temperatura, presión), y el tiempo en el que la estructura
estará sujeta a las cargas (flujo plástico).
141
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.5.1.4 TIPOS DE ANÁLISIS
4.3.1 Materiales
1 Elasticidad
2 Plasticidad
3 Visco-elasticidad
4 Visco-plasticidad
5 Flujo plástico
6 Daño
Las propiedades plásticas de los materiales varían de acuerdo a las condiciones externas
del problema i.e., temperatura, humedad, tiempo de aplicación de las cargas, entre otros
factores. Entonces, la deformación plástica de diferentes materiales, en condiciones de
temperatura controlada, es independiente del tiempo; mientras que algunos otros que
trabajan a altos niveles de temperatura, experimentan deformaciones plásticas que se
incrementan con el tiempo (flujo plástico). El comportamiento de estos últimos
materiales se asemeja más al de un fluido viscoso. De acuerdo a esto, se tienen modelos
constitutivos de plasticidad dependientes del tiempo y otros que son independientes. El
flujo plástico corresponde a los modelos de plasticidad que dependen del tiempo.
142
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
f Punto B
Punto A
f0
E E
Punto C
e p
donde
𝜎𝑓𝑜 esfuerzo de fluencia inicial, en MPa
𝜎𝑓 esfuerzo de fluencia, en MPa
𝜀𝑝 deformación unitaria plástica, adimensional
𝜀𝑒 deformación unitaria elástica, adimensional
𝐸 módulo de elasticidad, en Mpa
1 Perfectamente plásticos
2 Elasto—plásticos con endurecimiento
3 Elasto—plásticos con ablandamiento
4 Elasto—plásticos ideales
143
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.5.1.4 TIPOS DE ANÁLISIS
Ante este hecho, los modelos de plasticidad no son consistentes con la física del
problema; debido a que en su formulación no se toma en cuenta la disipación de energía
que se da en el proceso de plastificación. Esto se observa fácilmente en el modelo elasto-
plástico ideal. La energía disipada corresponde al área bajo la curva. Aquí es posible
incrementar la carga indefinidamente; lo que conduce a una energía infinita. Este hecho
no corresponde con la física real de problema.
144
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
f o f o f o f o
f o f o
fo fo
ƒ ƒ
a) Curva de (a)
comportamiento
(a) (b)(b)
b) ƒ Curva de comportamiento elastoplástico
ƒ
plástico perfecto plástico
(a)(a) (b) perfecto
(b)
max
max
max
omax
f fo
f o
fo
f
o fo
f o f o
ƒ ƒ max max
(c)(c) (d)(d)
ƒ ƒ max
max
c) curva de comportamiento
(c)(c) con d) curva de comportamiento
(d)
(d) con
endurecimiento por deformación ablandamiento
donde
𝜎𝑓𝑜 esfuerzo de fluencia inicial, en MPa
𝜀𝑓 deformación unitaria asociada al esfuerzo de fluencia inicial, adimensional
𝜎𝑚𝑎𝑥 esfuerzo máximo, en MPa
𝜀𝑚𝑎𝑥 deformación unitaria asociada al esfuerzo máximo, adimensional
Los modelos de daño permiten simular problemas en los que es importante conocer los
mecanismos de falla que desarrollan las estructuras ante condiciones de carga que
provocan su colapso total.
145
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.5.1.4 TIPOS DE ANÁLISIS
Esfuerzo , en MPa
146
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
P
control de
desplazamiento
c
Punto c
Carga, en N
c P
0 Desplazamiento, en cm c
(b)
donde
𝛿𝑐 desplazamiento de un elemento arco, debido a una carga P, en m
Punto de
control de
desplazamiento
P
Esfuerzo en MPa
Módulo de
c eslasticidad, en MPa
P
Deformación adimensional
Se muestra el problema de una viga en voladizo en la Fig 4.19, en la que se asume una
rigidez infinita para la viga; es decir, no se acorta ni alarga y tampoco experimenta
147
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.5.1.4 TIPOS DE ANÁLISIS
donde
𝜑′ rotación de un elemento infinitamente rígido, en rad
Para incluir la no linealidad geométrica en los análisis, es necesario definir dos sistemas
de referencia:
''
Punto a''
y Punto a y''
x x''
donde
𝛺 volumen de un sólido, en m3
𝛺′ ′ volumen de un sólido en un sistema de referencia espacial, en m3
En el contexto del método de los elementos finitos, existen tres tipos de descripciones
utilizadas para este tipo de no linealidad:
1 Lagrangiana total
2 Lagrangiana actualizable
3 Euleriana
Las dos primeras son las más comunes en la mecánica de sólidos; mientras que la tercera
se utiliza por lo general en la mecánica de fluidos o cuando las deformaciones de los
sólidos son muy grandes.
Existen casos, en la no linealidad del material, para los que la diferencia entre las dos
descripciones llega a ser importante. Por ejemplo, los problemas que involucran modelos
de plasticidad y/o daño. Para estos casos, es mejor utilizar una descripción Lagrangiana
actualizable.
La matriz geométrica se obtiene a través de considerar los términos no lineales entre los
desplazamientos y las deformaciones. De esta forma, la matriz total en coordenadas
globales de un elemento, se obtiene como, ec 4.27.
𝐊𝒕 = 𝐊 + 𝐊𝐺 (4.27)
149
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.5.1.4 TIPOS DE ANÁLISIS
donde
𝐊 matriz global de rigidez
𝐊𝐺 matriz global de rigidez geométrica
𝐊𝒕 matriz total en coordenadas globales
f v1
Posición deformada
Nodo 1 del elemento cable f v2
Tc
uv1 Nodo 2
Tc
uv2
Nodo 1 Nodo 2
Tc Tc
donde
u𝑣1 desplazamiento vertical del nodo 1, en m
u𝑣2 desplazamiento vertical del nodo 2, en m
𝑓𝑣1 fuerza vertical generada debido a la tensión del cable en el nodo 1, en N
𝑓𝑣2 fuerza vertical generada debido a la tensión del cable en el nodo 2, en N
𝑇𝑐 tensión inicial en el cable, en N
𝑇𝑐
𝑓𝑣1 = (u − u𝑣2 ) (4.28)
𝑙 𝑣1
150
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
Combinando las ecs 4.28 y 4.29, se pueden expresar las fuerzas verticales en función
de los desplazamientos mediante la relación matricial, ec 4.30.
𝑓𝑣1 𝑇𝑐 1 −1 u𝑣1
[ ]= [ ][ ] (4.30)
𝑓𝑣2 𝑙 −1 1 u𝑣2
La matriz de rigidez local del elemento mostrada en la ec 4.31, se observa que la matriz
geométrica no depende directamente de las propiedades mecánicas del material, y sólo
es función de la longitud del elemento y de la fuerza de tensión.
′ ′
𝑓𝑥1 1 0 −1 0 𝑑𝑥1 0 0 0 0 𝑑𝑥1
′ ′
𝑓𝑦1 𝐴𝐸 𝑑𝑦1 𝑇𝑐 0 −1] 𝑑𝑦1
= [ 0 0 0 0]
′ + [ 1 0
′ (4.31)
𝑓𝑥2 𝑙 −1 0 1 0 𝑑𝑥2 𝑙 0 0 0 0 𝑑𝑥2
[𝑓𝑦2 ] 0 0 0 0 ′
[𝑑𝑦2 ] 0 −1 0 1 [𝑑𝑦2 ]
′
donde
𝐴 área de la sección transversal, en m2
𝑓𝑥1 , 𝑓𝑦1 , 𝑓𝑥2 , 𝑓𝑦2 coeficiente del vector local de fuerza de un elemento cable,
en N
′
𝑑𝑥1 ′
, 𝑑𝑦1 ′
, 𝑑𝑥2 ′
, 𝑑𝑦2 coeficiente del vector local de desplazamiento de un elemento
cable, en m
La matriz geométrica 𝐊 𝐺 está dada como, ecs 4.32. La matriz global de rigidez 𝐊, o
mecánica, del elemento. (ec 4.33)
0 0 0 0
𝑇𝑐 0 1 0 −1 ]
𝐊𝐺 = [ (4.32)
𝑙 0 0 0 0
0 1 0 −1
1 0 −1 0
𝐴𝐸 0 0 0 0]
𝐊= [ (4.33)
𝑙 −1 0 1 0
0 0 0 0
151
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.5.1.4 TIPOS DE ANÁLISIS
z, w
uz1 uz2
ux1 ux2 x, u
donde
𝑢𝑥1 , 𝑢𝑥2 desplazamientos nodales en dirección x de un elemento estructural, en el
nodo 1 y 2, en m
𝑢𝑧1 , 𝑢𝑧2 desplazamientos nodales en dirección z de un elemento estructural, en el
nodo 1 y 2, en m
θ𝑧1 , θ𝑧2 rotaciones nodales sobre el eje 𝑧, en rad y/o grados
′ ′
𝑓𝑦1 12 6𝑙 −12 6𝑙 𝑑𝑦1 36 3𝑙 −36 3𝑙 𝑑𝑦1
𝐸𝐼 ′ 𝑇𝑐 ′
𝑀𝑧1 4𝑙 2 2𝑙 2 ] 𝜃𝑧1 4𝑙 2 −𝑙 2 ] 𝜃𝑧1 (4.34)
= 3 [ 6𝑙 −6𝑙
′ + [ 3𝑙 −3𝑙
′
𝑓𝑦2 𝑙 −12 −6𝑙 12 −6𝑙 𝑑𝑦2 30𝑙 −36 −3𝐿 36 −3𝑙 𝑑𝑦2
[𝑀𝑧2 ] 6𝑙 2𝑙 2 −6𝑙 4𝑙 2 ′
[ 𝜃𝑧2 ] 3𝑙 −𝑙 2 −3𝑙 2
4𝑙 [ 𝜃𝑧2 ]
′
donde
𝑓𝑦1 , 𝑀𝑧1 , 𝑓𝑦2 , 𝑀𝑧2 coeficientes del vector local de fuerzas para un elemento
viga, N y N-m
′
𝑑𝑦1 ′
, 𝜃𝑧1 ′
, 𝑑𝑦2 ′
, 𝜃𝑧2 coeficientes del vector local de desplazamientos de un
elemento viga, en m
36 3𝑙 −36 3𝑙
𝑇𝑐 4𝑙 2 −𝑙 2 ]
𝐤𝐺 = [ 3𝑙 −3𝑙 (4.35)
30𝑙 −36 −3𝐿 36 −3𝑙
3𝑙 −𝑙 2 −3𝑙 4𝑙 2
donde
𝐤𝐺 matriz local de rigidez geométrica
152
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
12 6𝑙 −12 6𝑙
𝐸𝐼 4𝑙 2 2𝑙 2 ]
𝐤 = 3 [ 6𝑙 −6𝑙 (4.36)
𝑙 −12 −6𝑙 12 −6𝑙
6𝑙 2𝑙 2 −6𝑙 4𝑙 2
Por lo tanto, las fuerzas totales que actúan en elemento viga son, ec 4.37
𝐟𝑇 = 𝐟 + 𝐟𝐺 = [𝐤 + 𝐤 𝐺 ]𝐝′ = 𝐤 𝑇 𝐝′ (4.37)
donde
𝐟 vector local de fuerzas
𝐟𝐺 vector local de fuerzas asociadas a la matriz de rigidez geométrica
𝐟𝑇 vector total de fuerzas en coordenadas globales
𝐤𝑇 matriz local de rigidez total en coordenadas globales
𝐝′ vector local de desplazamientos
El uso de la matriz de rigidez es una aproximación general que permite incluir los efectos
secundarios en los análisis estáticos y dinámicos de todo tipo de sistemas estructurales.
Sin embargo, en la ingeniería estructural a este tipo de análisis, comúnmente se le
conoce con el nombre de análisis P-Delta.
Los análisis P-Delta toman en cuenta los efectos secundarios debidos a la no linealidad
geométrica, especialmente para edificios flexibles sujetos a la acción de grandes cargas
laterales o estructuras flexibles que cubren claros largos, como son los puentes, y que
soportan cargas gravitacionales. En este tipo de estructuras, la geometría cambia debido
a las deformaciones, y los efectos de las cargas aplicadas a lo largo de la estructura
deformada. Los tipos de efectos P-Delta, comúnmente utilizados, son dos:
Cuando el peso de la estructura es alto de forma proporcional con su rigidez lateral, los
efectos P-Delta influyen en la amplificación de las fuerzas y los desplazamientos. Bajo
ciertas condiciones, las amplificaciones llegan a ser del orden del 25% o más.
En un modelo de múltiples grados de libertad para un edificio de varios niveles, las masas
están concentradas en cada uno de los entrepisos, como se muestra en la Fig 4.23a. Por
la acción de fuerzas y/o desplazamientos horizontales de las masas. La masa de cada
153
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.5.1.4 TIPOS DE ANÁLISIS
uei Mv
uei Mv
W ei F ei
Piso i Piso i F ei
Piso i W ei Piso i
hei
hei
Piso i -1 Piso i -1
Piso i -1 Piso i -1 F ei-1
F ei-1
W u
Fei F=ei =eiW eieiuei
heihei
(a) (b)
donde
𝑢ei desplazamiento en el entrepiso i, en m
𝑊ei peso del entrepiso i, en N
ℎei altura del entrepiso i, en m
𝐹ei fuerza del entrepiso i debido a los efectos de segundo orden, en N
𝐹ei−1 fuerza del entrepiso i-1 debido a los efectos de segundo orden, en N
𝑀𝑣 momento de volteo adicional debido a los efectos P-Delta, en N-m
Las fuerzas laterales que se muestran en la Fig 4.23b, se obtienen para cada nivel y se
suman al vector de cargas externas, y se modifica el sistema de ecuaciones para el
equilibrio horizontal; quedando (ec 4.38)
𝐊 𝒍 = 𝐅𝒉 + 𝐋𝒌 𝐮 (4.38)
154
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
donde
𝐊𝒍 matriz de rigidez lateral, adimensional
𝐮 vector de desplazamientos, adimensional
𝐅𝒉 vector de fuerzas laterales externas, adimensional
𝐋𝒌 matriz que contiene la relación entre la masa y la altura del entrepiso
Pc
Deformación e, adimensional
donde
𝑃𝑎 carga axial, en N
𝑃𝐶 carga crítica, en N
𝜋 2 𝐸𝐼𝑍
𝑃𝐸 = (4.39)
𝑙2
donde
𝐸𝐼𝑍 rigidez elástica a flexión, en N-m2
𝑃𝐸 carga de Euler, en N
155
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.5.1.4 TIPOS DE ANÁLISIS
a a
Pa Pa
(a) (b)
donde
u desplazamiento, en m
𝜋 2 𝐸𝐼𝑍
𝑃𝐶 = (4.40)
(KL)2
donde
KL longitud efectiva de la columna, en m
156
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
Pa Pa
l /4
l l/2 = KL
l /4
Pa Pa
(a) (b)
𝜋 2 𝐸𝐼𝑧 4𝜋 2 𝐸𝐼𝑧
𝑃𝑐𝑟 = = (4.41)
(𝑙 ⁄2)2 𝑙2
donde
𝑃𝑐𝑟 carga crítica de pandeo, en N
Pa Pa Pa Pa Pa Pa Pa Pa Pa
Pa Pa Pa Pa Pa Pa
157
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.5.1.4 TIPOS DE ANÁLISIS
Se muestra la gráfica correspondiente a cada uno de estos casos, el valor teórico para
el coeficiente crítico de pandeo y valores de diseño recomendados para 𝑘𝑝 , cuando se
tienen condiciones cercanas a las ideales, mostradas en la Fig 4.28.
donde
𝑘𝑝 coeficiente teórico crítico de pandeo
158
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
Este procedimiento se aplica con fines de diseño por lo que se establecen ecuaciones
que especifican las resistencias necesarias para flexión y carga axial, de cada uno de los
elementos que forman la estructura. La resistencia a flexión, ante efectos de segundo
orden y la resistencia axial deben de calcularse por medio de las ecs 4.42 y 4.43.
𝐶𝑚
𝐵1 = ≥1 (4.44)
𝑃𝑟
1 − 𝛼𝑝
𝑃𝑒1
1
𝐵2 = ≥1 (4.45)
𝛼 𝑝 ∑ 𝑃𝑛𝑡
1−
∑ 𝑃𝑒2
donde
𝑀𝑢 resistencia a flexión de segundo orden requerida, en N-m
𝑀𝑛𝑡 momento de primer orden, considerando que no hay traslación
lateral del marco, en N-m
𝑀𝑙𝑡 momento de primer orden, causado únicamente por traslación
lateral del marco, en N-m
𝑃𝑟 resistencia axial de primer orden requerida, en N
𝑃𝑢 resistencia axial de segundo orden requerida, en N
𝑃𝑒1 resistencia elástica crítica de pandeo del elemento en el plano de
flexión, calculada bajo la consideración de que no existe
desplazamiento lateral, en N
𝑃𝑒2 resistencia elástica crítica de pandeo para el entrepiso, calculada a
partir de un análisis de pandeo; considerando que no existen
desplazamientos laterales, en N
𝑃𝑛𝑡 fuerza axial de primer orden, considerando que no hay traslación
lateral del marco, en N
𝑃𝑙𝑡 fuerza axial de primer orden, causada únicamente por traslación
lateral del marco, en N
𝐵1 𝑦 𝐵2 factores de amplificación, adimensional
𝐶𝑚 coeficiente que se determina considerando que no hay traslación
lateral del marco, adimensional
𝛼𝑝 factor de resistencia del método de análisis AISC, adimensional
Para elementos sujetos a cargas transversales entre los apoyos se puede utilizar 𝐶𝑚 = 1.
Para miembros con momentos en los extremos 𝐶𝑚 se estima con la ec 4.46.
159
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.5.1.4 TIPOS DE ANÁLISIS
𝑀1
𝐶𝑚 = 0.60 − 0.40 ( ) (4.46)
𝑀2
donde
𝑀1 𝑦 𝑀2 momentos mínimo y máximo, respectivamente, en los extremos
del elemento y en el plano de flexión considerado, en N-m
Los factores 𝛼 𝑝 se toma igual a 1 para diseño por resistencia última y 1.6 para esfuerzos
permisibles. Para calcular 𝐵1 puede realizarse un análisis de primer orden, por lo lanto.
𝜋 2 𝐸𝐼𝑧
𝑃𝑒1 = (4.48)
(KL)2
Para marcos con carga gravitacional la sumatoria ∑ 𝑃𝑒1 se calcula con la ec 4.49. para
todos los tipos de sistemas resistentes a cargas laterales ∑ 𝑃𝑒2 se define mediante la
ec 4.50.
𝜋 2 𝐸𝐼𝑧
∑ 𝑃𝑒1 = ∑ (4.49)
(KL)2
∑ ℎ𝑖 𝑙
∑ 𝑃𝑒2 = 𝑅𝑀 (4.50)
Δ𝐻
donde
Δ𝐻 distorsión de entrepiso de primer orden debido a las fuerzas laterales.
Su valor varía sobre el área en planta de la estructura por lo que su
valor será la deriva media ponderada en proporción a la carga vertical
o, alternativamente, la deriva máxima, adimensional
𝑅𝑀 factor de arriostramiento, 1.0 para sistemas con marcos arriostrados,
0.85 para marcos rígidos y sistemas combinados, adimensional
160
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
hei
hei
hei
l l
Fig 4.29 Marco plano característico de una estructura sujeta a carga de empuje lateral
Vigas
Columnas
Muros
Constantes
Variables
161
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.5.1.4 TIPOS DE ANÁLISIS
EQ3
Nodo monitoreado
EQ2
EQ1
Detección de la primera
articulación plástica
(a)Primer
(a) Primermecanismo
mecanismo
EQ3
EQ2
Articulaciones plásticas
alcanzadas en el segundo
incremento de carga
EQ1
6. Una vez que se obtiene el primer punto con comportamiento inelástico, los
elementos con secciones dañadas, comúnmente idealizados como una rótula
plástica, se reduce su rigidez por flexión de la estructura, por lo cual se define una
nueva matriz de rigidez del sistema estructural.
162
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
7. Se aplica un segundo incremento del vector de fuerzas hasta que aparecen los
segundos puntos con comportamiento inelástico, Fig 4.30b.
𝐸𝑄1 = Δ𝐸𝑄1
2
𝑑𝑎𝑧 1
= 𝑑𝑎𝑧 1
+ Δ𝑑𝑎𝑧 (4.54)
donde
𝑀𝑜1 momento flexionante de los elementos estructurales, para el primer
incremento de carga de un análisis de empuje lateral, en N-m
𝑀𝑜2 momento flexionante de los elementos estructurales, para el segundo
incremento de carga de un análisis de empuje lateral, en N-m
Δ𝑀𝑜1 primer incremento de momento flexionante, debido al primer
incremento de carga lateral, en N-m
Δ𝑉𝑜1 primer incremento de fuerza cortante, debido al primer incremento de
carga lateral, en N
𝑉𝑜1 fuerza cortante de los elementos estructurales, para el primer
incremento de carga, de un análisis de empuje lateral, en N
𝑉𝑜2 fuerza cortante de los elementos estructurales, para el segundo
incremento de carga, de un análisis de empuje lateral, en N
Δ𝑑𝑎𝑧1
primer incremento de desplazamiento de azotea, debido al primer
incremento de carga lateral, en m
1
𝑑𝑎𝑧 desplazamiento de azotea, debido al primer incremento de carga
lateral, en m
2
𝑑𝑎𝑧 desplazamiento de azotea, debido al segundo incremento de carga
lateral, en m
9. Se aplica un tercer incremento del vector de fuerzas hasta que aparecen los
siguientes puntos con comportamiento inelástico, Fig 4.30b.
163
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.5.1.4 TIPOS DE ANÁLISIS
Desplazamiento de azotea, en m
(a) (b)
Los procedimientos dinámicos no tienen una respuesta única, como sucede para el caso
estático; en vez de ello se tiene un conjunto de soluciones que corresponden a cada uno
de los tiempos de interés. Considere una viga simplemente apoyada sometida a una
carga estática, Fig 4.32a, su respuesta depende únicamente de la carga y esta se puede
obtener mediante la aplicación de la ecuación de equilibrio. Por otro lado, si se le aplica
la misma carga pero dinámicamente, como muestra en la Fig 4.32b, la respuesta de la
viga no solo depende de esta carga sino también, de las fuerzas de inercia que se oponen
al movimiento.
164
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
P P(t)
Fuerzas de inercia
(a) (b)
donde
P(t) fuerza en función del tiempo, en N
P fuerza, en N
Fig 4.32 Viga simplemente apoyada, a) bajo una carga estática, b) bajo una carga
dinámica
La masa total del sistema 𝑚 se concentra en el bloque rígido que está restringido en su
movimiento vertical, y sólo se puede desplazar horizontalmente; entonces, el
movimiento del sistema completo se define mediante el desplazamiento u(𝑡) . La
resistencia al movimiento horizontal se da mediante el resorte elástico 𝑘𝑟𝑜 , y la disipación
de energía se representa mediante el mecanismo de amortiguamiento 𝑐 . La carga
dinámica externa, con variación en el tiempo, es 𝑃(𝑡).
165
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.5.1.4 TIPOS DE ANÁLISIS
u(t)
c
(a)
m P(t )
kro
u(t)
ƒ D(t )
(b) ƒ I (t ) P(t )
ƒ S(t )
donde
u(t) desplazamiento en función del tiempo, en m
𝑐 mecanismo de amortiguamiento, adimensional
𝑚 masa total del sistema, en N
𝑘𝑟𝑜 rigidez de un resorte elástico de un oscilador de un grado de libertad
𝑓𝐼 (𝑡) fuerza de inercia en función del tiempo, en N/m
𝑓𝐷 (𝑡) fuerza de amortiguamiento en función del tiempo, N
𝑓𝑆 (𝑡) fuerza del resorte en función del tiempo, en N/m
𝑃(𝑡) carga dinámica externa en función del tiempo, en N
Cada una de las fuerzas del lado izquierdo de la ec 4.55 están en función del
desplazamiento u(𝑡) o de sus derivadas con respecto al tiempo.
166
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
donde
ü (𝑡) aceleración en función del tiempo de un sistema de un grado de libertad
en, en m/s²
u̇ (𝑡) velocidad en función del tiempo de un sistema de un grado de libertad
en, en m/s
u(𝑡) desplazamiento en función del tiempo de un sistema de un grado de
libertad, en m
167
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.5.1.4 TIPOS DE ANÁLISIS
donde
𝑴 matriz de masa
𝑪 matriz de amortiguamiento
𝑲 matriz de rigidez de la ecuación de equilibrio dinámico
𝑭𝒂 vector de fuerzas externas dinámicas aplicadas a la estructura
𝐔̈ vector de aceleración de la ecuación de equilibrio dinámico
𝐔 ̇ vector de velocidad de la ecuación de equilibrio dinámico
𝐔 vector de desplazamiento de la ecuación de equilibrio dinámico
donde
𝑭𝐼 (𝑡) = 𝑴Ü vectores fuerzas de inercia para un sistema de uno o varios
vectores grados de libertad
𝑭𝐷 (𝑡) = 𝑪U̇ vectores fuerzas de amortiguamiento para un sistema de uno o
varios grados de libertad
𝑭𝐸 (𝑡) = 𝑲U vectores fuerzas elásticas o internas para un sistema de uno o
varios grados de libertad
Existen diferentes métodos paso a paso, pero todos ellos se basan en una división de la
carga y la respuesta de la estructura en el tiempo; dada como una secuencia de
intervalos ∆𝑡, pasos. De esta forma, la solución de la ecuación de movimiento se hace
únicamente para cada intervalo de tiempo, y no para el tiempo total 𝑡.
Los métodos paso a paso más comunes de amplía aplicación en la ingeniería estructural
son:
1 Método de Newmark
2 Método de las diferencias centrales
3 Método de Houbolt
4 Método 𝜃 de Wilson
5 Método de la aceleración promedio
6 Método 𝛼 de Hilber, Hughes y Taylor
168
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
Δ𝑡 2
𝐮𝑡 = 𝐮𝑡−Δ𝑡 + Δ𝑡𝐮̇ 𝑡−Δ𝑡 + 𝐮̈ + βΔ𝑡 3 𝐮̈ (4.62)
2 𝑡−Δ𝑡
(𝐮̈ 𝑡 − 𝐮̈ 𝑡−Δ𝑡 )
𝐮̈ = (4.64)
Δ𝑡
Mediante la sustitución de la ec 4.64 en las ecs 4.63 y 4.62, se obtienen las ecuaciones
estándar de Newmark:
1
𝐮𝑡 = 𝐮𝑡−Δ𝑡 + Δt𝐮̇ 𝑡−Δ𝑡 + ( − β) Δ𝑡 2 𝐮̈ 𝑡−Δ𝑡 + βΔ𝑡 2 𝐮̈ 𝒕 (4.65)
2
donde
𝐮𝑡 vector de desplazamientos con respecto al tiempo
𝐮𝑡−Δ𝑡 vector de desplazamientos con respecto a un incremento de tiempo
anterior
𝐮̇ 𝑡 vector de velocidades con respecto al tiempo
𝐮̇ 𝑡−Δ𝑡 vector de velocidades con respecto a un incremento de tiempo
anterior
𝐮̈ 𝒕 vector de aceleraciones con respecto al tiempo
𝐮̈ 𝑡−Δ𝑡 vector de aceleraciones con respecto a un incremento de tiempo
anterior
Δ𝑡 incremento de tiempo
β, γ parámetros para obtener una solución que sea estable y que
converja.
169
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.5.1.4 TIPOS DE ANÁLISIS
Cáculo inicial
1 1 1
𝑏1 = 𝑏2 = 𝑏3 = β −
β∆𝑡 2 β∆𝑡 2
(4.67)
𝑏4 = γ∆𝑡𝑏1 𝑏5 = 1 + γ∆𝑡𝑏2 𝑏6 = ∆𝑡(1 + γ𝑏3 − γ)
̅ = 𝑲 + 𝑏1 𝑴 + 𝑏4 𝑪
𝐾 (4.68)
̅𝑡 = 𝐋𝐁𝐋𝑇
𝐾 (4.69)
donde
u0 desplazamiento inicial, en s
u̇ 0 velocidad inicial. en m/s
ü 0 aceleración inicial de la gravedad. en m/s²
̅
𝐾 matriz de rigidez efectiva
̅𝑡
𝐾 matriz de rigidez efectiva triangular
𝐋 matriz que contiene operadores diferenciales
𝐁 matriz de transformación
̅𝑡
𝐋𝐮𝐋𝑇 = 𝑭 (4.71)
donde
̅𝑡
𝑭 vector de desplazamiento en el tiempo
170
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
donde
𝑡0 tiempo inicial, en s
𝑡1 tiempo 1, en s
u1 desplazamiento en un tiempo 1, en m
u̇ 1 velocidad inicial en un tiempo 1, en m/s
ü 1 aceleración inicial en un tiempo 1, en m/s²
Δ𝑡
u̇ 1 = u0 + (ü + ü 1 ) (4.73)
2 0
171
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.5.1.4 TIPOS DE ANÁLISIS
Δ𝑡 2 Δ𝑡 2
u1 = u0 + u̇ 0 Δ𝑡 + ü 0 + ü (4.74)
3 6 1
Si se toma a β = 1/6, con γ = 1/2, se obtienen una variación para la velocidad cuadrática
y en los desplazamientos cúbica, Fig 4.35. Este método se le conoce como de aceleración
lineal.
1 1 1
𝛾≥ , 𝛽≤ 𝑦 Δ𝑡 ≤ (4.75)
2 2 𝜔𝑚𝑎𝑥 √𝛾⁄2 − 𝛽
172
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
donde
𝜔𝑛𝑚𝑎𝑥 frecuencia natural máxima de la estructura
1
2𝛽 ≥ 𝛾 ≥ (4.76)
2
Existen otros métodos a partir de utilizar diferentes valores para los parámetros β y γ,
que pertenecen a la familia de métodos de Newmark (ver Tabla 4.2).
El método de análisis modal espectral, es una aproximación que permite calcular los
desplazamientos y esfuerzos máximos inducidos a las estructuras por un sismo. La
principal limitación del método es que no es aplicable para análisis no lineal.
El método modal espectral requiere como dato de partida para su aplicación, conocer los
modos y frecuencias naturales del sistema de múltiples grados de libertad; es decir, que
se conocen los valores de las frecuencias naturales y de los modos de vibrar.
donde
𝑢̈ 𝑆 histograma de las aceleraciones en el tiempo, en m/s2
𝐔 = ∑ 𝐮𝑖 𝚽𝑖 (4.78)
𝑖=1
donde
𝐮𝑖 vector de desplazamiento asociado al modo 𝑖.
𝚽𝑖 matriz de formas modales asociado al modo 𝑖
𝚽𝑖𝑇 𝑲𝚽𝑖 𝐮𝑖 (𝑡) + 𝚽𝑖𝑇 𝑴𝚽𝑖 𝐮̈ 𝑖 (𝑡) + 𝚽𝑖𝑇 𝑪𝚽𝑖 𝐮̇ 𝑖 (𝑡) = −𝑢̈ 𝑆 (𝑡)𝚽𝑖𝑇 𝑴 (4.79)
2
𝜔𝑛𝑎𝑡𝑖 u𝑖 (𝑡) + ü 𝑖 (𝑡) + 2𝜉𝑖 𝜔𝑛𝑎𝑡𝑖 u̇ 𝑖 (𝑡) = −𝑢̈ 𝑆 (𝑡)(Φ𝑖𝑇 𝑀𝑎𝑖 )/𝑀 (4.80)
donde
𝜔𝑛𝑎𝑡𝑖 frecuencia natural de la estructura asociada al modo 𝑖, en rad/s
𝜉𝑖 coeficiente de amortiguamiento asociado al modo 𝑖, adimensional
ü 𝑖 (𝑡) aceleración del modo 𝑖, en m/s2
u̇ 𝑖 (𝑡) velocidad del modo 𝑖, en m/s
𝑀𝑎𝑖 masa asociada al modo 𝑖, N
donde
𝜔𝑜 frecuencia natural para un sistema de un grado de libertad dinámico,
en rad/s
𝜉𝑜 coeficiente de amortiguamiento para un sistema de un grado de
libertad dinámico, adimensional
ü 𝑖 (𝑡) aceleración del modo 𝑖, en m/s2
u̇ 𝑖 (𝑡) velocidad del modo 𝑖, en m/s
uos = Γ𝑖 𝑆𝑑 (4.82)
donde
uos desplazamiento máximo de la coordenada modal, en m
𝑆𝑑 desplazamiento espectral, en m
Γ𝑖 factor de participación del modo 𝑖, adimensional
174
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
Por lo tanto, los desplazamientos relativos máximos asociados con el modo 𝑖 están dados
por la ecuación
donde
Φ𝑖 forma modal para un sistema de un grado de libertad dinámico, en
m
Üos = ü os Φ𝑖 = Γ𝑖 𝑆𝑎 Φ𝑖 (4.84)
donde
Üos vector de pseudo—aceleraciones, adimensional
ü os aceleración máxima de la coordenada modal, en m/s2
𝑆𝑎 seudoaceleración, en m/s²
donde
𝑷𝑒𝑞 vector de cargas equivalente, adimensional
Los valores máximos modales no se suman en forma algebraica sino a través de la suma
cuadrática; es decir que el vector de desplazamientos máximos de todos los modos
puede aproximarse, para cada componente 𝑗, por la ecuación
𝑁
𝑗 2
𝑢max ≈ √∑(𝑢𝑖,𝑚𝑎𝑥 ) (4.86)
𝑖=1
donde
𝑗
umax desplazamientos máximos combinados de todos los modos, en m
𝑢𝑖,𝑚𝑎𝑥 desplazamientos máximos asociado al modo 𝑖, en m
175
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.5.1.4 TIPOS DE ANÁLISIS
𝑁
2
𝐸max ≈ √∑(𝐸𝑖,max ) (4.87)
𝑖=1
donde
𝐸max elementos mecánicos máximos combinados de todos los modos
𝐸𝑖,max elementos mecánicos máximos del modo 𝑖
donde
1
𝑃𝑖,𝑒𝑞 fuerza equivalente del nivel 1, en N
2
𝑃𝑖,𝑒𝑞 fuerza equivalente del nivel 2, en N
3
𝑃𝑖,𝑒𝑞 fuerza equivalente del nivel 3, en N
4
𝑃𝑖,𝑒𝑞 fuerza equivalente nivel 4, en N
𝑅𝑒𝑖 suma de las fuerzas equivalentes de los entrepisos, en N
La fuerza total de cortante acumulada en la base del edificio está dada por la suma
algebraica de las componentes individuales del vector 𝑃𝑖,𝑒𝑞 . La reacción en la base para
el modo 𝑖, expresada en función de los desplazamientos y de la masa de cada piso,
resulta
176
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
𝑁 𝑁 𝑁
𝑗 𝑗 𝑗 𝑗 𝑗
𝑅𝑟𝑖 = ∑ 𝑃𝑖,𝑒𝑞 = ∑ 𝑚𝑒 𝑆𝑎,𝑖 = Γ𝑖 𝑆𝑎,𝑖 ∑ 𝑚𝑒 Φ𝑖 (4.88)
𝑗=1 𝑗=1 𝑗=1
donde
𝑚𝑒 masa del entrepiso, en N
𝑅𝑟𝑖 fuerza total acumulada del edificio
Φ𝑖 forma modal para un sistema de un grado de libertad dinámico, en m
(Φ𝑖𝑇 𝑴𝐁)2
̅𝑖 =
𝑀 (4.89)
(Φ𝑖𝑇 𝑴Φ𝑖 )
donde
̅𝑖
𝑀 masa modal de modo i
Se observa en la ec 4.89 que la masa modal es siempre positiva, y su valor es tal que
la suma de la 𝑀̅𝑖 de todos los modos es igual al 100% de la masa de toda la estructura.
Sobre la base de este concepto, es posible decidir si el número de modos considerados
es suficiente. Normalmente se considera que es suficiente cuando la suma de la masa
modal de los modos incluidos en el análisis alcanza o supera el 90% de la masa de la
estructura.
177
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.5.1.4 TIPOS DE ANÁLISIS
El calor generado
internamente (por
ejemplo, debido a la
corriente eléctrica,
hidratación del cemento)
Temperatura prescrita
Fig 4.37 Sólido continuo con cargas térmicas y condiciones de frontera térmicas
La temperatura que se transfiere a través de las fronteras del continuo, es análoga a las
cargas aplicadas en la superficie; es equivalente a las fuerzas de cuerpo en los análisis
de esfuerzos.
𝐊 𝑇𝐸 𝐓𝑇 = 𝐐𝑻 (4.90)
donde
𝐊 𝑇𝐸 matriz depende de la conductividad del material
𝐓𝑇 vector de temperaturas nodales
𝐐𝑻 vector de cargas térmicas
Los elementos finitos utilizados para análisis de esfuerzos se pueden utilizar para los
análisis térmicos y por tanto el mismo tipo de malla. Una vez conocida la distribución de
temperatura se puede utilizar el mismo modelo para conocer el estado de esfuerzos
dentro del sólido dependientes de la temperatura.
178
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
𝜕𝑻
𝑓ℎ = −𝒌𝑡𝑒 (4.91)
𝜕𝑥
donde
𝑻 temperatura, °C
𝑓ℎ flujo de calor por unidad de área, en J
𝒌𝑡𝑒 conductividad térmica del material, en W/(m-K)
El signo negativo significa que el calor se transfiere en una dirección opuesta a la del
incremento de temperatura.
ƒ hx T ƒ hx T
x x
ƒ hy
(a) Material isotrópico (b) Material anisotrópico
donde
𝑓ℎ𝑦 , 𝑓ℎ𝑥 flujo de calor por unidad de área paralelo al gradiente de temperatura
179
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.5.1.4 TIPOS DE ANÁLISIS
𝐓𝑇1
𝐓
𝑇(𝑥, 𝑦, 𝑧) = [𝑁1 𝑁2 ⋯ 𝑁𝑛 ] { 𝑇2 } o 𝑇 = 𝐍𝐓𝑇 (4.92)
⋮
𝐓𝑇𝑛
donde
𝐍 matriz que contiene las funciones de forma, adimensional
𝑁1 𝑁2 ⋯ 𝑁𝑛 componentes de la matriz de funciones de forma,
adimensional
𝑇 temperatura, en grados centígrados
𝐓𝑇 vector de temperaturas nodales
𝐓𝑇1 , 𝐓𝑇2 , 𝐓𝑇𝑛 coeficientes del vector de temperaturas nodales
𝜕/𝜕𝑥
𝐁={ }𝐍 (4.94)
𝜕/𝜕𝑦
donde
𝐁 matriz de transformación, adimensional
𝐊 𝑇𝐸 = ∫ 𝐁 𝑇 𝒌𝑡𝑒 𝐁 𝑑𝛺 (4.95)
𝛺
donde
𝐊 𝑇𝐸𝐹 matriz de conductividad del elemento finito
180
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
P
P
P
P
181
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.5.1.4 TIPOS DE ANÁLISIS
182
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
5.1 Introducción
183
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
En este subcapítulo se presentan algunas ayudas para que los usuarios del software de
análisis estructural puedan adquirir las bases teóricas de los métodos y puedan modelar
una estructura de forma adecuada. Se ofrecen las bases necesarias para que el analista
realice un correcto mallado de los elementos estructurales, se dan recomendaciones de
como considerar las discontinuidades de un elemento estructural en el mallado y
finalmente la forma de interpretar los resultados de un análisis estructural utilizando un
software comercial de elementos finitos.
Para establecer los objetivos del análisis, con respecto a la utilidad de los resultados, se
deben considerar los siguientes puntos:
185
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.6 ASPECTOS PRÁCTICOS DEL MODELADO POR EL MÉTODO DE LOS ELEMENTOS FINITOS
Para recopilar la información necesaria para construir el modelo numérico, que ofrezca
resultados congruentes con el fenómeno físico, se debe tener presente las
incertidumbres del problema.
Posterior a la solución del modelo numérico definido, se debe verificar y validar la calidad
de los resultados obtenidos.
Para la construcción del modelo numérico se tienen los siguientes cuatro pasos a seguir:
a) Los elementos finitos de orden mayor ofrecen mejores resultados que elementos
finitos de orden menor.
b) Es preferible utilizar mallas finas con elementos finitos simples, evitando el uso
mallas gruesas con elementos finitos de orden mayor.
c) Para el análisis de esfuerzos de un elemento estructural se sugiere utilizar
elementos finitos rectangulares.
186
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
En los siguientes puntos se describen varios factores significativos para definir el modelo
matemático de un problema físico.
y
x y
x
(a) Malla 1 (565 grados de libertad) (b) Malla 2 (3155 grados de libertad)
187
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.6 ASPECTOS PRÁCTICOS DEL MODELADO POR EL MÉTODO DE LOS ELEMENTOS FINITOS
aef
bef
oblicuo
< 120°
> 60°
bmef < 5a
donde
𝑏𝑚𝑒𝑓 base mayor de elemento finito trapezoidal, en m
𝑎𝑚𝑒𝑓 base menor de elemento finito trapezoidal, en m
188
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
< 120°
5. Posición del nodo medio. Cuando se utilizan elementos con nodos medios es
necesario localizarlos en el centro. Se debe evitar colocar nodos en los cuartos,
como se muestra en la Fig 6.6, ya que provocan esfuerzos singulares.
189
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.6 ASPECTOS PRÁCTICOS DEL MODELADO POR EL MÉTODO DE LOS ELEMENTOS FINITOS
aef
nodo cuarto
no válido
> aef /4
> bef /4
bef
En la Fig 6.7 se muestra la transición entre diferentes elementos finitos, llevando a una
incompatibilidad en la cinemática de desplazamientos del problema. En este ejemplo se
utilizaron elementos finitos lineales y cuadráticos, cuya variación de los desplazamientos
es lineal y cuadrática, respectivamente, lo que provoca que entre los nodos se tengan
huecos.
cuadrático lineal
Es recomendable construir un modelo con elementos finitos simples, y del mismo orden,
que puedan representar el problema físico que se desea resolver. Sin embargo, existen
190
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
Cuando se requiere realizar transiciones de malla como las mostradas en la Fig 6.8, es
preferible utilizar elementos lineales simples.
Simetría del problema. Muchas estructuras se caracterizan por poseer algún tipo de
simetría o anti—simetría debido a su geometría, condiciones de apoyo, cargas y
materiales. En la Fig 6.9 se muestran diferentes tipos de simetría en las estructuras. Un
análisis previo del problema, puede ayudar a identificar estas condiciones; permitiendo
que sólo una parte del problema se analice. El objetivo de esto, es eliminar costo
computacional innecesario en su análisis.
191
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.6 ASPECTOS PRÁCTICOS DEL MODELADO POR EL MÉTODO DE LOS ELEMENTOS FINITOS
Los resultados obtenidos mediante estás simplificaciones, son exactamente los mismos
a los que se obtendrían si se analiza el problema completo. Sin embargo, se debe tener
presente que, para poder realizar un análisis de la estructura bajo estos términos, es
necesario que se tenga simetría en todos y cada uno de los aspectos que conforman el
problema completo. En los siguientes puntos se describen los diferentes tipos de
simetría:
a) Simetría con respecto a un plano: Este es el caso más frecuente que se puede
encontrar en los problemas de la ingeniería estructural, se obtiene haciendo un efecto
de espejo a una parte del problema. En la Fig 6.10 se muestra un elemento que posee
este tipo de simetría.
192
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
u1
y
u2
u3
De manera general, si se tiene una malla en la que los nodos en el eje de simetría son
los que se muestran en la Fig 6.11, los apoyos que se deben imponer en estos puntos
deben de cumplir que el desplazamiento horizontal sea cero, es decir
u1 = u2 = u3 = 0 (6.1)
donde
u1 desplzamiento en el nodo 1, en m
u2 desplzamiento en el nodo 2, en m
u3 desplzamiento en el nodo 3, en m
193
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.6 ASPECTOS PRÁCTICOS DEL MODELADO POR EL MÉTODO DE LOS ELEMENTOS FINITOS
La viga de la Fig 6.12, presenta simetría con respecto a todas las variables del problema,
por lo que a la mitad de su longitud se coloca un apoyo como el que se muestra, en
donde se impide el desplazamiento horizontal y el giro en ese punto es cero.
P P
3 1 1 3
8 l 8 l 8 l 8 l
P
Desplazamiento = libre
Rotación = 0
donde
P fuerza, en m
𝑙 longitud de la viga
Para el caso de la viga de la Fig 6.13, se tiene anti—simetría con respecto a las cargas
pero simetría de material, apoyos y geometría. Debido a esto, se coloca un apoyo en el
punto medio de la viga en el que se impide el desplazamiento vertical y horizontal, pero
se permite el giro.
194
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
P P
3 1 1 3
8 l 8 l 8 l 8 l
P
Desplazamiento = libre
Rotación = 0
A continuación, se dan algunas reglas que pueden ayudar a definir las condiciones de
frontera, apoyos, en problemas con simetría y anti—simetría de cargas (Ver
Tablas 6.1 y 6.2).
b) Simetría axial: Se dice que una estructura o un sólido tiene simetría axial, cuando
el problema completo se puede generar mediante una rotación de 360 grados de un
plano, alrededor de un eje axial. Entonces, un problema 3D se puede simplificar a uno
en 1D o 2D. En la Fig 6.14 se muestran dos problemas de sólidos en los que se simplifican
los modelos de 3D a 2D.
195
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.6 ASPECTOS PRÁCTICOS DEL MODELADO POR EL MÉTODO DE LOS ELEMENTOS FINITOS
Las cargas aplicadas a un problema con simetría axial, no necesariamente deben ser
simétricas o anti—simétricas. Cualquier tipo de carga con asimetría axial se puede
expresar por medio de una superposición de cargas con ayuda de las series de Fourier.
c) Simetría cíclica: Estructuras en las que las cargas y la geometría se repiten con
un cierto patrón, presentan una simetría cíclica. En estos casos, el comportamiento de
las estructuras se obtiene analizando únicamente un sector de ellas; imponiendo las
condiciones apropiadas en las fronteras del sector a analizar, Fig 6.15.
P
unA utA
Lado A
P unB P
utB
P
Lado B
Célula representativa
donde
u𝑛𝐴 desplazamiento normal del lado A, en m
u𝑛𝐵 desplazamiento normal del lado B, en m
u𝑡𝐴 desplazamiento tangencial del lado A, en m
u𝑡𝐵 desplazamiento tangencial del lado B, en m
196
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
Las condiciones de frontera para la estructura mostrada en la Fig 6.15 deben de ser tal
que, todas las variables en el lado “A”, del sector a analizar, sean cinemáticamente
compatibles con las del lado “B”.
Las variables involucradas en las ecs 6.2 y 6.3 son incógnitas en el problema; por lo que
estás ecuaciones son restricciones que se deben incluir en su solución.
donde
u𝐴𝑥 desplazamiento normal del lado A en dirección x, en m
u𝐵𝑥 desplazamiento normal del lado B en dirección x, en m
197
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.6 ASPECTOS PRÁCTICOS DEL MODELADO POR EL MÉTODO DE LOS ELEMENTOS FINITOS
P
u Ax
Punto a
Punto b
uBx
Célula representativa
Otro ejemplo de una estructura, con simetría repetitiva, se muestra en la Fig. 6.17.
198
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
a) Geométricas
b) De cargas
c) Por las condiciones de apoyo
d) Debidas al material
La barra mostrada en la Fig 6.18a, sujeta a fuerza axial tiene un cambio en su sección
transversal; dando origen a una discontinuidad. De igual forma sucede con la placa que
se muestra en la Fig. 6.18b, que experimenta un cambio en su espesor. La malla de
elementos finitos debe incluir nodos en estos puntos para delimitar este cambio de
geometría.
(a)
Linea de nodos
t p1
t p2
(b)
donde
𝑡𝑝1 espesor de la placa 1, en m
𝑡𝑝2 espesor de la placa 2, en m
Fig 6.18 a) Barra con variación de área, b) Placa con diferentes espesores
199
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.6 ASPECTOS PRÁCTICOS DEL MODELADO POR EL MÉTODO DE LOS ELEMENTOS FINITOS
200
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
donde
P fuerza, en N
201
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.6 ASPECTOS PRÁCTICOS DEL MODELADO POR EL MÉTODO DE LOS ELEMENTOS FINITOS
Discontinuidad cuando existe un cambio de las propiedades del material. Ver Fig 6.22.
La definición de los objetivos del análisis se debe hacer con base al conocimiento general
del problema a resolver; tomando en cuenta los siguientes puntos:
202
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
1. Tipo de análisis
Estático
Modal
Armónico
Dinámico transitorio
Espectral
De pandeo
Dinámico explícito
Sólidos 1D
Sólidos en 2D
o Esfuerzos planos
o Deformaciones planas
o Axisimétricos
Flexión
o Vigas
o Placas
o Membranas
o Cascarones
Sólidos en 3D
Como regla general se sugiere utilizar modelos en dos dimensiones; si las condiciones
del problema lo permiten. Si se tiene una geometría compleja, se puede utilizar modelos
en tres dimensiones.
203
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.6 ASPECTOS PRÁCTICOS DEL MODELADO POR EL MÉTODO DE LOS ELEMENTOS FINITOS
204
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
6.3.2 Verificación
El proceso de verificación consiste en revisar que los datos introducidos al modelo sean los
correctos, que el programa de computadora está trabajando adecuadamente y que los
errores en los resultados del análisis estén dentro de los límites permitidos.
Por medio de este proceso se evalúa si el modelo numérico definido es congruente con el
problema físico a analizar. Esto se puede hacer en tres etapas:
1. Discretización de la malla
2. Selección del tipo de elemento a utilizar
3. Sistema de unidades a utilizar
4. Definición de las condiciones de frontera
Etapa 3. Revisar la congruencia de los resultados, los cuales deben verificarse conforme
a:
1. Condiciones de equilibrio
2. Exactitud de los resultados obtenidos de acuerdo al modelo numérico utilizado. Tipo
de elemento, condiciones de frontera, localización de nodos y cargas.
6.3.3 Validación
Es el proceso por medio del cual se establece si el modelo numérico adoptado para
analizar un fenómeno físico reproduce satisfactoriamente todas las características de
éste.
1 Resultados experimentales
2 Soluciones exactas del problema
3 Análisis simplificados
4 Modelos similares, previamente validados
205
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.6 ASPECTOS PRÁCTICOS DEL MODELADO POR EL MÉTODO DE LOS ELEMENTOS FINITOS
5 Datos históricos
No ¿Se reunieron Si
Falla Pasa
los criterios?
Los principales resultados de un análisis mediante elemento finito son dados en los
puntos nodales de la malla, estos pueden ser: desplazamientos, temperaturas. u otras
cantidades basadas en las variables nodales. Los esfuerzos obtenidos por el programa,
generalmente se dan en los puntos de integración de Gauss del elemento; sin embargo,
algunos programas estiman la proyección de estos esfuerzos en los nodos.
Existen dos formas mediante las que el programa de computadora reporta los resultados
obtenidos del análisis numérico y gráfico. Los resultados gráficos son más fáciles de
interpretar que los numéricos; pero no tienen mucha precisión. Las gráficas de los
desplazamientos y esfuerzos indican el comportamiento de la estructura analizada.
Los criterios para graficar los resultados obtenidos de un modelo numérico, ayudan a
interpretar de manera visual los resultados obtenidos de un programa de elementos
finitos. Se muestra esquemáticamente un método para visualizar los resultados en la
Fig 6.25.
Si
Si
Si
207
C.1.5 ANÁLISIS DE ESTRUCTURAS
C.1.5.6 ASPECTOS PRÁCTICOS DEL MODELADO POR EL MÉTODO DE LOS ELEMENTOS FINITOS
Debido a que los desplazamientos son de magnitud muy pequeña, en comparación con
las dimensiones del problema, es necesario dar un factor de amplificación para las
gráficas deformadas que se utilizarán en la animación.
208
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
son, estos campos son continuos dentro de los elementos, pero discontinuos en las
fronteras de los mismos.
El graficar el estado de esfuerzos de todo el problema permite ver las partes donde se
espera tener mayor concentración de estos. Sin embargo, debe tenerse en cuenta que
el programa presenta una gráfica alisada y no necesariamente corresponde a la
distribución de esfuerzos reales. En la Fig 6.27a se muestra una gráfica en la que el
programa dibuja el estado de esfuerzos en el sólido como si fuera un campo continuo;
mientras que en la Fig 6.27b se muestra la distribución de esfuerzos como son en
realidad, discontinuos entre elementos.
a) b)
209
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
REFERENCIAS
Adams, V y Askenazi, A (1998) Building Better Products with Finite Element Analysis,
Ed. OnWord Press.
Atkinson, J (2007) The Mechanics of Soils and Foundations, Ed. Taylor & Francis.
Bazant, Z P y Planas, J (1998) Fracture and Size Effect in Concrete and Other Quasibrittle
Materials, Ed. CRC Press.
Belytschko, T y Black, T (1999) “Elastic crack growth in finite elements with minimal
remeshing”, International Journal for Numerical Methods in Engineering, vol. 45, pág.
601—620.
Bhavikatti, S S (2005) Finite Element Analysis, Ed. New age international publishers.
Broughton, P y Ndumbaro P (1994) The Analysis of Cable and Catenary structures, Ed.
Thomas Telford.
211
C.1.5 ANÁLISIS DE ESTRUCTURAS
REFERENCIAS
Carroll, W F (1999) A Primer for Finite Elements in Elastic Structures, Ed. John Wiley &
Sons.
Celigueta J T (2008) Método de los Elementos Finitos para Análisis Estructural, Escuela
Superior de Ingenieros de San Sebastián, Universidad de Navarra, España.
Cook R D (1995) Finite Element Modeling for Stress Analysis, Ed. John Wiley & Sons.
Crisfield, M A (1991) Non—Linear Finite Element Analysis of Solids and Structures. Vol.
1: Essentials, Ed. John Wiley & Sons.
Dolbow, J E (1999) An Extended Finite Element Method with Discontinuous Enrichment for
Applied Mechanics, Tesis de doctorado, Northwestern University, EUA.
Doyle, J F (1997) Wave Propagation in Structures: Spectral Analysis Using Fast Discrete
Fourier Transforms, Ed. Springer, Segunda edición.
212
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
Fish, J y Belytschko, T (2007) A first Course in Finite Elements, Ed. John Wiley & Sons.
Galambos, T V (1998) Stability Design Criteria for Metals Structures, Ed. John Wiley &
Sons.
Hartmann, F y Katz, C (2007) Structural Analysis with Finite Elements, Ed. Springer—
Verlag.
Heyman, J (1996) Elements of the Theory of Structures, Ed. Cambridge University Press.
IStructE (2002) Guidelines for the Use of Computers for Engineering Calculations,
Institute of Structural Engineers.
Jirásek, M y Bazant, Z P (2002) Inelastic Analysis of Structures, Ed. John Wiley & Sons.
213
C.1.5 ANÁLISIS DE ESTRUCTURAS
REFERENCIAS
Langhaar, H L (1932) Energy Methods in Applied Mechanics, Ed. John Wiley & Sons.
Lee, U; Kim, J. y Leung A Y (2000) The Spectral Element Method in Structural Dynamics,
The Shock and Vibration Digest, vol. 32, pág. 451—465.
Lee, U (2009) Spectral element method in structural dynamics, Ed. John Wiley & Sons.
Lepi, S M (1998) Practical Guide to Finite Elements. A solid mechanics approach, Ed.
Marcel Dekker.
Liu, G R (2003) Mesh Free Methods. Moving beyond the finite element method, Ed. CRC
Press.
Liu, G R y Quek, S S (2003) The Finite Element Method. A practical course, Ed. Elsevier
BH.
Logan, D L (2007) A First Course in the Finite Element Method, Ed. Thomson.
MacLeod, I A (1995) A Strategy for the Use of Computers in Structural Engineering, The
structural Engineer, 73:1, 366—370.
MacLeod, I A (2005) Modern Structural Analysis. Modelling process and guidance, Ed.
Thomas Telford Publishing.
Manual de usuario, DIANA (2008) Analysis Procedures, Ed. TNO DIANA BV.
214
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
Melenk, J M y Babuska, I (1996) “The partition of unity finite element method: basic theory
and applications”, Computer Methods in Applied Mechanics and Engineering, vol. 39, pág.
289—314.
Mohammadi, S (2008) Extenden Finite Element Method for Fracture Analysis of Structures,
Ed. Blackwell Publishing.
Morris, A (2008) A Practical Guide to Reliable Finite Element Modelling, Ed. John Wiley
& Sons.
NAFEMS, http://www.nafems.org.
Oñate, E (2009 y 2013) Structural Analysis with the Finite Element Method. Vol. 1 y 2, Ed.
Springer.
Schueller, W (2008) Building Supports Structures. Analysis and Design with SAP2000
software, Ed. Computer and Structures.
Segerlind, L J (1984) Applied Finite Element Analysis, Ed. John Wiley & Sons.
215
C.1.5 ANÁLISIS DE ESTRUCTURAS
REFERENCIAS
Ugural, A C y Fenster, S K (2003) Advanced Strength and Applied Elasticity, Ed. Prentice—
Hall.
Ventsel, E y Krauthammer, T (2001) Thin Plates and Shells, Ed. Marcel Dekker.
Washizu, K (1975) Variational Methods in Elasticity and Plasticity, Ed. Pergamon Press.
Wallerstein, D V (2002) A Variational Approach to Structural Analysis, Ed. John Wiley &
Sons.
Wilson, E (1911) CAL—91. Computer Assisted Learning of Static and Dynamic Analysis
of Structural Systems, Report no. UCB/SEMM—91/01, University of California at
Berkeley, C.A., U.S.A.
Ye, J (2008) Structural and Stress Analysis. Theories, tutorials and examples, Ed. Taylor
& Francis Group.
Zienkiewicz, O C; Taylor, R L y Zhu, J Z (2005) The Finite Element Method. Its basis and
fundamentals, Ed. Elsevier BH.
Zienkiewicz, O C y Taylor, R L. (2005) The finite element method for solid and structural
mechanics, Ed. Elsevier BH.
Zienkiewicz, O C y Taylor R L (2005) The finite element method. Its basis and
fundamentals, Ed. Elsevier BH.
216
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
AYUDAS DE DISEÑO
Una viga de sección variable sujeto a una carga distribuida se muestra en Fig Ad.1. Por
simplicidad, sin limitar la validez del procedimiento mostrado, se considera que la sección
transversal es rectangular y que la única dimensión que varía a lo largo de su longitud es
el peralte, es decir la base se mantiene constante.
y
q
x
Nodo a
Nodo b hv(x
bv
l1 l2 l3
donde
𝑏𝑣 ancho de la sección transversal de una viga, en m
ℎ𝑣 (𝑥) altura de la sección transversal de una viga de sección variable a una distancia
𝑥, en m
𝑙1 , 𝑙2 , 𝑙3 diferentes longitudes de tramos de una viga de sección variable, en m
𝑞 carga uniformente distribuida, en N/m
1 𝑙 N2 1 𝑙 M𝑧 2 1 𝑙 V𝑦 2
Π= ∫ 𝑑𝑥 + ∫ 𝑑𝑥 + ∫ 𝑑𝑥 (Ad.1)
2 0 𝐸𝐴 2 0 𝐸𝐼𝑧 2 0 𝐺𝐴𝑐
donde
Π energía potencial de cuerpo, en N-m
N fuerza normal, en N
M𝑧 momento flexionante alrededor del eje z, en N-m
V𝑦 cortante en dirección 𝑦, en N
𝐸 módulo de elasticidad, en Mpa
𝐺 módulo de cortante, en N/m
𝐴 área de la sección transversal del elemento, en m2
𝐴𝑐 área efectiva de cortante, en m2
𝐼𝑧 momento de inercia de la sección transversal de un elemento, en m4
𝑙 longitud de un elemento barra, viga marco, cable, en m
𝑙
1 𝑙1 N 2 1 𝑙1 M𝑧 2 1 𝑙1 V𝑦 2
∫ 𝑑𝑥 = ∫ 𝑑𝑥 + ∫ 𝑑𝑥 + ∫ 𝑑𝑥 +
0 2 0 𝐸𝐴 2 0 𝐸𝐼𝑧 2 0 𝐺𝐴𝑐
Para mostrar el procedimiento de cálculo se aplica a una viga doblemente empotrada con
una carga uniformemente distribuida. La sección transversal es rectangular con ancho
constante y peralte variable, como se observa en la Fig Ad.2.
218
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
y
q
x
hva Nodo a
bv
donde
ℎ𝑣𝑎 altura de la sección transversal de una viga de sección variable en el nodo a,
en m
ℎ𝑣𝑏 altura de la sección transversal de una viga de sección variable en el nodo b,
en m
ℎ𝑣 (𝑥) altura de la sección transversal de una viga de sección variable a una distancia
x, en m
𝐴𝑣a
𝐴(𝑥) = (𝑙 + 𝑥) (Ad.3)
𝑙
𝑥 3
𝐼𝑧 (𝑥) = 𝐼𝑧 (1 + ) (Ad.4)
𝑙
donde
𝐴(𝑥) área de la sección transversal de una viga de sección variable a una
distancia x, en m2
𝐴𝑣a área de la sección transversal de una viga de sección variable en el
nodo a, en m2
𝐼𝑧 (𝑥) momento de inercia de la sección transversal de una viga a una
distancia x, en m4
Para aplicar el método de las fuerzas es necesario definir una estructura isostática en la
que se han eliminado las fuerzas redundantes. Para este caso, se elimina el apoyo A con lo
que la estructura isostática es la que se muestra en la Fig Ad.3.
219
C.1.5 ANÁLISIS DE ESTRUCTURAS
AYUDAS DE DISEÑO
Nodo a Nodo b
l
M
donde
N fuerza normal, en N
V cortante, en N
M momento flexionante, en t-m
1 𝑙 N2 𝑙 𝑙
N2
Π𝑁 = ∫ 𝑑𝑥 = ∫ 𝑑𝑥 (Ad.5)
2 0 𝐸𝐴 2𝐸𝐴 0 (𝑙 + 𝑥)
donde
Π𝑁 energía de deformación complementaria, debido al efecto de la fuerza
normal, en N-m
N(𝑥) = Na (Ad.6)
donde
N(𝑥) fuerza normal de un elemento estructural a una distancia x, en N
Na fuerza normal de un elemento estructural en el nodo a, en N
𝑙 𝛼ℎ 𝑙
𝐟𝐟 = log (1 + ) (Ad.7)
𝐸𝐴(𝑥) ℎ𝑣 (𝑥)
donde
𝐟𝐟 matriz de flexibilidades
𝛼ℎ factor que relaciona los peraltes de una vida de sección transversal
variable, adimensional
220
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
𝐟𝐟 −1 𝐟𝐟 −1 𝐁 T
𝐤=[ ] (Ad.8)
𝐁𝐟𝐟 −1 𝐁𝐟𝐟 −1 𝐁 T
donde
𝐁 matriz de transformación
k matriz local de rigidez
𝐁 = −1 (Ad.9)
𝐸𝐴(𝑥) 1 1 −1
𝐤= [ ] (Ad.10)
𝑙 𝛼ℎ 𝑙 −1 1
log (1 + )
ℎ𝑣 (𝑥)
ℎ𝑣𝑎 − ℎ𝑣𝑏
𝛼ℎ = (Ad.11)
𝑙
𝛼ℎ 𝑙
log [1 + ]
ℎ𝑣 (𝑥)
𝐟𝐟 = (Ad.12)
𝑏𝑣 𝐸𝛼ℎ
𝛼ℎ 𝑙
log [1 + ] 𝑙
ℎ𝑣𝑎
lim = (Ad.13)
𝛼ℎ →0 𝑏𝑣 𝐸𝛼ℎ 𝐸𝐴𝑣a
𝑓11 0 0 0 0 0
0 𝑓22 0 0 0 𝑓26
0 0 𝑓33 0 −𝑓35 0
[𝐟𝐟 ] = (Ad.14)
0 0 0 𝑓44 0 0
0 0 −𝑓53 0 𝑓55 0
[0 𝑓62 0 0 0 𝑓66 ]
𝑙
𝑑𝑧
𝑓11 = ∫ (Ad.15)
0 𝐸 𝐴𝑣a
221
C.1.5 ANÁLISIS DE ESTRUCTURAS
AYUDAS DE DISEÑO
𝑙 2 𝑙
𝑧 𝑑𝑧 𝑑𝑧
𝑓22 = ∫ +∫ (Ad.16)
0 𝐸𝐼𝑧 0 𝐺 𝐴𝑐
𝑙
𝑧𝑑𝑧
𝑓26 = ∫ (Ad.17)
0 𝐸𝐼𝑧
𝑙 2 𝑙
𝑧 𝑑𝑧 𝑑𝑧
𝑓33 = ∫ +∫ (Ad.18)
0 𝐸𝐼𝑧 𝐺
0 𝐴𝑐
𝑙
𝑧𝑑𝑧
𝑓35 = ∫ (Ad.19)
0 𝐸𝐼𝑧
𝑙
𝑑𝑧
𝑓44 = ∫ (Ad.20)
0 𝐸𝐼𝑧
𝑙
𝑑𝑧
𝑓55 = ∫ (Ad.21)
0 𝐸𝐼𝑧
𝑙
𝑑𝑧
𝑓66 = ∫ (Ad.22)
0 𝐸𝐼𝑧
donde
𝑓11 , 𝑓12 , 𝑓13 , 𝑓14 , 𝑓15 … 𝑓66 coeficientes de la matriz de flexibilidades
donde
ΠM,V energía de deformación complementaria, debido al efecto del
momento flexionante y de la fuerza cortante
𝐺 módulo de cortante
M𝑧 momento flexionante alrededor del eje “z”, en N-m
V𝑦 cortante en dirección “y”, en N
222
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
V𝑦a
M𝑧 (𝑥) = [𝑥 − 1] [ ] (Ad.26)
M𝑧a
donde
M𝑧 (𝑥) momento flexionante alrededor del eje z, en N-m
V𝑦 (𝑥) cortante en y a una distancia x, en N
V𝑦a fuerza cortante en y en el nodo a, en N
M𝑧a momento flexionante alrededor del eje z en el nodo a, en N-m
Ad.2 Ejemplo
Calcular los momentos en los extremos de las barras de la viga mostrada en la Fig Ad.4
empleando a) el método de Cross, b) Comparar los resultados con un programa
comercial. Considérese EI constante.
1 2 3 4
6m 6m 6m
Para elementos 1 y 3
P𝑎𝑙 𝑏𝑙2
𝑀𝐸11 = − (Ad.28)
𝑙2
P𝑎𝑙2 𝑏𝑙
𝑀𝐸12 = + (Ad.29)
𝑙2
223
C.1.5 ANÁLISIS DE ESTRUCTURAS
AYUDAS DE DISEÑO
donde
𝑀𝐸11 momento de empotramiento perfecto en el nodo 1 del elemento 1, en N-m
𝑀𝐸21 momento de empotramiento perfecto en el nodo 2 del elemento 1, en N-m
𝑎𝑙 longitud medida del origen a la posición de la fuerza P, en m
𝑏𝑙 longitud medida de la posición de la carga P al nodo final del elemento, en m
P fuerza, en N
1
(5 )(3)(3)2
𝑀𝐸1,2 =− = −3.75 𝑚2
(6)2
1
(5)(3)2 (3)
𝑀𝐸2,1 =+ = +3.75 𝑚2
(6)2
3
(5)(3)(3)2
𝑀𝐸3,4 =− = −3.75 𝑚2
(6)2
3
(5)(3)2 (3)
𝑀𝐸4,3 =+ = +3.75 𝑚2
(6)2
Para el elemento 2
2
𝑞𝑙 2
𝑀𝐸2,3 =−
12
3
𝑞𝑙 2
𝑀𝐸3,2 =+
12
2
(2)(6)2
𝑀𝐸2,3 =− = −6 𝑚2
12
3
(2)(6)2
𝑀𝐸3,2 =+ = +6 𝑚2
12
b) Momento de distribución.
𝑘𝑟𝑛
𝑀𝐷𝑖𝑛 = 𝑀𝐸𝑖𝑛
∑ 𝑘𝑟𝑛
donde
𝑀𝐸𝑖𝑛 momento de empotramiento del elemento 𝑛 en el nodo 𝑖, en N-m
𝑛
𝑘𝑟 rigidez rotacional del elemento n, en N/m.
224
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
c) Momento de transporte.
1
𝑀𝑇𝑖𝑛 = 𝑀𝑒𝑖
2
donde
𝑀𝑇𝑖𝑛 momento de transporte del nodo 𝑖 para el elemento 𝑛, en N-m
𝑀𝐷𝑖𝑛 momento de distribución en el nodo 𝑖 debido a la rigidez rotacional
para el elemento n, en N-m
Primera liberación
Para calcular los momentos de distribución es necesario determinar cuáles son los nodos
que tienen posibilidad de rotar, para este caso son los nodos 2 y 3.
2
∑ 𝑀𝐸1,2 = 3.75 + 6 = +2.25 t-m
Los momentos se consideran de signo contrario ya que en esta etapa se liberan para
distribuirlos a los extremos de los elementos que concurren sobre el nodo en estudio.
𝐼
𝑀𝐷21 = 6 (2.25) = +1.125 t-m
𝐼 𝐼
( + )
6 6
1
𝑀𝑇11 = 1.125 = +0.56 t-m
2
𝐼
𝑀𝐷22 = 6 (+2.25) = +1.125 t-m
𝐼 𝐼
( + )
6 6
1
𝑀𝑇32 = 1.125 = +0.56 t-m
2
225
C.1.5 ANÁLISIS DE ESTRUCTURAS
AYUDAS DE DISEÑO
Segunda liberación
2
∑ 𝑀𝐸2,3 + 𝑀𝑇23 = (+3.75 − 6) − 0.56 = −2.81 t-m
𝐼
𝑀𝐷32 = 6 (−2.81) = −1.41 t-m
𝐼 𝐼
( + )
6 6
1
𝑀𝑇22 = (−1.41) = −0.71 t-m
2
𝐼
𝑀𝐷33 = 6 (−2.81) = −1.41 t-m
𝐼 𝐼
( + )
6 6
1
𝑀𝑇43 = (−1.41) = −0.71 t-m
2
M=M21 + 𝑀𝐷
M=3.75+1.125= 4.875 t—m
M=M23 + 𝑀𝐷- 𝑀𝑇
M=-6+1.125-0.71=-5.585 t—m;
+4.875-5.585= -0.71≠0.
Se observa que no hay equilibrio, el proceso aún no converge, se requiere realizar más
iteraciones.
226
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
Se presenta la geometría y las cargas que actúan sobre la trabe en la Fig Ad.5.
227
C.1.5 ANÁLISIS DE ESTRUCTURAS
AYUDAS DE DISEÑO
Fig Ad.5 Geometría del modelo estructural y cargas que actúan sobre la estructura
Se observa que los resultados obtenidos tanto por el método de Cross como por el
programa comercial son similares.
Fig Ad.7 Elementos mecánicos en los extremos de las barras que conforman a la
estructura
Ad.3 Ejemplo
Caso 1
a) Desplazamiento en el punto b
b) Momento flexionante al centro del claro.
c) Fuerza cortante en el punto a
d) Esfuerzo debido al momento flexionante obtenido en b)
e) Esfuerzo debido a la fuerza cortante obtenido en c)
f) Reacción en el punto a
Caso 2
a) Desplazamiento en el punto b.
b) Momento flexionante al centro del claro.
c) Fuerza cortante en el punto a.
e) Esfuerzo debido al momento flexionante obtenido en b)
f) Esfuerzo debido a la fuerza cortante obtenido en c)
g) Reacción en el punto a
229
C.1.5 ANÁLISIS DE ESTRUCTURAS
AYUDAS DE DISEÑO
a b c
Caso 1
M𝑧
𝜎𝑏𝑥 = 𝑦 (Ad.30)
𝐼𝑧 ℎ
donde
𝜎𝑏𝑥 esfuerzo debido al momento flexionante, en kg/cm2
M𝑧 momento flexionante alrededor del eje z, en kg-cm
𝐼𝑧 momento de inercia de la sección transversal, en cm4
𝑦ℎ distancia medida a partir del eje neutro, en cm
𝑏𝑣 ℎ𝑣3 (20)(40)3
𝐼𝑧 = = = 106,666.67 cm4
12 12
donde
ℎ𝑣 peralte de la viga en cm
ℎ𝑣
𝑦ℎ = = 20 cm
2
(200,000) kg
𝜎𝑏𝑥 = (20) = 37.5
106,666.67 cm2
231
C.1.5 ANÁLISIS DE ESTRUCTURAS
AYUDAS DE DISEÑO
𝑉𝑦 𝜎𝑏𝑝 2
𝜏= [( ) − 𝑦ℎ2 ] (Ad.31)
2𝐼𝑧 2
donde
𝜏 esfuerzo cortante, en kg/cm2
V𝑦 cortante en dirección y, en kg
𝜎𝑏𝑝 esfuerzo principal, en kg/cm2
2050 40 2 kg
𝜏= [( ) − 0] = 3.84
2(106,666.67) 2 cm2
Caso 2
232
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
2. El esfuerzo promedio debido a la flexión en centro del claro, en la fibra más alejada
del eje neutro (nodos 535, 514 y 513) es de 37.537 kg/cm2, ver Tabla Ad.2. Se
muestran los elementos finitos en los cuales se estimó los esfuerzos normales en
la Fig Ad 17.
Nodo 𝜎(kg/cm2)
535 37.437
514 37.737
513 37.437
Promedio 37.537
233
C.1.5 ANÁLISIS DE ESTRUCTURAS
AYUDAS DE DISEÑO
3. El esfuerzo promedio debido cortante en los apoyos, sobre el eje neutro (nodos
247, 249 y 329) es de 3.831 kg/cm2, ver Tabla Ad.3. Se muestran los elementos
finitos en los cuales se estimó los esfuerzos cortantes en la Fig Ad. 18, localizados
sobre el apoyo “a”.
Nodo 𝜏(kg/cm2)
247 3.851
249 3.792
329 3.851
Promedio 3.831
234
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
37.537
M𝑧 = 106,666.67 = 2 𝑡𝑜𝑛 − 𝑚
20
(2)(3.831)(106,666.67)
V𝑦 = = 2043.2 𝑘𝑔
40 2
[( ) − 02 ]
2
Fy
Nodo
(t)
247 0.523
249 1.005
329 0.523
Suma 2.051
235
C.1.5 ANÁLISIS DE ESTRUCTURAS
AYUDAS DE DISEÑO
Elemento Elemento
Cálculo de: viga finito
Desplazamiento en b (mm) -1.403 -1.438
Momento flector en b (t-m) -2 -2
Fuerza cortante en a (t) 2.05 2.04
Esfuerzo debido a momento en b
(kg/cm2) 37.5 37.537
Esfuerzo debido a cortante en a
(kg/cm2) 3.84 3.831
Reacción en a (t) 2.05 2.05
236
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
ANEXOS
La viga es un elemento estructural que se utiliza con mayor frecuencia dentro de las
edificaciones. Puede ser viga simple o continua. Se trata de un elemento estructural que
se utiliza para transmitir las cargas, provenientes de las losas a las columnas. En la Fig An.1
se muestra una superficie de rodamiento de un puente, la cual se puede idealizar como
una viga simplemente apoyada.
En la Fig An.2 se muestra la calzada de un puente apoyada sobre columnas. Esta condición
hace que todos los tramos del puente estén unidos a lo largo de su longitud, dando lugar
a que las vigas se comporten como vigas continúas.
237
C.1.5 ANÁLISIS DE ESTRUCTURAS
ANEXOS
Las armaduras son estructuras que están formadas por barras, éstas pueden trabajar en
el plano (armaduras en 2D) o en el espacio (armaduras en 3D). Las barras que componen
a una armadura trabajan exclusivamente a fuerza axial, ya sea a tensión o compresión. En
la Fig An.3 se muestra la estructuración de un puente a base de armaduras.
238
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
Las armaduras empezaron a ser utilizadas en el siglo XIX, un ejemplo es la torre Eiffel
construida en Paris, Francia en 1889, esta estructura está formada por diferentes
armaduras espaciales, tal y como se muestra en la Fig An.4.
An.1.3 Retículas
Las retículas son estructuras que tienen la misma configuración de un marco plano, pero a
diferencia de este último, las cargas se aplican en dirección perpendicular al plano. Este
tipo de estructuras están formadas por elementos los cuales trabajan esencialmente a
flexión y a torsión. Se emplean en parrillas de cimentación, voladizos, etc. La Fig An.5
Muestra una retícula plana típica de acero.
239
C.1.5 ANÁLISIS DE ESTRUCTURAS
ANEXOS
An.1.4 Arcos
En la Fig An.6 se muestran diferentes estructuras de arcos y las partes que lo componen.
En el subcap 3 se hace énfasis en que un elemento arco trabaja básicamente a compresión;
despreciando los efectos de flexión y de cortante.
240
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
Clave
Dovela
Trasdós
Salmer
Flecha
Imposta
Intradós
Claro
Clave
Contraclave
Dovela
Riñon
Trasdós
Flecha
Intradós
Linea de
Salmer arranque
Imposta
Claro
Jamba Pilar
Flecha
Salmer
Claro
o
di
Ra
241
C.1.5 ANÁLISIS DE ESTRUCTURAS
ANEXOS
El uso del arco como elemento estructural data de tiempos ancestrales. Comúnmente se
ha utilizado para construir puentes, acueductos y edificios Fig An.7.
Puente
Acueducto
Edificio
242
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
Los marcos son estructuras utilizadas principalmente en las edificaciones. Su uso más
extendido es en la estructuración de edificios de distinta índole, aunque también existen
aplicaciones en naves industriales y puentes, entre muchas otras. Los elementos que
componen a un marco se clasifican en dos grupos: a) vigas, que son los elementos
horizontales del marco y trabajan esencialmente a flexión y a cortante y b) columnas que
trabajan principalmente a flexión, carga axial y cortante. En la Fig An.10 se muestra una
estructura industrial formada por marcos planos de acero
Los contravientos son elementos estructurales que se utilizan en los marcos que componen
a las edificaciones, con la finalidad de aumentar su rigidez y con ello soportar cargas
laterales importantes sin experimentar grandes desplazamientos y deformaciones. Esto se
presenta regularmente por cargas accidentales como sismo o viento. En la Fig An.11 se
muestra un edificio formado por marcos con contravientos de acero.
243
C.1.5 ANÁLISIS DE ESTRUCTURAS
ANEXOS
Los marcos tridimensionales son estructuras formadas principalmente por marcos planos.
La Fig An.12 muestra la estructuración, la cual está formada por marcos de acero
tridimensionales.
Este tipo de edificios son los que se construyen con mayor frecuencia. En general están
estructurados por marcos de concreto reforzado y acero. En la Fig An.13 se muestra la
estructuración de un edificio de concreto reforzado.
244
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
Una nave industrial es una estructura propia para realizar actividades de producción,
trasformación, manufactura, ensamble, almacenaje, distribución entre otros. En muchos
casos están formadas por elementos estructurales prefabricados, los cuales pueden librar
grandes claros sin existencia de apoyos intermedios, de tal manera que la nave pueda
operar sin obstáculos y restricciones. Se muestra una nave industrial en la Fig An.14.
245
C.1.5 ANÁLISIS DE ESTRUCTURAS
ANEXOS
Algunas estructuras utilizan cables como elemento estructural, estos sirven para cubrir
distancias largas y soportan cargas a tensión. En la Fig An.15 se muestra una columna, la
cual utiliza un sistema de cables para su estabilidad.
246
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
Las estructuras geodésicas pueden sufrir pandeo global sin que ninguna de las barras
comprimidas que la forman haya sufrido pandeo local. Eso implica que un cálculo como
estructura lineal convencional, y comprobación posterior de pandeo local, puede no ser
adecuado en muchos casos y para grandes claros se requiere un cálculo no lineal para
determinar sus cargas críticas y asegurarse de que no se producen fenómenos de
inestabilidad estática, Fig An.17.
Este tipo de estructuras son capaces de transmitir cargas a los apoyos mediante
solicitaciones y deformaciones de sus componentes, en este caso una disposición especial
de barras, con o sin cables, donde cada componente se va a poder mover. En la Fig An.18
muestra un puente desplegable.
247
C.1.5 ANÁLISIS DE ESTRUCTURAS
ANEXOS
Este tipo de estructuras están formadas por armaduras espaciales hechas principalmente
de acero, se caracterizan por tener una base ancha, luego estrechándose desde el fondo
hasta el punto medio de la estructura, Fig An.19.
An.1.15 Postes
Un Poste es un elemento estructural tipo columna formada por concreto, acero o madera,
se utiliza para la construcción de tendidos de líneas, redes eléctricas, telefónicos, etc.
Fig An.20.
Las estructuras atirantadas utilizan cables como elemento estructural. Este es el caso de
los puentes atirantados. Ver Fig An.21.
Son estructuras que se localizan fuera de la costa. Por ejemplo, plataformas marinas, las
cuales están formadas por distintos elementos estructurales como: armaduras espaciales,
vigas, columnas (Fig An.22).
249
C.1.5 ANÁLISIS DE ESTRUCTURAS
ANEXOS
Estos muros tienen básicamente la función de absorber una parte de la fuerza cortante
horizontal, con la finalidad de controlar los desplazamientos laterales (sismo y viento). Ver
Fig An 23.
An.1.19 Placas
Las placas son estructuras rígidas, las cuales dispersan las cargas aplicadas según un
patrón multidireccional, siguiendo generalmente las rutas más cortas y más rígidas hasta
los apoyos. Un ejemplo de uso de una placa se muestra en la Fig An.24.
250
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
An.1.21 Losas
Las losas son placas de concreto armado, las cuales están apoyados generalmente sobre
el terreno o vigas y muros. La función principal de estos elementos estructurales es
distribuir las cargas sobre los apoyos. En la Fig An.26 se muestra una losa de cimentación.
251
C.1.5 ANÁLISIS DE ESTRUCTURAS
ANEXOS
An.1.22 Membranas
Una membrana es un elemento estructural de espesor pequeño, sin rigidez flexional que
soporta tensiones y esfuerzos normales. Ver Fig An.27.
An.1.23 Cascarones
Los cascarones son estructuras resistentes por forma, lo suficientemente finas como para
no desarrollar momentos apreciables, pero con el suficiente espesor para soportar cargas
axiales y cortantes; por lo que son ideales para la construcción con concreto reforzado. Ver
Fig An.28.
252
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
An.1.24 Tanques
Los tanques son estructuras fundamentales para el control del abastecimiento de algunos
recursos tales como agua potable, gas, etc. (Fig An.29).
An.1.25 Chimeneas
Las chimeneas son estructuras que están construidas de concreto reforzado u otro material,
que sirven para dejar salir a la atmosfera gases resultantes de una combustión o de una
reacción química, tal como se muestra en la Fig An.30.
253
C.1.5 ANÁLISIS DE ESTRUCTURAS
ANEXOS
An.1.26 Silos
Un silo es una estructura diseñada para almacenar grano y otros materiales a granel. Los
más comunes tienen forma cilíndrica, asemejándose a una torre, construida de madera,
metal o concreto reforzado, tal como se muestra en la Fig An.31.
254
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
donde
𝐔̈ vector de aceleración de la ecuación de equilibrio dinámico
𝐔 vector de desplazamiento de la ecuación de equilibrio dinámico
𝑴 Matriz de masa
𝑲 matriz de rigidez de la ecuación de equilibrio dinámico
𝐮𝑡 = 𝑞𝑛 (𝑡)𝚽 (An.2)
donde
𝑞𝑛 (𝑡) desplazamiento con respecto al tiempo de un oscilador de un grado
de libertad, en m
𝚽 matriz de formas modales
𝐮𝑡 vector de desplazamientos con respecto al tiempo
𝐴𝑛 , 𝐵𝑛 constantes de integración, en m
𝜔𝑛𝑎𝑡 frecuencia natural, rad/s
donde
𝜹𝑎𝑣 vector formado por las amplitudes de los desplazamientos normales
(autovectores)
Para que esta ecuación tenga una solución distinta de la trivial es necesario que su
determinante sea nulo,
255
C.1.5 ANÁLISIS DE ESTRUCTURAS
ANEXOS
2
[𝑲 − 𝜔𝑛𝑎𝑡 𝑴] = 0 (An.5)
𝑲𝐀 = 𝑴𝐀𝛀 (An.6)
donde
𝐀 matriz de autovectores
𝛀 matriz diagonal de autovalores,
2
𝜔𝑛𝑎𝑡 0 ⋯ 0
2
𝜔𝑛𝑎𝑡 ⋯ 0
𝛀=[ 0 ⋱ ⋮ ] (An.8)
⋮ ⋮
⋯ 2
𝜔𝑛𝑎𝑡
0 0
𝐀𝐓 𝑲𝐀 = 𝛀 (An.9)
𝐀𝐓 𝑴𝐀 = 𝑰 (An.10)
donde
𝑰 matriz identidad
𝜹𝑎𝑣1 𝜹𝑎𝑣2 ⋯ 𝜹𝑎𝑣𝑛 coeficientes de la matriz de autovectores
El método de Jacobi generalizado comienza utilizando una matriz anulación (𝐐) que anula
simultáneamente el mismo término situado fuera de la diagonal en las matrices. De la
misma forma una nueva matriz de transformación anula un nuevo término situado fuera
de la diagonal. Mediante este proceso se logra diagonalizar las matrices 𝑲 y 𝐌 que se
transforman, respectivamente, en lo siguiente.
̂ 𝑘+1 = 𝐐𝑇𝑘 𝐊 𝑘 𝐐𝑘
𝐊 (An.11)
256
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
̂ 𝑘+1 = 𝑷𝑇𝑘 𝐌𝑘 𝑷𝑘
𝐌 (An.12)
donde
̂
𝐌 matriz diagonal de masas
̂
𝐊 matriz diagonal de rigideces para una iteración
𝐐 matriz de anulación
𝐐𝑇𝑘 𝐐𝑘 = 𝑰 (An.13)
𝟏
̂ −𝟐
𝐀 = 𝐐𝐌 (An.14)
̂ −𝟏 𝐊
𝛀=𝑴 ̂ (An.15
siendo
1 0 0
−
𝑴112 1 ⋯ 0
− ⋯
̂ −𝟏/𝟐
𝐌 = 0 𝑴222 ⋱ ⋮ (An.16)
1
⋮ ⋮ ⋯ −
2
[ 0 0 𝑴𝑛𝑚 ]
1 0 0
−
𝑴112 1 ⋯ 0
− ⋯
̂ −𝟏
𝐌 = 0 𝑴222 ⋱ ⋮ (An.17)
1
⋮ ⋮ ⋯ −
[ 0 0 𝑴𝑛𝑛2 ]
donde
1 1 1
− − −
𝑴112 , 𝑴222 , 𝑴𝑛𝑛2 coeficientes de la matriz diagonal de masas
2
𝑲𝜹𝑎𝑣 = 𝜔𝑛𝑎𝑡 𝑴𝜹𝑎𝑣 (An.18)
257
C.1.5 ANÁLISIS DE ESTRUCTURAS
ANEXOS
𝜹𝑻𝑎𝑣 𝑲𝜹𝑎𝑣
𝑅(𝜹𝑎𝑣 ) = (An.20)
𝜹𝑻𝑎𝑣 𝐌𝜹𝑎𝑣
donde
𝑅(𝜹𝑎𝑣 ) cociente de Rayleigh
donde
𝜔𝑛𝑎𝑡 𝑚𝑖𝑛 frecuencia natural mínima
𝜔𝑛𝑎𝑡 𝑚𝑎𝑥 frecuencia natural máxima
2
𝑲𝜹𝑎𝑣1 = 𝜔𝑛𝑎𝑡 𝑴𝜹𝑎𝑣1 (An.22)
donde
𝜹𝑎𝑣𝐼 cociente de Rayleigh para el modo 1
Para ello se procede de forma iterativa, partiendo con un vector inicial que suponemos
es una estimación de un múltiplo el primer modo. La primera aproximación del vector
inicial que satisface la ec An.19, se obtiene resolviendo la ecuación
donde
𝐯 vector inicial para estimar el vector de formas modales, adimensional
𝐯1 vector para estimar el vector de formas modales asociado a la
iteración 1
𝐯2 vector para estimar el vector de formas modales asociado a la
iteración 2
Continuando las iteraciones de esta forma, se llega a un vector final con el que finaliza
el proceso iterativo. Se estima que el primer modo normalizado, respecto a la matriz de
masa 𝑴, es
𝐯𝑓
𝚽𝑚 = 1 (An.25)
[𝐯𝑓𝑇 𝑴𝐯𝑓 ]2
258
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
donde
𝚽𝑚 modo normalizado, respecto a la matriz de masa
𝐯𝑓 vector final para estimar el vector de formas modales
𝐯𝑓𝑇 𝑲𝐯𝑓
𝜔𝑛𝑎𝑡 = (An.26)
𝐯𝑓𝑇 𝑴𝐯𝑓
Este método es útil para problemas de grandes dimensiones en los que se desea conocer
solamente unos modos y frecuencias de vibrar de la estructura.
El objetivo fundamental del método es resolver el problema para los valores propios más
pequeños y sus correspondientes vectores característicos.
2
[𝑲 − 𝜔𝑛𝑎𝑡 𝑴]𝜹𝑎𝑣1 = 0 (An.27)
donde
𝚿0 matriz inicial para estimar el vector de formas modales
𝚿1 matriz de la primera iteración
donde
𝜿1 , 𝝁1 matrices asociadas al problema de valores característicos
259
C.1.5 ANÁLISIS DE ESTRUCTURAS
ANEXOS
Una línea de influencia es un diagrama que define la variación de una fuerza interna,
esfuerzo, reacción o deformación en un punto fijo de una estructura elástica lineal a
medida que una carga unitaria se mueve sobre ella.
El método de las líneas de influencia permite definir la colocación de una carga viva o
variable de tal manera en la que se producen efectos máximos de cortante, flexión,
reacciones y deflexiones tanto para cargas puntuales como para cargas distribuidas.
Este método se utiliza mucho para cargas vivas sobre puentes, puentes grúas, bandas
transportadoras y especialmente en aquellas estructuras con cargas móviles.
Para estimar una línea de influencia se utiliza una carga unitaria ya que de acuerdo con
los conceptos de linealidad, proporcionalidad y superposición es posible determinar la
función específica simplemente multiplicando el valor de la línea de influencia por el
valor de la carga real.
Para encontrar los valores máximos de este efecto la carga puntual se debe colocar en
el punto de máxima ordenada.
Para determinar una línea de influencia de una carga distribuida se puede usar la línea
de la carga puntual para determinar los tramos donde colocar la carga distribuida que
produzca los valores máximos en un punto.
260
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
Se empieza a mover la carga puntual a diferentes distancias, para cada una se calcula
𝑅𝐴 .
Para encontrar el valor de la reacción en el apoyo, para cada posición de una carga
unitaria aplicada, se utiliza la ec An.44.
𝑃∗ (𝑙 − 𝑥 ∗ )
𝑅𝐴 = (An.31)
𝑙
donde
𝑅𝐴 reacción en el apoyo A, en N
𝑃∗ carga unitaria aplicada a una distancia x, en N
𝑥 ∗
distancia donde se aplica carga unitaria, en m
𝑙 longitud del elemento estructural, en m
261
C.1.5 ANÁLISIS DE ESTRUCTURAS
ANEXOS
Internacional Técnico
MASA m
MKS CGS FPS MKS CGS PLS
Internacional Técnico
FUERZA F
MKS CGS PLS MKS CGS PLS
Internacional Técnico
MOMENTO M
MKS CGS PLS MKS CGS PLS
262
MANUAL DE DISEÑO DE OBRAS CIVILES – CFE
C.1.5 ANÁLISIS DE ESTRUCTURAS
Internacional Técnico
ESFUERZO σ
MKS CGS PLS MKS CGS PLS
Internacional Técnico
RIGIDEZ K
MKS CGS PLS MKS CGS PLS
263
CONTENIDO GENERAL DEL MANUAL
265
SECCIÓN A - HIDROTECNIA
Tema Capítulo Descripción
1 Hidrología
1 Consideraciones generales de hidrología
2 Precipitación
3 Escurrimiento
4 Pérdidas
5 Relación entre precipitación y escurrimiento
6 Análisis estadístico
7 Tormentas de diseño
8 Avenidas de diseño
9 Tránsito de avenidas y simulación del funcionamiento de vasos de almacenamiento
10 Pronóstico de avenidas
11 Geohidrología
12 Sistemas de información geográfica, con aplicación en hidrología e hidráulica
13 Hidrología e hidráulica en cuencas pequeñas
2 Hidráulica
1 Consideraciones generales de hidráulica
2 Planeación de sistemas de aprovechamiento hidroeléctrico
3 Evaluación del impacto ambiental para la ejecución de proyectos hidroeléctricos
4 Conducciones a presión
5 Obra de toma para plantas hidroeléctricas
6 Obra de toma y sistemas de enfriamiento para plantas termoeléctricas
7 Máquinas hidráulicas
8 Transitorios hidráulicos en conductos a presión
9 Turbinas y análisis dinámico
10 Escurrimiento a superficie libre
11 Tránsito hidrológico e hidráulico, en cauces
12 Obra de desvío
13 Obra de excedencias
14 Diseño de aireadores en estructuras vertedoras
15 Hidráulica fluvial
16 Sedimentación en embalses y desarenadores
17 Obras marítimas
18 Terminales marítimas para recibo y manejo de combustibles
19 Interacción mar - tierra
20 Dragado
21 Instrumentación hidráulica
22 Técnicas experimentales de hidráulica
23 Métodos numéricos de hidráulica
24 Análisis de riesgo y confiabilidad en hidrología
25 Evaluación de la seguridad hidráulica de las estructuras
26 Rehabilitación hidráulica de las estructuras
267
SECCIÓN B - GEOTECNIA
Tema Capítulo Descripción
1 Geología
1 Información geológica existente
2 Obtención y análisis de información geológica
3 Presentación de información geológica.
4 Estudios de peligro sísmico
2 Mecánica de Suelos
1 Clasificación y caracterización de los depósitos de suelos
2 Pruebas de laboratorio para determinar las propiedades de los suelos y enrocamientos
3 Pruebas de campo para determinar propiedades de los suelos y enrocamientos
4 Comportamiento de suelos parcialmente saturados y aplicaciones
5 Cimentaciones en suelos
6 Estructuras de retención
7 Estabilidad de taludes en suelos, excavaciones y laderas
8 Mejoramiento de suelos
9 Flujo de agua en suelos
3 Mecánica de rocas
1 Pruebas de laboratorio de mecánica de rocas
2 Pruebas de campo de mecánica de rocas
3 Caracterización de macizos rocosos
4 Cimentaciones en rocas
5 Estabilidad y tratamientos de taludes rocosos
6 Inyecciones en macizos rocosos
7 Métodos de excavación en macizos rocosos
8 Caracterización de geomateriales frontera
4 Diseño de Presas
1 Consideraciones generales de presas
2 Presas de enrocamiento con núcleo de arcilla
3 Presas de enrocamiento con cara de concreto
4 Presas de enrocamiento con núcleo de asfalto
5 Presas de jales
6 Presas de Hard Fill
5 Túneles y lumbreras
1 Túneles y lumbreras en suelos
2 Túneles y lumbreras en rocas
3 Cavernas para casa de máquinas
6 Control de calidad
1 Control estadístico geotécnico de la construcción
2 Análisis de riesgos y confiabilidad en geotecnia
7 Instrumentación
1 Consideraciones generales de instrumentación geotécnica
2 Instrumentación geotécnica
3 Instrumentación estructural
4 Evaluación de la seguridad geotécnica de las estructuras
5 Rehabilitación geotécnica de las estructuras
8 Aplicaciones generales
1 Criterios para la exploración y pruebas de campo y laboratorio en geotecnia
2 Criterios para el empleo de métodos numéricos en la geotecnia
268
SECCIÓN C - ESTRUCTURAS
Tema Capítulo Descripción
1 Criterios generales de análisis y diseño
1 Métodos de diseño para estructuras
2 Acciones
3 Diseño por sismo
4 Diseño por viento
5 Análisis de estructuras
2 Diseño de Estructuras Especiales
1 Diseño estructural de cimentaciones
2 Diseño estructural de líneas de transmisión
3 Diseño estructural de torres de enfriamiento
4 Diseño estructural de tanques y depósitos
5 Diseño estructural de tuberías
6 Diseño estructural de chimeneas
7 Diseño estructural de puentes
8 Tópicos estructurales en obras hidráulicas
9 Diseño de estructuras industriales
10 Diseño estructural de subestaciones
11 Diseño estructural de compuertas
12 Diseño estructural de casa de máquinas
13 Diseño estructural de túneles y lumbreras
3 Monitoreo y control de respuestas estructurales
1 Monitoreo de estructuras
2 Diseño de estructuras con aislamiento de base
3 Diseño de estructuras con disipadores de energía sísmica
4 Análisis térmico en estructuras
4 Evaluación de la seguridad y rehabilitación de estructuras
1 Evaluación estructural de la seguridad de las estructuras existentes
2 Rehabilitación estructural
5 Materiales de construcción
1 Fibras de acero, polipropileno, carbón y vidrio
2 Concreto masivo
3 Concreto compactado con rodillo
4 Concreto autocompactable
5 Concreto asfáltico
6 Tecnología del concreto
6 Diseño de presas de concreto o mampostería
1 Consideraciones generales de presas rígidas
2 Presas de gravedad (Convencional y concreto compactado con rodillo)
3 Presas de arco
4 Presas de machones
269