Universidad Tecnológica de Panamá
Trabajo Final de Modelación
de Estructuras
Presentado por:
Kevin Chirú
Guillermo Alvarado
Luis Jaén
Andony Torres
Reina Pinzón
Presentado a:
Profesor Daniel Ulloa
Tabla de Contenido
01 02 03
Resolución del Resolución del Resolución del
Primer Parcial Segundo Parcial Tercer Parcial
• Problema 1 • Problema 2 • Problema 1
• Problema 2 • Problema 2 (parcial • Problema 2
alterno) • Problema 3 (adicional)
Examen Parcial 1
Problema 1
La siguiente estructura tiene un modulo de elasticidad e y momento de inercia I. ignore el efecto
de esfuerzos axiales.
DATOS: I=4.5e4 cm4, E=210 Gpa, Hb=20kN, Mb=115.0 kNm
Empleando el principio de superposición:
a) Dibuje la(s) estructura(s) primaria(s) y la(s) estructura(s) con la(s) reacción(es) redundante(s)
u u u
dB1 dB2 dB3
HB MB DB/A
B´ B´ B´
dA3
dA1 dA2
HA HA
H0
v v v
dA1+dA2-dA3=0 HA
dB1+dB2-dB3=dB
Fórmulas de Deflexiones y Pendientes
L=11.0m; x=6.0m
L=11.0m; a=6.0m
L=11.0m; x=6.0m
REEMPLZANDO CADA UNA DE LAS ECUACIONES
DE PENDIENTE DEFLEXIÓN PARA CADA CASO
𝐻𝐵𝑎2 3ℎ−𝑎
d𝐴1 = 6𝐸𝐼
= 0.034286𝑚
𝑀𝐵 𝑎2
d𝐴2 = = 0.021905𝑚
2𝐸𝐼
−𝑀𝐻𝐴𝑎3
d𝐴3 = 3𝐸𝐼 = 0.000762𝐻𝐴
𝐻𝐴𝑎2 3𝐿−𝑎
d𝐵3 = 6𝐸𝐼
= −0.12643𝑚
𝐻𝐴𝑎3 𝐻𝐴𝑎2 𝐿 − 𝑎
d𝐵3 = dA3+D𝐵/A = + = −0.12643𝑚
3𝐸𝐼 2𝐸𝐼
𝐻𝐵ℎ3 𝐻𝐵ℎ2 𝐻𝐴𝑎3 𝐻𝐴𝑎2 𝐿 − 𝑎
dB =dB1 +dB2 +dA3 + A3* (h-a) =
3𝐸𝐼
+
2𝐸𝐼
+
3𝐸𝐼
+
2𝐸𝐼
= 0.041093𝑚
RESOLVIENDO LOS SISTEMAS DE ECUACIONES DE COMPATIBILIDAD
dA1+dA2-dA3=0 HA = −73.75𝑘𝑁
𝑀𝑂 = 0 𝐻𝐵ℎ + 𝑀𝐵 + 𝐻𝐴𝑎 + 𝑀𝑂 𝑀𝑂 = 170.5𝑘𝑁𝑚
=0
𝐻𝑖 = 0 𝐻𝑂 = 53.75𝑘𝑁
Problema 2 parcial 1
método de la doble integral
Reacciones:
𝑀𝐶 = 7200𝑙𝑏 ∗ 𝑓𝑡 Planteamiento para la formula de la doble
𝑅𝐶 = 1200 𝑙𝑏 integral
Transformación de unidades 𝑑 2 𝑦 𝑀𝑥
=
𝐸 = 10𝑥106 𝑙𝑏ൗ 2 ≈ 1.44𝑥109 𝑙𝑏ൗ 2 𝑑𝑥 2 𝐸𝐼
𝑖𝑛 𝑓𝑡
4 −3
𝐼 = 40𝑖𝑛 ≈ 1.93𝑥10 𝑓𝑡 4
-Forma de la ecuación a utilizar
(𝐸𝐼)𝑑 2 𝑦 = ඵ 𝑀𝑥 𝑑𝑥 2
Sección a-B de la viga
valido de 0 ≤ 𝑥 ≤ 6
• Aplicando el método de la doble integral
𝐸𝐼𝑑 2 𝑦 = ඵ −100𝑥 2 𝑑𝑥 2
𝑥3
𝐸𝐼𝑑𝑦 = න −100 + 𝐶1 𝑑𝑥
3
𝑥4
𝐸𝐼𝑦 = −25 + 𝐶1 𝑥 + 𝐶2
3
Se realiza sumatoria de momentos en el punto
indicado:
𝑥
𝑀𝑥 + 200 𝑥 =0
2
𝑀𝑥 = −100𝑥 2
Sección b-c
valido de 6 ≤ 𝑥 ≤ 10
Aplicando doble integral
2𝐸𝐼𝑑 2 𝑦 = ඵ −1200𝑥 + 3600 𝑑𝑥 2
2𝐸𝐼𝑑𝑦 = න(−600 𝑥 2 + 3600𝑥 + 𝐶3 )𝑑𝑥
2𝐸𝐼𝑦 = −200𝑥 3 + 1800𝑥 2 + 𝐶3 𝑥 + 𝐶4
𝑀𝑥 + 200 6 𝑥 − 3 = 0
𝑀𝑥 = −1200 𝑥 − 3
𝑀𝑥 = −1200𝑥 + 3600
Determinación de constantes de integración
● Para encontrar C3
Cuando x=10 ; φ=0 (ángulo de rotación)
‘dy=φ ; representa la rotación del elemento como la derivada de la deformación o desplazamiento’
𝐸𝐼𝑑𝑦 = −600𝑥 2 + 3600𝑥 + 𝐶3
0 = −600 102 + 3600 10 + 𝐶3
𝐶3 = 24000
Debido a que es un empotramiento en el extremo derecho podemos determinar que tanto el
ángulo de rotación como el desplazamiento en este punto son iguales a 0
DETERMINACIÓN DE CONSTANTES DE INTEGRACIÓN
● Para determinar C4
Cuando x=10; y=0
𝐸𝐼𝑦 = −200𝑥 3 + 1800𝑥 2 + 2400𝑥 + 𝐶4
Remplazando los valores
0 = −200 103 + 1800 102 + 2400 10 + 𝐶4
𝐶4 = −220000
Determinación de constantes de integración
● Para determinar C1
En el punto x=6 existe una compatibilidad donde el desplazamiento y el ángulo de rotación
utilizando cualquier ecuación de los segmentos debe ser el mismo ya que representan un punto
en común
En x=6 𝑑𝑦𝐴−𝐵 = 𝑑𝑦𝐵−𝐶
𝑥3 −600𝑥 2 + 3600𝑥 + 24000
−100 + 𝐶1 = ( )
3 2
2
100 63
𝐶1 = −300 6 + 1800 6 + 12000 +
3
𝐶1 = 19200
Determinación de constantes de integración
● Para determinar C2
Cuando x=6 ; 𝑦𝐴−𝐵 = 𝑦𝐵−𝐶
25 4
− 𝑥 + 19200𝑥 + 𝐶2 = −100𝑥 3 + 900𝑥 2 + 12000𝑥 − 110000
3
8400 + 𝐶2 = −27200
𝐶2 = −35600
Ecuaciones de elástica de cada sección
● Para 0 ≤ 𝑥 ≤ 6
−25𝑥 4 19200𝑥 35600
𝑦𝐴−𝐵 = + −
3𝐸𝐼 𝐸𝐼 𝐸𝐼
● Para 6 ≤ 𝑥 ≤ 10 (denominadores simplificados entre 2 debido al 2i de la segunda
sección)
−100𝑥 3 900𝑥 2 12000𝑥 110000
𝑦𝐵−𝐶 = + + −
𝐸𝐼 𝐸𝐼 𝐸𝐼 𝐸𝐼
Desplazamiento máximo de la viga
● El desplazamiento máximo de la viga al tratarse de un canto libre podemos deducir que
ocurre en el extremo opuesto al empotrado por ende el desplazamiento máximo ocurre
en el punto A cuando x=0, tomando esto en cuenta se debe utilizar la ecuación que
cumpla para esta sección de la viga. En este caso la ecuación perteneciente a esa
sección de la viga es 𝑦𝐴−𝐵
Cuando x=0
25 04 19200 0 35600
𝑦𝐴−𝐵 =− + −
3𝐸𝐼 𝐸𝐼 𝐸𝐼
35600
𝑦𝐴−𝐵 =− = 0.013𝑓𝑡 ≈ 0.15𝑖𝑛
2779200
Parcial 2-Problema 2
La figura muestra una barra uniforme cargada con su propio peso (asumido
concentrado en el medio de la barra) y una fuerza P aplicada en el extremo
inferior de la barra. La barra tiene un módulo de elasticidad igual a E y una
sección transversal A, su densidad es ρ y la fuerza de gravedad está
definida por el valor g.
2𝐸𝐴
La barra está unida a un resorte cuya rigidez es k. Se nos indica que 𝑘 =
𝐿
Nota: se asume que todos los valores están dados en unidades consistentes
a) Cree un modelo, para aproximar el comportamiento 𝑤 =𝑚∗𝑔
de la estructura real, tomando en consideración la
𝑉 = 𝐴 ∗ 2𝐿
información suministrada. 𝑚 = 𝜌∗𝑉
𝑤 = 2𝜌𝑔𝐴𝐿
b) Calcule las matrices de rigidez individuales de cada
elemento de su estructura.
1 1 −1 2 1 −1 3 1 −1
𝐾 = 𝑘1 𝐾 = 𝑘2 𝐾 = 𝑘3
−1 1 −1 1 −1 1
Ensamble de Matriz Global
Desplazamientos
conocidos 𝑢𝑠 Cargas desconocidas
𝑃𝑓
𝑘1 𝐾 0 −𝑘1 𝐾𝑠𝑠
0 𝑑1 𝑅1
𝑠𝑓
0 𝑘3 0 −𝑘3 𝑑2
= −𝑅2
−𝑘1 𝐾𝑓𝑓 0 𝑘1 + 𝑘2 𝐾𝑓𝑠 −𝑘2 𝑑3 𝑤
0 −𝑘3 −𝑘2 𝑘2 + 𝑘3 𝑑4 𝑃
Desplazamientos Cargas conocidas 𝑃𝑠
desconocidos 𝑢𝑓
𝑢 𝑃
𝐾
c) Escriba la ecuación matricial de equilibrio.
𝑘1 0 −𝑘1 0 𝑑1 𝑅1
0 𝑘3 0 −𝑘3 𝑑2
= 𝑅2
−𝑘1 0 𝑘1 + 𝑘2 −𝑘2 𝑑3 −𝑤
0 −𝑘3 −𝑘2 𝑘2 + 𝑘3 𝑑4 −𝑃
Condiciones de borde 𝑑1 = 𝑑2 = 0
d) Considerando las condiciones de borde, obtenga y escriba la ecuación
matricial reducida del sistema estructural.
𝑘1 + 𝑘2 −𝑘2 𝑑3 𝑤
=
−𝑘2 𝑘 2 + 𝑘 3 𝑑4 𝑃
e) Obtenga los desplazamientos nodales del modelo estructural. −𝐸𝐴 −𝐸𝐴 𝐸𝐴
1+2
2𝐸𝐴 𝐸𝐴 𝐴𝑑𝑗 = −1 =
− 𝐿 𝐿 𝐿
𝐿 𝐿 𝑑3 = 𝑤
𝐸𝐴 3𝐸𝐴 𝑑4 𝑃
− −𝐸𝐴 2+1
−𝐸𝐴 𝐸𝐴
𝐿 𝐿 𝐴𝑑𝑗 = −1 =
Cálculo de la matriz inversa 𝐿 𝐿 𝐿
−1 −1
𝐾𝑓𝑓 𝐾𝑓𝑓 𝑢𝑓 = 𝑃𝑓 𝐾𝑓𝑓
3𝐸𝐴 2+2
2𝐸𝐴 2𝐸𝐴
𝐴𝑑𝑗 = −1 =
𝐿 𝐿 𝐿
−1 1 𝑡
𝐾𝑓𝑓 = ∗ 𝐴𝑑𝑗 𝐾𝑓𝑓
𝐾𝑓𝑓 𝑡
3𝐸𝐴 𝐸𝐴
Cálculo de la determinante −1 1
𝐾𝑓𝑓 = ∗ 𝐿 𝐿
2𝐸𝐴 3𝐸𝐴 𝐸𝐴 𝐸𝐴 5𝐸 2 𝐴2 5𝐸 2 𝐴2 𝐸𝐴 2𝐸𝐴
𝐾𝑓𝑓 = − − − = 𝐿2
𝐿 𝐿 𝐿 𝐿 𝐿2 𝐿 𝐿
Cálculo de las adjuntas 3𝐸𝐴 𝐸𝐴
𝑑3 1 𝑤
𝐴𝑑𝑗
2𝐸𝐴
= −1 1+1
3𝐸𝐴
=
3𝐸𝐴 = ∗ 𝐿 𝐿
𝑑4 5𝐸 2 𝐴2 𝐸𝐴 2𝐸𝐴 𝑃
𝐿 𝐿 𝐿
𝐿2 𝐿 𝐿
3𝐸𝐴 𝐸𝐴 𝐸𝐴 2𝐸𝐴
𝑤 +𝑃
𝑤
𝐿
+𝑃
𝐿 = 5𝐿 3𝑊 + 𝑃 𝑑4 = 𝐿 𝐿 = 𝐿 𝑊 + 2𝑃
𝑑3 = 5𝐸 2 𝐴2 5𝐸𝐴
5𝐸 2 𝐴2 𝐸𝐴
𝐿2
𝐿2
f) Adicional: calcule las reacciones en los apoyos.
−𝐸𝐴
∗ 𝑑3 = 𝑅1
𝐿
−𝐸𝐴 30𝐿2 𝜌𝑔 5𝑃𝐿
∗ + = 𝑅1
𝐿 𝐸 𝐸𝐴
𝑅1 = −30𝐿𝐴𝜌𝑔 − 5𝑃
−2𝐸𝐴
∗ 𝑑4 = −𝑅2
𝐿
−2𝐸𝐴 2𝐿2 𝜌𝑔 4𝑃𝐿
∗ + = 𝑅2
𝐿 𝐸 5𝐸𝐴
8
−4𝐿𝐴𝜌𝑔 − 𝑃 = 𝑅2
5
Examen Parcial 2
Problema 2- Alterno
Paso 1
Calculamos el valor de k para cada uno de los elementos. desconocidas mediante
la ecuación.
Paso 2
Se determina la matriz de rigidez de cada elemento
Examen Parcial 3
Problema 1
La cercha mostrada en la figura consiste de dos barras (elementos de cercha) y un elemento de
resorte, cuya constante de rigidez es k= 2EA/L.. La cercha está cargada con una fuerza puntual
aplicada en el nudo 1 como se indica en la figura. Se indica la longitud, área transversal y
módulo de elasticidad de los elementos de cercha.
Observe que el sistema global de coordenadas está rotado 45°, como se indica.
Datos obtenidos para la resolución
Rigidez da cada elemento Matriz de rigidez global Vector de fuerzas de la
reducida estructura
0.5 0.5 2 −1 1 𝑃
𝐾 (1) = 𝑘 (1) 𝐾 (𝐺) = 𝑘 𝑃
0.5 0.5 −1 1 2 −𝑃 −𝑃
Una vez 0
3 0 organizado
−𝑃 0
0.5 −0.5 4 −𝑃
𝐾 (2) = 𝑘 (2) 5 𝑃
−0.5 0.5 𝑂 𝑃
6
1 −1
𝐾 (3) = 𝑘 (3)
−1 1
1. La matriz de rigidez, en coordenadas globales, del elemento de cercha se define como:
𝑒
(𝑒)
𝐸 𝐴(𝑒) 𝑎 −𝑎 𝑎 −𝑎
[𝐾 ]= = 𝑘 (𝑒)
𝐿(𝑒) −𝑎 𝑎 −𝑎 𝑎
¿Cuál es el valor de a?
1 −1 2 1 −1 𝑐 𝑠
a= a= 𝑐 𝑠𝑐 a= 𝑘 (𝑒) a=
−𝑠 𝑐
−1 1 𝑠𝑐 𝑠2 −1 1
A B C D
2. La matriz de rigidez, en coordenadas globales, de cualquiera de los elementos de la cercha se
define como:
¿De cuál de los 3 elementos se trata?
𝑒
𝐸 𝐴(𝑒) 𝑎 −𝑎 𝑎 −𝑎 𝐸𝐴 1 1 1
[𝐾 (𝑒) ] = = 𝑘 (𝑒) , 𝑆𝑖 𝑘 (𝑒) = 𝑦𝑎=
𝐿(𝑒) −𝑎 𝑎 −𝑎 𝑎 𝐿 2 1 1
1 2 3
3. ¿Cuál es la matriz de rigidez reducida de la estructura? (en coordenadas globales)
𝐸𝐴 2 −1 𝐸𝐴 5 −1 2𝐸𝐴 2 0 𝐸𝐴 5 1
2𝐿 −1 1 𝐿 −1 1 𝐿 0 2 2𝐿 1 1
A B C D
4. El vector {f} de fuerzas en los nudos de la estructura (en coordenadas globales), que
incluye las fuerzas aplicadas en los nudos y las reacciones en los apoyos, es:
𝑃 𝑃 −𝑃 𝑃 𝑃
𝑃 −𝑃 −𝑃 −𝑃 𝑃
𝑃 0 0 0 −𝑃
−𝑃 0 0 0 𝑃
−𝑃 −𝑃 𝑃 −𝑃 𝑃
0 −𝑃 −𝑃 𝑃 0
A B C D E
Problema 02
Paso 1 Paso 2
Se identifica la matriz de rigidez reducida Aplicamos condiciones de borde
768 -192
-192 432
Paso 3 PASO 5
Calculamos los desplazamientos desconocidos Cálculo de las fuerzas
𝑃𝑠 = 𝑘𝑠𝑓 * 𝑈𝑓
-256 -192
0.017578
-192 -144
𝑃𝑠 =
-512 384 * 0.007812
384 -288
𝑅1𝑋 -6
PASO 4 𝑅1𝑦 -4.5
=
Aplico condiciones de borde para calcular las fuerzas 𝑅3𝑥 -6
𝑅3𝑦 4.5
Paso 1
-Encontrar los desplazamientos desconocidos mediante la
ecuación.
Paso 2
Se calculan las fuerzas desconocidas mediante la ecuación
¡Gracias!
CREDITS: This presentation template was created by
Slidesgo, including icons by Flaticon and infographics
& images by Freepik
