Tarea 2 - Rodríguez - Gutiérrez

Integrantes: Sofía Rodríguez y Natalia Gutiérrez
Departamento de Ingeniería Civil y Ambiental

Facultad de Ingeniería
Análisis Avanzado de Estructuras ICYA 4422, 2023‐10
Tarea 2
Profesor Juan Carlos Reyes
jureyes@uniandes.edu.co
Objetivos y parámetros
Objetivo : Analizar estructuras usando análisis estático lineal. Codificar el método

de análisis estático lineal en Matlab. Aplicar los conceptos de matriz de
rigidez y cargas debidas a deformaciones mecánicas y no mecánicas
como presfuerzo, temperatura, etc.
Parámetros : La tarea debe ser desarrollada en grupos de máximo dos estudiantes y
debe ser entregada antes del día domingo 25 de febrero antes de las 11:30
p.m. en Bloque Neón. La calificación máxima es 5.0 y la mínima es 0.0
Problema 1 (implementación manual)

Considere el puente continuo postensado de la Figura P1 sometido a una carga vertical
distribuida de 190 kN/m (esta carga equivalente incluye aproximadamente las cargas muerta y
viva). Simplifique el modelo estructural despreciando los desplazamientos horizontales de todos
los nodos y cambiando adecuadamente los apoyos (y tipos de elemento). La carga total de
presfuerzo es de 20 MN aplicada con excentricidades de 0, ‐1.161, 0.675, ‐1.161 y 0 en las
secciones a, b, c, d y e, respectivamente; la carga de 20MN ya incluye el efecto de algunas
pérdidas de presfuerzo; no considere pérdidas de presfuerzo adicionales. Suponga que la
trayectoria de todos los cables es parabólica en cada luz. El punto f sufre un asentamiento de 3
cm hacia abajo con respecto a los puntos a y e.
a) Calcule y dibuje el diagrama de momentos de la viga debido a la carga vertical distribuida.
b) Calcule y dibuje el diagrama de momentos de la viga debido al efecto del asentamiento en
el nodo f.
c) Calcule y dibuje el diagrama de momentos de la viga debido a la carga de presfuerzo.
d) Sume los diagramas de momento obtenidos en a), b) y c) y concluya sobre el efecto del
presfuerzo.
e) Verifique los resultados con SAP2000 o con Opensees.
ICYA 4422 – Primer semestre de 2023 – Pag. 1


Figura P1. Vista isométrica, longitudinal y sección transversal del puente

a) Calcule y dibuje el diagrama de momentos de la viga debido a la carga vertical
distribuida
A continuación, se realizan los pasos para obtener el diagrama de momentos en la viga

mediante el método de matriz de rigidez
Modelo estructural simplificado
Paso 1. Dado que el enunciado señala que los desplazamientos horizontales se desprecian se
realizan los cambios de apoyo requeridos para cumplir esta condición, además, teniendo en
cuenta que los giros en los apoyos se consideran triviales se realizan los siguientes cambios en
los elementos


Ilustración 1. Grados de libertad de la estructura
Luego de realizar dichas transformaciones se obtiene un grado de libertad en el nodo central

(nodo c) de la estructura.
Pasos 2 y 3
U1=1
Ilustración 2. U1=1


Ilustración 3. Fuerzas al hacer U1=1
3𝐸𝐼 3𝐸𝐼
𝐾11 =
+
𝐿 𝐿
3 ⋅ 28 × 106 𝑘𝑃𝑎 ⋅ 4.02 𝑚4 3 ⋅ 28 × 106 𝑘𝑃𝑎 ⋅ 4.02 𝑚4
𝐾𝑙𝑙 = +
34 𝑚 42 𝑚
𝐾11 = 9931764.71 + 8040000
𝐾11 = 17971764.7 𝑘𝑁 − 𝑚
[𝐾𝑙𝑙 ] = [17971764. 𝑘𝑁 − 𝑚]
Paso 4.
{𝐹𝑙 } = {𝐹𝑛 } − {𝐹0 }
{𝐹𝑛 } = {0}
Elemento ac
Ilustración 4. Momento en el extremo del elemento ac por acción de la carga distribuida


𝑤𝐿2 190 ∙ 342

{𝑃0 }𝑎𝑐 = − =− = −27455 𝑘𝑁 − 𝑚
8 8
Elemento ce
Ilustración 5. Momento en el extremo del elemento ce por acción de la carga distribuida
𝑤𝐿2 190 ∙ 422

{𝑃0 }𝑐𝑒 = = = 41895 𝑘𝑁 − 𝑚
8 8
Vector de fuerzas de la estructura
{𝐹0 } = {𝑃0 𝑎𝑐 + 𝑃0 𝑐𝑒 } = {−27455 + 41895} = {14440 𝑘𝑁 − 𝑚}
{𝐹𝑙 } = {𝐹𝑛 } − {𝐹0 }
{𝐹𝑙 } = {0} − {14440} = {−14440 𝑘𝑁 − 𝑚}
Paso 5.
{𝑈𝑙 } = [𝐾𝑙𝑙 ]−1 {𝐹𝑙 }

1
{𝑈𝑙 } = ⋅ −14440
17971764.6
{𝑈𝑙 } = −0.000803 𝑟𝑎𝑑
Paso 6.
Para calcualr las fuerzas internas por el efecto de las cargas aplicadas sobre el elemento, es
necesario sumar al hacer 𝑈1 = 1 para calcular la matriz de rigidez y las fuerzas por la carga


distibuida. Asimismo, como el diagrama de momento depende del diagrama de cortante, es
necesario calcular también el diagrama de cortante.
3𝐸𝐼
3𝐸𝐼 ⋅ 𝑈1
⋅ 𝑈1 𝐿
𝐿 e
a c
3𝐸𝐼
3𝐸𝐼 3𝐸𝐼 3𝐸𝐼 ⋅ 𝑈1
⋅ 𝑈1 ⋅ 𝑈1 ⋅ 𝑈1 𝐿2
𝐿2 𝐿2 𝐿 2
Ilustración 6. Efecto de U1=1
190 kN/m 190 kN/m

c e
a c e
3𝑤𝐿 5𝑤𝐿 5𝑤𝐿 3𝑤𝐿

8 8 8 8
Ilustración 7. Efecto de las cargas sobre el elemento
3𝐸𝐼 3𝑤𝐿
𝑉𝑎𝑐 = 2
⋅ 𝑈1 +
𝐿 8
𝑘𝑁
3 ⋅ 28000 𝑀𝑝𝑎 ⋅ 4.02 𝑚4 3 ⋅ 190 𝑚 ⋅ 34 𝑚
𝑉𝑎𝑐 = ⋅ −0.000803 +
34 𝑚 8
𝑉𝑎𝑐 = 2187.79 𝑘𝑁


