Departamento de Ingeniería Civil y Ambiental

Facultad de Ingeniería 
  Análisis Avanzado de Estructuras ICYA 4422, 2020‐10 

Tarea 2 
Profesor Juan Carlos Reyes 
[email protected] 
 
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 miércoles 19 de febrero a las 2:00 p.m. en el buzón de la oficina 
ML‐330  ubicado  cerca  a  la  cafetera  de  las  oficinas  del  tercer  piso.  La  calificación 
máxima es 5.0 y la mínima es 0.0 
 
Problema 1 (implementación manual) 
Una excelente opción para rehabilitar vigas de puentes con problemas de deflexiones es adicionar cables postensados 
externos  como  se  muestra  en  la  Figura  P1a.  Por  ejemplo,  en  Bogotá,  el  puente  de  la  Avenida  68  con  Calle  26  fue 
rehabilitado  con  ese  sistema.  Suponga  que  la  viga  “abc”  de  la  Figura  P1b  esta  postensada  externamente;  los  cables 
postensados se pueden representar mediante los elementos “ad” y “dc”, los cuales tienen una carga inicial de presfuerzo 
igual  a  qo.  La  viga  “abc”  tiene  una  carga  distribuida  vertical  de  32  kN/m  y  EcIc=7814000  kN‐m2.  En  una  primera  fase 
constructiva,  la  carga  vertical  actúa  sobre  la  viga  “abc”; posteriormente,  en una  segunda  fase,  se  instalan  los  cables 
presforzados  “ad”  y  “dc”,  los  cuales  tienen  un  EsAs=315840  kN.  Los  grados  de  libertad  de  los  nodos  “b”  y  “d”  esta 
relacionados por un “constraint” que hace que los desplazamientos verticales y el giro de estos dos nodos sean iguales. 
Por simplicidad se pueden despreciar las deformaciones axiales de la viga “abc”. 
a) Calcule la deflexión vertical del punto “b” de la estructura en función de la carga de presfuerzo qo al final de la fase 
2. Use el mínimo número de grados de libertad posible; el problema se puede resolver solo usando solo un grado de 
libertad. 
b) Calcule el valor de qo requerido para que la deflexión de la viga por carga vertical sea 0.5cm hacia arriba. 
c) Calcule la fuerza de postensado que se debería hacer en el gato durante el proceso constructivo para que la deflexión 
de la viga por carga vertical sea 1 cm hacia arriba; esta es la fuerza que el diseñador estructural debería consignar en 
los planos. Verifique su respuesta usando Opensees. 
 

 
(a) 
a  c
b

ab=bc=17.0m 
d bd=1.55m 

 
(b) 
Figura P1. (a) Aplicaciones de postensado externo. (b) Viga del puente a analizar. 
 
ICYA 4422 – Primer semestre de 2020 – Pag. 1 
 Departamento de Ingeniería Civil y Ambiental
Facultad de Ingeniería 
  Análisis Avanzado de Estructuras ICYA 4422, 2020‐10 

Problema 2 
Calcule  el  vector  de  cargas  𝑃   para  el  elemento  de  longitud  𝐿  mostrado  en  la  Figura  P2  en  función  de  los  datos 
mostrados. 
 
EI, EA

e 

Presfuerzo fo 
 
Figura P2. Viga. 
 
Problema 3 (implementación manual) 
Considere el puente continuo postensado de la Figura P3 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 
la deformación axial de los elementos y cambiando adecuadamente los apoyos. 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 resorte en el nodo e tiene una constante 
kL=1000MN/m y representa la rigidez lateral de la cabeza del estribo y el terraplen. El punto f sufre un asentamiento de 
3cm 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.  
a 

b  280kN
190kN/m 190kN/m 
c 
2400kN kL 
a c e 
d b d 

e Elemento fc 
A=7.02m2  13 m
I=1.23m4 
E=28000 MPa 
f 

34 m 42 m 

12000
2000 1000 400 500 4200

200
300
E=28000 MPa 
E=28000 MPa 
I=4.02 m44; A=5.64 m
I=4.02 m ; A=5.64 m22  
1200 Centroide= 1.3553 m 
Centroide= 1.3553 m 
(Medido desde la fibra inferior) 
(Medido desde la fibra inferior) 
300
200 *Dimensiones en mm 
 
Figura P3. Vista isométrica, longitudinal y sección transversal del puente 

ICYA 4422 – Primer semestre de 2020 – Pag. 2 
 Departamento de Ingeniería Civil y Ambiental
Facultad de Ingeniería 
  Análisis Avanzado de Estructuras ICYA 4422, 2020‐10 

Problema 4 (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 contenidos en la variable Elemento{i}. 
Adicional a la información de entrada de la función Po de la tarea 1, Elemento{i} contiene 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  Elemento{i}.lambda.  Cuando  se  llame  la  función,  su  resultado  se  puede  asignar  a 
Elemento{i}.Po 
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 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. 
 
Problema 5 (implementación computacional) 
La estructura de la Figura P5 tiene elementos viga‐columna de acero estructural ASTM A992 (𝐸 =200 GPa, 𝛼 =1.17E‐
5/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 4. 
b) ¿Por qué la deformada del elemento “bc” muestra curvaturas negativas mientras que el momento genera tensiones 
en la cara interna? Justifique su respuesta. 
 
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 
X  h= 4.50E2 mm 
 
Figura P5. Estructura a analizar 
 
 

ICYA 4422 – Primer semestre de 2020 – Pag. 3 
 Departamento de Ingeniería Civil y Ambiental
Facultad de Ingeniería 
  Análisis Avanzado de Estructuras ICYA 4422, 2020‐10 

Problema 6 (implementación computacional) 
Usted ha sido comisionado para revisar el puente peatonal atirantado de la Figura P6 utilizando el código desarrollado 
en el problema 4. 
a) Dibuje la deformada y los diagramas de fuerza cortante, momento flector y carga axial de todos los elementos, 
utilizando el modelo estructural plano de la Figura P3.  
b) Calcule la carga de presfuerzo necesaria en el cable “bd“ para reducir el desplazamiento vertical del nodo “b” a cero. 
Desprecie el efecto de la catenaria en el cable. No resuelva el problema iterativamente sino plantee alguna ecuación 
similar a la utilizada en el documento “Puentes peatonales atirantados ‐ Garcia y Perez (1996)” disponible en la 
carpeta “Material Adicional” en SicuaPlus. 
c) ¿Qué fuerza de presfuerzo debe aplicar el gato durante el proceso constructivo para para reducir el desplazamiento 
vertical del nodo “b” a cero?  
d) Verifique con SAP2000 su respuesta del literal c) y d). 
e) ¿Qué piensa usted que puede pasar si alguien corta los cables? ¿Por qué? Justifíquelo numéricamente. 
 

d
Modelo 
Estructural 

a b c 

 
Figura P6. Vista isométrica y detalles del puente peatonal 
 
 

ICYA 4422 – Primer semestre de 2020 – Pag. 4