Libro Res Mat 2
Libro Res Mat 2
9 de septiembre de 2022
“libroResMat2” — 2022/9/9 — 11:53 — page i — #2
El contenido de este documento es publicado bajo una licencia Creative Commons Attribution-ShareAlike
4.0 International License. Ver detalles en creativecommons.org/licenses/by-sa/4.0.
This work is published under a CC BY-SA license (Creative Commons Attribution-ShareAlike 4.0 Interna-
tional License), which means that you can copy, redistribute, remix, transform, and build upon the content
for any purpose, even commercially, as long as you give appropriate credit, provide a link to the license,
and indicate if changes were made. If you remix, transform, or build upon the material, you must distribute
your contributions under the same license as the original. License details: creativecommons.org/licenses/by-
sa/4.0.
The use of general descripitive names, resgistered names, trademarks, etc. in this publication does not
imply, even in the absence of a specific statement, that such names are exempt from the relevant protective
laws and regulations and therefore free for general use.
Documento producido usando software libre: LATEX, TeXstudio, Geany, Paraview, Inkscape y GNU-Octave.
“libroResMat2” — 2022/9/9 — 11:53 — page ii — #3
Lista de contenidos
1. Introducción 1
1.1. Análisis y modelado estructural . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.1.1. Definición y clasificación de estructuras . . . . . . . . . . . . . . . . 1
1.1.2. Análisis de estructuras . . . . . . . . . . . . . . . . . . . . . . . . . . 5
1.1.3. Modelado estructural . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.1.4. Ejemplo: Análisis simplificado de puente metálico . . . . . . . . . . . 8
1.2. Ecuaciones de la teoría de vigas . . . . . . . . . . . . . . . . . . . . . . . . . 10
1.2.1. Ecuaciones de la teoría de vigas en flexión pura . . . . . . . . . . . . 11
1.3. Conceptos básicos de Elasticidad Lineal . . . . . . . . . . . . . . . . . . . . . 13
1.3.1. Problema de Elasticidad Lineal . . . . . . . . . . . . . . . . . . . . . 13
1.3.2. Principios energéticos en elasticidad lineal . . . . . . . . . . . . . . . 15
1.3.3. Aplicación al desarrollo de teorías estructurales . . . . . . . . . . . . 17
1.4. Ejercicios . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
Unidad temática 1
Introducción
Definición: Estructura
barra: elemento con una dimensión (eje) considerablemente mayor que las otras. So-
porta tensiones de compresión o tracción según su eje.
(a) Ejemplo de estructura civil: vía férrea de tren de Cor- (b) Ejemplo de estructura mecánica: atracción mecánica
covado en Río de Janeiro, Brasil. en Montevideo, Uruguay.
viga (de eje recto): elemento con geometría de barra en su configuración natural (sin
cargas), que puede ser sometido a fuerzas según su eje (normales) o transversales (cor-
tantes) así como también a momentos según su eje (torsores) o transversales (flectores).
pilar: elemento con geometría de barra que se encuentra sometido a flexión y compre-
sión según su eje, siendo la compresión preponderante.
cable: también utilizado como tensor, son elementos de barra sometidos principal-
mente a tracción que no soportan considerable flexión o compresión según su eje.
losa y cáscara: elementos con una geometría tal que una dimensión es considerable-
mente menor que las otras dos. La dimensión menor define el espesor, mientras que las
otras dos definen el plano medio (o superficie media) de la losa (o cáscara). Soportan
flexión o tensión con vector contenido en el plano medio. Las losas tienen un plano
medio plano, mientras que en las cáscaras, este se encuentra dado por una superficie
con curvatura no nula.
membrana: elementos superficiales muy flexibles de pequeño espesor que no son ca-
paces de soportar compresión o flexión. En la Figura 1.2 se observa la cubierta exterior,
de la Terminal Madero de la Ciudad de Buenos Aires, formada y soportada por com-
ponentes de membrana, barras, vigas y cables.
otros: existen otros elementos estructurales con geometrías más complejas no asocia-
dos directamente a un estado tensional como los anteriores. Algunos ejemplos de estas
estructuras son: zapatas, cabezales, muros de contención (trabajando como losa o bajo
estados planos de deformación) y arcos, o vigas de eje curvo.
Figura 1.2: Foto exterior de las cubiertas de la Terminal Madero Ciudad de Buenos Aires, Argentina.
Estructura aporticada: formada por barras unidas donde al menos dos barras están
conectadas por un nudo rígido o un vínculo que transmite momentos.
1. todos los ejes de sus barras (barras, vigas o pilares) pertenecen a un plano y este
plano es plano de simetría de todas las barras de la estructura
2. todos los vectores de fuerza aplicados a la estructura son vectores contenidos en
el plano de simetría de la estructura
3. todos los momentos aplicados son aplicados en puntos que pertenecen al plano
medio de la estructura y están dados por un vector perpendicular al plano de
simetría de la estructura
4. las restricciones cinemáticas son tales que las reacciones correspondientes cum-
plan con las condiciones anteriores de fuerza y momento.
3. todos los momentos aplicados son aplicados en puntos que pertenecen al plano de
las barras de la estructura y están dados por vectores contenidos en dicho plano
4. las restricciones cinemáticas son tales que las reacciones correspondientes cum-
plan con las condiciones anteriores de fuerza y momento.
Actividad
Considere un edificio de viviendas de tres pisos de altura con planta rectan-
gular. Identificar los componentes estructurales de la estructura. Clasifique
la estructura de acuerdo a los criterios vistos.
El proceso de analizar una estructura y obtener todas las solicitaciones de sus elementos
estructurales, es usualmente llamado resolver una estructura.
Las fuerzas aplicadas a los elementos estructurales y la vinculación existente entre estos
elementos determinan las deformaciones y solicitaciones de la estructura. Para definir el tipo
de análisis a realizar se deberá categorizar las cargas y vínculos entre componentes.
A continuación se definen distintos tipos de cargas que suelen ser aplicados a estructuras.
Luego se enumeran algunos de los vínculos entre componentes estructurales más importan-
tes.
Cargas
Las cargas o fuerzas externas que pueden ser aplicadas a estructuras pueden tener diver-
sas causas, por lo que tendrán diferentes características respecto a: frecuencia de ocurrencia,
magnitudes, confiabilidad, etc. Las normas o códigos establecen metodologías para conside-
rar la distinta naturaleza de las acciones al diseñar estructuras.
En este material se considera que toda carga es de tipo estática o permanente, indepen-
dientemente de su origen. Estas cargas son provocadas por acciones cuya magnitud o sentido
no varía en el tiempo, como por ejemplo el peso propio de cada elemento estructural.
Las cargas o esfuerzos producidos por otras acciones como viento, sismos, impactos, etc.,
deben ser analizadas de forma diferente aplicando modelos numéricos o criterios de diseño
específicos definidos por normativas correspondientes. Como ejemplo, se puede ver la forma
en que se definen cargas en la norma europea1 .
En este texto se asume que, sea cual fuere su origen, toda carga podrá ser considerada
como estática. Existe una importante rama de la Ingeniería Estructural, enfocada al análisis
dinámico de estructuras, en la cual se destacan referencias bibliográficas como (Clough and
Penzien, 1993).
libre: un nodo o un extremo de una viga que no está vinculado a tierra o que no tiene
desplazamiento ni giro impuesto ni restringido en ninguna dirección o ángulo,
giro impedido: vínculo que impide el giro de un nodo según un cierto vector,
vínculo mixto: vínculo con combinaciones que no pertenecen a ninguna de las cate-
gorías anteriores.
Existen otros tipos de vínculos que pueden ser considerados entre elementos estructura-
les, pero no serán vistos en el curso.
1 Eurocódigo 1: UNE-EN 1991-1-2:2019 https://www.une.org/encuentra-tu-norma/busca-tu-
norma/norma?c=N0061459
(a) Elementos de transmisión de cargas y so- (b) Elementos de transmisión de cargas y soporte.
porte.
Figura 1.3: Fotos de interior de terminal Madero, Ciudad de Buenos Aires. (izq.), vínculo a tierra articulado (der.).
Actividad
Enumere las cargas presentes durante la vida útil del edificio de tres pisos
considerado en la actividad anterior. Analice qué tipo de apoyo puede tener
la estructura sobre el suelo así como también qué vínculos pueden existir
entre los elementos de la estructura.
Estructura Real El concepto de Estructura Real (ER), podrá ser utilizado para referirse a
una estructura existente (ya construida), o a aquella que se desea construir y que es imaginada
por los diseñadores. En toda ER se cuenta con una geometría con imperfecciones debidas al
3. a través del uso de un MN obtener las solicitaciones a las que cada componente de
la estructura estará sometido
4. diseñar los componentes de las estructuras para soportar las solicitaciones (acero,
secciones, tipo de hormigón, etc.)
𝓁
𝑧
2 3 4 5 6 7 8 9
1 𝑥
𝑃 𝑃 𝑃 𝑃 𝑃 𝑃 𝑃
(c) Diagrama de deformada con factor de escala 10 y escala de colores de tensiones axiales.
Figura 1.4: Ejemplo simplificado de las principales componentes en el proceso de diseño o verificación de una es-
tructura.
En la Figura 1.4a se muestra una foto de un puente metálico, que será considerado como
la estructura real.
Para pasar al EB se consideran diversas hipótesis, entre las que se destacan:
se consideran dos reticulados planos con todos sus nodos desplazandose en el plano que
los contiene, en lugar de una estructura tridimensional (teniendo en cuenta el estado
de cargas simétrico),
se asume que el material es elástico lineal y que la estructura sufre pequeñas deforma-
ciones y desplazamientos para las cargas aplicadas.
En la Figura 1.4b se muestra el EB obtenido para uno de los reticulados planos considerando
las hipótesis mencionadas. Es importante destacar que otras hipótesis pueden ser considera-
das, por ejemplo, el vínculo entre los perfiles puede no ser articulado.
Finalmente se utiliza el MEF para obtener las magnitudes deseadas resolviendo numéri-
camente las ecuaciones planteadas por el MN. Se modela cada barra con un elemento de barra
plano de dos nodos, con funciones de interpolación lineales. Luego de resolver el sistema de
ecuaciones lineales correspondiente, se obtienen los desplazamientos y luego de utilizar un
factor de escala adecuado (10 en este caso) se obtiene una estructura deformada dada por la
Figura 1.4c. La escala de colores representa los valores de tensiones axiales.
En este caso para resolver el Modelo Numérico se utilizó la herramienta de elementos
finitos ONSAS 2 . El archivo de entrada usado para resolver el reticulado, con la geometría y
conectividad de la estructura, está disponible en el repositorio de este libro.
En este texto nos concentraremos en el estudio de métodos para la resolución del Esquema
Básico de Cálculo, a través de la resolución del Modelo Numérico o del uso de soluciones
analíticas. De aquí en adelante se podrá utilizar la palabra estructura para hacer referencia al
Esquema Básico de Cálculo.
𝐞𝑦
𝐞𝑧 G
𝐞𝑥
el eje 𝐞𝑥 pasa por los baricentros de las secciones transversales, denotados por 𝐺,
2 Código abierto disponible en: www.onsas.org
𝐞𝑦
𝜃(𝑥)
𝐞𝑥 𝑣(𝑥)
𝑑𝑣
tan(𝜃(𝑥)) = (𝑥) (1.7)
𝑑𝑥
por lo que se puede obtener la función 𝜃(𝑥) como:
𝑑𝑣
𝜃(𝑥) = arctan (𝑥) . (1.8)
( 𝑑𝑥 )
Usando la identidad:
𝑑 arctan(𝑢) 1
= , (1.9)
𝑑𝑢 1 + 𝑢2
y la regla de la cadena se obtiene:
𝑑𝜃 1 𝑑2𝑣
(𝑥) = 2 2
(𝑥). (1.10)
𝑑𝑥 𝑑𝑣
1 + ( 𝑑𝑥 (𝑥)) 𝑑𝑥
𝑑𝑥 1
(𝜃(𝑥)) = ,
𝑑𝜃 𝑑𝜃
(𝑥)
𝑑𝑥
por lo tanto se obtiene
2
𝑑𝑥 𝑑𝑣 1
(𝜃(𝑥)) = 1 + (𝑥) 𝑑2 𝑣
(1.11)
𝑑𝜃 ( ( 𝑑𝑥 ))
𝑑𝑥 2 (𝑥)
𝑑2𝑣 1
𝜅(𝑥) = . (1.12)
𝑑𝑥 2 𝑑𝑣 2 3/2
(1 + ( 𝑑𝑥 ) )
En el caso de que la viga desarrolle pequeños giros 𝜃 ≪ 1, o que esa sea la hipótesis
𝑑𝜃 𝑑2𝑣
𝜅(𝑥) = (𝑥) = 2 . (1.14)
𝑑𝑥 𝑑𝑥
La deformación axial en la sección 𝑥 y en la fibra ubicada en la posición 𝑦 puede ser
también aproximada como:
𝑑2𝑣 𝑑𝜃
𝜀(𝑥, 𝑦) = −𝑦 2
(𝑥) = −𝑦 (𝑥), (1.15)
𝑑𝑥 𝑑𝑥
aproximación que será útil para el desarrollo del elemento finito de viga.
Ecuaciones de equilibrio
Sea un sólido ocupando una región Ω del espacio, como el mostrado en la Figura 1.7,
sometido a una fuerza de volumen 𝐛, cuya inercia puede ser despreciada, las ecuaciones pun-
𝐭̂
Γ𝐭
Γ𝐮
Ecuación constitutiva
Para materiales hiperelásticos, las tensiones están relacionadas de forma directa con el
tensor de deformaciones 𝜺 a través de la ecuación
𝜕Ψ
𝝈= (𝜺), (1.17)
𝜕𝜺
donde Ψ es la función de densidad de energía de deformación y 𝜺 es el tensor de deformaciones
infinitesimales.
Para materiales con comportamiento elástico lineal se tiene
1
Ψ(𝜺) = 𝜺 ∶ ℂ[𝜺], (1.18)
2
donde ℂ es el tensor constitutivo. Usando la Ecuación (1.17) se obtiene la ecuación constitu-
tiva en la forma:
𝜎 = ℂ[𝜺], (1.19)
es decir una relación lineal entre tensión y deformación.
En este texto se consideran materiales isótropos, los cuales tienen igual comportamiento
constitutivo en todas las direcciones. Se puede mostrar que el tensor constitutivo está dado
por dos parámetros (parámetros de Lamé), permitiendo obtener la siguiente expresión de la
ecuación constitutiva:
𝝈(𝜺) = 𝜆 Tr(𝜺)𝐈 + 2𝜇𝜺, (1.20)
donde Tr es la función traza y 𝜆 y 𝜇 son los parámetros de Lamé. Estos últimos están relacio-
nados con el módulo de Young 𝐸 y el coeficiente de Poisson 𝜈 a través de las expresiones
𝐸𝜈 𝐸
𝜆= y 𝜇= . (1.21)
(1 + 𝜈)(1 − 2𝜈) 2(1 + 𝜈)
Relaciones cinemáticas
El tensor de deformaciones está relacionado con el campo de desplazamientos 𝐮 a través
de la expresión:
1
𝜀 = (∇u + ∇u𝑇 ). (1.22)
2
Condiciones de contorno
Finalmente se presentan las condiciones de contorno en desplazamientos sobre Γ𝐮 y en
tensiones sobre Γ𝐭 {
𝐮 = 𝐮̂ en Γ𝐮
(1.23)
𝝈[𝐧] = 𝐭 en Γ𝐭
̂
donde el miembro izquierdo representa el trabajo de las fuerzas internas y el miembro derecho
el trabajo de las fuerzas externas.
Esta relación permite plantear el problema de encontrar 𝝈 como el problema de encontrar
cuál es el campo tensorial que verifica la igualdad anterior para cualquier campo de despla-
zamientos virtual.
Por otra parte, la energía potencial de deformación está dada por la integral de la densidad
de energía de deformación:
1
Π𝑖𝑛𝑡 (𝐮) = ∫ Ψ(𝐮)d𝑉 = ℂ[𝜺(𝐮)] ∶ 𝜺(𝐮)𝑑𝑉 . (1.28)
Ω 2 ∫Ω
Se define la energía potencial total del sistema para el campo de desplazamientos 𝐮, de-
notada por Π(𝐮) y dada por la expresión:
∫ 𝐛𝐴 ⋅ 𝐮𝐵 d𝑉 + ∫ 𝐭𝐴 ⋅ 𝐮𝐵 d𝑆 = ∫ 𝐛𝐵 ⋅ 𝐮𝐴 d𝑉 + ∫ 𝐭𝐵 ⋅ 𝐮𝐴 d𝑆. (1.31)
Ω Γ𝐭 Ω Γ𝐭
⎡𝜀𝑥 0 0 ⎤
𝜺(𝐱) = ⎢ 0 −𝜈𝜀𝑥 0 ⎥ ∀𝐱 ∈ Ω (1.32)
⎢ ⎥
⎣0 0 −𝜈𝜀𝑥 ⎦
⎡𝐸𝜀𝑥 0 0⎤
𝝈=⎢ 0 0 0⎥ , (1.33)
⎢ ⎥
⎣ 0 0 0⎦
por lo que la única ecuación constitutiva relevante en los principios energéticos es 𝜎𝑥 = 𝐸𝜀𝑥 .
A partir de la ecuación de equilibrio y considerando que no hay fuerzas de volumen
aplicadas en el interior de la barra se tiene que
𝜕𝜎𝑥 𝜕𝜀𝑥
=0=𝐸 , (1.34)
𝜕𝑥 𝜕𝑥
por lo tanto la deformación es uniforme en la barra y además puede ser escrita en términos
de los desplazamientos de los nodos de los extremos. Si 𝑢2 es el desplazamiento del nodo del
extremo derecho y 𝑢1 el del extremo izquierdo, se tiene:
𝑢2 − 𝑢1
𝜀𝑥 = (1.35)
𝓁
donde 𝓁 es el largo de la barra.
Sustituyendo estas expresiones en la definición de energía potencial se obtiene una ex-
presión de Π en función de los desplazamientos nodales:
1 1
Π(𝑢1 , 𝑢2 ) = 𝐸 (𝑢2 − 𝑢1 )2 𝑑𝑉 − 𝐹1 𝑢1 − 𝐹2 𝑢2 (1.36)
2 ∫Ω 𝓁2
calculando la integral y derivando respecto a 𝑢1 y 𝑢2 se obtiene:
𝐸𝐴 1 −1 𝑢1 𝐹
∇Π(𝑢1 , 𝑢2 ) = − 1 (1.37)
𝓁 [−1 1 ] [𝑢2 ] [𝐹2 ]
por lo que la condición de punto crítico es
𝐸𝐴 1 −1 𝑢1 𝐹
= 1 (1.38)
𝓁 [−1 1 ] [𝑢2 ] [𝐹2 ]
1.4 Ejercicios
Se pide:
a) trazar los diagramas de solicitaciones (cortante y momento) de la estructura,
b) dimensionar la mínima sección admisible, considerando que 𝑏 debe ser múltiplo de 5 cm,
c) bosquejar la deformada de la estructura y calcular los giros en los apoyos A y B, junto
con el desplazamiento vertical en C, suponiendo que la rigidez flexional 𝐸𝐼 es uniforme y que
𝐸 = 10 GPa.
d) indicar el valor de la directa de todas las barras.
Se pide:
a) calcular las reacciones y trazar los diagramas de solicitaciones,
b) bosquejar la deformada de la estructura,
c) calcular los giros sobre los apoyos y los desplazamientos de C y H en función de 𝐸𝐼 (valor
uniforme) y despreciando la deformación por directa.
Se pide:
a) obtener analíticamente los valores de la carga 𝑃 para los cuales las directas en las barras
𝑎 y 𝑏 son iguales en módulo,
b) verificar el resultado de la parte anterior utilizando alguna herramienta informática,
c) para el valor de 𝑃 positivo obtenido en la parte a), dimensionar una sección transversal for-
mada por 2 perfiles PNC a colocar en todas las barras de la estructura por igual, considerando
que 𝜎𝑎𝑑𝑚 = 140 MPa.
Ejercicio 1.4 Considere la estructura mostrada en la figura, donde la carga dibujada con
un círculo representa una carga movil unitaria.
Se pide:
a) dibujar las líneas de influencia de la reacción en el apoyo C, del cortante en P y del mo-
mento flector en S para la estructura,
b) observando los diagramas obtenidos, indicar que ubicaciones de cargas uniformes 𝑞 = 20
kN/m provocarían máximas y mínimas RC, VP y MS.
Se pide:
a) calcular las reacciones y trazar diagramas de solicitaciones aplicando el método de las
ecuaciones angulares,
b) verificar el resultado de la parte a) utilizando alguna herramienta informática.
Unidad temática 2
Métodos energéticos aplicados a reticulados
En esta unidad temática se presentan los conceptos básicos necesarios para abordar el Mé-
todo de los Desplazamientos (MD) y el Método de las Fuerzas (MF), en ambos casos, aplicados
al análisis de estructuras reticuladas. Se comienza presentando los conceptos de indetermina-
ción cinemática e indeterminación estática para reticulados planos. Luego se presenta breve-
mente el Método de los desplazamientos y finalmente, para realizar el desarrollo del Método
de las Fuerzas, se toma como referencias principal (Reddy, 2002).
𝑔𝑘 = 𝑔𝑙 𝑛𝑁 − 𝑐𝑎 , (2.1)
Por otra parte, si se considera cada nodo de la estructura como un componente aislado se
pueden plantear dos ecuaciones de equilibrio de fuerzas a cada componente, por lo tanto se
tiene:
𝑔ℎ = 4 + 6 − 2 × 4 = 2. (2.3)
Tener una estructura con grado de hiperestacidad 2, es equivalente a decir que si se aplican
las ecuaciones de equilibrio únicamente, tendríamos un sistema compatible indetermina-
do. El subespacio de soluciones del sistema será de dimensión 2.
𝑔ℎ = 0
𝑔ℎ > 0.
En este caso deben ser utilizadas oras ecuaciones como por ejemplo compatibilidad de
desplazamientos.
𝑔ℎ < 0.
𝑔ℎ𝐸 = 𝑛𝑅 − 3 (2.4)
de la estructura 𝑢𝑖 :
𝑛𝑑
𝑢(𝑃) = ∑ 𝑢𝑖 𝜙𝑖 (𝑃), (2.5)
𝑖=1
𝓁𝑒 − 𝑥 𝑒 𝑥𝑒
𝑁1𝑒 (𝑥 𝑒 ) = y 𝑁2𝑒 (𝑥 𝑒 ) = , (2.8)
𝓁𝑒 𝓁𝑒
siendo 𝓁𝑒 el largo de la barra 𝑒.
Usando la relación de deformación-desplazamiento y derivando (respecto a la coordenada
local) se obtiene la expresión de la deformación en el dominio del elemento:
𝜕𝑢𝑒𝐿 𝑒
𝜀 𝑒 (𝑥 𝑒 ) = 𝑒
(𝑥 ) = 𝐵1,𝐿 (𝑥 𝑒 )𝑢𝑒1,𝐿 + 𝐵2,𝐿
𝑒
(𝑥 𝑒 )𝑢𝑒2,𝐿 (2.9)
𝜕𝑥 𝑒
𝑦
𝑣2,𝐿
𝑢2,𝐿
𝑣1,𝐿
𝑢1,𝐿
donde
−1 1
𝑒
𝐵1,𝐿 y 𝐵2,𝐿
(𝑥 𝑒 ) = 𝑒
(𝑥 𝑒 ) = 𝑒 , (2.10)
𝓁𝑒 𝓁
son las derivadas de las funciones de interpolación respecto a la coordenada local. Se observa
que, para esta interpolación, la deformación axial es uniforme en todo el elemento (es decir
que no depende de 𝑥 𝑒 ) y vale
𝑢𝑒 − 𝑢𝑒
𝜀 𝑒 = 2,𝐿 𝑒 1,𝐿 . (2.11)
𝓁
Esto puede ser escrito en forma matricial como:
donde,
1 𝑇
[−1, 0, 1, 0] y 𝐮𝑒𝐿 = [𝑢1,𝐿 , 𝑣1,𝐿 , 𝑢2,𝐿 , 𝑣2,𝐿 ] .
𝐁𝑒𝐿 = (2.13)
𝓁𝑒
Usando la ecuación constitutiva lineal se obtiene la tensión de cada elemento como:
𝜎𝑒 = 𝐸𝑒 𝜀 𝑒 , (2.14)
nente axial. Usando la expresión general dada por la Ecuación (1.28) y la forma del tensor de
deformación dada por la Ecuación (1.32), se obtiene:
1 2
Π𝑒𝑖𝑛𝑡 (𝐮𝑒𝐿 ) = 𝐸𝑒 (𝜀 𝑒 (𝐮𝑒𝐿 )) d𝑉 , (2.15)
2 ∫Ω𝑒
donde Ω𝑒 representa el volumen ocupado por la barra 𝑒. Sustituyendo la relación de la Ecua-
ción (2.12) y operando se obtiene:
1
Π𝑒𝑖𝑛𝑡 (𝐮𝑒𝐿 ) = (𝐮𝑒𝐿 )𝑇 𝐊𝑒𝐿 𝐮𝑒𝐿 , (2.16)
2
donde 𝐊𝑒𝐿 es una matriz dada por:
𝓁𝑒
𝐊𝑒𝐿 = ∫ (𝐁𝑒𝐿 )𝑇 𝐸𝑒 𝐴𝑒 𝐁𝑒𝐿 d𝑥. (2.17)
0
⎡ 1 0 −1 0⎤
𝑒 𝑒 ⎢
𝐸 𝐴 0 0 0 0⎥⎥
𝐊𝑒𝐿 = 𝑒 ⎢ . (2.18)
𝓁 ⎢−1 0 1 0⎥
⎢0 0 0 0⎥⎦
⎣
Para poder resolver problemas con varios elementos de barra, con sus ejes orientados en
diferentes direcciones, es necesario trabajar con un sistema de coordenadas global. Conside-
rando que el sistema de coordenadas local del elemento forma un ángulo 𝛼 𝑒 como se muestra
en la Figura 2.2 (medido antihorario desde el sistema global), la matriz de cambio de base 𝐑𝑒
está dada por:
⎡ 𝑐 −𝑠 0 0 ⎤
⎢ ⎥
𝑠 𝑐 0 0⎥
𝐑𝑒 = ⎢ 𝐮𝑒 = 𝐑𝑒 𝐮𝑒𝐿 (2.19)
⎢0 0 𝑐 −𝑠 ⎥
⎢0 0 𝑠 𝑐 ⎥
⎣ ⎦
donde 𝑐 = cos(𝛼 𝑒 ) y 𝑠 = sin(𝛼 𝑒 ) y 𝐮𝑒 es el vector de desplazamientos del elemento en coor-
denadas globales.
Actividad
Demostrar la identidad dada por la Ecuación 2.19 aplicando la relación de
cambio de base vista en cursos de Álgebra Lineal:
formación se tiene
1
Π𝑒𝑖𝑛𝑡 (𝐮𝑒 ) = (𝐮𝑒 )𝑇 𝐊𝑒 𝐮𝑒 , (2.20)
2
donde la matriz 𝐊𝑒 en coordenadas globales está dada por:
⎡ 𝑐2 𝑐𝑠 −𝑐 2 −𝑐𝑠 ⎤
𝑒 𝑒 ⎢
𝐸 𝐴 𝑐𝑠 𝑠2 −𝑐𝑠 −𝑠 2 ⎥⎥
𝐊𝑒 = 𝐑𝑒 𝐊𝑒𝐿 (𝐑𝑒 )𝑇 = 𝑒 ⎢ 2 . (2.21)
𝓁 ⎢−𝑐 −𝑐𝑠 𝑐2 𝑐𝑠 ⎥
⎢−𝑐𝑠 −𝑠 2 𝑐𝑠 𝑠 2 ⎥⎦
⎣
La energía potencial de deformación de toda la estructura puede ser calculada sumando
las energías de las 𝑛𝑒 barras que forman la estructura:
𝑒 𝑛 𝑒 𝑛
1 1
Π𝑖𝑛𝑡 (𝐮) = ∑ Π𝑒𝑖𝑛𝑡 (𝐮𝑒 ) = ∑ (𝐮𝑒 )𝑇 𝐊𝑒 𝐮𝑒 = 𝐮𝑇 𝐊𝐮 (2.22)
𝑒=1 𝑒=1
2 2
La energía potencial total estará dada por lo tanto, por la suma de ambas energías:
1
Π(𝐮) = 𝐮𝑇 𝐊𝐮 − 𝐮𝑇 𝐟. (2.24)
2
Ecuación de rigidez
Sustituyendo la expresión matricial Π𝑖𝑛𝑡 de la Ecuación (2.22) en las condiciones del Teore-
ma de Castigliano dadas por la Ecuación (2.26), se obtiene el siguiente sistema de ecuaciones
lineales:
𝐊𝐮 = 𝐟 , (2.27)
el cual acompañado de las condiciones de contorno representan el sistema a formar y resolver
para aplicar el Método de los Desplazamientos en su forma matricial.
La matriz 𝐊 es llamada matriz de rigidez global de la estructura ya que vincula despla-
zamientos con fuerzas aplicadas. Es importante destacar que la cantidad de incógnitas del
sistema es igual al grado de indeterminación cinemática 𝑔𝑘 introducido en la Unidad Temá-
tica 1.
Cálculo de solicitaciones
Las solicitaciones de interés en el análisis de reticulados son las fuerzas directa 𝑁 𝑒 en
cada elemento de barra
𝑁 𝑒 = 𝜎 𝑒 𝐴𝑒 . (2.28)
Para ello se pueden utilizar las ecuaciones (1.32), (2.14) y (2.19):
Finalmente se menciona que también es posible analizar reticulados con resortes elásticos
en sus nodos. Para esto se puede adicionar a las ecuaciones los correspondientes términos de
energía potencial elástica. En el caso de un resorte lineal de constante elástica 𝑘 colocado en
el grado de libertad 𝑗, la energía potencial está dada por
1
Π𝑟𝑒𝑠 (𝐮) = 𝑘𝑢2𝑗 , (2.30)
2
que puede ser sumada a la energía potencial de deformación Π𝑖𝑛𝑡 .
encontrar una carpeta UT2 dentro de la carpeta octave. Allí se puede explorar el archivo
FEMTrusS.m el cual hay una implementación del método para reticulados planos
El objetivo del código es mostrar a los estudiantes la aplicación directa de los conceptos
vistos, así como también permitir un espacio para el desarrollo propio incluyendo nuevas
funcionalidades como por ejemplo la capacidad de resolver reticulados tridimensionales.
Este código está en proceso continuo de revisión por lo que, para la resolución de los
ejercicios prácticos, se recomienda usar los otros programas recomendados en el curso. En
caso de detectar errores en el FEMTrusS se agradece notificarlo creando un issue en el sitio
del repositorio.
𝓁 1
3
2
4
1 3
𝓁 𝓁
⎡1 1 −1 −1⎤
9
210 × 10 ⋅ 1 × 10 −3
1 ⎢⎢ 1 1 −1 −1⎥⎥
𝐊1 = √ N/m (2.31)
2 ⋅ 10 2 ⎢−1 −1 1 1⎥
⎢−1 −1 1 1 ⎥⎦
⎣
Operando se tiene
⎡1 1 −1 −1⎤
⎢1 1 −1 −1⎥⎥
𝐊1 = 7.425 × 106 ⎢ N/m (2.32)
⎢−1 −1 1 1⎥
⎢−1 −1 1 1 ⎥⎦
⎣
Para el elemento 2 se tiene
⎡0 0 0 0⎤
⎢0 1 0 −1⎥⎥
𝐊2 = 210 × 106 ⎢ N/m (2.33)
⎢0 0 0 0⎥
⎢0 −1 0 1 ⎥⎦
⎣
Para el elemento 3 se tiene
⎡1 −1 −1 1 ⎤
⎢−1 1 1 −1⎥⎥
𝐊3 = 7.425 × 106 ⎢ N/m (2.34)
⎢−1 1 1 −1⎥
⎢1 −1 ⎥
−1 1 ⎦
⎣
Se observa que el elemento 2 tiene una rigidez mayor (dada por su sección transversal)
aunque esta rigidez solamente vincula los desplazamientos verticales de los nodos 2 y 3.
Considerando que los grados de libertad fijos o de desplazamiento nulo son 1, 2, 5, 6, 7 y
8 y ensamblando las matrices en los grados de libertad restantes (3 y 4) se tiene.
1 1 0 0 1 −1
𝐊 = 7.425 × 106 N/m + 210 × 106 N/m + 7.425 × 106 N/m (2.35)
[1 1] [0 1] [−1 1]
14.84 0 104
106 N/m 𝐮 = N. (2.36)
[ 0 224.84] [ 0 ]
Ahora que están determinados los desplazamientos de todos los nodos de la estructura es
posible determinar la directa 𝑁 𝑒 en todas las barras usando la Ecuación (2.29):
𝑁 1 = −𝑁 3 = 7.07 kN y 𝑁 2 = 0. (2.38)
Las directas también pueden ser obtenidas usando la ecuación de rigidez para cada ele-
mento de barra, dada por:
𝐊𝑒 𝐮𝑒 = 𝐟 𝑒 .
Los resultados obtenidos pueden ser verificados utilizando el código mencionado en la
sección anterior. En la Figura 2.4 se muestra el diagrama de deformada con factor de escala
1000.
Figura 2.4: Gráficos de deformada generado por el código FEMTrusS para el ejemplo presentado usando un factor
de escala 1000.
Figura 2.5: Gráficos de directas generado por el código FEMTrusS para el ejemplo presentado.
14.84 × 106 + 𝑘 0
𝐊= (2.39)
[ 0 224.84 × 106 ]
Actividad
Enumerar las 9 incógnitas de fuerzas (solicitaciones internas y reacciones)
del Ejemplo de la Sección 2.2.5, aplicar las ecuaciones de equilibrio y despe-
jar todas las incógnitas en función de la directa de la barra 3 (𝑁 3 ). Escribir
el conjunto de los vectores de ℝ9 que cumplen el equilibrio como una suma
de vectores dependientes de 𝑃 y 𝑁 3 .
𝑃 𝑃 𝐸0 𝐸1 𝐸2
1
𝑋1 𝑋2
1
= + +
Figura 2.6: Ejemplo de descomposición de estados canónicos para aplicación del MF.
donde el subíndice 𝑗 representa los valores del estado 𝑗, 𝑋𝑗 representa un factor a determinar,
𝐍0 corresponde a las directas para el estado de cargas 𝐟0 y 𝐍𝑗 corresponde a las directas para
el estado de cargas unitarias 𝐟𝑗 .
Esto mismo también es realizado para calcular las reacciones de la estructura:
𝑔ℎ
R(𝐗) = R0 + ∑ 𝑋𝑗 R𝑗 (2.42)
𝑗=1
Si se considera que las directas 𝐍 y las reacciones 𝐑 son las incógnitas principales del
método entonces, se puede decir que los estados canónicos permiten interpolar los posibles
estados tensionales. Los estádos canónicos en el MF cumplen el rol de las funciones de inter-
polación de desplazamientos en el MD, mientras que los valores de las incógnitas 𝐗 cumplen
el rol de los desplazamientos nodales.
1 (𝑁 𝑒 (𝐟))2 𝑒
(Π∗𝑖𝑛𝑡 )𝑒 (𝐟) = 𝓁. (2.45)
2 𝐸𝑒 𝐴𝑒
donde las fuerzas externas 𝐟 están relacionadas con la directas a través de las ecuaciones de
equilibrio. La energía de deformación complementaria de la estructura completa está dada
por la suma:
𝑛
1 𝑒 (𝑁 𝑒 (𝐟))2 𝑒
Π∗𝑖𝑛𝑡 (𝐟) = ∑ 𝓁. (2.46)
2 𝑒=1 𝐸𝑒 𝐴𝑒
Para estructuras hiperestáticas, se puede considerar la relación entre 𝐟 y 𝐗 y considerar
el vector 𝐗 como variable:
1 (𝑁 𝑒 (𝐗))2 𝑒
(Π∗𝑖𝑛𝑡 )𝑒 (𝐗) = 𝓁. (2.47)
2 𝐸𝑒 𝐴𝑒
𝑛 𝑔ℎ 𝑔ℎ
1 𝑒 (𝑁0𝑒 + ∑𝑗=1 𝑋𝑗 𝑁𝑗𝑒 )2 𝑒
Π∗ (𝐗) = ∑ 𝓁 − 𝐮𝑇 𝐟0 − 𝐮𝑇 ∑ 𝑋𝑗 𝐟𝑗 . (2.49)
2 𝑒=1 𝐸 𝑒 𝐴𝑒 ( 𝑗=1 )
𝜕Π∗ 𝜕Π∗𝑖𝑛𝑡
(𝐗) = 0 𝑖 = 1, … , 𝑔ℎ ⇒ (𝐗) = 𝐮𝑇 𝐟𝑖 𝑖 = 1, … , 𝑔ℎ. (2.51)
𝜕𝑋𝑖 𝜕𝑋𝑖
Considerando que los valores de fuerzas externas de los estados canónicos distintos al fun-
damental son unitarios se tiene:
𝜕Π∗𝑖𝑛𝑡
(𝐗) = 𝛿𝑖 𝑖 = 1, … , 𝑔ℎ (2.52)
𝜕𝑋𝑖
𝜕Π∗𝑖𝑛𝑡
(𝐗) = 0 𝑖 = 1, … , 𝑔ℎ. (2.53)
𝜕𝑋𝑖
𝐌𝑓 𝐗 = 𝐛𝑓 (2.56)
3. calcular las directas y reacciones para cada uno de los estados canónicos,
Cálculo de desplazamiento
El segundo teorema de Castigliano puede ser aplicado también para calcular el despla-
zamiento de cualquier punto de la estructura. Dado un vector de fuerza 𝐗𝑠𝑜𝑙 obtenido como
solución de la ecuación de flexibilidad, se considera un estado canónico adicional 𝐸𝑔ℎ+1 con
una fuerza externa unitaria según el desplazamiento deseado. Se aplica Castigliano y a través
de la derivada de la energía respecto el valor correspondiente de 𝑋 y considerando que este
valor 𝑋 debe ser nulo se obtiene el desplazamiento deseado.
𝜕Π∗𝑖𝑛𝑡
𝛿𝑔ℎ+1 = (𝐗 = 𝐗𝑠𝑜𝑙 , 𝑋𝑔ℎ+1 = 0) (2.59)
𝜕𝑋𝑔ℎ+1
Existe un caso particular destacable. Si en el estado canónico 𝐸0 existe una única fuerza
aplicada de magnitud 𝑃, y además el estado canónico 𝐸𝑔ℎ+1 , necesario para calcular el des-
plazamiento, cumple que 𝐸0 = 𝑃𝐸𝑔ℎ+1 , entonces se puede probar que
𝜕Π∗𝑖𝑛𝑡
𝛿𝑔ℎ+1 = (𝐗 = 𝐗𝑠𝑜𝑙 , 𝑋𝑔ℎ+1 = 0). (2.60)
𝜕𝑃
Elemento 𝑁 𝑁1 𝓁/𝐸𝐴
√0
1 2𝑃 1
√ 6.73 × 10−8
2 -𝑃 - 2 4.76 × 10−9
3 0 1 6.73 × 10−8
𝑃 2 𝑃
𝑋1 1
=
1 3
+
1
2
1 3 4
Las directas en las barras son calculadas usando la ecuación Ecuación (2.41):
𝑁 1 = −𝑁 3 = 7.07 kN y 𝑁 2 = 0. (2.62)
2.4 Ejercicios
Ejercicio 2.1
Considere el reticulado de la figura donde todas las barras están construidas con barras
de acero (𝐸 = 210 GPa) de sección tubular con 𝜙𝑒𝑥𝑡 = 100 mm y espesor 𝑡 = 5 mm.
Ejercicio 2.2
En la estructura reticulada mostrada en la figura cada barra está formada por 2 perfiles
PNC10 de acero (𝐸 = 210 GPa).
Se pide:
a) Calcular las reacciones, los esfuerzos en todas las barras y el desplazamiento del nodo D
mediante el método de las fuerzas.
Ejercicio 2.3
En la estructura reticulada mostrada en la figura todas las barras tienen la misma sección
y están formadas por el mismo material.
Ejercicio 2.4
En la estructura reticulada mostrada en la figura, todas las barras tienen la misma sección
y el mismo material.
Ejercicio 2.5 Sea el reticulado y las cargas que se muestran en la figura, donde todos los
triángulos son equiláteros y las barras tienen todas mismo módulo de elasticidad y área.
Se pide:
a) Calcular la solicitación de cada barra utilizando el método de las fuerzas.
b) Dimensionar la estructura con perfiles PNI, suponiendo que todas las barras son iguales,
para 𝑃 = 300 kN, 𝑎 = 3 m y 𝜎𝑎𝑑𝑚 = 140 MPa.
Unidad temática 3
Métodos de Desplazamientos para pórticos
Hipótesis
Se considera un sólido con una geometría como la mostrada en la Figura 1.5, es decir una
viga. Se considera que el campo de desplazamientos de los puntos de la viga, están dados por
el vector formado por las funciones: 𝑢(𝑥, 𝑦, 𝑧), 𝑣(𝑥, 𝑦, 𝑧) y 𝑤(𝑥, 𝑦, 𝑧), representando desplaza-
mientos en 𝑥, 𝑦 y 𝑧 respectivamente. Asumiendo que la flexión se produce en el plano 𝑥 − 𝑦,
se consideran las siguientes hipótesis:
1. Los desplazamientos transversales (flecha) de todos los puntos en una sección trans-
versal (ubicada en la coordenada 𝑥) son pequeños e iguales al desplazamiento del eje
de la viga.
𝑣(𝑥, 𝑦, 𝑧) = 𝑣(𝑥).
Tensor de deformaciones
Se comienza calculando las componentes del tensor de deformaciones aplicando la re-
lación desplazamientos-deformación al campo de desplazamientos considerado (ver Ecua-
𝜕𝜃
𝜀𝑥 (𝑥, 𝑦) = 𝜀𝐺 (𝑥) − 𝑦 (𝑥). (3.2)
𝜕𝑥
donde 𝜀𝐺 es la deformación axial del eje de la viga, y está dada por
𝜕𝑢𝐺
𝜀𝐺 (𝑥) = (𝑥). (3.3)
𝜕𝑥
La distorsión angular 𝛾𝑥𝑦 está dada por:
𝜕𝑣
𝛾𝑥𝑦 (𝑥, 𝑦) = −𝜃(𝑥) + (𝑥). (3.4)
𝜕𝑥
La no dependencia de 𝛾𝑥𝑦 respecto a 𝑦 implica que la cara permanece plana. Al imponer que la
distorsión angular (asociada a la deformación por cortante) sea nula, se aporta la condición:
𝜕𝑣
𝜃= , (3.5)
𝜕𝑥
𝜕2 𝑣
𝜀𝑥 (𝑥, 𝑦) = 𝜀𝐺 (𝑥) − 𝑦 (𝑥). (3.6)
𝜕𝑥 2
lo que representa una extensión de la expresión obtenida en la Ecuación (1.15) para el caso
con deformación axial.
Actividad
Demostrar que, considerando las hipótesis de desplazamiento mencionadas,
la única componente del tensor de deformaciones no nula, es 𝜀𝑥𝑥 .
Tensor de tensiones
Se considera que el material es elástico lineal e isótropo. Se utiliza la correspondiente
ecuación constitutiva y el tensor de deformaciones obtenido, teniendo la relación:
𝜎𝑥 = 𝐸𝜀𝑥 , (3.7)
o también
𝜕2 𝑣
𝜎𝑥 (𝑥, 𝑦) = 𝐸𝜀𝐺 (𝑥) − 𝐸𝑦
(𝑥). (3.8)
𝜕𝑥 2
Esta componente de 𝝈 es la única que produce trabajo interno en la expresión de la energía
de deformación interna dada por la Ecuación (1.28).
𝜕𝑢𝐺 𝜕2 𝑣
𝑁 (𝑥) = ∫ 𝐸(𝑥) − 𝐸𝑦 2 (𝑥) d𝐴. (3.10)
𝐴(𝑥) ( 𝜕𝑥 𝜕𝑥 )
Usando que el vector 𝐞𝑥 pasa por el punto 𝐺, baricentro de la sección transversal, se tiene
Se puede destacar que la convención de signo de directa positiva a tracción es coherente con
el resultado obtenido ya que en dicho caso 𝜀𝐺 > 0.
Para el cálculo del momento según el eje 𝐞𝑧 se considera la suma (integral) del momento
de diferenciales de área por su brazo respectivo:
Sustituyendo las expresiones de las tensiones y la deformación dada por la Ecuación (3.2) se
obtiene:
𝜕2 𝑣
𝑀𝑧 (𝑥) = ∫ 𝑦𝐸𝜀𝐺 (𝐞𝑦 ∧ 𝐞𝑥 ) ⋅ 𝐞𝑧 d𝐴 + ∫ −𝑦 2 (𝑥) (𝐞𝑦 ∧ 𝐞𝑥 ) ⋅ 𝐞𝑧 d𝐴. (3.13)
𝐴(𝑥) 𝐴(𝑥) 𝜕𝑥 2
𝜕2 𝑣
𝑀𝑧 (𝑥) = 𝐸 ∫ 𝑦 2 d𝐴 (𝑥) (3.14)
𝐴(𝑥) 𝜕𝑥 2
usando ahora como hipótesis que la barra tiene sección transversal uniforme 𝐴(𝑥) = 𝐴 y
definiendo el segundo momento de inercia respecto de 𝑧 como 𝐼𝑧 (𝑥) = ∫𝐴 𝑦 2 𝑑𝐴 se obtiene:
𝜕2 𝑣
𝑀𝑧 (𝑥) = 𝐸𝐼𝑧 (𝑥). (3.15)
𝜕𝑥 2
𝑁 (𝑥) 𝜕2 𝑣 𝑀𝑧 (𝑥)
𝜀𝐺 (𝑥) = 2
(𝑥) = (3.16)
𝐸𝐴 𝜕𝑥 𝐸𝐼𝑧
𝑁 (𝑥) 𝑀𝑧 (𝑥)
𝜀𝑥 (𝑥, 𝑦) = −𝑦 , (3.17)
𝐸𝐴 𝐸𝐼𝑧
y multiplicando por 𝐸 ambos miembros y usando la ecuación constitutiva se tiene
𝑁 (𝑥) 𝑀𝑧 (𝑥)
𝜎(𝑥, 𝑦) = −𝑦 (3.18)
𝐴 𝐼𝑧
Convenciones de signo
Para el desarrollo de métodos de análisis de pórticos es útil y necesario definir una con-
vención de signos diferente a la usada para las solicitaciones internas. En la Figura 3.1 se
muestran dos convenciones de signo a ser utilizadas en este documento y en todo el curso.
Convención de signos 1: solicitaciones internas
𝑀1 𝑉1 𝑞 𝑀2
𝑁1 𝑉2
𝑁2
Figura 3.1: Convenciones de signo para momentos nodales y cargas externas en coordenadas locales.
Las solicitaciones internas definidas con la convención de signos usual para 𝜎 (positivo
tracción) corresponden a la convención de signos 1 de la figura. Esta convención corresponde
a aquella en la cual momentos positivos representan una tracción de fibras inferiores.
Por otra parte, la convención 2 es útil para el desarrollo de métodos matriciales o compu-
tacionales, asociada a las fuerzas externas aplicadas.
Ecuaciones de equilibrio
Las ecuaciones para vigas sometidas a cargas transversales 𝑞 y axiales 𝑏 distribuidas por
unidad de longitud están dadas por:
𝑑𝑁
(𝑥) = −𝑏(𝑥) (3.19)
𝑑𝑥
𝑑𝑉
(𝑥) = 𝑞(𝑥) (3.20)
𝑑𝑥
𝜕𝑀
(𝑥) = 𝑉 (𝑥) (3.21)
𝜕𝑥
La deducción de estas ecuaciones a partir del equilibrio de un segmento diferencial fue
realizado en cursos anteriores. Por otra parte, este desarrollo también puede ser realizado a
partir del teorema de trabajo virtual, de forma similar a como es hecho en la sección 5.4.5 de
(Hughes, 1987).
Figura 3.2: Esquemas de viga: (a) eje de viga en configuración de referencia (línea punteada) y eje de viga deformada
(línea continua), (b) eje de viga de referencia con fuerzas aplicadas.
Por ahora se asume que no hay cargas aplicadas en el tramo intermedio de la viga, ni
distribuidas ni puntuales. Se omitirá el sub-índice 𝑧 en 𝑀𝑧 y 𝜃𝑧 .
El desplazamiento en cada punto de la viga está dado por las funciones 𝑢(𝑥, 𝑦) y 𝑣(𝑥, 𝑦).
Para deducir las relaciones esfuerzo desplazamiento se utilizará el primer teorema de Cas-
tigliano y se comienza planteando la expresión de la energía potencial de deformación en
función de las funciones de desplazamiento.
Desarrollando se tiene:
1
Π𝑖𝑛𝑡 (𝑢, 𝑣) = 𝐸 (𝜀𝐺 (𝑥))2 d𝐴d𝑥
2 ∫𝓁 ∫𝐴
𝜕2 𝑣
… − ∫ ∫ 𝐸𝑦𝜀𝐺 (𝑥) 𝜕𝑥 2 (𝑥)d𝐴d𝑥
𝓁 𝐴
2
1 𝜕2 𝑣
… + ∫ ∫ 𝐸𝑦 2 (𝑥) d𝐴d𝑥. (3.23)
2 𝓁 𝐴 ( 𝜕𝑥 2 )
Usando que el origen del sistema de coordenadas es el baricentro de la sección y siendo que
la sección transversal es uniforme, se tiene que:
2
1 1 𝜕2 𝑣
Π𝑖𝑛𝑡 (𝑢, 𝑣) = 𝐸 (𝜀 (𝑥)) 2
d𝐴d𝑥 + 𝐸𝑦 2
(𝑥) d𝐴d𝑥 . (3.24)
2 ∫𝓁 ∫𝐴 2 ∫𝓁 ∫𝐴
𝐺
( 𝜕𝑥 2 )
⏟⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏟⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏟ ⏟⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏟⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏞⏟
axial flexión
En esta expresión fueron definidos los términos de energía de deformación por esfuerzo
axial (dependiente de 𝑢)y por esfuerzo de flexión (dependiente de 𝑣). Esto permite desacoplar
los efectos y hacer desarrollos independientes, por lo que a continuación se presentan expre-
siones de las energías de deformación axial y flexional, asumiendo cada efecto por separado.
1 𝐸𝐴
Π𝑖𝑛𝑡,axial (𝑢1 , 𝑢2 ) = (𝑢2 − 𝑢1 )2 . (3.26)
2 𝓁
𝑣(𝑥) = 𝜑𝑣1 (𝑥)𝑣1 + 𝜑𝜃1 (𝑥)𝜃1 + 𝜑𝑣2 (𝑥)𝑣2 + 𝜑𝜃2 (𝑥)𝜃2 , (3.29)
3𝑥 2 2𝑥 3
𝜑𝑣1 (𝑥) = 1− + 3 (3.30)
( 𝓁2 𝓁 )
2
2𝑥 𝑥3
𝜑𝜃1 (𝑥) = 𝑥− + 2 (3.31)
( 𝓁 𝓁 )
2 3
3𝑥 2𝑥
𝜑𝑣2 (𝑥) = − 3 (3.32)
( 𝓁2 𝓁 )
2
𝑥 𝑥3
𝜑𝜃2 (𝑥) = − + 2 . (3.33)
( 𝓁 𝓁 )
En la Figura 3.3 se muestran los gráficos de las funciones de interpolación, los cuales
corresponden a las funciones de las elásticas para desplazamientos o giros unitarios de cada
grado de libertad correspondiente.
Cualquier elástica 𝑣 dada por fuerzas y momentos nodales puede ser escrita como com-
binación de estas cuatro funciones de interpolación considerando como factores los valores
de flecha y giro como se muestra en la Ecuación (3.29).
Considerando esta interpolación en el término de energía de deformación a flexión de la
Ecuación (3.24), y utilizando que la sección transversal es uniforme en el elemento de viga se
obtiene
2
1 𝜕2 𝑣
Π𝑖𝑛𝑡,flexión (𝑣1 , 𝜃1 , 𝑣2 , 𝜃2 ) = 𝐸𝐼 ∫ 2
(𝑥, 𝑣1 , 𝜃1 , 𝑣2 , 𝜃2 ) d𝑥. (3.34)
2 𝓁 ( 𝜕𝑥 )
(a) Funciones de interpolación 𝜑𝑣1 y 𝜑𝑣2 . (b) Funciones de interpolación 𝜑𝜃1 y 𝜑𝜃2 .
𝐮 = [𝑢1 , 𝑣1 , 𝜃1 , 𝑢2 , 𝑣2 , 𝜃2 ]𝑇 . (3.37)
y aplicar las seis condiciones dadas por el primer teorema de Castigliano. Estas condiciones
se pueden agrupar en, por un lado las asociadas a los grados de libertad 𝑢1 y 𝑢2 :
𝜕Π𝑖𝑛𝑡,axial
(𝑢1 , 𝑢2 ) = 𝐹𝑥,1 (3.39)
𝜕𝑢1
𝜕Π𝑖𝑛𝑡,axial
(𝑢1 , 𝑢2 ) = 𝐹𝑥,2 (3.40)
𝜕𝑢2
𝜕Π𝑖𝑛𝑡,flexión
(𝑣1 , 𝜃1 , 𝑣2 , 𝜃2 ) = 𝐹𝑦,1 (3.41)
𝜕𝑣1
𝜕Π𝑖𝑛𝑡,flexión
(𝑣1 , 𝜃1 , 𝑣2 , 𝜃2 ) = 𝑀1 (3.42)
𝜕𝜃1
𝜕Π𝑖𝑛𝑡,flexión
(𝑣1 , 𝜃1 , 𝑣2 , 𝜃2 ) = 𝐹𝑦,2 (3.43)
𝜕𝑣2
𝜕Π𝑖𝑛𝑡,flexión
(𝑣1 , 𝜃1 , 𝑣2 , 𝜃2 ) = 𝑀2 . (3.44)
𝜕𝜃2
Las condiciones asociadas a 𝑢1 y 𝑢2 son equivalentes a relaciones ya vistas, por ejemplo
en la Ecuación (1.38).
Las condiciones asociadas a las otras variables requieren un desarrollo. A modo de ejem-
plo se presenta el desarrollo para la condición del momento del nodo izquierdo, es decir, la
Ecuación (3.42). Usando la interpolación de la flecha y calculando la derivada, se tiene:
𝓁
𝜕Π𝑖𝑛𝑡 1 𝜕2 𝑣
= 𝐸𝐼 ∫ 2 (𝑣 , 𝜃 , 𝑣 , 𝜃 )𝜑′′ d𝑥 (3.45)
𝜕𝜃1 2 0 ( 𝜕𝑥 2 1 1 2 2 𝜃1 )
para evaluar la derivada de 𝑣 se utilizan las expresiones de las derivadas de las funciones de
interpolación dadas por:
6 12𝑥 4 6𝑥
𝜑𝑣′′1 (𝑥) = − + 3 𝜑𝜃′′1 (𝑥) =− + 2 (3.46)
( 𝓁2 𝓁 ) ( 𝓁 𝓁 )
6 12𝑥 2 6𝑥
𝜑𝑣′′2 (𝑥) = − 3 𝜑𝜃′′2 (𝑥) = − + 2 . (3.47)
( 𝓁2 𝓁 ) ( 𝓁 𝓁 )
𝜕Π𝑖𝑛𝑡
= 𝐾𝑣1 ,𝜃1 𝑣1 + 𝐾𝜃1 ,𝜃1 𝜃1 + 𝐾𝑣2 ,𝜃1 𝑣2 + 𝐾𝜃2 ,𝜃1 𝜃2 (3.48)
𝜕𝜃1
donde los coeficientes 𝐾 están dados por las integrales:
𝓁 𝓁
𝐾𝑣1 ,𝜃1 = 𝐸𝐼 ∫ 𝜑𝑣′′1 𝜑𝜃′′1 d𝑥 𝐾𝜃1 ,𝜃1 = 𝐸𝐼 ∫ 𝜑𝜃′′1 𝜑𝜃′′1 d𝑥 (3.49)
0 0
𝓁 𝓁
𝐾𝑣2 ,𝜃1 = 𝐸𝐼 ∫ 𝜑𝑣′′2 𝜑𝜃′′1 d𝑥 𝐾𝜃2 ,𝜃1 = 𝐸𝐼 ∫ 𝜑𝜃′′2 𝜑𝜃′′1 d𝑥 (3.50)
0 0
24 84 1 72 1 6
𝐾𝑣1 ,𝜃1 = 𝐸𝐼 − + = 𝐸𝐼 2 (3.51)
( 𝓁2 2 𝓁 2 3 𝓁 2 ) 𝓁
16 48 1 36 1 4
𝐾𝜃1 ,𝜃1 = 𝐸𝐼 2
− 2
+ 2
𝓁 = 𝐸𝐼 (3.52)
(𝓁 2 𝓁 3 𝓁 ) 𝓁
24 84 1 72 1 6
𝐾𝑣2 ,𝜃1 = 𝐸𝐼 − 2 + 2
− 2
= −𝐸𝐼 2 (3.53)
( 𝓁 2 𝓁 3 𝓁 ) 𝓁
8 36 1 36 1 2
𝐾𝜃2 ,𝜃1 = 𝐸𝐼 − + 𝓁 = 𝐸𝐼 (3.54)
( 𝓁2 2 𝓁2 3 𝓁2 ) 𝓁
Se obtiene por lo tanto que la ecuación de Castigliano para el giro del primer nodo, dada
por la Ecuación (3.42), puede ser escrita como:
6 4 6 2
𝐸𝐼 2
𝑣1 + 𝐸𝐼 𝜃1 − 𝐸𝐼 2 𝑣2 + 𝐸𝐼 𝜃2 = 𝑀1 . (3.55)
𝓁 𝓁 𝓁 𝓁
Repitiendo el procedimiento para las otras condiciones de Castigliano para la flexión se
llega a:
12 6 12 6
𝐸𝐼 𝑣 + 𝜃 − 𝑣 + 𝜃 = 𝐹𝑦,1 (3.56)
( 𝓁3 1 𝓁2 1 𝓁3 2 𝓁2 2 )
6 4 6 2
𝐸𝐼 𝑣 + 𝜃 − 𝑣 + 𝜃 = 𝑀1 (3.57)
( 𝓁2 1 𝓁 1 𝓁2 2 𝓁 2 )
12 6 12 6
𝐸𝐼 − 3 𝑣1 − 2 𝜃1 + 3 𝑣2 − 2 𝜃2 = 𝐹𝑦,2 (3.58)
( 𝓁 𝓁 𝓁 𝓁 )
6 2 6 4
𝐸𝐼 𝑣 + 𝜃 − 𝑣 + 𝜃 = 𝑀2 , (3.59)
( 𝓁2 1 𝓁 1 𝓁2 2 𝓁 2 )
Las seis condiciones de Castigliano establecen una relación lineal entre las fuerzas 𝐟 y los
desplazamientos 𝐮, la cual puede ser escrita en forma matricial como:
𝐊𝐮 = 𝐟, (3.60)
⎡ 𝐸𝐴 0 0 −
𝐸𝐴
0 0 ⎤
⎢ 𝓁 𝓁 ⎥
⎢ 𝐸𝐼 𝐸𝐼 𝐸𝐼 𝐸𝐼 ⎥
⎢ 0 12 6 0 −12 6 2 ⎥
⎢ 𝓁3 𝓁2 𝓁3 𝓁 ⎥
⎢ 𝐸𝐼 𝐸𝐼 𝐸𝐼 𝐸𝐼 ⎥
⎢ 0 6 2 4 0 −6 2 2 ⎥
⎢ 𝓁 𝓁 𝓁 𝓁 ⎥
𝐊 = ⎢ 𝐸𝐴 𝐸𝐴 ⎥. (3.61)
⎢− 0 0 0 0 ⎥
⎢ 𝓁 𝓁 ⎥
⎢ 𝐸𝐼 𝐸𝐼 𝐸𝐼 𝐸𝐼 ⎥
⎢ 0 −12 −6 0 12 −6 2 ⎥
⎢ 𝓁3 𝓁2 𝓁3 𝓁 ⎥
⎢ 𝐸𝐼 𝐸𝐼 𝐸𝐼 𝐸𝐼 ⎥
⎢ 0 6 2 2 0 −6 2 4 ⎥
⎣ 𝓁 𝓁 𝓁 𝓁 ⎦
Es importante destacar que la matriz de rigidez es simétrica, lo que puede ser justificado
por el hecho de que la función de energía potencial total es una función de tipo 𝐶 2 y que la
matríz 𝐊 es la hesiana de dicha función.
Por otra parte también se destaca que la matriz 𝐊 tiene valores propios nulos con sus
correspondientes vectores propios asociados. (pertenecientes al subespacio asociado al va-
lor propio cero 𝑆0 ). Esto quiere decir que existen movimientos (dados por los vectores de
desplazamientos en dicho subespacio) que se corresponden con fuerzas nodales nulas, por
lo tanto estos movimientos no introducen deformaciones o tensiones en el elemento. Dicho
de otra forma, estos movimientos serán los movimientos de cuerpo rígido de la viga. Estos
movimientos son 3 por lo que es necesario eliminar tres grados de libertad (a través de las
condiciones de contorno) para obtener una matriz reducida con valores propios positivos y
por lo tanto invertible. De esta forma la solución del sistema es única.
Este método permite obtener valores exactos de giros y flechas nodales, sin embargo, es
necesario integrar las ecuaciones de equilibrio para obtener los diagramas de solicitaciones
y elástica en el interior del tramo.
A continuación se consideran dos casos importantes de cargas en el tramo, con un enfo-
que como el presentado en (Oñate, 2013) donde se puede encontrar un desarrollo más com-
pleto.
Carga distribuida
En el caso de carga distribuida 𝑞(𝑥) (de acuerdo a la convención 2), se tiene que la energía
potencial de esta fuerza distribuida está dada por:
𝓁
𝑞
Π𝑒𝑥𝑡 = − ∫ 𝑞(𝑥)𝑣(𝑥)d𝑥. (3.62)
0
𝑒𝑞 𝑞𝓁
𝐹𝑦,1 = (3.64)
2
𝑒𝑞 𝑞𝓁2
𝑀1 = (3.65)
12
𝑒𝑞 𝑞𝓁
𝐹𝑦,2 = (3.66)
2
𝑒𝑞 𝑞𝓁2
𝑀2 = − (3.67)
12
Carga puntual
Se considera ahora que la carga puntual 𝑃 es aplicada (según la convención 2) en el punto
de coordenada 𝑥 = 𝑥𝑃 con 𝑥𝑃 ∈ (0, 𝓁) y que 𝑞(𝑥) = 0.
La energía potencial de la fuerza en el tramo está dada por:
Π𝑃𝑒𝑥𝑡 = −𝑃𝑣(𝑥𝑃 ) = −𝑃𝜑𝑣1 (𝑥𝑃 )𝑣1 − 𝑃𝜑𝜃1 (𝑥𝑃 )𝜃1 − 𝑃𝜑𝑣2 (𝑥𝑃 )𝑣2 − 𝑃𝜑𝜃2 (𝑥𝑃 )𝜃2 , (3.68)
por lo tanto las fuerzas nodales equivalentes a la carga puntual en 𝑥𝑃 son calculadas evaluando
las funciones 𝜑. Por ejemplo el momento nodal en 1 es
𝑒𝑞 2𝑥𝑃2 𝑥𝑃3
𝑀1 = 𝑃𝜑𝜃1 (𝑥𝑃 ) = 𝑃 𝑥𝑃 − + 2 (3.69)
( 𝓁 𝓁 )
𝑒𝑞 1
𝑀1 = 𝑃𝑥𝑃 (𝓁 − 𝑥𝑃 )2 . (3.70)
𝓁2
𝑒𝑞 3𝑥 2 2𝑥 3 𝑒𝑞 3𝑥 2 2𝑥 3
𝐹𝑦,1 = 𝑃 1− + 3 𝐹𝑦,2 = 𝑃 − 3 . (3.72)
( 𝓁2 𝓁 ) ( 𝓁2 𝓁 )
Las fuerzas externas equivalentes, tanto para el caso de carga distribuida como el de carga
puntual, son sumadas al término de fuerzas externas nodales, obteniendo una ecuación de
rigidez más general
𝐊𝐮 = 𝐟 + 𝐟 𝑒𝑞 . (3.73)
𝐊𝐿 𝐮𝐿 = 𝐅𝐿 . (3.75)
𝐊𝐺 𝐮 = 𝐅𝐺 . (3.77)
Apoyos elásticos
Para la resolución mediante análisis matricial o de forma computacional, las fuerzas in-
troducidas por los resortes pueden ser consideradas a través de una matriz de rigidez para
cada elemento (en coordenadas globales):
⎡𝑘𝑢1 0 0 0 0 0⎤
⎢ ⎥
⎢0 𝑘𝑣1 0 0 0 0⎥
⎢0 0 𝑘𝜃1 0 0 0⎥
𝑒
𝐊𝑆 = ⎢ (3.78)
⎢0 0 0 𝑘𝑢2 0 0 ⎥⎥
⎢0 0 0 0 𝑘𝑣2 0⎥
⎢0 0 0 0 0 𝑘𝜃2 ⎥
⎣ ⎦
que puede ser ensamblada y sumada a la matriz de rigidez global, obteniendo un sistema dado
por
(𝐊𝐺 + 𝐊𝑆 ) 𝐮 = 𝐅𝐺 . (3.79)
Actividad
Calcular la flecha máxima de una viga de largo 𝓁, biempotrada, con rigi-
dez flexional uniforme 𝐸𝐼 y una carga puntual 𝑃 aplicada en la mitad de su
longitud, usando el método de los desplazamientos en su forma matricial,
utilizando dos elementos.
lo que es equivalente a:
2
1 𝓁 1 𝓁 𝜕2 𝑣
𝐸𝐴𝜀 2
d𝑥 ≪ 𝐸𝐼 d𝑥. (3.81)
2 ∫0 𝐺
2 ∫0 ( 𝜕𝑥 2 )
𝑁 (𝑥) 𝜕2 𝑣 𝑀(𝑥)
𝜀𝐺 (𝑥) = y 2
(𝑥) = , (3.82)
𝐸𝐴 𝜕𝑥 𝐸𝐼
se tiene que esa condición es equivalente a
𝓁 𝓁
𝑁 (𝑥)2 𝑀(𝑥)2
∫ d𝑥 ≪ ∫ d𝑥.. (3.83)
0 𝐸𝐴 0 𝐸𝐼
Esta condición puede ser asumida al inicio del análisis de un pórtico, y verificada luego
de haber calculado los diagramas de directa y momentos flectores. Se puede decir que para
cierto tipo de diagramas, esta condición se cumple si 𝐴 ≫ 𝐼 .
En caso de que esta condición sea asumida, la energía de deformación axial sería despre-
ciada del desarrollo que llevó a la Ecuación (3.60), por lo que solo se tendrían los términos de
flexión y la barra no tendría desplazamiento relativo, es decir 𝑢2 = 𝑢1 .
𝜃2
𝜃1 𝜓
𝑣1 𝑣2
𝑀 𝑒𝑚𝑝 = −𝑀 𝑒𝑞
⎧
⎪ 2𝐸𝐼 𝑒𝑚𝑝
⎪ (2𝜃1 + 𝜃2 − 3𝜓) + 𝑀1 = 𝑀1
⎪ 𝓁
⎨ (3.86)
⎪
⎪ 2𝐸𝐼 𝑒𝑚𝑝
⎪ (𝜃1 + 2𝜃2 − 3𝜓) + 𝑀2 = 𝑀2
⎩ 𝓁
Empotramiento perfecto
Actividad
Sea el pórtico mostrado en la figura, donde cada barra tiene sección uniforme
de inercia 𝐼 , y está formado por un material con módulo de Young 𝐸.
3m 6m
B C
1350 kN
6m
Por otra parte a partir de combinaciones lineales de las ecuaciones de momento y cortante
se obtiene las relaciones:
⎧ 𝑒𝑞 𝑒𝑞
⎪ 𝑒𝑞 𝑀1 + 𝑀2 𝑀1 + 𝑀2
⎪
⎪ 𝐹𝑦,1 + 𝐹𝑦,1 = +
⎪ 𝓁 𝓁
⎨ (3.87)
⎪
⎪ 𝑒𝑞 𝑒𝑞
⎪ 𝑒𝑞 𝑀1 + 𝑀2 𝑀1 + 𝑀2
⎪ 𝐹 + 𝐹𝑦,2 =− −
⎩ 𝑦,2 𝓁 𝓁
Agrupando se puede obtener que esto es equivalente a las ecuaciones:
⎧
⎪ 𝑀1 + 𝑀2 𝑡𝑟−𝑖𝑠𝑜
⎪
⎪ 𝐹𝑦,1 = + 𝐹𝑦,1
⎪ 𝓁
⎨ (3.88)
⎪
⎪ 𝑀1 + 𝑀2
⎪ 𝑡𝑟−𝑖𝑠𝑜
⎪ 𝐹𝑦,2 =− + 𝐹𝑦,2
⎩ 𝓁
donde
𝑒𝑞 𝑒𝑞
𝑡𝑟−𝑖𝑠𝑜 𝑀1 + 𝑀2
𝑒𝑞
𝐹𝑦,1 = −𝐹𝑦,1 +
𝓁
𝑒𝑞 𝑒𝑞
𝑒𝑞 𝑀 + 𝑀2
𝑡𝑟−𝑖𝑠𝑜
𝐹𝑦,2 = −𝐹𝑦,2 − 1
𝓁
donde 𝐹𝑦,1
𝑡𝑟−𝑖𝑠𝑜
y 𝐹𝑦,2
𝑡𝑟−𝑖𝑠𝑜
son las cargas equivalentes correspondientes a las reacciones de una viga
simplemente apoyada (cargas de tramo isostáticas). Esto se desarrollará con mayor detalle en
el momento de aplicación en práctico.
2𝐸𝐼 𝑒𝑞
(𝜃1 + 2𝜃2 − 3𝜓) = 0 + 𝑀2 (3.89)
𝓁
por lo que despejando 𝜃2 se tiene
𝑒𝑞
𝜃1 3 𝑀 𝓁
𝜃2 = − + 𝜓+ 2 . (3.90)
2 2 4𝐸𝐼
Sustituyendo en la ecuación del momento 𝑀1 se tiene
𝑒𝑞
2𝐸𝐼 𝜃1 3 𝑀 𝓁 𝑒𝑞
2𝜃1 + − + 𝜓 + 2 − 3𝜓 = 𝑀1 + 𝑀1 . (3.91)
𝓁 ( ( 2 2 4𝐸𝐼 ) )
3𝐸𝐼 𝑒𝑚𝑝′
(𝜃1 − 𝜓) + 𝑀1 = 𝑀1 (3.93)
𝓁
𝑒𝑞
𝑒𝑚𝑝′ 𝑒𝑞
donde 𝑀1 = −𝑀1 + 𝑀22 es el momento de empotramiento perfecto (reacción) de una viga
empotrada-articulada. Fácilmente obtenible utilizando tablas.
De forma análoga, si se tiene una articulación en el nodo 1, se impone 𝑀1 = 0 y se tiene
por lo que despejando 𝜃2 se tiene
𝑒𝑞
𝜃2 3 𝑀 𝓁
𝜃1 = − + 𝜓+ 1 . (3.94)
2 2 4𝐸𝐼
obteniendo la expresión de momento:
3𝐸𝐼 𝑒𝑚𝑝′
(𝜃2 − 𝜓) + 𝑀2 = 𝑀2 (3.95)
𝓁
1. Ecuaciones de momentos nodales por barra: dadas por las Ecuaciones (3.86) para barras
con nodos rígidos en ambos extremos y la Ecuación (3.92) para barras con un extremo
articulado.
4. Equilibrios de equilibrio de pisos. Para cada piso del pórtico se debe verificar el equili-
brio de fuerzas horizontales o cortantes. Esto se realiza aislando el elemento o conjunto
de elementos horizontales correspondiente de la estructura, considerando las acciones
realizadas por las otras barras y planteando el equilibrio de fuerzas.
3.3 Ejemplo
Sea la estructura mostrada en la Figura 3.5 donde las barras tienen módulo de Young 𝐸
y sección transversal de Inercia 𝐼 . Se desea determinar desplazamientos y giros nodales. El
procedimiento de resolución es presentado de forma esquemática a continuación.
La estructura es desplazable y el desplazamiento horizontal de 𝐵 es una incógnita a de-
terminar. Se utilizará como variable el desplazamiento del piso 𝐵 − 𝐶, dado por Δ𝐵𝐶 = 𝜓𝐴𝐵 𝓁𝐴𝐵 ,
por lo que de acuerdo con la convención utilizada Δ𝐵𝐶 es positivo para desplazamientos del
piso hacia la izquierda. Las otras variables a determinar son los giros 𝜃𝐵 y 𝜃𝐶 . Se utilizarán en
todas las ecuaciones las condiciones de contorno: desplazamientos y giros nulos en A y D.
Usando la ecuación de momento en B en la barra BA se tiene:
2𝐸𝐼 Δ𝐵𝐶
𝑀𝐵𝐴 = (2𝜃𝐵 − 3 ). (3.97)
6 6
3m 6m
B C
1350 kN
6m
9m
1 1 1 3 4
2𝐸𝐼 2 + 𝜃 + 𝜃 − Δ = − 𝑃. (3.99)
( ( 6 9 ) 𝐵 9 𝐶 62 𝐵𝐶 ) 3
Esta consiste en la primer ecuación que relaciona las tres incógnitas a determinar.
Por otra parte se realiza el mismo procedimiento para el nodo C. Para la barra BC se tiene
2𝐸𝐼 𝑃 ⋅ 32 ⋅ 6
𝑀𝐶𝐵 = (2𝜃𝐶 + 𝜃𝐵 ) − (3.100)
9 92
y para la barra CD
2𝐸𝐼 3
𝑀𝐶𝐷 =
(2𝜃𝐶 − Δ𝐵𝐶 ) (3.101)
9 9
Usando el equilibrio de momentos en C y simplificando se tiene:
1 4 3 2
2𝐸𝐼 𝜃 + 𝜃 − Δ = 𝑃. (3.102)
( 9 𝐵 9 𝐶 92 𝐵𝐶 ) 3
⎡ 𝜃𝐵 ⎤ ⎡−4/3⎤
𝐊 ⎢ 𝜃𝐶 ⎥ = 𝑃 ⎢ 2/3 ⎥ (3.105)
⎢ ⎥ ⎢ ⎥
⎣Δ𝐵𝐶 ⎦ ⎣ 0 ⎦
donde 𝐊 es la matriz del sistema dada por:
1 1 1
⎡2 ( 6 + 9 ) 9 − 632 ⎤
𝐊 = 2𝐸𝐼 ⎢ 1
9
4
9 − 932 ⎥ (3.106)
⎢ 3 1 ⎥
⎣ − 62 − 932 6 ( 613 + 93 )⎦
Esta matriz es la llamada matriz de rigidez. Previo a la resolución del sistema es impor-
tante verificar que la misma es simétrica. Para resolver se pueden ejecutar los siguientes
comandos en Octave:
Obteniendo la solución:
−1.90385𝑃 0.93930𝑃 −3.43990𝑃
𝜃𝐵 = , 𝜃𝐶 = y Δ𝐵𝐶 = (3.107)
𝐸𝐼 𝐸𝐼 𝐸𝐼
Usando Ftool se verifican los resultados. En la Figura 3.6 se muestran los diagramas de
momento y cortante.
En la Figura 3.7 se muestran los diagramas de deformada y directas.
3.4 Ejercicios
En cada ejercicio o estructura descrita a continuación, se pide:
Ejercicio 3.1
Considere la viga mostrada en la figura, compuesta por madera (𝐸 = 10 GPa) y con una
sección transversal rectangular de 20 cm de ancho y 40 cm de altura.
50 kN
10 kN/m
B C
A
4m 2m 3m
Ejercicio 3.2
La estructura de acero de la figura (𝐸 = 210 GPa), compuesta por una viga (ABC) formada
por 1 PNI28 y un pilar (BD) formada por un perfil PNI22.
20 kN/m
B C
A
3m
6m 6m
Ejercicio 3.3
La estructura de acero (𝐸 = 210 GPa) compuesta por una viga (ABC), dos pilares (AD y
BE) y un tensor (BD), cuya sección es 1PNI18.
Ejercicio 3.4
En la estructura de hormigón de la figura (E=25 GPa) los pilares son de 20x30 cm mientras
que la viga es de 15x40 cm.
Ejercicio 3.5
En la estructura de acero de la figura (E=210 GPa), las barras BD, CF y DE están confor-
madas por 2PNC24 soldados ([]), mientras que la barra AB está compuesta por 1PNI20.
Ejercicio 3.6
En la estructura de hormigón de la figura (E=30 GPa) las secciones son rectangulares de
12x30 cm.
Ejercicio 3.7
Sean las siguientes estructuras conformadas por barras de EI=cte y de longitud L so-
metidas a un desplazamiento impuesto Δ. Se pide calcular los momentos flectores en los
empotramientos.
Δ Δ
a) b)
Ejercicio 3.8
Sea la viga continua que se muestra en la siguiente figura, compuesta por un material
con módulo de elasticidad 𝐄=29 GPa y sección rectangular 20x60 cm2 . La viga se encuentra
sometida a una carga uniformemente distribuida 𝑞= 30 kN/m y los apoyos 𝐁, 𝐂 y 𝐃 tienen un
asentamiento 𝛿𝐵 =1.5 cm, 𝛿𝐶 =4.0 cm y 𝛿𝐷 =1.9 cm respectivamente. Se pide calcular los giros
correspondientes a la mínima cantidad de incógnitas del problema.
Ejercicio 3.9
Sea el pórtico que se muestra en la siguiente figura, compuesto por barras de rigidez
flexional EI uniforme. ´
C
D
30 kN
4.8 m
A
B
3.6 m 6m
Se pide:
Unidad temática 4
Simetría y Líneas de Influencia en pórticos
En esta unidad se presentan dos conceptos o métodos específicos del análisis de pórticos
planos. Por una parte simetría y antisimetría, con su aplicación al análisis simplificado de
pórticos planos. Por otra parte se presenta una introducción a la aplicación de líneas de in-
fluencia en pórticos hiperestáticos. Para la sección de simetría se utiliza un enfoque basado en
las ecuaciones de equilibrio para obtener las condiciones mecánicas y cinemáticas dadas por
la simetría, aplicándolas luego en pórticos de forma similar a como es realizado en (Cervera
Ruiz and Blanco Díaz, 2002b). La sección de Líneas de Influencia está basada en (Celigüeta,
2003).
𝑞
𝐹𝑦 𝐹𝑦
𝑀 𝑀
𝐹𝑥 𝐹𝑥
𝑥
Figura 4.1: Esquema de viga sometiga a cargas simétricas.
Las cargas nodales pueden ser provocadas por agentes externos o bien deberse a reaccio-
nes asociadas a vínculos que restringen desplazamientos nodales. Si las cargas y las restric-
ciones en desplazamiento (vínculos a tierra) son simétricos entonces las reacciones asociadas
a dichos vínculos también lo serán.
Se desea obtener cuales son las propiedades que cumplen las solicitaciones y los despla-
zamientos para este caso. Para esto se plantean las expresiones matemáticas de cada una de
las hipótesis y ecuaciones disponibles y se procede a realizar el desarrollo.
𝓁 𝓁
𝑉 −= 𝐹𝑦 𝑉 = −𝐹𝑦 (4.1)
( 2) (2)
𝓁 𝓁
𝑀 − = 𝑀𝑛 𝑀 = 𝑀𝑛 . (4.2)
( 2) (2)
Por otra parte, la simetría de la carga distribuida es equivalente a decir que la función de
la carga 𝑞(𝑥) ∶ [− 2𝓁 , 2𝓁 ] → ℝ es una función par. La definición de función par establece que:
𝓁 𝓁
𝑞(𝑥) = 𝑞(−𝑥) ∀𝑥 ∈ − , . (4.3)
[ 2 2]
𝑑𝑉
= 𝑞 (4.4)
𝑑𝑥
𝑑𝑀
= 𝑉 (4.5)
𝑑𝑥
donde nuevamente se recuerda que se utiliza la convención de signos 1 de la UT3. También se
cuenta con las relaciones fuerza desplazamiento obtenidas en la UT3 a partir de la ecuación
constitutiva y las hipótesis cinemáticas:
𝑑𝜃 𝑀
= (4.6)
𝑑𝑥 𝐸𝐼
𝑑𝑣
= 𝜃 (4.7)
𝑑𝑥
Equilibrio global Finalmente se establece la condición del equilibrio global para cualquier
segmento de viga considerado. En particular se consideran las dos mitades de la viga por lo
que se realiza un corte en 𝑥 = 0. El equilibrio de fuerzas verticales de la mitad derecha está
dado por:
𝓁/2
𝑉 (0) + ∫ 𝑞(𝑢)𝑑𝑢 + 𝐹𝑦 = 0 (4.8)
0
donde 𝑉 (0) es el cortante en la sección del plano de simetría., Por otra parte también se tiene
la ecuación del equilibrio de fuerzas de la mitad izquierda, dado por:
0
−𝑉 (0) + ∫ 𝑞(𝑢)𝑑𝑢 + 𝐹𝑦 = 0. (4.9)
−𝓁/2
Desarrollo de cortante Se pasa a buscar cual es la forma del cortante 𝑉 (𝑥) bajo las hipó-
tesis consideradas. El cortante puede ser calculado integrando la primer ecuación puntual de
equilibrio:
𝑥
𝑉 (𝑥) = ∫ 𝑞(𝑢)d𝑢 + 𝑉 (0). (4.13)
0
donde se puede nuevamente ver que esto es equivalente a una ecuación de equilibrio del
segmento [−𝑥, 0]. Usando que 𝑞 es par se tiene
−𝑥
𝑉 (𝑥) = − ∫ 𝑞(𝑠)d𝑠. (4.16)
0
por lo que queda demostrado que 𝑉 es una función impar y el cortante es nulo en el eje de
simetría.
Desarrollo de momento El momento interno en una sección cualquiera 𝑥 puede ser ob-
tenido a partir de la segunda ecuación puntual de equilibrio:
𝑥
𝑀(𝑥) = ∫ 𝑉 (𝑢)d𝑢 + 𝑀(0). (4.18)
0
Desarrollo de giro El giro puede ser calculado integrando la ecuación de momento cur-
vatura: 𝑥
𝑀(𝑢)
𝜃(𝑥) = ∫ d𝑢 + 𝜃(0). (4.21)
0 𝐸𝐼
Haciendo cambio de variable 𝑠 = −𝑢 y usando que el momento es par se obtiene:
−𝑥
𝑀(𝑠)
𝜃(𝑥) = − ∫ d𝑠 + 𝜃(0). (4.22)
0 𝐸𝐼
Si se considera cada mitad de la viga por separado se puede ver que el giro del punto en
el eje de simetría debería tener valores opuestos debido a la simetría de las cargas. Por otra
parte por continuidad de la viga en dicho punto se cumple que los giros de ambos extremos
deben ser iguales por lo que se concluye que dicho giro debe de ser nulo, por lo tanto 𝜃(0) = 0.
Se tiene por lo tanto que:
−𝑥
𝑀(𝑠)
𝜃(𝑥) = − ∫ d𝑠 = −𝜃(−𝑥) ⇒ 𝜃(𝑥) = −𝜃(−𝑥) . (4.23)
0 𝐸𝐼
por lo que se demostró que el giro es impar y tiene valor nulo en el eje de simetría.
𝐸𝐴 ⇒
𝐸𝐴/2
𝑞 𝑞 𝑞
𝐸𝐼 ⇒
𝐸𝐼 /2
𝑞 𝑞 𝑞
2 1 4
3m
2m 2m
La flecha con un círculo representa la carga aplicada considerada móvil, en todo el tramo
horizontal de la estructura y de módulo unitario.
I 1 II
Δ𝑦,𝑃
Δ𝑦,3 Δ𝑦,3
La condición de apoyo fijo en 3 es Δ𝑦,3 = 0, por lo que imponiendo esto, se puede obtener
el valor de la reacción deseada:
Δ𝐼𝑦,3
𝑅𝑦,3 = − 𝐼 𝐼 . (4.28)
Δ𝑦,3
Para poder obtener un método simple de cálculo, se debe poder calcular el desplazamiento
del estado I de otra forma. Para esto, se utiliza el Teorema de Betti, el cual permite establecer
una relación entre el trabajo de las cargas del estado I con los desplazamientos del estado II
y las cargas del estado II y los desplazamientos del estado I.
La relación establecida por el Teorema de Betti es:
1Δ𝐼𝑦,𝑃
𝐼
= 1Δ𝐼𝑦,3 , (4.29)
expresión que permite calcular la línea de influencia realizando el análisis de un único EBC
con un único estado de cargas estáticas. En los materiales complementarios audiovisuales se
mostrarán los diagramas utilizando herramientas computacionales y los desarrollos analíti-
cos correspondientes.
I Δ𝑦,𝑀𝑖 II
1 1
1
Δ𝑦,𝑃
Δ𝑦,𝑀𝑑
Por otra parte, utilizando superposición, los desplazamientos en la estructura real (con el
valor de cortante 𝑉𝑀 a determinar), pueden ser escritos como:
Δ𝐼𝑦,𝑀𝑖 + 𝑉𝑀 Δ𝐼𝑦,𝑀𝑖
𝐼
= − (Δ𝐼𝑦,𝑀𝑑 + 𝑉𝑀 Δ𝐼𝑦,𝑀𝑑
𝐼
) (4.34)
Δ𝐼𝑦,𝑀𝑖 + Δ𝐼𝑦,𝑀𝑑
𝑉𝑀 = − (4.35)
Δ𝐼𝑦,𝑀𝑖
𝐼 + Δ𝐼𝑦,𝑀𝑑
𝐼
Δ𝐼𝑦,𝑃
𝐼
𝑉𝑀 = − (4.37)
Δ𝐼𝑦,𝑀𝑖
𝐼 + Δ𝐼𝑦,𝑀𝑑
𝐼
4.3 Ejercicios
En todos los ejercicios presentados a continuación se deberá trabajar con estructuras
simplificadas a través del uso de la simetría o antisimetría cuando corresponda, así como
también, resortes equivalentes para algunas componentes estructurales.
Ejercicio 4.1
En la estructura de acero (𝐸 = 210 GPa) de la figura, la viga continua AE está conformada
por 2 PNC 18 soldados ([]) y cada uno de los tensores BF y DG tiene una sección transversal
constituida por 2 barras 𝜙16.
Ejercicio 4.2
Para el marco (formado por cuatro nudos con uniones
rígidas) de hormigón (𝐸 = 30 GPa) mostrado en la
figura, conformado por secciones de dimensiones 15 ×
30 cm (inercia 𝐼0 ) y 30 × 30 cm (inercia 2𝐼0 ):
a) Calcular las reacciones y trazar diagramas de soli-
citaciones.
b) Bosquejar la deformada de la estructura.
Ejercicio 4.3
Sea la estructura de la figura, conformada por barras de acero (E=210 GPa) de sección
tubular de diámetro exterior 𝜙 = 125 mm y espesor 𝑡 = 5 mm. Considerando 𝑃 = 10 kN, se
pide:
Ejercicio 4.4
El pórtico de acero (𝐸 = 210 GPa) mostrado en la figura está construido con perfiles con
sección dada por un perfil PNI 24.
Ejercicio 4.5
La estructura de acero (𝐸 = 210 GPa) mostrada en la figura recibe la descarga distri-
buida de un entrepiso sobre el travesaño BC. El marco ABCD consiste en 2 perfiles PNC 20
apareados (][), mientras que el tensor AD se compone de 2 𝜙16.
Se pide:
a) Hallar los valores de 𝛼 para los cuales el máximo momento que tracciona las fibras supe-
riores y el máximo momento que tracciona fibras inferiores en el travesaño BC, son iguales
en módulo.
b) Para el mayor valor de 𝛼, trazar los diagramas de solicitaciones.
Ejercicio 4.6
Considere la estructura de la Figura 4.7, con 𝐸𝐼 uniforme, sometida a las cargas que se
indican. Se sabe que 𝐿 = 2 m y 𝐻 = 100 kN.
Los pilares (barras BE y CF) están formados por tubulares de acero de diámetro exterior
𝜙𝑒𝑥𝑡 = 45 cm. La sección transversal del piso superior (barras AB, BC y CD) es una sección
compuesta mostrada en la Figura 4.8. Se considera 𝐸𝐴𝐶𝐸𝑅𝑂 = 210 GPa y 𝐸𝑀𝐴𝐷𝐸𝑅𝐴 = 10.5 GPa.
Se pide:
a) Determinar el espesor 𝑡 de los tubulares de forma tal que se cumpla la hipótesis de EI=cte.
b) Determinar el valor de la constante k de los resortes de forma que el desplazamiento
horizontal del piso superior sea 𝛿 = 2.6 mm hacia la derecha.
c) Trazar diagramas de solicitaciones y bosquejar la deformada de la estructura.
Ejercicio 4.7
La figura muestra una viga atensorada. La viga DAEBF consiste en una escuadría de ma-
dera de sección 20 × 60 cm, y el tensor ACB en una barra redonda de acero (𝐸 = 210 GPa) de
𝜙16. El puntal CE se considera infinitamente rígido. Se considera que 𝐸𝑀𝐴𝐷𝐸𝑅𝐴 = 𝐸𝐴𝐶𝐸𝑅𝑂 301 .
acero (𝐸 = 210 GPa) que vincula la barra BC con AD en sus puntos medios, compuesta por
1𝜙 20.
Considerando Si 𝐻 = 8 kN y 𝐿 = 2 m, se pide:
a) Calcular el desplazamiento del punto C.
b) Obtener los diagramas de solicitaciones de la estructura.
c) Bosquejar la deformada de la estructura.
Ejercicio 4.10
Sea la viga continua de la figura con 𝐸𝐼 =cte y con vanos de luz 𝑙.
Unidad temática 5
Estructuras tridimensionales de barras
5.1 Emparrillados
En esta sección se presentan de forma esquemática los conceptos más importantes para la
determinación de desplazamientos y solicitaciones en estructuras planas de barras, sometidas
a fuerzas perpendiculares al plano de la estructura. Estas estructuras, tal como fue visto en
la Sección 1.1.1, son llamadas estructuras planoespaciales o emparrillados.
5.1.1. Torsión
En cursos anteriores se abordó la teoría de torsión de Saint-Venant para vigas de sección
transversal con doble eje de simetría. A continuación se enumeran algunas de las magnitu-
des y conceptos más importantes, aplicando el enfoque utilizado en (Wunderlich and Pilkey,
2002).
Sea una barra con un extremo empotrado sometida a un momento torsor 𝑀𝑡 se verifica:
𝑀𝑡 = 𝐺𝐽 Θ (5.1)
𝐽 = ∫ 𝜓𝑃 d𝐴. (5.2)
𝐴
u obtenida de fórmulas para algunos casos particulares (Young and Budynas, 2001).
Si se considera un elemento con giros nodales 𝜃1 y 𝜃2 se tiene:
𝜃2 − 𝜃1
𝑀𝑡 = 𝐺𝐽 (5.4)
𝓁
1 𝐺𝐽
Π𝑖𝑛𝑡,torsión (𝜃1 , 𝜃2 ) = (𝜃2 − 𝜃1 )2 . (5.5)
2 𝓁
A partir de la relación entre 𝜓𝑃 y las tensiones rasantes se puede calcular también el valor
de la tensión rasante máxima debida a torsión, dado por:
𝑀𝑡
𝜏𝑚𝑎𝑥 = (5.6)
𝑊𝑡
donde 𝑊𝑡 es el módulo resistente torsional.
𝑣𝑦,1 , 𝐹𝑦,1
𝜃𝑥,1 , 𝑀𝑥,1
𝜃𝑥,1,𝐿 , 𝑀𝑥,1,𝐿
𝜃𝑧,1 , 𝑀𝑧,1 𝑣𝑦,2 , 𝐹𝑦,2
𝜃𝑧,1,𝐿 , 𝑀𝑧,1,𝐿
𝛼𝑦 𝜃𝑧,2 , 𝑀𝑧,2 𝜃𝑥,2 , 𝑀𝑥,2
𝜃𝑧,2,𝐿 , 𝑀𝑧,2,𝐿
𝜃𝑥,2,𝐿 , 𝑀𝑥,2,𝐿
Figura 5.1: Esquema de elemento de barra a flexo-torsión. Sistemas de coordenadas locales (L) y coordenadas globales
representados.
como:
Π𝑖𝑛𝑡 (𝐮𝐿 ) = Π𝑖𝑛𝑡,torsión (𝜃𝑥,1,𝐿 , 𝜃𝑥,2,𝐿 ) + Π𝑖𝑛𝑡,flexión (𝑣𝑦,1,𝐿 , 𝜃𝑧,1,𝐿 , 𝑣𝑦,2,𝐿 , 𝜃𝑧,2,𝐿 ) (5.19)
donde Π𝑖𝑛𝑡,flexión está dada por la Ecuación (3.34). La energía potencial de las fuerzas externas
está dada por Π𝑒𝑥𝑡 (𝐮𝐿 ) = −𝐮𝑇𝐿 𝐟𝐿 donde 𝐟𝐿 está dado por:
𝑇
𝐟𝐿 = [𝑀𝑥,1,𝐿 𝐹𝑦,1,𝐿 𝑀𝑧,1,𝐿 𝑀𝑥,2,𝐿 𝐹𝑦,2,𝐿 𝑀𝑧,2,𝐿 ] . (5.20)
Utilizando el principio de mínima energía potencial total y aplicando las condiciones del
primer teorema de Castigliano se tienen las siguientes cuatro condiciones para flexión:
𝜕Π𝑖𝑛𝑡,flexión 𝜕Π𝑖𝑛𝑡,flexión
= 𝐹𝑦,1,𝐿 = 𝑀𝑧,1,𝐿 (5.21)
𝜕𝑣𝑦,1,𝐿 𝜕𝜃𝑧,1,𝐿
𝜕Π𝑖𝑛𝑡,flexión 𝜕Π𝑖𝑛𝑡,flexión
= 𝐹𝑦,2,𝐿 = 𝑀𝑧,2,𝐿 (5.22)
𝜕𝑣𝑦,2,𝐿 𝜕𝜃𝑧,2,𝐿
𝜕Π𝑖𝑛𝑡,torsión 𝜕Π𝑖𝑛𝑡,torsión
= 𝑀𝑥,1,𝐿 = 𝑀𝑥,2,𝐿 (5.23)
𝜕𝜃𝑥,1,𝐿 𝜕𝜃𝑥,2,𝐿
Sustituyendo las expresiones de las energías de deformación se obtiene que las condicio-
nes de Castigliano son equivalentes al sistema de ecuaciones:
𝐊𝐿 𝐮𝐿 = 𝐟𝐿 (5.24)
⎡ 𝐺𝐽 0 0 −
𝐺𝐽
0 0 ⎤
⎢ 𝓁 𝓁 ⎥
⎢ 𝐸𝐼𝑧 𝐸𝐼𝑧 𝐸𝐼𝑧 𝐸𝐼𝑧 ⎥
⎢ 0 12 6 0 −12 6 2 ⎥
⎢ 𝓁3 𝓁2 𝓁3 𝓁 ⎥
⎢ 𝐸𝐼𝑧 𝐸𝐼𝑧 𝐸𝐼𝑧 𝐸𝐼𝑧 ⎥
⎢ 0 6 2 4 0 −6 2 2 ⎥
⎢ 𝓁 𝓁 𝓁 𝓁 ⎥
𝐊𝐿 = ⎢ 𝐺𝐽 𝐺𝐽 ⎥. (5.25)
⎢− 0 0 0 0 ⎥
⎢ 𝓁 𝓁 ⎥
⎢ 𝐸𝐼𝑧 𝐸𝐼𝑧 𝐸𝐼𝑧 𝐸𝐼𝑧 ⎥
⎢ 0 −12 −6 0 12 −6 2 ⎥
⎢ 𝓁3 𝓁2 𝓁3 𝓁 ⎥
⎢ 𝐸𝐼𝑧 𝐸𝐼𝑧 𝐸𝐼𝑧 𝐸𝐼𝑧 ⎥
⎢ 0 6 2 2 0 −6 2 4 ⎥
⎣ 𝓁 𝓁 𝓁 𝓁 ⎦
𝐊𝐮 = 𝐟 (5.26)
donde la matriz de rigidez 𝐊 es la matriz de rigidez en coordenadas globales dada por
𝐊 = 𝐑𝐊𝐿 𝐑𝑇 . (5.27)
Para un emparrillado formado por 𝑛 elementos de barra se debe calcular cada matriz en
coordenadas globales y realizar el ensamblado de forma similar a como es realizado para
elementos de reticulados o pórticos.
5.1.3. Ejemplo
Sea la estructura mostrada en la Figura 5.2, donde todas las barras tienen inercia flexional
𝐼 e inercia torsional 𝐽 a determinar y están formadas por un material de módulo de Young 𝐸.
Se considera que 𝐸 = 210 GPa, 𝜈 = 0.3, 𝓁 = 2 m, la sección transversal es cuadrada de lado 5
cm. Se desea calcular los desplazamientos nodales para una carga puntual 𝑃 = 4 kN.
𝑧 𝑥
𝓁 𝓁/2
𝓁/2
La barra con la carga puntual aplicada tiene una geometría y condiciones de vínculo
simétricos por lo que el giro 𝜃𝑥 debe ser nulo. Considerando está condición de simetría y
dividiendo la estructura en dos mitades, se obtiene el esquema básico simplificado que se
muestra en la Figura 5.3.
1 𝑧 𝑥
𝑦𝐿 𝑃/2
𝑧𝐿 𝑥𝐿
𝑦𝐿
𝑥𝐿 3
𝓁 2
𝑧𝐿
𝓁/2
⎡0 0 1 0 0 0⎤
⎢ ⎥
⎢0 1 0 0 0 0⎥
⎢−1 0 0 0 0 0⎥
𝐑=⎢ . (5.28)
⎢0 0 0 0 0 1⎥⎥
⎢0 0 0 0 1 0⎥
⎢0 0 0 −1 0 0⎥
⎣ ⎦
El módulo de corte es 𝐺 = 2(1+𝜈)
𝐸
= 80.7 GPa. La inercia flexional es 𝐼 = 5.208 × 10−7 m4 y
la inercia torsional, para una sección cuadrada, está dada por 𝐽 = 0.141(0.05)4 = 8.81 × 10−7
m4 .
Los desplazamientos generalizados para el nodo 2 son:
Para validar los resultados obtenidos con el código se realiza otro modelo usando la he-
rramienta ONSAS. Se realiza un primer modelo con análisis lineal y la misma geometría y el
resultado es verificado.
Finalmente para comprender la deformación de la estructura se realizada un análisis no
lineal para una carga máxima 20 veces superior a 𝑃. La deformada para el estado de cargas
máximo con el análisis no lineal es mostrado en la Figura 5.4, donde el factor de escala es 1.
Figura 5.4: Deformada emparrillado, con análisis no lineal y estado de carga 20 𝑃 (factor de esacala: 1).
5.3 Ejercicios
Ejercicio 5.1
Sea la estructura mostrada en la figura, construida con caños de acero (𝐸 = 210 GPa y
𝜈 = 0.3) de 𝜙𝑒𝑥𝑡 = 100 mm y 𝑡 = 8 mm.
𝑃
𝐿
𝐿
2. Análisis matricial.
Ejercicio 5.2
Se construye un emparrillado con caños de acero (E=210 GPa, 𝜈 = 0.3) de 𝜙𝑒𝑥𝑡 = 100 mm
y 𝑡 = 8 mm, sometido a dos estados de carga mostrados en las figuras.
Considerando que todas las barras están compuestas por el mismo tubular que el empa-
rrillado, se pide:
c) Modelar la estructura con un programa computacional y comparar los resultados con los
obtenidos en las partes anteriores en la zona del emparrillado.
Ejercicio 5.3
Sea un emparrillado modificado con respecto al del Ejercicio 2, en el cual se agrega una
barra adicional, que está sometida a una carga 𝑞 = 2 kN/m, como se muestra en la figura.
Se pide:
a) Resolver mediante análisis matricial. Trazar diagramas de solicitaciones.
b) Comparar resultados obtenidos en a) con los obtenidos utilizando alguna herramienta
computacional.
Ejercicio 5.4
Sea la estructura mostrada en la figura, construida con barras de acero (E=210 GPa, 𝜈 =
0.3) de sección circular maciza de diámetro 𝜙𝑒𝑥𝑡 = 105 mm.
Ejercicio 5.5
La estructura tridimensional de hormigón (𝐸 = 30 GPa y 𝜈 = 0.2) mostrada en las está
conformada por pilares de sección cuadrada 30 × 30 cm y vigas de 𝑏 = 20 cm y ℎ = 50 cm.
Ejercicio 5.6
Se desea estudiar un aerogenerador cuyas dimensiones se muestran en la figura a la iz-
quierda. Todas las palas (o aspas) tienen el mismo largo y tienen mismo ángulo de apertura
120◦ . En la figura a la derecha se muestra un estado de cargas estáticas idealizado, de interés
para el diseño. Las fuerzas actuantes son:
Una fuerza del viento sobre la torre que se ha simplificado como una fuerza en B hori-
zontal de 20 kN.
El peso de las aspas de 30 kN que se aplica como una carga puntual en el punto E.
80 kN F
q=6.5 kN/m
15 m
C
D
2m
7m
E
12.5 m
30 kN 80 kN
H G
80 kN
20 kN
B
12.5 m
Ejercicio 5.7
Se desea analizar la estructura de un pasamanos considerando un correspondiente Es-
quema Básico de Cálculo como se muestra en las figuras.
La estructura está compuesta por tubulares de acero (𝐸 = 210 GPa y 𝜈 = 0.3) de 𝜙𝑒𝑥𝑡 = 80
mm y 𝑡 = 6 mm. Sobre las barras transversales actúa una carga q=3 kN/m y L=1.0 m. Se pide:
a) Considerar que las barras no presentan rigidez a torsión y calcular analíticamente las
solicitaciones en cada una de las barras.
b) Considerar la rigidez a torsión en las barras y obtener resultados de modelar computacio-
nalmente la estructura.
Ejercicio 5.8
Sea una estructura formada por barras de sección cuadrada de 40 cm de lado, compuestas
por un material elástico lineal. La estructura tiene aplicadas dos cargas puntuales: en su nodo
A, una carga 𝑃𝐴 = 500 kN y en su nodo B, otra carga 𝑃𝐵 = 250 kN. Por otra parte, la estructura
está sometida a una carga distribuida uniforme 𝑞 = 250 N/m con sentido según 𝑦𝐺 . No se
considera el peso propio.
A la izquierda muestra el esquema de la estructura en el plano 𝑥 − 𝑧 y a la derecha el
esquema básico estructural con ejes locales propuestos para cada barra.
10 m 20 m
𝑧𝐿
𝐴 𝑦𝐿 𝑥𝐿
𝐵 𝑧𝐿
𝑦𝐿 𝑥𝐿 𝑃𝐵
𝑃𝐴
𝑞
40 m 𝑃𝐵 𝑃𝐴
𝑦𝐿 𝑥𝐿
𝑧𝐿
𝑥𝐺 𝑦𝐺
𝐶 𝑥𝐺
𝑦𝐺
𝑧𝐺
𝑧𝐺
Trazar los diagramas de todas las solicitaciones de la estructura, de acuerdo con los ejes
locales indicados para cada barra. Verificar que el momento flector 𝑀𝑦 es nulo en la sección
C.
Considerar que todas las vigas son de hormigón (E=30 GPa, 𝜈 = 0.2) y su sección es 30x80
cm. Las mismas están cargadas con una carga q=15 kN/m y L=5.0 m.
a) Modelar la estructura computacionalmente. Analizar el comportamiento de las solicita-
ciones y la deformada.
b) Si se desprecia la rigidez a torsión de las vigas, ¿cómo se ven afectadas las solicitaciones
y deformaciones de las vigas? Analizar.
Unidad temática 6
Análisis Seccional
La definición en términos de la función tensión puede ser enunciada como, que dada una
sección transversal ubicada en un punto de coordenada 𝑥, la Línea Neutra es el conjunto de
puntos (𝑦, 𝑧) que verifican: 𝜎(𝑥, 𝑦, 𝑧) = 0.
Existe una correspondencia biunívoca entre cada línea neutra y un punto de aplicación
de la carga. A continuación se muestran los desarrollos para obtener la LN a partir del punto
de aplicación y luego el desarrollo inverso.
𝑦 torsor equiv. 𝑦
𝐴
𝑁
𝑀𝑦
𝑧 𝑧 𝑀𝑧
𝑁
𝑥 𝑥
Figura 6.1: Esquema de sección, sistema de coordenadas considerado y equivalencia de torsores de carga y solicita-
ciones internas.
El punto A tiene coordenadas (𝑦𝐴 , 𝑧𝐴 ) las cuales pueden ser calculadas como
𝑀𝑧 (𝑥) 𝑀𝑦 (𝑥)
𝑦𝐴 = − y 𝑧𝐴 = (6.1)
𝑁 (𝑥) 𝑁 (𝑥)
𝑦𝐴 𝑧𝐴
1+ 2
𝑦 + 2𝑧 = 0 (6.4)
𝜌𝑧 𝜌𝑦
Se observa que la recta está contenida en el 𝑝𝑙𝑎𝑛𝑜 𝑦𝑧 y no pasa por el origen, esto se debe
a la existencia de una fuerza de directa aplicada.
𝑎 𝑏
𝑦 + 𝑧 + 1 = 0. (6.5)
𝑐 𝑐
Igualando los coeficientes de los polinomios de las expresión de las Ecuaciones (6.4) y
(6.5) se tiene:
𝑎 𝑏
𝑦𝐴 = 𝜌𝑧2 𝑧𝐴 = 𝜌𝑦2 (6.6)
𝑐 𝑐
Dado que hay directa aplicada la recta que define la LN no pasa por el origen, por lo tanto
se cumple que 𝑐 ≠ 0.
Esto es equivalente a decir que, para cualquier punto de aplicación perteneciente al NC,
la línea neutra correspondiente no corta a la sección.
𝑁 𝑦𝐴 (𝑡) 𝑧𝐴 (𝑡)
𝜎(𝑥, 𝑦(𝑡̂), 𝑧(𝑡̂)) = 1 + 2 𝑦(𝑡̂) + 2 𝑧(𝑡̂) = 0 (6.9)
𝐴( 𝜌𝑧 𝜌𝑦 )
𝜕𝜎 𝑦𝐴 (𝑡) 𝑧𝐴 (𝑡)
(𝑥, 𝑦(𝑡̂ = 𝑡), 𝑧(𝑡̂ = 𝑡)) = 0 = 𝑦(𝑡)
̇ + 2 𝑧(𝑡) ̇ (6.10)
𝜕𝑡̂ 𝜌𝑧2 𝜌𝑦
Combinando la Ecuación (6.9) con la Ecuación (6.10) y usando que es regular y no pasa
por el origen, se obtiene la curva paramétrica del borde del núcleo central:
𝑧(𝑡)
̇
𝑦𝐴 (𝑡) = − 𝜌𝑧2 (6.11)
𝑧(𝑡)𝑦(𝑡)
̇ − 𝑦(𝑡)𝑧(𝑡)
̇
𝑦(𝑡)
̇
𝑧𝐴 (𝑡) = 𝜌𝑦2 (6.12)
𝑧(𝑡)𝑦(𝑡)
̇ − 𝑦(𝑡)𝑧(𝑡)
̇
Se conoce que los puntos de aplicación (𝑦1 , 𝑧1 ) y (𝑦2 , 𝑧2 ) pertenecen al NC y tienen aso-
ciados las líneas neutras 𝐿𝑁1 y 𝐿𝑁2 respectivamente. La unión de los puntos de aplicación de
las fuerzas de directa 𝑁1 y 𝑁2 forman una recta 𝑟, que se puede expresar como:
(𝑧1 − 𝑧2 )
𝑟 ∶ 𝑧 − 𝑧1 = (𝑦 − 𝑦1 ).
(𝑦1 − 𝑦2 )
𝑦(𝑧1 − 𝑧2 ) − 𝑧(𝑦1 − 𝑦2 )
=1 (6.13)
𝑧1 𝑦2 − 𝑦1 𝑧2
𝑦𝐴 (𝑧1 − 𝑧2 ) − 𝑧𝐴 (𝑦1 − 𝑦2 )
= 1, (6.14)
𝑧 1 𝑦2 − 𝑦1 𝑧 2
o también
𝑧𝐴 𝑦2 − 𝑦𝐴 𝑧2 𝑦𝐴 𝑧1 − 𝑧𝐴 𝑦1
+ = 1, (6.15)
𝑧 1 𝑦2 − 𝑦1 𝑧 2 𝑧 1 𝑦2 − 𝑦1 𝑧 2
Por otra parte, la fuerza de directa 𝑁 actuando en este punto puede ser descompuesta
en dos fuerzas de directa aplicadas simultáneamente en (𝑦1 , 𝑧1 ) y (𝑦2 , 𝑧2 ) respectivamente, se
proponen las siguientes expresiones para dichas fuerzas,
𝑧 𝐴 𝑦2 − 𝑦𝐴 𝑧 2 𝑦𝐴 𝑧 1 − 𝑧 𝐴 𝑦1
𝑁1 = 𝑁 𝑁2 = 𝑁 (6.16)
𝑧 1 𝑦2 − 𝑦1 𝑧 2 𝑧1 𝑦2 − 𝑦1 𝑧2
Se demuestra que estas dos fuerzas de directa actuando simultáneamente tienen un torsor
equivalente al de la fuerza 𝑁 aplicada en (𝑦𝐴 , 𝑧𝐴 ).
La equivalencia en el momento 𝑀𝑦 es
𝑀𝑦 = 𝑁1 𝑧1 + 𝑁2 𝑧2 = 𝑁 𝑧𝐴 , (6.17)
𝑀𝑧 = −(𝑁1 𝑦1 + 𝑁2 𝑦2 ) = −𝑁 𝑦𝐴 , (6.18)
𝑦𝐴 (𝑧1 − 𝑧2 ) − 𝑧𝐴 (𝑦1 − 𝑦2 )
𝑁1 + 𝑁2 = 𝑁 = 𝑁, (6.19)
𝑧 1 𝑦2 − 𝑦1 𝑧 2
𝑁 = 𝑁1 + 𝑁2 = 𝑁 (𝜆1 + 𝜆2 )
6.1.3. Ejemplos
Sección elíptica
En la Figura 6.5a se tiene una sección regular de forma elíptica cuyo contorno responde
a la siguiente curva paramétrica:
1 https://en.wikipedia.org/wiki/Convex_cone
{
𝑦(𝑡) = 𝐴𝑐𝑜𝑠(𝑡) 𝑦(𝑡)
̇ = −𝐴𝑠𝑒𝑛(𝑡)
(6.23)
𝑧(𝑡) = 𝐵𝑠𝑒𝑛(𝑡) 𝑧(𝑡)
̇ = 𝐵𝑐𝑜𝑠(𝑡)
𝐴
𝑦𝐴 (𝑡) = − 𝑐𝑜𝑠(𝑡) (6.24)
4
𝐵
𝑧𝐴 (𝑡) = − 𝑠𝑒𝑛(𝑡) (6.25)
4
Para el caso particular de una sección circular de radio 𝑅 se tiene que 𝐴 = 𝐵 = 𝑅, por lo tanto
la curva paramétrica que describe al NC de una sección circular es:
𝑅
𝑦𝐴 (𝑡) = − 𝑐𝑜𝑠(𝑡) (6.26)
4
𝑅
𝑧𝐴 (𝑡) = − 𝑠𝑒𝑛(𝑡) (6.27)
4
Sección rectangular
En la Figura 6.6a se tiene una sección de forma rectangular de ancho 𝑎 y alto 𝑏. Se plantea
la curva paramétrica correspondiente al lado superior de la sección como:
{
𝑦(𝑡) = 𝑏/2 𝑦(𝑡)
̇ =0
(6.28)
𝑧(𝑡) = 𝑡 𝑧(𝑡)
̇ =1
Utilizando la Ecuación (6.11) y la Ecuación (6.12) se obtiene la curva paramétrica que
relevancia.
𝑁 = ∫ 𝜎(𝑧)d𝐴, (6.31)
𝐴
donde 𝐴 representa la región de la sección que está comprimida, es decir donde 𝜎 es no nula.
En la Figura 6.7a se observa que la distribución de tensiones es lineal y que 𝜎 se puede
escribir como,
𝜎𝑚𝑎𝑥
𝜎(𝑧) = 𝑘𝑧 = − 𝑧 (6.32)
𝑐
Sustituyendo en la Ecuación (6.31) se obtiene,
𝜎𝑚𝑎𝑥
𝑁 = ∫ 𝜎(𝑧)d𝐴 = − 𝑧d𝐴 (6.33)
𝐴 𝑐 ∫𝐴
Por otra parte se plantea el equilibrio de momentos respecto del origen de coordenadas
definido,
𝜎𝑚𝑎𝑥
𝑁 (𝑐 − 𝑑) = ∫ 𝜎(𝑧)𝑧d𝐴 = − 𝑧 2 d𝐴 (6.34)
𝐴 𝑐 ∫𝐴
Se define la inercia de primer orden (o momento estático) referida a la línea neutra, que
será función de la posición de la línea neutra 𝑐, como:
6.2.2. Ejemplo
Sección rectangular
𝑐
𝑏𝑐3
𝐼𝐿𝑁 (𝑐) = ∫ 𝑧 2 d𝐴 = 𝑏 ∫ 𝑧 2 d𝑧 = (6.40)
𝐴 0 3
𝑐
𝑏𝑐 2
𝜇𝐿𝑁 (𝑐) = ∫ 𝑧d𝐴 = 𝑏 ∫ 𝑧d𝑧 = (6.41)
𝐴 0 2
Se tiene entonces que la relación entre el momento de segundo orden y el momento de
primer orden es,
𝐼𝐿𝑁 (𝑐) 2
= 𝑐 (6.42)
𝜇𝐿𝑁 (𝑐) 3
Utilizando la Ecuación (6.39) se tiene que,
2
𝑐−𝑑 = 𝑐 (6.43)
3
Obteniendo que,
𝑐 = 3𝑑 (6.44)
Por otra parte utilizando la Ecuación (6.37) se tiene que,
𝜎𝑚𝑎𝑥 𝑏𝑐 2
𝑁 =− (6.45)
𝑐 2
Obteniendo que la tensión máxima de compresión es,
2𝑁
𝜎𝑚𝑎𝑥 = − (6.46)
𝑏𝑐
6.3 Ejercicios
Ejercicio 6.1
La pieza de la figura soporta las cargas 𝑃 = 30 kN indicadas. Obtener el diagrama de
tensiones normales en la sección 𝑥 − 𝑥.
Ejercicio 6.2
Ejercicio 6.3
Una barra recta, formada por un perfil normalizado L 100𝑥100𝑥10, de 3𝑚 de longitud,
está empotrada en uno de sus extremos y sometida a su peso propio (𝑔 = 10𝑚/𝑠 2 ), como se
muestra en la figura.
Ejercicio 6.4
Para la generación de columnas de iluminación se suelen emplear elementos prefabrica-
dos como los de la Figura 6.9a. Se considera una de esas columnas de sección tubular de radio
externo 𝑅 y radio interno 𝑟 = 𝛼𝑅, como indica la Figura 6.9b.
a) Determinar 𝛼 para que el perímetro interior sea el contorno del núcleo central de la sec-
ción.
b) Si se aplica una directa de compresión 𝑃 que varía de ubicación entre A y B, y el material es
tal que 𝜎𝑎𝑑𝑚,𝑡𝑟𝑎𝑐 = 9𝑀𝑃𝑎 y 𝜎𝑎𝑑𝑚,𝑐𝑜𝑚𝑝 = 40𝑀𝑃𝑎, calcular 𝑃𝑎𝑑𝑚 para 𝑅 = 50𝑐𝑚 y el 𝛼 determinado.
Ejercicio 6.5
En la Figura 6.10a se presenta una viga pretensada de largo 𝐿 = 25𝑚. El esquema básico
de cálculo de la viga se presenta en la Figura 6.10b y su sección se muestra en la Figura 6.11.
Pretensar un elemento estructural implica introducirle esfuerzos previamente a su puesta
en funcionamiento con el propósito de contrarrestar aquellos que serán ocasionados por la
Figura 6.9
aplicación de las cargas que actuarán cuando ella entre en servicio. El pretensado de esta viga
se puede representar mediante una fuerza concentrada 𝐹 de compresión en el eje vertical, a
una cierta distancia 𝑒 del borde inferior.
Figura 6.10
Ejercicio 6.6
Ejercicio 6.7
Sobre la sección de la figura se aplican las cargas 𝑃 y 𝑄 de igual valor y signo. El punto de
aplicación de P es fijo y coincide con el baricentro 𝐺, mientras que el de 𝑄 se desplaza sobre
el perímetro de la sección.
Ejercicio 6.8
Ejercicio 6.9
La zapata de la figura trasmite al terreno la descarga vertical P cuyo punto de aplicación
es el indicado en la figura.
Calcular 𝑃𝑎𝑑𝑚 si 𝜎𝑡𝑒𝑟𝑟𝑒𝑛𝑜 = 0, 6𝑀𝑃𝑎.
Ejercicio 6.12
Figura 6.12
Figura 6.13
Ejercicio 6.13
En algunos casos de estructuras de transmisión eléctrica es usual utilizar secciones circu-
lares de hormigón armado para transmitir las cargas de la torre a la fundación. Sea la sección
circular de diámetro Φ=3.2 m que se muestra en la siguiente figura.
𝑀 Φ
𝐺
𝑁
Se sabe que la tensión normal máxima es 𝜎=0.171 MPa y que la longitud comprimida
medida desde un extremo de la sección es 𝑐=1.836 m. Asumiendo que el hormigón es un
material que no soporta tracciones, que se tiene un diagrama de tensiones de forma triangular
se pide:
a) Obtener la pareja de solicitaciones (𝑁 , 𝑀) que debe actuar en el baricentro de la sección
para satisfacer las condiciones dadas.
Unidad temática 7
Estabilidad Estructural
En esta unidad temática se consideran nuevas referencias, dado que se debe introducir un
concepto fundamentalmente diferente al resto de las unidades temáticas. La hipótesis funda-
mental modificada: es que el equilibrio deja de realizarse en la configuración de referencia
para realizarlo en una configuración deformada (bajo hipótesis de pequeños giros), permi-
tiendo modelar el fenómeno de inestabilidad estructural. Las principales referencias utilizadas
son (Yoo and Lee, 2011; Bazzano and Pérez Zerpa, 2017).
una tipología estructural de puente que se estaba popularizando en ese momento. Los puentes
con sección cajón de acero de Milford Haven de 1970 en Gales, Koblenz de 1971 en Alemania
y West Gate de 1970 en Australia, son algunos de los casos en dicha serie de colapsos. En
todos ellos, la inestabilidad de las placas de acero, que conforman los cajones, fue un factor
determinante en los colapsos. Consultar el artículo de Firth (2010)1 para conocer la historia
de dicha tipología, así como la serie de colapsos en la década de los 70.
Estos colapsos fueron el disparador de un gran esfuerzo de investigación en estructuras
con el objetivo de generar reglas de diseño para evitar nuevos accidentes en esta clase de
estructuras. En Inglaterra, dicho impulso resultó en la serie de reglas de diseño conocidas co-
mo IDWR, que luego pasarían a formar parte de BS 5400-3, la antigua norma británica para
diseño de puentes de acero, hoy reemplazada por el Eurocódigo EN 1993. Dicha investigación
avanzó enormemente el entendimiento de cómo se comportan los cajones de acero, en par-
ticular la estabilidad de las chapas rigidizadas que los componen. El lector interesado puede
leer el reporte de la comisión investigadora del colapso de West Gate2 , en el cual se puede
apreciar la magnitud y complejidad de la estructura que colapsó.
Las normas modernas exigen, o permiten, que los problemas comunes de inestabilidad
sean considerados mediante reglas relativamente sencillas. Todos los materiales incluyen re-
glas en sus respectivas normas para la evaluación de la estabilidad estructural, tanto a nivel de
componentes así como a nivel global de la estructura. Sin embargo, el diseño de estructuras
de acero es el que hace mayor hincapié en este tema, dado que en general dichas estructuras
se forman a partir de un gran número de componentes esbeltos.
A modo de ejemplo, se refiere a la sección E3 de la norma AISC 360-16 o sección 6.3.1
de la norma EN 1993-1-1 para ver el tratamiento de resistencia axial de columnas esbeltas de
acero que fallan por inestabilidad. El lector encontrará que las normas hacen referencia a 𝐹𝑒
(AISC 360-16) y 𝑁𝑐𝑟 (EN 1993-1-1) que son esencialmente la carga crítica que Euler descubrió
en 1757.
Las normas técnicas exigen también la evaluación de la estabilidad global de la estructura.
El capítulo C de la norma AISC 360-16 trata sobre este tema, mientras que la norma EN 1992-1-
1 lo hace en el Anexo H. En pocas palabras, se requiere que el profesional evalúe la estabilidad
de la estructura en su conjunto y que determine los efectos que puedan surgir de la esbeltez
global de la estructura. Estas evaluaciones son parte básica del trabajo del profesional que
diseña una estructura.
Finalmente, como motivación se muestran diagramas de deformadas de dos modelos
computacionales con inestabilidad. En la Figura 7.1a se muestra un modelo de una chapa
sometida a cargas de compresión meridional, mientras que en la Figura 7.1b se observa una
viga sometida a cargas verticales mostrando pandeo lateral.
Ambos elementos y cargas estructurales están presenten en puentes del tipo de los cita-
dos en la revisión histórica. Se cuenta entonces con la motivación de entender los modelos
matemáticos que permiten predecir y evitar este tipo de fenómenos y las consecuentes fallas.
1 Disponible en istructe.org/webtest/files/48/488ca532-a956-4929-9f50-aee0d5317afc.pdf, últ. acceso 15/Nov/2020
2 Disponible en parliament.vic.gov.au/papers/govpub/VPARL1971-72No2.pdf, ult. acceso 15/Nov/2020
(a) Pandeo de Chapa a compresión (b) Pandeo lateral de viga por flexión
En esta Unidad Temática se introduce del concepto de Estabilidad dado por la siguiente
definición.
La condición de mínima energía potencial total no solo garantiza el equilibrio de las es-
tructuras, sino que también implica la estabilidad de las mismas.
Caso lineal
Para el modelo estructural lineal (i.e. equilibrio en la configuración de referencia, peque-
ños desplazamientos y material elástico lineal) se vio que la energía potencial total de la
estructura está dada por la expresión general:
1
Π(𝐮) = 𝐮 ⋅ 𝐊𝐮 − 𝐅𝐺 ⋅ 𝐮 (7.2)
2
donde 𝐅𝐺 es el vector de fuerzas generalizadas y 𝐊 la matriz de rigidez. Se puede verificar que
para estructuras con condiciones de apoyo suficientes y que no tengan mecanismos internos
𝐮 ⋅ 𝐊𝐮 > 0 ∀𝐮 ≠ 0 ∈ . (7.3)
Desarrollo lineal
En este desarrollo, se considera como hipótesis que la energía de deformación axial es des-
preciable respecto a la energía de deformación por flexión. Asimismo, también se desprecia
la energía de deformación por corte o distorsión angular.
En la Figura 7.2 se muestra una columna donde cada sección transversal ubicada en la
posición 𝑥 está formada por un material de módulo de Young 𝐸(𝑥) y tiene inercia 𝐼 (𝑥). La
columna está sometida a una carga de compresión 𝑃 y una carga distribuida transversal 𝑞(𝑥).
1 1
Π𝑖𝑛𝑡 = ∫ 𝝈 ∶ 𝜺d𝑉 = ∫ 𝜎𝑥 𝜀𝑥 d𝑉 , (7.4)
2 Ω 2 Ω
donde sustituyendo las Ecuaciones (1.20) y (3.2) y usando que se desprecia energía de defor-
mación axial 𝜀𝐺 , se tiene:
2 2
1 𝓁 𝜕2 𝑣 1 𝓁 𝜕2 𝑣
Π𝑖𝑛𝑡 = 𝐸(𝑥)𝑦 2
d𝐴d𝑥 = 𝐸𝐼 (𝑥) d𝑥. (7.5)
2 ∫0 ∫𝐴(𝑥) ( 𝜕𝑥 2 ) 2 ∫0 ( 𝜕𝑥 2 )
Se debe ahora obtener una expresión de la energía potencial de las fuerzas externas, como:
𝓁
Π𝑒𝑥𝑡 = −𝑃Δ𝑏 − ∫ 𝑞(𝑥)𝑣(𝑥)d𝑥. (7.6)
0
El descenso Δ𝑏 puede ser visto como una función no lineal de la función de flecha 𝑣:
√
2
𝜕𝑣
Δ𝑏 (𝑣) = 𝓁 − ∫ 1− d𝑠.
C𝑑 ( 𝜕𝑠 )
Es posible por lo tanto, establecer y aplicar un orden de aproximación para esta función,
en el caso de aplicar primer orden se puede mostrar que Δ𝑏 = 0, lo que sería coherente con
la teoría lineal.
En el caso de establecer una aproximación de segundo orden se puede aplicar un
desarrollo de Taylor respecto al origen, obteniendo:
2 2
1 𝜕𝑣 1 𝓁 𝜕𝑣
Δ𝑏 ≈ 𝓁 − ∫ 1− d𝑠 = ∫ d𝑥. (7.8)
C𝑑 [ 2 ( 𝜕𝑠 ) ] 2 0 ( 𝜕𝑥 )
Se tiene entonces que la energía potencial de las fuerzas Π𝑒𝑥𝑡 puede ser aproximada por
la siguiente expresión cuadrática en 𝑣:
𝓁 𝓁 2
1 𝜕𝑣
Π𝑒𝑥𝑡 = − 𝑃 ∫ d𝑥 − ∫ 𝑞(𝑥)𝑣(𝑥)d𝑥 (7.9)
2 0 ( 𝜕𝑥 ) 0
2 2
1 𝓁 𝜕2 𝑣 1 𝓁
𝜕𝑣 𝓁
Π(𝑣) = ∫ 𝐸𝐼 (𝑥) 2
d𝑥 − 𝑃 ∫ d𝑥 − ∫ 𝑞(𝑥)𝑣(𝑥)d𝑥, (7.10)
2 0 ( 𝜕𝑥 ) 2 0 ( 𝜕𝑥 ) 0
donde se puede observar que esta función (u operador) es cuadrático en 𝑣 con un nuevo
término respecto al desarrollo de la teórica lineal.
donde 𝜀 es un real pequeño arbitrario. Se puede mostrar que 𝜂 debe cumplir condiciones de
contorno homogéneas.
La expresión anterior representa una pequeña perturbación arbitraria de la deformación
de la columna respecto de la configuración de equilibrio representada en la Figura 7.3.
𝓁 2 2
𝐸𝐼 𝜕2 𝑣 𝜕2 𝜂 𝑃 𝜕𝑣 𝜕𝜂
Π=∫ 2
+𝜀 2 − +𝜀 − 𝑞(𝑥) (𝑣(𝑥) + 𝜀𝜂(𝑥)) d𝑥 (7.12)
0 [ 2 ( 𝜕𝑥 𝜕𝑥 ) 2 ( 𝜕𝑥 𝜕𝑥 ) ]
𝑑Π || 𝓁
𝜕2 𝑣 𝜕2 𝜂 𝜕𝑣 𝜕𝜂
| =∫ 𝐸𝐼 −𝑃 − 𝑞(𝑥)𝜂(𝑥) d𝑥 = 0 ∀𝜂 (7.15)
𝑑𝜀 || 0 [ 𝜕𝑥 2 𝜕𝑥 2 𝜕𝑥 𝜕𝑥 ]
𝜀=0
Para simplificar la Ecuación (7.15) se aplica integración por partes y las condiciones de
contorno homogéneas de 𝜂 𝜂(0) = 𝜂(𝓁) = 0. Se tiene entonces que para el segundo término
de la Ecuación (7.15),
𝓁
𝓁
𝜕𝑣 𝜕𝜂 𝜕𝑣 || 𝓁 2
𝜕𝑣 𝓁 2
𝜕𝑣
∫ d𝑥 = 𝜂| − ∫ 𝜂d𝑥 = − ∫ 𝜂d𝑥 (7.16)
0 𝜕𝑥 𝜕𝑥 𝜕𝑥 || 0 𝜕𝑥 2
0 𝜕𝑥 2
0
Por otra parte para el primer término de la Ecuación (7.15) integrando por parte una vez
se tiene,
𝓁
𝓁
𝜕2 𝑣 𝜕2 𝜂 𝜕2 𝑣 𝜕𝜂 || 𝓁
𝜕 𝜕2 𝑣 𝜕𝜂
∫ 𝐸𝐼 𝜕𝑥 2 𝜕𝑥 2 d𝑥 = 𝐸𝐼 𝜕𝑥 2 𝜕𝑥 || − ∫ 𝜕𝑥 (𝐸𝐼 𝜕𝑥 2 ) 𝜕𝑥 d𝑥
0 | 0
0
La condición anterior debe cumplirse para cualquier función 𝜂 arbitraria, por lo que ne-
𝜕2 𝜕2 𝑣 𝜕2 𝑣
𝐸𝐼 (𝑥) + 𝑃 (𝑥) − 𝑞(𝑥) = 0 Ecuación Euler-Lagrange (7.19)
𝜕𝑥 2 ( 𝜕𝑥 2 ) 𝜕2 𝑥
𝜕2 𝑣 ||
𝐸𝐼 2 | = 0 Condición de contorno (7.20)
𝜕𝑥 ||
𝑥=0
𝜕2 𝑣 ||
𝐸𝐼 2 | = 0 Condición de contorno (7.21)
𝜕𝑥 ||
𝑥=𝓁
𝜕 𝜕2 𝑣 𝜕𝑣
𝑉 (𝑥) = 𝐸𝐼 2 (𝑥) + 𝑃 (𝑥) (7.22)
𝜕𝑥 ( 𝜕𝑥 ) 𝜕𝑥
𝜕2 𝑣
𝑀𝑧 (𝑥) = 𝐸𝐼 (𝑥) (7.23)
𝜕𝑥 2
Donde 𝑣𝑝 (𝑥) será una solución particular que dependerá de la carga distribuida 𝑞(𝑥) y
𝑣ℎ (𝑥) será la solución de la ecuación homogénea, esto es la solución de la ecuación diferencial
cuando 𝑞(𝑥) = 0.
En este paso se establece una hipótesis frecuente en análisis de pórticos 𝐸𝐼 uniforme en
todo el elemento de viga. Además se define un parámetro 𝑘 como:
√
𝑃
𝑘= (7.25)
𝐸𝐼
De esta manera la Ecuación (7.19) se puede reescribir como,
𝜕4 𝑣 2
2𝜕 𝑣
(𝑥) + 𝑘 (𝑥) = 0 (7.26)
𝜕𝑥 4 𝜕2 𝑥
𝜕2 𝑣
Utilizando un cambio de variable 𝑣̂ = 𝜕𝑥 2 y sustityendo en la Ecuación (7.26) se tiene
𝜕2 𝑣̂
(𝑥) + 𝑘 2 𝑣(𝑥)
̂ =0
𝜕𝑥 2
por lo que la solución homogénea sería
𝑣(𝑥)
̂ = 𝐴̂ cos(𝑘𝑥) + 𝐵̂ sin(𝑘𝑥).
donde las constantes multiplicando cos y sin fueron renombradas y las constantes 𝐴, 𝐵, 𝐶 y
𝐷 deben determinarse empleando las condiciones de contorno sobre la solución particular
en cada problema como se verá en la Sección 7.3.
𝑣(0) = 0 ⇒ 𝐴 + 𝐷 = 0, (7.32)
𝑣(𝓁) =
0 ⇒ 𝐴 cos(𝑘𝓁) + 𝐵 sin(𝑘𝓁) + 𝐶𝓁 + 𝐷 = 0, (7.33)
𝜕2 𝑣
𝑀(0) = 𝐸𝐼 2 (0) = 0 ⇒ −𝐴𝑘 2 = 0, (7.34)
𝜕𝑥
𝜕2 𝑣
𝑀(𝓁) = 𝐸𝐼 2 (𝓁) = 0 ⇒ −𝐴𝑘 2 cos(𝑘𝓁) − 𝐵𝑘 2 sin(𝑘𝓁) = 0, (7.35)
𝜕𝑥
De las ecuaciones anteriores y teniendo en cuenta que 𝑘 > 0 se obtiene que se deben
cumplir las siguientes condiciones
𝐴 = 𝐵 sin(𝑘𝓁) = 𝐶 = 𝐷 = 0 (7.36)
𝑘𝓁 = 𝑛𝜋 con 𝑛 ∈ ℕ (7.37)
En ese caso cualquier valor de 𝐵 cumple la condición por lo que existen infinitas solucio-
𝜋 2 𝐸𝐼
𝑃𝑐𝑟 = (7.38)
𝓁2
𝐿𝑃 = 𝛽𝓁 (7.39)
Donde 𝛽 depende de los vínculos que tenga el elemento. De esta manera la expresión de
Carga Crítica (7.38) se puede generalizar como,
𝜋 2 𝐸𝐼
𝑃𝑐𝑟 = (7.40)
𝐿2𝑃
𝜋 2 𝐸𝐼
𝑃𝑐𝑟 =
𝓁2
mientras que para un elemento en ménsula, con un extremo empotrado y otro libre, se verá
mas adelante que la Carga Crítica es,
𝜋 2 𝐸𝐼
𝑃𝑐𝑟 =
4𝓁2
Se observa que para el caso del elemento simplemente apoyado la luz de pandeo es 𝐿𝑃 = 𝓁 y
por lo tanto 𝛽 = 1, mientras que para el caso del elemento en ménsula la luz de pandeo es
𝐿𝑃 = 2𝓁 y por lo tanto 𝛽 = 2.
Se puede observar en la Figura 7.5 que el caso del elemento en ménsula puede ser asimila-
do al caso de un elemento simplemente apoyado de longitud 2𝓁. Resolviendo analíticamente
otros casos y realizando un razonamiento similar se puede ver que la luz de pandeo es tam-
bién la distancia entre dos puntos de momento nulo. Esto último implica que la longitud de
pandeo también se puede interpretar como aquella que separa puntos con curvatura acotada,
debido a que los puntos de momento nulo coinciden con los puntos de inflexión de la elástica.
7.3.3. Esbeltez
Se define la esbeltez 𝜆 de un elemento comprimido de la siguiente manera,
𝐿𝑃
𝜆= (7.41)
𝜌
𝐼
𝜌2 = , (7.42)
𝐴
tal como se vió en la Sección 6.1.
Si se sustituye en la Ecuación (7.40) se tiene
𝜋 2 𝐸𝐴𝜌2 𝜋 2 𝐸𝐴
𝑃𝑐𝑟 = = 2
𝐿2𝑃 𝐿𝑃
(𝜌)
Obteniendo que,
𝜋 2 𝐸𝐴
𝑃𝑐𝑟 = (7.43)
𝜆2
De igual manera y usando que 𝜎 = 𝑃/𝐴 se puede definir una tensión de compresión
crítica como
𝜋2𝐸
𝜎𝑐𝑟 = 2 (7.44)
𝜆
Cuyo valor tiene como cota superior la tensión admisible 𝜎𝑎𝑑𝑚 del material.
Hasta ahora se ha estudiado el fenómeno de pandeo con propiedades geométricas 𝐼 y 𝜌
genéricas. Sin embargo, si los vínculos del elemento son iguales en ambas direcciones la carga
crítica quedará determinada por el menor momento de inercia. Por otra parte, para analizar
el fenómeno de pandeo en elementos con vínculos distintos en cada dirección, será nece-
sario estudiar ambas direcciones por separado con el momento de inercia que corresponda,
quedando en general determinado por la dirección de mayor esbeltez.
𝑒 𝑁
𝜕4 𝑣 𝜕2 𝑣
𝐸𝐼 (𝑥) + 𝑃 (𝑥) = 0
𝜕𝑥 4 𝜕𝑥 2
La solución y sus derivadas están dadas por:
𝑉 (𝓁) = 0, (7.50)
𝑣′ (0) = 0 ⇔ 𝐵𝑘 + 𝐶 = 0, (7.54)
𝐵=0 (7.55)
𝜕𝑣2
𝑀(𝓁) = −𝑃𝑒 ⇔ (𝓁) = −𝑘 2 𝑒 (7.56)
𝜕𝑥 2
por lo tanto
−𝑘 2 𝐴 cos(𝑘𝓁) − 𝑘 2 𝐵 cos(𝑘𝓁) = −𝑘 2 𝑒 (7.57)
Esto es equivalente a
𝐴 cos(𝑘𝓁) = 𝑒 (7.58)
Finalmente usando la condición de desplazamiento nulo en el empotramiento se tiene
𝑣(0) = 0 ⇔ 𝐴 + 𝐷 = 0 (7.59)
Es necesario estudiar de forma independiente las soluciones obtenidas para los casos con
y sin imperfecciones.
Si se verifica
𝜋
𝑘𝓁 = + 𝑛𝜋 con 𝑛 ∈ ℕ, (7.61)
2
entonces cualquier valor de 𝐴 cumple la condición por lo que existen infinitas soluciones. En
particular la menor directa asociada a estos valores de 𝑘 es
𝜋 2 𝐸𝐼
𝑃𝑐𝑟 = (7.62)
4𝓁2
1
𝑣(𝓁) = 𝑒 1 − (7.64)
( cos(𝑘𝓁) )
𝜕2 𝑣 cos(𝑘𝑥)
2
(𝑥) = −𝑘 2 𝑒 . (7.67)
𝜕𝑥 cos(𝑘𝓁)
1.5
1
𝑘𝓁
= 0.5 𝑒
0.5
= 0.2 𝑒
= 0.1 𝑒
0 =0 𝑒
0 5 10 15 20 25 30 35
Desplazamiento 𝑣(𝓁)
7.4 Ejercicios
Nota: Salvo aclaración contraria, se considerará que las piezas pueden pandear en torno
a cualquiera de sus direcciones principales.
Ejercicio 7.1
Determinar los coeficientes 𝛽 tal que la luz de pandeo es 𝐿𝑃 = 𝛽𝓁 para las condiciones de
borde mostradas en la Figura 7.8.
Ejercicio 7.2
La barra horizontal de la figura está soportada por dos columnas. Cada columna está
articulada en su parte superior a la barra horizontal, mientras que el apoyo inferior de la
columna izquierda es un empotramiento y el de la derecha es un apoyo fijo. Ambas columnas
son barras de acero 𝐸 = 200𝐺𝑃𝑎 de sección transversal cuadrada con un ancho de 15𝑚𝑚.
a) Si 𝑎 = 0.4𝑚, ¿Cuál es el valor crítico de la carga 𝑄?
b) Si 𝑎 varía entre 0 y 1 m. ¿Cuál es el valor máximo de 𝑄𝑐𝑟𝑖𝑡𝑖𝑐𝑜 , y cuál es el valor correspon-
diente de 𝑎?
Nota: Estudiar la inestabilidad sólo en el plano de la figura.
Ejercicio 7.3
Una columna rectangular, con dimensiones de sección
transversal 𝑏 × ℎ, está articulada en los extremos A y C,
como se indica en la Figura.
La columna está restringida en el plano de la figura en
la mitad de su altura, pero puede deformarse libremente
en el plano perpendicular al de la figura.
Determinar la relación ℎ/𝑏 tal que la carga crítica sea la
misma para el pandeo en los dos planos principales de
la columna.
Nota: Analizar el problema a partir de los posibles mo-
dos de pandeo obtenidos en un pilar biarticulado.
Ejercicio 7.4
Se tiene la columna de la Figura (con 𝐸𝐼 uniforme) de largo 𝐿. La misma está soportada
en su parte superior por un resorte elástico lineal de constante 𝑟 = 2𝐸𝐼 /𝐿3 contenido en
el 𝑝𝑙𝑎𝑛𝑜 𝑥𝑦, mientras que en el plano perpendicular se encuentra en ménsula (libre en su
extremo superior). Se pretende que la pieza sea capaz de soportar una carga 𝑄 = 200𝑘𝑁 con
una altura 𝐿 = 3𝑚.
Hallar la ecuación que permite determinar la luz de pandeo en el plano del resorte y de-
mostrar que 𝐿𝑝 = 1, 56𝐿. Se sugiere plantear el origen de coordenadas en el extremo superior
del pilar.
Bibliografía
Brady, S. (2014). The Quebec bridge collapse: A preventable failure (part 1). The Structural
Engineer, 92(1):20–21.
Bui, N. N., Ngo, M., Nikolic, M., Brancherie, D., and Ibrahimbegovic, A. (2014). Enriched
Timoshenko beam finite element for modeling bending and shear failure of reinforced
concrete frames. Computers and Structures, 143:9–18.
Cervera Ruiz, M. and Blanco Díaz, E. (2002a). Mecánica de estructuras Libro 1 Resistencia de
materiales. Edicions UPC, 2 edition.
Cervera Ruiz, M. and Blanco Díaz, E. (2002b). Mecánica de estructuras Libro 2 Métodos de
Análisis. Edicions UPC, Barcelona, 2 edition.
Eaton, J. W., Bateman, D., Hauberg, S., and Wehbring, R. (2015). GNU Octave: a high-level
interactive language for numerical computations. CreateSpace Independent Publishing
Platform, 4 edition.
Firth, I. (2010). Lessons from history - the steel box girder story. The Structural Engineer,
88(3):18–24.
Forets, M., Freire Caporale, D., and Pérez Zerpa, J. M. (2022). Combining set propagation
with finite element methods for time integration in transient solid mechanics problems.
Computers & Structures, 259:106699.
Frischkorn, J. and Reese, S. (2015). Solid-beam finite element analysis of Nitinol stents. Com-
puter Methods in Applied Mechanics and Engineering, 291:42–63.
Fuchs, M. B. (2016). The Force Method. In Structures and Their Analysis, pages 195–233.
Springer International Publishing, Cham.
Bibliografía 138
Holzapfel, G. A., Gasser, T. C., and Ogden, R. W. (2000). A new constitutive framework for
arterial wall mechanics and a comparative study of material models. Journal of Elasticity,
61(1-3):1–48.
Krenk, S. and Høgsberg, J. (2013). The Force Method. In Statics and Mechanics of Structures,
pages 227–265. Springer Netherlands, Dordrecht.
Mofid, M., Shadnam, M. R., and McCabe, S. L. (2001). On a new approach for finding the kern
of a beam section. Engineering Structures, 23(9):1080–1084.
Oñate, E. (2013). Structural Analysis with the Finite Element Method Linear Statics, volume 2 of
Lecture Notes on Numerical Methods in Engineering and Sciences. Springer Netherlands,
Dordrecht.
Pereira, J., Xavier, J., Morais, J., and Lousada, J. (2014). Assessing wood quality by spatial varia-
tion of elastic properties within the stem: Case study of Pinus pinaster in the transverse
plane 1. Canadian Journal of Forest Research, 44(2):107–117.
Pérez Zerpa, J., Castrillo, P., and Baño, V. (2017). Development of a method for the identifica-
tion of elastoplastic properties of timber and its application to the mechanical characte-
risation of Pinus taeda. Construction and Building Materials, 139:308–319.
Pilkey, W. D. (2002). Analysis and Design of Elastic Beams. John Wiley & Sons, Inc., New York.
Reddy, J. N. (2002). Energy Principles and Variational Methods in Applied Mechanics. Wiley.
Taroco, E. O., Blanco, P. J., and Feijóo, R. A. (2019). Introduction to the Variational Formulation
in Mechanics. Wiley.
Timoshenko, S. (1940). Strength of Materials, Part I: Elementary theory and problems. D. Van
Nostrand, 2 edition.
Yoo, C. H. and Lee, S. C. (2011). Stability of Structures. In Stability of Structures, pages 521–
523. Elsevier.
Yoshiaki, G., Yasuhito, W., Toshihiro, K., and Makoto, O. (1992). Elastic buckling phenomenon
applicable to deployable rings. International Journal of Solids and Structures, 29(7):893–
909.
Young, W. C. and Budynas, R. G. (2001). Roark’s Formulas for Stress and Strain. McGraw-Hill,
7 edition.
Zienkiewicz, O. C. (1972). Introductory Lectures on the Finite Element Method. Springer Vienna,
Vienna.
Apéndice A
Soluciones de los ejercicios
En este capítulo se presentan los resultados numéricos solución de cada ejercicio práctico.
No se presenta un desarrollo extensivo de la solución, se motiva a los estudiantes a asistir a
práctico para evacuar dudas del desarrollo.
a) Diagrama de cortantes en kN
b)
16.2 × 103 2𝑏
. ⩽ 8 × 106 ⇒ 𝑏 ⩾ (.1448)1/3 ⇒ 𝑏 = 15 cm
𝑏(2𝑏)3 /12 2
d) cero.
Ejercicio 1.2
a) Directas en kN
cortantes en kN
momentos en kN m
Ejercicio 1.3
b) 2 PNC 6.5.
Ejercicio 1.4
a)
b)
RC
VP
MS
• Máximo: De A a C y de E a F → MC,MAX = 50 kN m.
• Mínimo: De C a E → MC,MIN = -160 kN m.
Ejercicio 1.5
Reacciones,
MA = −5.2 kN m.
RA = 5.04 kN.
Momentos en kN m
Reacciones
RB = 50 kN hacia abajo.
Directas en barras
Barra N [kN]
AB 50
BC 0
AC −70.7
Ejercicio 2.2
Reacciones
Directas en barras
Barra N [kN]
AB −18.0
AE 22.5
BE 46.5
BC −22.5
EC −62.0
EF 100.0
CF −60.0
CD −100.0
FD 80.0
Ejercicio 2.3
Reacciones
RA = 20 kN hacia la izquierda.
Barra N [kN]
AB 15.68
BC 4.42
AD −4.42
BE 0.00
CF 4.42
AE 6.25
DB −7.89
BF 7.89
CE −6.25
DE −14.42
EF −5.58
Ejercicio 2.4
Reacciones
RF = 9.52 kN.
RG = 5.24 kN.
Directas en barras
Barra N [kN]
AB −7.86
BC −7.86
AD 9.45
DE 0.72
EC 9.45
BD 8.58
BE 8.58
FB 9.52
GC 5.24
Ejercicio 2.5
a) Directas en barras
Barra N [kN]
AB 0.086 P
AF −0.172 P
BC 0.259 P
BF 0.172 P
BG −0.172 P
CD 0.313 P
CG −0.982 P
CH −0.928 P
CI −0.094 P
DE 0.086 P
DH −0.226 P
DI 0.226 P
EI −0.172 P
FG −0.172 P
GH 0.233 P
HI −0.118 P
b) PNI 16
Ejercicio 3.2
Ejercicio 3.3
(kN)
(kNm)
Ejercicio 3.4
(kN)
(kNm)
Ejercicio 3.5
(kN)
(kNm)
Ejercicio 3.6
(kN)
(kNm)
Ejercicio 3.7
𝑎)
𝑏)
Ejercicio 3.8
Los ángulos solución son 𝜃𝐵 = −3.89 × 10−3 rad y 𝜃𝐶 = −1.33 × 10−3 rad.
Ejercicio 3.9
1. aplicando que cada barra restringe su deformación axial, se puede determinar como
incógnitas, los giros nodales y el desplazamiento axial del piso C-D.
(kN)
(kN)
(kN.m)
Deformada
Ejercicio 4.2
(kN) (kN)
(kN.m) Deformada
Ejercicio 4.3
(kN) (kN)
(kN.m) Deformada
Ejercicio 4.4
a) 𝑃 = 17.5 kN
b)
Giro en C: 𝜃𝐶 = 5.6 × 10−4 rad.
Giro en D: 𝜃𝐷 = −5.6 × 10−4 rad.
Desplazamiento de M: Δ𝑀 = 2.99 × 10−3 m (hacia abajo).
Ejercicio 4.5
(kN.m)
Ejercicio 4.6
a) 𝑡 = 2.5 cm
b) 𝑘 = 2000 kN/m
(kN)
(kN)
(kNmm)
Deformada
Ejercicio 4.7
a) 𝑘 = 1200 kN/m
b) 𝑎 = 2.82 m
c) 𝜎𝑚𝑎𝑥 = 6.85 MPa
(kN)
(kN)
(kN.m)
Ejercicio 4.8
a)
Giro en C: 𝜃𝐶 = 3.16 × 10−3 .
Desplazamiento de C: Δ𝐶 = 2.53 × 10−2 m (hacia la izquierda).
(kN) (kN)
(kN.m) deformada
Ejercicio 4.9
(kN) (kN)
(kN.m) deformada
Ejercicio 4.10
Incógnitas cinemáticas,
Sist. de coord
⎡𝜃𝑥,2 ⎤ ⎡−1.00 × 10−2 rad ⎤
⎢ 𝑣𝑦,2 ⎥ = ⎢ −2.56 × 10−3 m ⎥
⎢ ⎥ ⎢ ⎥
⎣𝜃𝑧,2 ⎦ ⎣−3.86 × 10 rad ⎦
−3
Ejercicio 5.2
Cortante 𝑉𝑦 (N)
Incógnitas cinemáticas,
Ejercicio 5.3
Sistema de coordenadas
Cortante 𝑉𝑦 (N)
Incógnitas cinemáticas,
Ejercicio 5.4
Nodos, sistema de ejes globales y locales iguales a los del Ejercicio 5.2.
Cortante 𝑉𝑦 (N)
Incógnitas cinemáticas,
Ejercicio 5.5
b) Directa (kN)
d)
Directa (kN)
Ejercicio 5.6
Diagramas en la torre,
Directa (kN) Cortante - Vy,local (kN) Flexión - Mz,local (kN.m)
Ejercicio 5.7
Ejercicio 5.8
5000.0
12.5
50.0
0.0
−750.0
500.0
Directa (kN) Momento 𝑀𝑦 (kNm) Momento 𝑀𝑧 (kNm)
250.0
2.5
5.0
7.5
500.0
37, 5 0.0
17.5
Torsor 𝑀𝑥 (kNm) Cortante 𝑉𝑦 (kN) Cortante 𝑉𝑧 (kN
Ejercicio 5.9
Ejercicio 5.10
Incógnitas cinemáticas,
Cortante 𝑉𝑦 (N)
⎡𝜃𝑥,2 ⎤ ⎡−4.26 × 10−3 rad ⎤
⎢ 𝑣𝑦,2 ⎥ = ⎢ 5.35 × 10−3 m ⎥
⎢ ⎥ ⎢ ⎥
⎣𝜃𝑧,2 ⎦ ⎣ 8.02 × 10 rad ⎦
−3
Ejercicio 6.1
Ejercicio 6.2
Ejercicio 6.3
Ejercicio 6.4
Ejercicio 6.5
a) Vértices del Núcleo Central: (-40, 0), (0, 8.6), (0, −8.6), (12, 10.5), (12,−10.5) y (31,0).
Ejercicio 6.6
a) Distancia de G a vértice del Núcleo Central: 17 cm. b) 𝐹 = 1.53 kN. c) 𝐹 = 1.08 kN. d)
𝐹 = 2.55 kN.
Ejercicio 6.7
b) No existen posiciones de 𝑄 tal que toda la sección esté comprimida (la resultante cae
siempre fuera del Núcleo Central).
c) 𝛼 ⩽ 0.27.
Ejercicio 6.8
Ejercicio 6.9
La línea neutra queda ubicada a 61.1 cm del borde izquierdo. 𝑃𝑎𝑑𝑚 = −81.6 kN.
Ejercicio 6.10
Ejercicio 6.11
a) Reacción Vertical 28.81 kN, Reacción Horizontal: 4.55 kN, Reacción Momento: 3.86
kNm.
b) 𝜎𝑚𝑎𝑥,𝑡𝑒𝑟𝑟𝑒𝑛𝑜 = 86.7 kPa.
Ejercicio 6.12
−𝑀𝑧 −500𝑘𝑁 𝑚 2
𝑁 = −750𝑘𝑁 =𝑒𝑦 = = 𝑚 𝑒𝑧 = 0
𝑁 −750𝑘𝑁 3
Finalmente, se plantea la ecuación de la línea neutra.
𝑒𝑦 𝑒𝑧 𝜌2
1+ 2
𝑦 + 2𝑧 = 0 ⇒ 𝑦 = −1 ⇒ 𝑦 = −0.02𝑚
𝜌 𝜌 𝑒𝑦
c) Dado que se tiene flexión compuesta en solo un sentido, basta con determinar las
tensiones para los puntos con mayor y menor coordenada 𝑦 respecto del baricentro.
𝑁 𝑒𝑦 𝑒𝑧
1 + 2𝑦 + 2𝑧
𝐴( 𝜌 𝜌 )
−750𝑘𝑁 (2/3)𝑚
𝑦 = 0.2𝑚 ⇒ 1+ 0.2𝑚 = −51.56 MPa
(0.42 )𝑚2 ( (0.42 /12)𝑚2 )
−750𝑘𝑁 (2/3)𝑚
𝑦 = −0.2𝑚 ⇒ 1+ × (−0.2𝑚) = 42.19 MPa
(0.42 )𝑚2 ( (0.42 /12)𝑚2 )
𝐺 0.02𝑚
42.19𝑀𝑃𝑎
Ejercicio 6.13
a) 𝛽 = 1, b) 𝛽 = 0.5 c) 𝛽 = 0.7
Ejercicio 7.2
Ejercicio 7.3
ℎ/𝑏 = 2
Ejercicio 7.4
𝑟
𝑘3 + (tan(𝑘𝓁) − 𝑘𝓁) = 0
𝐸𝐼
Apéndice B
Fórmulas útiles
⎛ 𝑐 −𝑠 0 0 ⎞
⎜ 𝑠 𝑐 0 0 ⎟⎟
R𝑒 = ⎜
⎜ 0 0 𝑐 −𝑠 ⎟
⎜ 0 0 𝑠 𝑐 ⎟⎠
⎝
Matriz de rigidez: - Coordenadas globales - K𝑒𝐺 = R𝑒 K𝑒𝐿 R𝑇
⎛ 𝑐2 ⎞ 𝑐𝑠 −𝑐2 −𝑐𝑠
Ω E ⎜ 𝑐𝑠 ⎟ 𝑠2 −𝑐𝑠 −𝑠 2
𝑒 𝑒 ⎜ ⎟
K𝑒𝐺 =
𝐿𝑒 ⎜ −𝑐2 ⎟ −𝑐𝑠 𝑐2 𝑐𝑠
⎜ −𝑐𝑠 ⎟ −𝑠 2 𝑐𝑠 𝑠2
⎝ ⎠
𝑐 = cos(𝛼), 𝑠 = sen(𝛼) siendo 𝛼 el ángulo que forman los ejes globales con los locales
(sentido: eje global hacia eje local antihorario positivo).
176
𝑔ℎ 𝑔ℎ
𝑁 𝑒 (𝑥) = 𝑁0𝑒 + ∑ 𝑋𝑗 𝑁𝑗𝑒 ; 𝑅𝑖 (𝑥) = 𝑅0𝑖 + ∑ 𝑋𝑗 𝑅𝑗𝑖
𝑗=1 𝑗=1
𝑒 𝑔ℎ
𝑛𝑒
𝑁𝑔ℎ+1 (𝑁0𝑒 + ∑𝑗=1 𝑋𝑗 𝑁𝑗𝑒 )
𝛿𝑔ℎ+1 = ∑ 𝓁𝑒
𝑒=1
𝐸𝑒 𝐴𝑒
Método de Slope-Deflection
𝜃2
𝜃1
𝑣1 𝑣2
𝑀1 𝑉1 𝑀2
𝑉2
Desplazamiento relativo:
𝑣2 − 𝑣1
𝜓=
𝐿
Momentos nodales:
{ 2𝐸𝐼 𝑒𝑚𝑝
𝑀12 = 𝐿 (2𝜃1 + 𝜃2 − 3𝜓) + 𝑀12
2𝐸𝐼 𝑒𝑚𝑝
𝑀21 = 𝐿 (𝜃1 + 2𝜃2 − 3𝜓) + 𝑀21
Fuerzas nodales:
⎧
⎪ 𝑀12 +𝑀21 𝑡𝑟𝑎𝑚𝑜−𝑖𝑠𝑜
⎪ 𝐹12 = 𝐿 + 𝐹12
⎨
⎪
⎪
⎩ 𝐹21 = − (𝑀12 +𝑀
𝐿
21 ) 𝑡𝑟𝑎𝑚𝑜−𝑖𝑠𝑜
+ 𝐹21
177
Torsión
𝜃2 − 𝜃1 𝐸
𝑀𝑡 = 𝐺𝐽 Θ, Θ= , 𝐺=
𝑙 2(1 + 𝜈)
𝑀𝑡
𝜏𝑀𝐴𝑋 =
𝑊𝑡
𝑣𝑦,1 , 𝐹𝑦,1
𝜃𝑥,1,𝐿 , 𝑀𝑥,1,𝐿
𝐱𝐺
𝛼𝑦
⎡ 𝐺𝐽 0 0 − 𝐺𝐽𝑙 0 0 ⎤⎥
⎢ 𝑙
⎢ 0 12 𝐸𝐼𝑙3𝑧 6 𝐸𝐼𝑙2𝑧 0 −12 𝐸𝐼𝑙3𝑧 6 𝐸𝐼𝑙2𝑧 ⎥⎥
⎢
⎢ 0 6 𝐸𝐼𝑙2𝑧 4 𝐸𝐼𝑙 𝑧 0 −6 𝐸𝐼𝑙2𝑧 2 𝐸𝐼𝑙 𝑧 ⎥
K𝑒𝐿 = ⎢ 𝐺𝐽 ⎥
⎢− 0 0 𝐺𝐽
0 0 ⎥
⎢ 𝑙 𝑙 ⎥
⎢ 0 −12 𝐸𝐼𝑙3𝑧 −6 𝐸𝐼𝑙2𝑧 0 12 𝐸𝐼𝑙3𝑧 −6 𝐸𝐼𝑙2𝑧 ⎥
⎢ ⎥
⎢ 0 6 𝐸𝐼𝑙2𝑧 2 𝐸𝐼𝑙 𝑧 0 −6 𝐸𝐼𝑙2𝑧 4 𝐸𝐼𝑙 𝑧 ⎥
⎣ ⎦
178
Matriz de rotación:
⎡𝑐𝑜𝑠(𝛼𝑦 ) 0 −𝑠𝑖𝑛(𝛼𝑦 ) 0 0 0 ⎤
⎢ 0 1 0 0 0 0 ⎥⎥
⎢
𝑒 ⎢ 𝑠𝑖𝑛(𝛼𝑦 ) 0 𝑐𝑜𝑠(𝛼𝑦 ) 0 0 0 ⎥
R =⎢ ⎥
⎢ 0 0 0 𝑐𝑜𝑠(𝛼𝑦 ) 0 −𝑠𝑖𝑛(𝛼𝑦 )⎥
⎢ 0 0 0 0 1 0 ⎥
⎢ ⎥
⎣ 0 0 0 𝑠𝑖𝑛(𝛼𝑦 ) 0 𝑐𝑜𝑠(𝛼𝑦 ) ⎦
𝛼𝑦 es el ángulo que forman los ejes globales con los locales (sentido: eje global hacia eje
local como se muestra en la figura).
𝑁 (𝑥) 𝑦𝐴 𝑧𝐴 𝐼𝑦 𝐼𝑧
𝜎(𝑥, 𝑦, 𝑧) = 1+ 2𝑦 + 2𝑧 , 𝜌𝑦2 = ; 𝜌𝑧2 = .
𝐴 ( 𝜌𝑧 𝜌𝑦 ) 𝐴 𝐴
179
𝜎𝑚𝑎𝑥𝑐
𝑁 = ∫ 𝜎𝑑𝐴; 𝜎 = 𝑘.𝑧 = − 𝑧
𝐴 𝑐
Equilibrio de momentos en la sección:
𝑁 (𝑐 − 𝑑) = ∫ 𝜎𝑧𝑑𝐴;
𝐴
𝑐 = 3𝑑
𝜎𝑚𝑎𝑥𝑐 𝑏 𝑐2
𝑁 =−
𝑐 2
2𝑁
𝜎𝑚𝑎𝑥𝑐 = −
𝑏𝑐
Estabilidad Estructural
𝑑2 𝑑2𝑣 𝑑2𝑣 𝑁
𝐸𝐼 + 𝑁 + 𝑞 = 0; 𝑘2 =
𝑑𝑥 2 ( 𝑑𝑥 2 ) 𝑑𝑥 2 𝐸𝐼
𝑑 𝑑2𝑣 𝑑𝑣 𝑑2𝑣
𝑉 (𝑥) = 𝐸𝐼 2 + 𝑁 𝑀(𝑥) = 𝐸 𝐼
𝑑𝑥 ( 𝑑𝑥 ) 𝑑𝑥 𝑑𝑥 2
𝐸 𝐼 𝜋2
𝜎𝑐𝑟,𝑖 = ; 𝐿𝑝,𝑖 = 𝛽𝑖 𝐿; 𝑖 ∈ {𝑦, 𝑧}
𝐴 𝐿𝑝,𝑖 2
√
𝜋2𝐸 𝐿𝑝,𝑖 𝐼𝑖
𝜎𝑐𝑟,𝑖 = 2; 𝜆𝑖 = ; 𝜌𝑖 =
𝜆𝑖 𝜌𝑖 𝐴
𝑀 𝑁
+ ⩽ 𝑚𝑖𝑛{𝜎𝑓 𝑙𝑢𝑒𝑛𝑐𝑖𝑎 , 𝜎𝑐𝑟 }
𝑊 𝐴
180
𝐼𝑥 𝑊𝑥 𝐼𝑦 𝑊𝑦
𝑦
ℎ/2 𝑥 𝑏ℎ3 𝑏ℎ2 ℎ𝑏3 ℎ𝑏2
ℎ/2
𝑏
12 6 12 6
𝑦
181
Cargas M M’ M M’
𝑄
𝑄𝑎𝑏 𝑄𝑎𝑏 𝑄𝑎𝑏 𝑄𝑎𝑏
(𝑙 + 𝑏) (𝑙 + 𝑎) 𝑏 𝑎
2𝑙2 2𝑙2 𝑙2 𝑙2
𝑎 𝑏
𝑄
3 3 1 1
𝑄𝑙 𝑄𝑙 𝑄𝑙 𝑄𝑙
16 16 8 8
𝑙/2 𝑙/2
𝑄 𝑄
3 𝑎 3 𝑎 𝑎 𝑎
𝑄𝑎 1 − ) 𝑄𝑎 1 − ) 𝑄𝑎 (1 − ) 𝑄𝑎 (1 − )
2 ( 𝑙 2 ( 𝑙 𝑙 𝑙
𝑎 𝑎
𝑄 𝑄
1 1 2 2
𝑄𝑙 𝑄𝑙 𝑄𝑙 𝑄𝑙
3 3 9 9
𝑙/3 𝑙/3 𝑙/3
𝑄 𝑄 𝑄
15 15 5 5
𝑄𝑙 𝑄𝑙 𝑄𝑙 𝑄𝑙
32 32 16 16
𝑙/4 𝑙/4 𝑙/4 𝑙/4
𝑞
1 2 1 2 1 2 1 2
𝑞𝑙 𝑞𝑙 𝑞𝑙 𝑞𝑙
8 8 12 12
𝑙
𝑞
𝑞𝑎2 𝑎 2 𝑞𝑎2 𝑎2 𝑞𝑎2 𝑎 𝑎2 𝑞𝑎2 𝑎 𝑎2
2− ) 2− 2 6−8 +3 2 4 −3 2
8 ( 𝑙 8 ( 𝑙 ) 12 ( 𝑙 𝑙 ) 12 ( 𝑙 𝑙 )
𝑎 𝑏
𝑞
9 2 7 2 11 2 5 2
𝑞𝑙 𝑞𝑙 𝑞𝑙 𝑞𝑙
128 128 192 192
𝑙/2 𝑙/2
𝑞 𝑞
𝑞𝑎2 𝑎 𝑞𝑎2 𝑎 𝑞𝑎2 𝑎 𝑞𝑎2 𝑎
3−2 ) 3−2 ) 3−2 ) 3−2 )
4 ( 𝑙 4 ( 𝑙 6 ( 𝑙 6 ( 𝑙
𝑎 𝑎
𝑞
𝑞𝑎𝑏𝑐 𝑐2 𝑞𝑎𝑏𝑐 𝑐2 𝑞𝑐 𝑐2 𝑞𝑐 2 𝑐2
𝑙+𝑏− 𝑙+𝑎− 𝑎𝑏2 + (𝑙 − 3𝑏) 𝑎 𝑏+ (𝑙 − 3𝑎)
𝑐/2𝑐/2 2𝑙2 ( 4𝑎 ) 2𝑙2 ( 4𝑏 ) 𝑙2 [ 12 ] 𝑙2 [ 12 ]
𝑎 𝑏
𝑞
𝑞𝑙𝑎 𝑎2 𝑞𝑙𝑎 𝑎2 𝑞𝑙𝑎 𝑎2 𝑞𝑙𝑎 𝑎2
3− 2 3− 2 3− 2 3− 2
16 ( 𝑙 ) 16 ( 𝑙 ) 24 ( 𝑙 ) 24 ( 𝑙 )
𝑏 𝑎 𝑏
𝑞
𝑞𝑎2 𝑎 𝑎2 𝑞𝑎2 𝑎2 𝑞𝑎2 𝑎 𝑎2 𝑞𝑎2 𝑎 𝑎2
40 − 45 + 12 2 10 − 6 2 10 − 15 + 6 2 5 −4 2
120 ( 𝑙 𝑙 ) 60 ( 𝑙 ) 30 ( 𝑙 𝑙 ) 20 ( 𝑙 𝑙 )
𝑎 𝑏
𝑞
7 2 1 2 1 2 1 2
𝑞𝑙 𝑞𝑙 𝑞𝑙 𝑞𝑙
120 15 30 20
𝑙
182
𝑑 𝜋𝐷4 𝑑4 𝜋𝐷3 𝑑4
1− 4 1− 4
32 ( 𝐷 ) 16 ( 𝐷 )
𝐷
𝑎
𝜋𝑎3 𝑏3 𝜋𝑎𝑏2
(𝑎2 + 𝑏2 ) 2
𝑎
𝑏 𝑏
/2
√
2𝑎
𝑎/
3 4 𝑎3
2𝑎
𝑎
𝑎/
80 20
/2
𝑎/2 𝑎/2
𝑎/2
𝑏/2 𝑏/2
𝑡
𝑠 𝑠𝑡 3 𝑠𝑡 2
́
𝑙𝑖𝑛𝑒𝑎 𝑚𝑒𝑑𝑖𝑎 3 3
𝑠1
𝑡1
𝑡2 𝑠𝑖 𝑡𝑖3 𝐽
∑𝑛𝑖
𝑠2
𝑡3
3 𝑡𝑚𝑎𝑥
́
𝑡4
𝑠4 𝑠3
𝑡
Ω 4Ω2
2Ω𝑡𝑚𝑖𝑛
𝑑𝑠
∫𝑠
𝑠 𝑡
́
𝑙𝑖𝑛𝑒𝑎 𝑚𝑒𝑑𝑖𝑎
Tabla B.1: 𝐽 corresponde a la inercia torsional. 𝑊𝑡 correspondé al módulo resistente de la sección a torsión. (*): Ver
valores de parámetros 𝛼 y 𝛽 en la Tabla B.2.
Material en revisión. Se agradecerán los avisos de errores vía gitlab.fing.edu.uy/jorgepz/libroResMat2/-/issues/new
“libroResMat2” — 2022/9/9 — 11:53 — page 183 — #188
183
Tabla B.2: Parámetros 𝛼 y 𝛽 para el cálculo de propiedades seccionales de torsión de una sección rectangular.