3𝐸𝐼 5𝑤𝐿
𝑉𝑐𝑎 = − 2
⋅ 𝑈1 +
𝐿 8
𝑘𝑁
3 ⋅ 28000 𝑀𝑝𝑎 ⋅ 4.02 𝑚4 5 ⋅ 190 𝑚 ⋅ 34 𝑚
𝑉𝑎𝑐 = − ⋅ −0.000803 +
34 𝑚 8
𝑉𝑎𝑐 = 4272.21 𝑘𝑁
3𝐸𝐼 5𝑤𝐿
𝑉𝑐𝑒 = ⋅ 𝑈1 +
𝐿2 8
𝑘𝑁
3 ⋅ 28000 𝑀𝑝𝑎 ⋅ 4.02 𝑚4 5 ⋅ 190 𝑚 42 𝑚
𝑉𝑐𝑒 = ⋅ −0.000803 +
42 𝑚 8
𝑉𝑐𝑒 = 4833.78 𝑘𝑁
3𝐸𝐼 3𝑤𝐿
𝑉𝑒𝑐 = − ⋅ 𝑈1 +
𝐿2 8
𝑘𝑁
3 ⋅ 28000 𝑀𝑝𝑎 ⋅ 4.02 𝑚4 3 ⋅ 190 𝑚 ⋅ 42 𝑚
𝑉𝑒𝑐 = − ⋅ −0.000803 +
422 𝑚 8
𝑉𝑒𝑐 = 3146.22 𝑘𝑁
𝑀𝑎𝑐 = 0 𝑘𝑁 ⋅ 𝑚
3𝐸𝐼 𝑤𝐿2
𝑀𝑐𝑎 = ∙ 𝑈1 + (− )
𝐿 8
𝑘𝑁
3 ∙ 28000𝑀𝑃𝑎 ∙ 4.02 𝑚4 190 𝑚 ∙ 342 𝑚
𝑀𝑐𝑎 = ∙ (−0.000803) + (− )
34 𝑚 8
𝑀𝑐𝑎 = −35435 𝑘𝑁 ⋅ 𝑚
3𝐸𝐼 𝑤𝐿2
𝑀𝑐𝑒 = ∙ 𝑈1 + (− )
𝐿 8


𝑘𝑁
3 ∙ 28000𝑀𝑃𝑎 ∙ 4.02 𝑚4 190 𝑚 ∙ 422 𝑚
𝑀𝑐𝑒 = ∙ (−0.000803) + (− )
42 𝑚 8
𝑀𝑐𝑒 = 35435 𝑘𝑁 ⋅ 𝑚
𝑀𝑒𝑐 = 0 𝑘𝑁 ⋅ 𝑚
190 kN/m 190 kN/m

e
c
35435 𝑘𝑁 ⋅ 𝑚 35435 𝑘𝑁 ⋅ 𝑚
2187.79 𝑘𝑁 4272.21 𝑘𝑁 4833.78 𝑘𝑁 3146.22 𝑘𝑁
Ilustración 8. Fuerzas internas sobre los elementos de la viga
4833.78 𝑘𝑁
2187.79 𝑘𝑁
+ +
a
- -
11.51 m 25.44 m
4272.21 𝑘𝑁 2838.69 𝑘𝑁
Ilustración 9. Diagrama de cortante por efecto de la carga distribuida


26050.7𝑘𝑁 ⋅ 𝑚
12590.7 𝑘𝑁 ⋅ 𝑚
+ +
+
+
a
- -
11.51 m 25.54 m
35435 𝑘𝑁 ⋅ 𝑚
Ilustración 10. Diagrama de momento por acción de la carga distribuida
b) Calcule y dibuje el diagrama de momentos de la viga debido al efecto del

asentamiento en el nodo f.
En primer lugar, es necesario calcular el vector de fuerzas {𝐹𝑙 } debido al asentamiento.
3𝐸𝐼 3𝐸𝐼
⋅ −0.03 𝑚 ⋅ −0.03 𝑚
𝐿2 𝐿2
3𝐸𝐼 3𝐸𝐼
3𝐸𝐼 ⋅ −0.03 𝑚 3𝐸𝐼 ⋅ −0.03 𝑚
⋅ −0.03 𝑚 𝐿3 ⋅ −0.03 𝑚 𝐿3
𝐿3 𝐿3
Ilustración 11. Efecto del asentamiento

{𝐹𝑛 } = 0
3𝐸𝐼 3𝐸𝐼
{𝐹𝑜 } = − ⋅ −0.03 𝑚 + ⋅ −0.03 𝑚
𝐿 𝐿
3 ⋅ 28000 𝑀𝑝𝑎 ⋅ 4.02 𝑚2 3 ⋅ 28000 𝑀𝑝𝑎 ⋅ 4.02 𝑚

{𝐹𝑜 } = − ⋅ −0.03 𝑚 + ⋅ −0.03 𝑚
34 𝑚 42 𝑚


{𝐹𝑜 } = 3020.46 𝑘𝑁 ⋅ 𝑚
{𝐹𝑙 } = {𝐹𝑛 } − {𝐹𝑜 }
{𝐹𝑙 } = {0 − 3020.46 𝑘𝑁 ⋅ 𝑚}
{𝐹𝑙 } = −3020.46 𝑘𝑁 ⋅ 𝑚
Ahora se prosigue a calcular el desplazamiento debido a la acción del asentamiento.
{𝑈𝐿 } = [𝐾𝑙𝑙 ]−1 ⋅ {𝐹𝑙 }
1
{𝑈𝑙 } = ⋅ −2030.46
17971764.7
{𝑈𝐿 } = −0.000168 𝑟𝑎𝑑
Una vez obtenido el valor del desplazamiento, es posible calculas las fuerzas internas en el
extremo de los elementos, superponiendo el efecto de hacer 𝑈1 = 1 para calcular la matriz de
rigidez y el efecto del asentamiento. El procedimiento se muestra a continuación.
3𝐸𝐼
3𝐸𝐼 ⋅ 𝑈1
⋅ 𝑈1 𝐿
𝐿 e
a c
3𝐸𝐼
⋅ 𝑈1 2
⋅ 𝑈1 ⋅ 𝑈1 𝐿2
𝐿2 𝐿 𝐿2


3𝐸𝐼 3𝐸𝐼
⋅ −0.03 𝑚 ⋅ −0.03 𝑚
𝐿2 𝐿2
3𝐸𝐼 3𝐸𝐼
3𝐸𝐼 ⋅ −0.03 𝑚 3𝐸𝐼 ⋅ −0.03 𝑚
⋅ −0.03 𝑚 𝐿3 ⋅ −0.03 𝑚 𝐿3
𝐿3 𝐿3
3𝐸𝐼 3𝐸𝐼
𝑉𝑎𝑐 = 2
⋅ 𝑈1 + 3 − 0.03 𝑚
𝐿 𝐿
3 ⋅ 28000 𝑀𝑝𝑎 ⋅ 4.02 𝑚4 3 ⋅ 28000 𝑀𝑝𝑎 ⋅ 4.02 𝑚4

