Departamento de Ingeniería Civil y Ambiental
Facultad de Ingeniería
Análisis Avanzado de Estructuras ICYA 4422, 2023-10
Tarea 4
Profesor Juan Carlos Reyes
jureyes@[Link]
Objetivos y parámetros
Objetivos : (1) Analizar estructuras teniendo en cuenta no-linealidades geométricas. (2) Codificar
en Matlab las ecuaciones básicas de elementos con no-linealidades geométricas. (3)
Determinar modos de pandeo. (4) Analizar estructuras considerando la secuencia
constructiva.
Parámetros : La tarea debe ser desarrollada en grupos de máximo dos estudiantes y debe ser subida
a Bloque Neón antes del domingo 16 de abril a las 11:30 p.m. La calificación máxima es
5.0 y la mínima es 0.0.
Problema 1 (Matlab)
La arquitectura moderna propone formas cada vez más atrevidas para las estructuras. La Figura P1 presenta la
geometría propuesta para un pórtico estructural. El pórtico de la Figura P1 también podría ser un pórtico de acero
donde se ha modelado explícitamente el desplome de las columnas (ver F.[Link] del NSR-10).
a) Codifique las siguientes funciones en Matlab o Phyton:
Función vo: calcula el vector {𝑣 } para un elemento a partir de los datos del elemento, cargas distribuidas
locales (𝑤 y 𝑤 ), y carga 𝑓 colocada a los tercios en dirección del eje local y.
Función Ue: calcula el vector {𝑈 } del elemento 𝑖 a partir de los desplazamientos globales de la estructura
almacenados en el vector {𝑈ℓ }.
Función TG (transformaciones geométricas): calcula el vector {𝑣} y las matrices [𝜆] y [𝑇] del elemento 𝑖.
Función kbq: calcula la matriz [𝑘 ] y el vector {𝑞} del elemento 𝑖. En esta función, es necesario calcular
primero 𝑞 como 𝐸𝐴/𝐿 (𝑣 − 𝑣 ). A partir de esta carga axial se calculan los otros parámetros. Suponga que
el material es lineal-elástico pero que la geometría es no-lineal.
Función KeP: calcula la matriz de rigidez [𝐾 ] y el vector {𝑃} del elemento 𝑖. Esta función requiere el uso de
todas las funciones anteriores. Note que para obtener {𝑃} es necesario calcular {𝑝} y transformarlo a
coordenadas globales usando la matriz [𝜆]. {𝑝} se calcula por equilibrio usando las fuerzas y momentos
almacenados en {𝑞}.
Función EstadoCR (estado formulación co-rotacional): calcula el estado de la estructura (vector {𝑃ℓ } y matriz
de rigidez tangente [𝐾ℓℓ ]), conocidos los seis desplazamientos en sus extremos y las propiedades de los
elementos. Esta función usa todas las funciones desarrolladas previamente.
b) Usando las funciones creadas en a), calcule el estado de la estructura (vector {𝑃ℓ } y matriz [𝐾ℓℓ ]) para los
desplazamientos mostrados en la Figura P2. El nodo “a” tiene grados de libertad 𝑈 , 𝑈 y 𝑈 y el “b” 𝑈 ,𝑈 y 𝑈 .
c) Calcule el vector {𝐹 } − {𝑃ℓ }. ¿Se cumple el equilibrio? Usando el método de Newton, haga una nueva propuesta
de desplazamientos donde el error en {𝐹 } − {𝑃ℓ } se reduzca.
120 kN 120 kN E=200000 MPa
b d Sección W10x12
(flexión alrededor del eje fuerte)
U1=54.52 mm
4m
U2=-16.00 mm
U3=-0.00889 rad
a c U4=54.22 mm
U5=-14.90 mm
4m 1m U6=-0.00792 rad
Figura P1. Estructura y sus desplazamientos.
ICYA 4422 – Primer semestre de 2023 – Pag. 1
Departamento de Ingeniería Civil y Ambiental
Facultad de Ingeniería
Análisis Avanzado de Estructuras ICYA 4422, 2023-10
Tarea 4
Profesor Juan Carlos Reyes
jureyes@[Link]
Objetivos y parámetros
Objetivos : (1) Analizar estructuras teniendo en cuenta no-linealidades geométricas. (2) Codificar
en Matlab las ecuaciones básicas de elementos con no-linealidades geométricas. (3)
Determinar modos de pandeo. (4) Analizar estructuras considerando la secuencia
constructiva.
Parámetros : La tarea debe ser desarrollada en grupos de máximo dos estudiantes y debe ser subida
a Bloque Neón antes del domingo 16 de abril a las 11:30 p.m. La calificación máxima es
5.0 y la mínima es 0.0.
Problema 2 (Manual)
Las vigas del pórtico plano de la Figura P2 son infinitamente rígidas. Las columnas son secciones W14×38 trabajando
alrededor de su eje fuerte. El pórtico se encuentra arriostrado fuera del plano. Realice las siguientes actividades:
a) Usando Sap2000, calcule la carga axial sobre las columnas del pórtico cuando una carga vertical de 500 kN/m en el
sentido de la gravedad actúa sobre las vigas bf, cg y dh. Aplique las cargas verticales al modelo dentro de un “load
pattern” llamado “Vertical”.
b) Calcule el factor por el cual se deben multiplicar las cargas axiales calculadas en la parte a) para producir pandeo
elástico; dibuje los modos de pandeo. Utilice formulación 𝑃Δ y 𝑃𝛿. Calcule el valor de en todas las columnas
para el primer modo de pandeo usando una carga axial igual a 𝛼1 𝑞1 . Verifique los resultados usando SAP2000 y
comente sobre cómo cambian los resultados si las columnas se dividen en cuatro pedazos.
d h
W14x38 W14x38
c g
W14x38 W14x38
b f
ab=bc=cd=4.00 m
W14x38 W14x38
ae=8.00m
a E=200 GPa e
Figura P3. Pórtico de tres pisos.
ICYA 4422 – Primer semestre de 2023 – Pag. 1
Departamento de Ingeniería Civil y Ambiental
Facultad de Ingeniería
Análisis Avanzado de Estructuras ICYA 4422, 2023-10
Tarea 4
Profesor Juan Carlos Reyes
jureyes@[Link]
Objetivos y parámetros
Objetivos : (1) Analizar estructuras teniendo en cuenta no-linealidades geométricas. (2) Codificar
en Matlab las ecuaciones básicas de elementos con no-linealidades geométricas. (3)
Determinar modos de pandeo. (4) Analizar estructuras considerando la secuencia
constructiva.
Parámetros : La tarea debe ser desarrollada en grupos de máximo dos estudiantes y debe ser subida
a Bloque Neón antes del domingo 16 de abril a las 11:30 p.m. La calificación máxima es
5.0 y la mínima es 0.0.
Problema 3 (SAP2000)
En el pasado, era muy común diseñar edificios aporticados con vigas de transición. Los edificios de Ecopetrol ubicados
en la carrera 13 con calle 36 son un ejemplo de esto. En esta clase de edificios, la modelación de la secuencia
constructiva puede ser determinante en los resultados del análisis. El pórtico de la Figura P4 tiene tres luces de 5
metros y cuatro pisos con alturas de entrepiso de 3m. La viga de 30×65cm es una viga de transición. Todas las vigas
tienen una carga vertical (hacia abajo) de 75 kN/m. El peso propio de los elementos se puede despreciar.
Usando el programa SAP2000, calcule el desplazamiento vertical máximo y los diagramas de momento del pórtico
haciendo los siguientes análisis:
a) Modelación completa del pórtico sometido a las cargas verticales mencionadas previamente. Conduzca un análisis
lineal estático convencional.
b) Modelación automática de secuencia constructiva sin usar efectos P-Delta.
c) Modelación automática de secuencia constructiva usando efectos P-Delta.
d) Compare los resultados obtenidos en las partes a), b), y c) y concluya.
V30x30
V30x35
C45x45 C45x45
V30x40
V30x65
C45x45 E=25000 MPa C45x45
Figura P3. Pórtico con viga de transición.
ICYA 4422 – Primer semestre de 2023 – Pag. 1