𝑉𝑎𝑐 = ⋅ −0.000168 𝑟𝑎𝑑 + ⋅ −0.03 𝑚
342 𝑚 343 𝑚
𝑉𝑎𝑐 = 208.67 𝑘𝑁
3𝐸𝐼 3𝐸𝐼
𝑉𝑐𝑎 = − ⋅ 𝑈1 − − 0.03 𝑚
𝐿2 𝐿3
3 ⋅ 28000 𝑀𝑝𝑎 ⋅ 4.02 𝑚4 3 ⋅ 28000 𝑀𝑝𝑎 ⋅ 4.02 𝑚4

𝑉𝑐𝑎 = − ⋅ −0.000168 𝑟𝑎𝑑 − ⋅ −0.03 𝑚
342 𝑚 343 𝑚
𝑉𝑎𝑐 = −208.67 𝑘𝑁
3𝐸𝐼 3𝐸𝐼
𝑉𝑐𝑒 = 2
⋅ 𝑈1 + 3 − 0.03 𝑚
𝐿 𝐿
3 ⋅ 28000 𝑀𝑝𝑎 ⋅ 4.02 𝑚4 3 ⋅ 28000 𝑀𝑝𝑎 ⋅ 4.02 𝑚4

𝑉𝑐𝑒 = ⋅ −0.000168 𝑟𝑎𝑑 + ⋅ −0.03 𝑚
422 𝑚 423 𝑚
𝑉𝑐𝑒 = −168.895 𝑘𝑁


3𝐸𝐼 3𝐸𝐼
𝑉𝑒𝑐 = − 2
⋅ 𝑈1 − 3 − 0.03 𝑚
𝐿 𝐿
3 ⋅ 28000 𝑀𝑝𝑎 ⋅ 4.02 𝑚4 3 ⋅ 28000 𝑀𝑝𝑎 ⋅ 4.02 𝑚4

𝑉𝑒𝑐 = − ⋅ −0.000168 𝑟𝑎𝑑 − ⋅ −0.03 𝑚
422 𝑚 423 𝑚
𝑉𝑐𝑐 = 168.895 𝑘𝑁
3𝐸𝐼 3𝐸𝐼
𝑀𝑐𝑎 = ⋅ 𝑈1 − 2 ⋅ −0.03 𝑚
𝐿 𝐿
3 ⋅ 2800 𝑀𝑝𝑎 ⋅ 4.02 𝑚4 3 ⋅ 28000 𝑀𝑝𝑎 ⋅ 4.02 𝑚4

𝑀𝑐𝑎 = ⋅ −0.000168 𝑟𝑎𝑑 − ⋅ −0.03 𝑚
34 𝑚 342
𝑀𝑐𝑎 = 7094.12 𝑘𝑁 ⋅ 𝑚
3𝐸𝐼 3𝐸𝐼
𝑀𝑒𝑐 = ⋅ 𝑈1 + 2 ⋅ −0.03 𝑚
𝐿 𝐿
3 ⋅ 2800 𝑀𝑝𝑎 ⋅ 4.02 𝑚4 3 ⋅ 28000 𝑀𝑝𝑎 ⋅ 4.02 𝑚4

𝑀𝑒𝑐 = ⋅ −0.000168 𝑟𝑎𝑑 + ⋅ −0.03 𝑚
42 𝑚 422
𝑀𝑒𝑐 = −7094.12 𝑘𝑁 ⋅ 𝑚
7094.12 𝑘𝑁 ⋅ 𝑚 7094.12 𝑘𝑁 ⋅ 𝑚
208.69 𝑘𝑁 208.69 𝑘𝑁 168.89 𝑘𝑁 168.89 𝑘𝑁
Ilustración 14. Fuerzas internas en los elementos de la viga por el efecto del asentamiento


208.69 𝑘𝑁
168.69 𝑘𝑁
Ilustración 15. Diagrama de fuerza cortante por el efecto del asentamiento
7094.12 𝑘𝑁
+ -
Ilustración 16. Diagrama de momentos por efecto del asentamiento
𝑀𝑎𝑐 = 0
3𝐸𝐼 3 ∙ 28000𝑀𝑃𝑎 ∙ 4.02

𝑀𝑐𝑎 = 2
∙ −0.03 = ∙ −0.03
𝐿 342
𝑀𝑐𝑎 = 292.11 ∙ −0.03 = −8763.32 𝑘𝑁 − 𝑚
3𝐸𝐼 3 ∙ 28000𝑀𝑃𝑎 ∙ 4.02

𝑀𝑐𝑒 = − ∙ −0.03 = − ∙ −0.03
𝐿2 422
𝑀𝑐𝑒 = −191429 ∙ −0.03 = 5742.86 𝑘𝑁 − 𝑚
𝑀𝑒𝑐 = 0


c) Calcule y dibuje el diagrama de momentos de la viga debido a la carga de
presfuerzo.
Diagrama viga presforzada
Para el caso del preesfuerzo, es necesario calcular el vector de fuerzas generado por este.
Como se tiene un único desplazamiento de giro, sólo interesan los momentos en los extremos
de los elementos.
Para el elemento ac se tiene lo siguiente.
𝑐=0
𝑏 = 0.15417
𝑎 = 0.005119
−20000
473.088
0
{𝑃𝑜 }𝑎𝑐 =
20000
−473.088
{ 16085 }
Para el elemento ce se tiene lo siguiente.
𝑐 = 0.675
𝑏 = −0.1569
𝑎 = 0.003354


−20000
−382.976
−168085
{𝑃𝑜 }𝑐𝑒 =
20000
382.976
{ 0 }
{𝐹𝑜 }𝑝𝑟𝑒𝑒𝑠𝑓𝑢𝑒𝑟𝑧𝑜 = 0 𝑘𝑁 ⋅ 𝑚
{𝐹𝑙 } = {𝐹𝑛 } − {𝐹𝑛 }

{𝐹𝑙 } = 0
Con este vector de fuerzas y la matriz de rigidez calculada previamente, es posible calcular el
desplazamiento por preesfuerzo como se muestra a continuación.
{𝑈𝐿 } = [𝐾𝐿𝐿 ]−1 {𝐹𝑙 }
{𝑈𝑙 } = 0 𝑟𝑎𝑑
Ahora bien, para el caso de preesfuerzo, las fuerzas internas de los elementos se calcula de la
siguiente manera.
{𝑃} = [𝐾𝑒 ] ⋅ {𝑈𝑒 } + {𝑃𝑜 } + {𝑃𝑝 }
3𝐸𝐼
3𝐸𝐼 ⋅ 𝑈1
⋅ 𝑈1 𝐿
𝐿 e
a c
3𝐸𝐼
⋅ 𝑈1 2
⋅ 𝑈1 ⋅ 𝑈1 𝐿2
𝐿2 𝐿 𝐿 2


16085 𝑘𝑁 ⋅ 𝑚 −16085 𝑘𝑁 ⋅ 𝑚
c e
a c e
473.088 𝑘𝑁 −473.008 𝑘𝑁 −382.976 𝑘𝑁 382.976 𝑘𝑁
Ilustración 18. Efecto del preesfuerzo
El vector de fuerzas {𝑃𝑝 } para el elemento ac se obtiene de la siguiente manera.
𝜃𝑖 = atan(−𝑒 ′ (0))
𝑒 ′ = 2 ⋅ 𝑎̃𝑥 2 + 𝑏̃
𝜃𝑖 = atan(−𝑏̃) = atan(−0.15417) = −8.76429°
𝑓𝑥𝑖 = 𝑓𝑜 ⋅ cos 𝜃𝑖 = 19766.5 𝑘𝑁
𝑓𝑦𝑖 = 𝑓𝑜 ⋅ sin 𝜃𝑖 = −3047.4 𝑘𝑁
𝜃𝑗 = atan(−𝑒 ′ (𝐿))
𝜃𝑖 = atan(𝑏̃)
𝑓𝑥𝑖 = 𝑓𝑜 ⋅ cos 𝜃𝑖 = −19634.4𝑘𝑁
𝑓𝑦𝑖 = 𝑓𝑜 ⋅ sin 𝜃𝑖 = 13253.2 𝑘𝑁


𝑓𝑥𝑖
19766.5 𝑘𝑁
−𝑓𝑦𝑖 3047.4 𝑘𝑁
𝑓𝑥𝑖 ⋅ 𝑒1 0
{𝑃𝑃 }𝑎𝑐 = −𝑓𝑥𝑗 =
−19634.4 𝑘𝑁
−𝑓𝑦𝑗 3806.76 𝑘𝑁
{ 13253.2 𝑘𝑁 ⋅ 𝑚}
{ 𝑓𝑥𝑗 ⋅ 𝑒3 }
El vector de fuerzas {𝑃𝑝 } para el elemento ce se obtiene de la misma manera, obteniendo lo

siguiente.
𝑓𝑥𝑖
19758.1 𝑘𝑁
−𝑓𝑦𝑖 −3101.08 𝑘𝑁
𝑓𝑥𝑖 ⋅ 𝑒1 −13336.7 𝑘𝑁 ⋅ 𝑚
{𝑃𝑃 }𝑎𝑐 = −𝑓𝑥𝑗 =
−19846 𝑘𝑁
−𝑓𝑦𝑗 2476.97 𝑘𝑁
{ 0 }
{ 𝑓𝑥𝑗 ⋅ 𝑒3 }
3𝐸𝐼
𝑀𝑐𝑎 = 16085 𝑘𝑁 ⋅ 𝑚 + ⋅ 𝑈1 + 13253.2 𝑘𝑁 ⋅ 𝑚
𝐿
3 ⋅ 28000 𝑀𝑝𝑎 ⋅ 4.02 𝑚4

𝑀𝑐𝑎 = 16085 𝑘𝑁 ⋅ 𝑚 + ⋅ 0 + 13253.2 𝑘𝑁 ⋅ 𝑚
34 𝑚
𝑀𝑐𝑎 = 29338.2 𝑘𝑁
𝑀𝑒𝑐 = 0 𝑘𝑁 ⋅ 𝑚
3𝐸𝐼
𝑀𝑐𝑒 = −16085 𝑘𝑁 ⋅ 𝑚 + ⋅ 𝑈1 + 19758.1 𝑘𝑁
𝐿
3 ⋅ 28000 𝑀𝑝𝑎 ⋅ 4.02 𝑚4

𝑀𝑐𝑒 = −16085 𝑘𝑁 ⋅ 𝑚 + ⋅ 0 − 13336.9 𝑘𝑁 ⋅ 𝑚
42 𝑚
𝑀𝑐𝑒 = −29421.7𝑘𝑁 ⋅ 𝑚


29338.2 𝑘𝑁 ⋅ 𝑚 29421.7 𝑘𝑁 ⋅ 𝑚
Ilustración 19. Fuerzas internas por el efecto del preesfuerzo
16636.66 𝑘𝑁 ⋅ 𝑚
16630.66 𝑘𝑁 ⋅ 𝑚
+ +
+
+
a
- -
29338.2 𝑘𝑁 ⋅ 𝑚 29421.7 𝑘𝑁 ⋅ 𝑚
Ilustración 20. Diagrama de momentos en la viga por efectos del preesfuerzo
f) Sume los diagramas de momento obtenidos en a), b) y c) y concluya sobre el efecto del
presfuerzo.
Al sumar los tres diagramas de momentos obtenidos anteriormente, se obtiene lo siguiente.
𝑀𝑐𝑎 = −35435 𝑘𝑁 ⋅ 𝑚 + 7094.12 𝑘𝑁 ⋅ 𝑚 + 29338.2 𝑘𝑁 ⋅ 𝑚

𝑀𝑐𝑎 = 997.32 𝑘𝑁 ⋅ 𝑚
𝑀𝑒𝑐 = 35435 𝑘𝑁 ⋅ 𝑚 − 7094.12 𝑘𝑁 ⋅ 𝑚 − 29421.7 𝑘𝑁 ⋅ 𝑚

𝑀𝑒𝑐 = −1080.82 𝑘𝑁 ⋅ 𝑚
𝑀𝑐𝑒 = 0 𝑘𝑁 ⋅ 𝑚


997.32 𝑘𝑁 ⋅ 𝑚 1080.82 𝑘𝑁 ⋅ 𝑚
Ilustración 21. Momentos en la viga al sumar el efecto de las tres cargas
1080.82 𝑘𝑁 ⋅ 𝑚 12940.69 𝑘𝑁 ⋅ 𝑚
997.32 𝑘𝑁 ⋅ 𝑚
+ +
+
a -
1554.35 𝑘𝑁 ⋅ 𝑚
Ilustración 22. Diagrama de momento en la viga por la combinación de las tres cargas
Con base en lo anterior, se puede evidenciar que los momentos generados por el preesfuerzo
contrarrestan en gran medida los momentos generados por la carga distribuida y el asentamiento.
Ya que la suma de los momentos generados por estas dos cargas en el nodo c, es casi de la
misma magnitud que el momento generado en este mismo nodo por efecto del preesfuerzo.
g) Verifique los resultados con SAP2000 o con Opensees.
Para el caso del efecto de las cargas distribuidas, se obtiene el siguiente diagrama de momentos, el
cual da muy similar al diagrama obtenido con el método manual.


Ilustración 23. Diagrama de momento para la viga en SAP2000 por efecto de la carga
distribuida
Ilustración 24. Diagrama de momento para el elemento ac en SAP2000 por efecto de la carga
distribuida
Ilustración 25. Diagrama de momento para el elemento cd en SAP2000 por efecto de la carga
distribuida
Para el caso del efecto del asentamiento, se tiene obtiene el siguiente diagrama de momentos
en SAP2000.
Ilustración 26. Diagrama de momento para la viga en SAP2000 por efecto del asentamiento


Ilustración 27. Diagrama de momento para el elemento ac en SAP2000
Ilustración 28. Diagrama de momento para el elemento cd en SAP2000
Para el caso del preesfuerzo se debe calcular la distancia a la cual se encuentra la excentricidad
de 1.161 de la siguiente manera.
𝑏 0.154176
𝑥𝑎𝑐 = − = = 15.0607 𝑚
2𝑎 2 ⋅ 0.005119
𝑏 0.156952
𝑥𝑎𝑐 = − = = 23.3956 𝑚
2𝑎 2 ⋅ 0.003354
Y se obtiene el siguiente diagrama de momentos.
Ilustración 29. Diagrama de momento en SAP2000 por acción del preesfuerzo
Ilustración 30. Diagrama de momento en el elemento ac en SAP2000


Ilustración 31. Diagrama de momento para el elemento ce en SAP2000
Finalmente, para el caso de la combinación de las tres cargas, se obtiene lo siguiente.
Ilustración 32. Diagrama de momento en la viga con SAP2000
Ilustración 33. Diagrama de momento del elemento ac con SAP2000
Ilustración 34. Diagrama de momento del elemento ce con SAP2000
Para cada uno de los casos, se realizó una tabla con la comparación entre los resultados obtenidos
con el método manual y los resultados obtenidos con SAP2000, las cuales se presentan a
continuación.
Tabla 1. Comparación entre el método manual y SAP2000 para el efecto de la carga

distribuida
Cargas distribuidas

Momento [kN*m] Manual SAP Error [%]
Mca -35435 -35435 0%
Mce 35435 35435 0%
Tabla 2. Comparación entre el método manual y SAP2000 para el efecto del asentamiento
Asentamiento
Mca 7094.12 7094.12 0%
Mce -7094.12 -7094.12 0%
Tabla 3. Comparación entre el método manual y SAP2000 para el efecto del preesfuerzo
Preesfuerzo
Mca 29338.2 29330.18 0.027%
Mce -29421.7 -29413.71 -0.027%
Tabla 4. Comparación entre el método manual y SAP2000 para el efecto de todas las cargas
juntas
Combinación
Mca 997.32 989.3 0.80%
Mce -1080.82 -1072.83 -0.74%
Con base en las Tabla 1, Tabla 2, Tabla 3 y Tabla 4, es posible evidenciar que los resultados
obtenidos con el método manual y con SAP2000 son muy similares, por lo que se puede concluir
que estos resultados son coherentes y correctos.
Problema 2 (implementación computacional)

Usando las funciones desarrolladas en la Tarea 1, codifique en Matlab (o Python) las siguientes
rutinas:
a) Función Po: calcula las fuerzas 𝑃 para un elemento a partir de los datos de elemento.
Adicional a la información de entrada de la función Po de la tarea 1, se deben incluir datos
de cambios de temperatura y presfuerzo axial y a flexión. Note que para obtener 𝑃 es
necesario calcular 𝑝 y transformarlo a coordenadas globales usando la matriz de
transformación.


b) Rutina AnálisisEstaticoLineal.m: calcula fuerzas desplazamientos, fuerzas internas y
reacciones de una estructura plana cualquiera. La entrada de datos puede ser mediante un
archivo de texto, Excel o similar que contiene la siguiente información:
▪ Matriz DatosElementos: cada fila representa un elemento y las columnas contienen los
siguientes datos: 𝐴, 𝐼, 𝐿, 𝐸, 𝜃, y Tipo. Dimensión: número de elementos × 6.
▪ Matriz GDL: cada fila representa un elemento y las columnas contienen la numeración
de desplazamientos del nodo inicial y del nodo final (seis columnas en total). Dimensión:
número de elementos × 6.
▪ Vector Ur: cada fila contiene el valor del desplazamiento conocido. Dimensión: número
de desplazamientos conocidos × 1.
▪ Matriz Resortes: cada fila contiene el número del grado de libertad y la constante del
resorte. Dimensión: número de resortes × 2.
▪ Vector Fn: cada fila contiene el valor de la fuerza en el grado de libertad. Dimensión:
número de desplazamientos desconocidos × 1.
▪ Matriz FuerzasElementos: cada fila representa un elemento y las columnas contienen
datos de todos los tipos de carga, por ejemplo: 𝑤 , 𝑤 , 𝑓, 𝜀 , 𝜙 , 𝑞 , 𝑓 , 𝑒 , 𝑒 , 𝑒 , etc.
Dimensión: número de elementos × 10. Al final del archivo de entrada de datos se puede
llamar la rutina AnálisisEstaticoLineal.m.
Tanto la función para calcular el vector Po como el programa de análisis estático lineal, se
mandan como anexo.

La estructura de la Figura P3 tiene elementos viga‐columna de acero estructural ASTM A992
(𝐸 =200 GPa, 𝛼 =1.17E5/oC). La temperatura al interior de la estructura es de ‐15°C mientras
que la temperatura en su exterior es de 45°C.
Suponga que la temperatura inicial es 25°C. Considere que todos los elementos se deforman
axialmente.
a) Dibuje la deformada y los diagramas de axial, cortante y momento usando análisis
estático lineal (implementación computacional) y el programa de computador
SAP2000. Compare los resultados obtenidos con las dos metodologías. Para la
implementación computacional utilice el programa desarrollado en el Problema 2.
Para usar el programa desarrollado en el punto 2 de la tarea, en primer lugar, es necesario

enumerar los grados de libertad como se muestra a continuación.


5 8
6 4 9 7
3
2 2
3 1
4
14
1
15 13
11
12 10
Ilustración 35. Etiqueta de los grados de liberta para la implementación computacional
Asimismo, es necesario calcular 𝜖𝑜 y 𝜙𝑜 para tener en cuenta los efectos por cambio de
temperatura.
𝑇𝑏 + 𝑇𝑡 −15 + 45
𝛥𝑇𝑎 = − 𝑇𝑜 = − 25 = −10°𝐶
2 2
𝜖0 = 𝛼 𝑇 ⋅ 𝛥𝑇𝑎 = 1.17 × 10−5 /°𝐶 ⋅ −10°𝐶 = −0.000117
𝛥𝑇𝑓 = 𝑇𝑏 − 𝑇𝑡 = −15 − 45 = −60 °𝐶
𝛥𝑇𝑓
∇Tf =
ℎ
El gradiente de temperatura varía para las vigas y las columnas debido a que tienen altura
diferente.
−60
∇Tfvigas = = −133.333
0.45 𝑚
𝜙0𝑣𝑖𝑔𝑎𝑠 = 𝛼 𝑇 ⋅ ∇Tfvigas = −0.00156
−60
∇Tfcolumnas = = −181.818
0.33 𝑚
𝜙0𝑐𝑜𝑙𝑢𝑚𝑛𝑎𝑠 = 𝛼𝑇 ⋅ ∇Tfcolumnas = −0.002127


Ahora bien, los datos ingresados como archivos de entrada son los siguientes.
Ilustración 36. Datos de cada elemento
Ilustración 37. Fuerzas en cada elemento
Ilustración 38. Etiqueta grados de libertad de cada elemento
Ilustración 39. Resortes


Ilustración 40. Vector Fn
Ilustración 41. Vector Ur
Una vez se corre el programa, se obtienen los siguientes valores de desplazamiento.


Ilustración 42. Vector Ul obtenido en Matlab

−0.0048 𝑚
−0.000449 𝑚
0.0024 𝑟𝑎𝑑
−0.0107 𝑚
{𝑈𝐿 } = 0.0084 𝑚
0.000705 𝑟𝑎𝑑
−0.0113 𝑚
−0.000486 𝑚
{ −0.0014 𝑟𝑎𝑑 }
Para cada uno de los elementos se obtienen las siguientes fuerzas internas.
Ilustración 43. Vector p elemento ab


Ilustración 44. Vector p elemento bc
Ilustración 45. Vector p elemento cd
Ilustración 46. Vector p elemento ed
Con base en los desplazamientos obtenidos, se obtiene la siguiente deformada.


c d
Ilustración 47. Deformada estructura implementación computacional
Asimismo, con base en las fuerzas internas obtenidas, se obtienen los siguientes diagramas de
fuerza axial cortante y momento.
c d
-
b 3.58 𝑘𝑁 +
7.166 𝑘𝑁
- e
7.166 𝑘𝑁
Ilustración 48. Diagrama de fuerza axial implementación computacional


c d
- 7.166 𝑘𝑁
b 6.209 𝑘𝑁
Ilustración 49. Diagrama de cortante implementación computacional
c d
- 35.83 𝑘𝑁 ⋅ 𝑚
- 35.83 𝑘𝑁 ⋅ 𝑚
b
- 69.96 𝑘𝑁 𝑚 e
a 69.96 𝑘𝑁 𝑚
Ilustración 50. Diagramas de momento obtenidos con la implementación computacional
Para modelar la estructura en SAP se ingresan los gradientes de temperatura calculados

previamente. Es importante tener en cuenta que los gradientes de temperatura a flexión se deben
ingresa con signo positivo, mientras que el gradiente axial si se ingresa con el signo calculado.
Una vez corrido el modelo de SAP, se obtiene la siguiente deformada.


Ilustración 51. Deformada SAP2000
Asimismo, se obtienen los siguientes desplazamientos en los nodos.


Ilustración 52. Desplazamiento en el nodo b obtenidos con SAP2000


Ilustración 53. Desplazamiento en el nodo c obtenidos con SAP2000
Ilustración 54. Desplazamientos en el nodo d obtenidos con SAP2000

Asimismo, se obtienen los siguientes diagramas de fuerzas internas con SAP2000.


Ilustración 55. Diagrama de carga axial obtenido con SAP2000
Ilustración 56. Diagrama de cortante obtenido con SAP2000


Ilustración 57. Diagrama de momento con SAP2000
Al comparar los valores obtenidos con ambos métodos se obtienen los siguientes porcentajes de
error.
Tabla 5. Comparación de los resultados obtenidos con el método computacional y SAP2000
GDL Computacional SAP2000 Error [%]

1 [m] -0.004847219 -0.00485 -0.057%
2 [m] -0.000449017 -0.000449 -0.004%
3 [rad] 0.002423609 0.00242 0.149%
4 [m] -0.010708345 -0.01071 -0.015%
5 [m] 0.008445363 0.00844 0.063%
6 [rad] 0.000705126 0.0007052 0.010%
7 [m] -0.011293345 -0.01129 -0.030%
8 [m] -0.000486983 -0.000487 -0.004%
9 [rad] -0.001430664 -0.00143 -0.046%
Con base en la Tabla 5 se puede evidenciar que los resultados obtenidos por ambos métodos son
muy cercanos, llegando a presentar errores de menor del 1%. Teniendo en cuenta lo anterior, es
posible concluir que los resultados obtenidos con el programa realizado en Matlab son
coherentes.


b) ¿Por qué la deformada del elemento “bc” muestra curvaturas negativas mientras que
el momento genera tensiones en la cara interna? Justifique su respuesta.
A pesar de que la deformada muestre curvaturas negativas, la temperatura es una acción no

mecánica. La cual genera tensión donde haya un aumento de la temperatura y compresión donde
haya una disminución de la temperatura. Asimismo, este efecto de la temperatura, se analiza
como si el elemento viniera con defectos de fabricación, es decir que ya viene deformado, y el
trata de volver a su posición original, haciendo que se genera tensión en la cara interna.
d
x c x
=30o
b x
Elementos ab y ed
A= 7.55E3 mm2
I= 1.28E8 mm4
x e
h= 3.30E2 mm
ab=4.00m Elementos bc y cd
cd=5.00m A= 6.65E3 mm2
a bc=5.50m Y
I= 2.12E8 mm4
ed=4.00m h= 4.50E2 mm
X
Figura P3. Estructura a analizar

Usted ha sido comisionado para revisar el puente peatonal atirantado de la Figura P4 utilizando
como base el código desarrollado en el problema 2. El puente está sometido a una carga vertical
distribuida de 30 kN/m sobre los elementos abcd y un asentamiento de la torre (nodos d, e y f)
de 1 cm hacia abajo.
a) Calcule los desplazamientos verticales de los nodos b y c por la acción de la carga
vertical y el asentamiento.
En primer lugar, es necesario realizar la enumeración de los grados de libertad para poder
ingresar la información al programa.


17
18 16
14
13
15
4
5
8 3 6 11
1 3 10
9 1 4 2 2 7 5 12
Ilustración 58. Enumeración de los grados de libertad
Los datos de entrada son los siguientes.
Ilustración 59. Input de grados de libertad de los elementos
Ilustración 60. Input datos elemento


Ilustración 61. Input datos resortes
Ilustración 62. Input vector Ur
Ilustración 63. Input vector Fn


Ilustración 64. Input fuerzas en los elementos
Al correr el programa, se obtienen los siguientes desplazamientos, donde la tercera fila es el

desplazamiento vertical del nodo b y la sexta fila es el desplazamiento vertical del nodo c. e
Ilustración 65. Desplazamientos por acción de la carga distribuida y el asentamiento usando

el método computacional
Al correr el programa SAP2000 se obtienen los siguientes desplazamientos.
Ilustración 66. Desplazamiento SAP2000 por efecto de la carga distribuida y el asentamiento

en el nodo b


Ilustración 67. Desplazamiento obtenido en SAP2000 por efecto de la carga distribuida y el

asentamiento en el nodo c
Al correr el modelo en OpenSees se obtienen los siguientes resultados.
Ilustración 68. Desplazamientos verticales en el nodo b y c con OpenSees
b) Someta el cable “ce” a un presfuerzo axial de 1 kN y calcule los desplazamientos.
Para el caso del preesfuerzo axial, se ingresan los mismos inputs del punto anterior teniendo en
cuenta que ya no hay desplazamientos y que en el input de fuerzas en los elementos sólo se tiene
el preesfuerzo axial del elemento ce.


Ilustración 69. Input de fuerzas en los elementos
Al correr el programa, se obtienen los siguientes desplazamientos, donde la tercera fila es el

desplazamiento vertical del nodo b y la sexta fila es el desplazamiento vertical del nodo c.
Ilustración 70. Desplazamientos por el preesfuerzo axial en el cable ce

Al correr el modelo de SAP2000 se obtienen los siguientes resultados
Ilustración 71. Desplazamiento en el nodo b en SAP2000 por el presfuerzo axial en el cable ce


Ilustración 72. Desplazamiento en el nodo c en SAP2000 por el preesfuerzo axial en el cable

ce
En OpenSees, se obtienen los siguientes resultados.
Ilustración 73. Desplazamiento vertical en el nodo b y c con OpenSees
c) Someta el cable “bf” a un presfuerzo axial de 1 kN y calcule los desplazamientos.
En este caso se realiza el mismo procedimiento que en el punto anterior, a diferencia de que la
carga de preesfuerzo axial va en el elemento bf. Al correr el programa se obtienen los siguientes
desplazamientos, donde la tercera fila es el desplazamiento vertical del nodo b y la sexta fila es
el desplazamiento vertical del nodo c.


Ilustración 74. Desplazamientos por el preesfuerzo axial del cable bf
Al correr el modelo en SAP2000 se obtienen los siguientes resultados.
Ilustración 75. Desplazamiento en el nodo b con SAP2000 para el preesfuerzo axial en el

cable bf


Ilustración 76. Desplazamiento en el nodo c en SAP2000 por el preesfuerzo axial en el

elemento ce
Al correr el código en OpenSes se obtienen los siguientes resultados.
Ilustración 77. Desplazamiento en el nodo b y c con OpenSees
d) Usando los resultados de a), b) y c), calcule la fuerza de presfuerzo necesaria para
disminuir los desplazamientos verticales de los puntos “b” y “c” a cero. Desprecie el
efecto de la catenaria en el cable. Plantee un sistema de ecuaciones similar al utilizado
en el documento “Puentes peatonales atirantados ‐ Garcia y Perez (1996)” disponible
en la carpeta “Material” en Bloque Neón. Reporte la fuerza de presfuerzo que debe
aplicar el gato durante el proceso constructivo.
Teniendo en cuenta que en los numerales anteriores se aplicó un preesfuerzo de 1 kN en los

cables, ahora se debe encontrar el valor del preesfuerzo escalando el desplazamiento obtenidos

cuando se realizó este procedimiento, para disminuir a cero los desplazamientos obtenidos por
el efecto de las cargas y el asentamiento. Con base en lo anterior, se obtienen las dos siguientes
ecuaciones, una para cada nodo.
0.2878 𝑚 = 0.0258 × 10−3 𝑚 ⋅ 𝑃𝑐𝑎𝑏𝑙𝑒 𝑐𝑒 + 0.242149 × 10−3 𝑚 ⋅ 𝑃𝑐𝑎𝑏𝑙𝑒 𝑏𝑓
0.0979 𝑚 = 0.11286 × 10−3 ⋅ 𝑃𝑐𝑎𝑏𝑙𝑒 𝑐𝑒 + 0.017219 × 10−3 𝑚 ⋅ 𝑃𝑐𝑎𝑏𝑙𝑒 𝑏𝑓
Al resolver el sistema de ecuaciones, se obtienen los siguientes valores de preesfuerzo para cada
uno de los cables.
𝑃𝑐𝑎𝑏𝑙𝑒 𝑐𝑒 = 697.451 𝑘𝑁
𝑃𝑐𝑎𝑏𝑙𝑒 𝑏𝑓 = 1114.21 𝑘𝑁
Al aplicar estos valores de preesfuerzo y correr el programa computacional, se obtienen los

siguientes resultados. Como e puede evidenciar, el desplazamiento en los nodos b y c son muy
cercanos a cero.
Ilustración 78. Desplazamientos en el nodo b y c al aplicar un preesfuerzo en los cables ce y

bf
Asimismo al aplicar estos valores de preesfuerzo en SAP2000, se puede evidenciar que para los
nodos b y c se obtiene un desplazamiento muy cercano a 0.


Ilustración 79. Desplazamiento en el nodo b
Ilustración 80. Desplazamiento en el nodo c
Finalmente, al correr el modelo en OpenSees también se obtienen valores muy cercanos a cero.


Ilustración 81. Desplazamiento en los nodos b y c en OpenSees por acción del preesfuerzo de
los cables bf y ce
Sin embargo, se debe tener en cuenta que esa no es la fuerza que se le debe aplicar al gato. Ya
que el cable ya ha tenido una deformación inicial debido a la acción de la carga distribuida, el
asentamiento y el preesfuerzo. Con base en lo anterior, se tiene lo siguiente.
−𝑘4 −𝑘4
−𝑘1 −𝑘4


𝑐𝑜𝑠
𝑠𝑒𝑛𝑜
1 𝑘𝑁
Ilustración 83. Efecto del preesfuerzo en el cable bf
𝑐𝑜𝑠
𝑠𝑒𝑛𝑜
1 𝑘𝑁
Ilustración 84. Efecto del preesfuerzo en el cable ce

Con base en lo anterior, la carga vertical y horizontal en el extremo de los cables bf y ce es el
siguiente.


𝐹𝑦 𝑏𝑓 = −𝑘4 ⋅ −0.01 𝑚 − cos 28.0725
𝐸𝐴
𝐹𝑦 𝑏𝑓 = − cos 28.0725 ⋅ sin 28.0725 ⋅ −0.01 𝑚 − cos 28.0725
𝐿
200000𝑀𝑝𝑎 ⋅ 0.001184 𝑚2
𝐹𝑦 𝑏𝑓 = − cos 28.0725 ⋅ sin 28.0725 ⋅ −0.01 𝑚 − cos 28.0725
√302 𝑚 + 162 𝑚
𝐹𝑦 𝑏𝑓 = 28.0368 𝑘𝑁
𝐹𝑥 𝑏𝑓 = −𝑘1 ⋅ −0.01 𝑚 − sin 28.0725
𝐸𝐴
𝐹𝑥 𝑏𝑓 = − ⋅ (cos 28.0725)^2 ⋅ −0.01 𝑚 − sin 28.0725
𝐿
200000𝑀𝑝𝑎 ⋅ 0.001184 𝑚2
𝐹𝑥 𝑏𝑓 = − ⋅ (cos 28.0725)^2 ⋅ −0.01 𝑚 − sin 28.0725
√302 𝑚 + 162 𝑚
𝐹𝑥 𝑏𝑓 = 53.753 𝑘𝑁
𝐹𝑦 𝑐𝑒 = −𝑘4 ⋅ −0.01 − cos 45
200000𝑀𝑝𝑎 ⋅ 0.001184 𝑚2
𝐹𝑦 𝑐𝑒 = − ⋅ cos 45 ⋅ sin 45 ⋅ −0.01 𝑚 − cos 45
15 ⋅ √2
𝐹𝑦 𝑐𝑒 = 55.107 𝑘𝑁
𝐹𝑥 𝑐𝑒 = −𝑘1 ⋅ −0.01 − sin 45
−𝐸𝐴
𝐹𝑥 𝑐𝑒 = ⋅ cos 452 ⋅ −0.01 𝑚 − sin 45
𝐿
𝐹𝑥 𝑐𝑒 = 55.1072 𝑘𝑁
Con estos valores de fuerzas verticales y horizontales en los cables, es posible calcular la carga
axial en estos por medio de Pitágoras.


𝑞01 𝑏𝑓 = √53.7532 + 28.03682 = 60.6263 𝑘𝑁
𝑞01 𝑐𝑒 = √55.1072 + 55.10722 = 77.933 𝑘𝑁
La fuerzas final en el gato es entonces, la fuerza de preesfuerxo calculada previamente menor

el valor de q01.
𝐺𝑎𝑡𝑜𝑏𝑓 = 1114.21 𝑘𝑁 − 60.6253 𝑘𝑁 = 1053.58 𝑘𝑁

𝐺𝑎𝑡𝑜𝑐𝑒 = 697.451 𝑘𝑁 − 77.933 𝑘𝑁 = 601.518 𝑘𝑁
e) Verifique sus respuestas con SAP2000 y Opensees. ¿Qué pasaría si alguien corta el
cable bf? Justifique su respuesta con SAP2000.
Al cortar el cable bf y dejar sólo el cable ce con el postensado calculado previamente, es posible
evidenciar que el desplazamiento del nodo b aumenta a 1.26 m. Lo cual quiere decir que el
puente colapsaría si esto ocurre, ya que todo el tablero tendría fallos estructurales graves.
Ilustración 85. Desplazamiento en el nodo b al cortante el cable bf


Elementos abcd:
I=72E‐4 m4; E=20000 MPa
Elementos bf y ce:
A=11.84E‐4 m2; E=200000 MPa
f 1m
15 m
b c d
15 m 15 m 15 m
Figura P4. Puente peatonal

Tarea 2 - Rodríguez - Gutiérrez

Cargado por

Copyright:

Formatos disponibles

Tarea 2 - Rodríguez - Gutiérrez

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Tarea 2 - Rodríguez - Gutiérrez

Cargado por

Copyright:

Formatos disponibles

Integrantes: Sofía Rodríguez y Natalia Gutiérrez

Departamento de Ingeniería Civil y Ambiental

Objetivo : Analizar estructuras usando análisis estático lineal. Codificar el método

Problema 1 (implementación manual)

ICYA 4422 – Primer semestre de 2023 – Pag. 1

Departamento de Ingeniería Civil y Ambiental

Figura P1. Vista isométrica, longitudinal y sección transversal del puente

A continuación, se realizan los pasos para obtener el diagrama de momentos en la viga

Modelo estructural simplificado

ICYA 4422 – Primer semestre de 2023 – Pag. 2

Departamento de Ingeniería Civil y Ambiental

Ilustración 1. Grados de libertad de la estructura

Luego de realizar dichas transformaciones se obtiene un grado de libertad en el nodo central

ICYA 4422 – Primer semestre de 2023 – Pag. 3

Departamento de Ingeniería Civil y Ambiental

Ilustración 3. Fuerzas al hacer U1=1

𝐾11 = 9931764.71 + 8040000

{𝐹𝑙 } = {𝐹𝑛 } − {𝐹0 }

Ilustración 4. Momento en el extremo del elemento ac por acción de la carga distribuida

Departamento de Ingeniería Civil y Ambiental

𝑤𝐿2 190 ∙ 342

Ilustración 5. Momento en el extremo del elemento ce por acción de la carga distribuida

𝑤𝐿2 190 ∙ 422

Vector de fuerzas de la estructura

{𝐹0 } = {𝑃0 𝑎𝑐 + 𝑃0 𝑐𝑒 } = {−27455 + 41895} = {14440 𝑘𝑁 − 𝑚}

{𝐹𝑙 } = {𝐹𝑛 } − {𝐹0 }

{𝐹𝑙 } = {0} − {14440} = {−14440 𝑘𝑁 − 𝑚}

{𝑈𝑙 } = [𝐾𝑙𝑙 ]−1 {𝐹𝑙 }

{𝑈𝑙 } = −0.000803 𝑟𝑎𝑑

ICYA 4422 – Primer semestre de 2023 – Pag. 5

Departamento de Ingeniería Civil y Ambiental

Ilustración 6. Efecto de U1=1

190 kN/m 190 kN/m

3𝑤𝐿 5𝑤𝐿 5𝑤𝐿 3𝑤𝐿

Ilustración 7. Efecto de las cargas sobre el elemento

ICYA 4422 – Primer semestre de 2023 – Pag. 6

Departamento de Ingeniería Civil y Ambiental

ICYA 4422 – Primer semestre de 2023 – Pag. 7

Departamento de Ingeniería Civil y Ambiental

190 kN/m 190 kN/m

2187.79 𝑘𝑁 4272.21 𝑘𝑁 4833.78 𝑘𝑁 3146.22 𝑘𝑁

Ilustración 8. Fuerzas internas sobre los elementos de la viga

Ilustración 9. Diagrama de cortante por efecto de la carga distribuida

ICYA 4422 – Primer semestre de 2023 – Pag. 8

Departamento de Ingeniería Civil y Ambiental

Ilustración 10. Diagrama de momento por acción de la carga distribuida

b) Calcule y dibuje el diagrama de momentos de la viga debido al efecto del

En primer lugar, es necesario calcular el vector de fuerzas {𝐹𝑙 } debido al asentamiento.

Ilustración 11. Efecto del asentamiento

3 ⋅ 28000 𝑀𝑝𝑎 ⋅ 4.02 𝑚2 3 ⋅ 28000 𝑀𝑝𝑎 ⋅ 4.02 𝑚

ICYA 4422 – Primer semestre de 2023 – Pag. 9

Departamento de Ingeniería Civil y Ambiental

{𝐹𝑙 } = {𝐹𝑛 } − {𝐹𝑜 }

Ahora se prosigue a calcular el desplazamiento debido a la acción del asentamiento.

{𝑈𝐿 } = [𝐾𝑙𝑙 ]−1 ⋅ {𝐹𝑙 }

{𝑈𝐿 } = −0.000168 𝑟𝑎𝑑

Ilustración 12. Efecto de U1=1

ICYA 4422 – Primer semestre de 2023 – Pag. 10