Análisis Matricial de Estructuras, VII Nivel, Carrera de Ingeniería Civil, ULEAM. Docente: Ing. Abel Zambrano Mejía, MSc.

37

UNIDAD 2: ANALISIS MATRICIAL DE EMPARRILLADOS PLANOS

Ejercicio de aplicación
Calcular el emparrillado plano por el método de rigidez usando MATLAB. La sección de todos los
elementos es de 30x50cm y el módulo de elasticidad del hormigón es de 1.8x106T/m2. Considere
las deformaciones por cortante. Trabajar en unidades T-m.

5T/m
5T/m

X
Z
Análisis Matricial de Estructuras, VII Nivel, Carrera de Ingeniería Civil, ULEAM. Docente: Ing. Abel Zambrano Mejía, MSc. 38

UNIDAD 2: ANALISIS MATRICIAL DE EMPARRILLADOS PLANOS

5m 5m 3.5m

x x
Z
5 z 7 z 4.5m
5.5m

2 x 3 x
1 x
x z x z
z
4 z z
6 4.5m

3.5m
Análisis Matricial de Estructuras, VII Nivel, Carrera de Ingeniería Civil, ULEAM. Docente: Ing. Abel Zambrano Mejía, MSc. 39

UNIDAD 2: ANALISIS MATRICIAL DE EMPARRILLADOS PLANOS

MODELO MATEMÁTICO DE LA ESTRUCTURA

Y
r20 r23
r19 r22
X r21 r24
Z 4 7 r18
r17 3 6 9
1
5 8 10
r12 r15
2
r11 r14
10 grados de libertad y 14 reacciones
r13 r16

Numeración de los grados de libertad y reacciones de la estructura

Análisis Matricial de Estructuras, VII Nivel, Carrera de Ingeniería Civil, ULEAM. Docente: Ing. Abel Zambrano Mejía, MSc. 40

MATRICES DE RIGIDEZ DE LOS ELEMENTOS EN COORDENADAS LOCALES

Área y momentos de inercia

2 𝑏 × ℎ3 0.3 × 0.53
𝐴 = 0.3 × 0.5 = 0.15𝑚 𝐼𝑧 = = = 0.003125𝑚4
12 12

La inercia alrededor del eje x (eje longitudinal) o inercia torsional para una sección rectangular, se calcula con la
siguiente fórmula:

3
1 0.3 0.34
𝐼𝑥 = 𝛽ℎ𝑏 𝛽 = − 0.21 1− 4
= 0.20869
3 0.5 12 0.5

1 𝑏 𝑏4
𝛽 = − 0.21 1 − 𝐼𝑥 = 0.20869 0.5 0.3 3 = 0.00281737𝑚4
3 ℎ 12ℎ4
Análisis Matricial de Estructuras, VII Nivel, Carrera de Ingeniería Civil, ULEAM. Docente: Ing. Abel Zambrano Mejía, MSc. 41

MATRICES DE RIGIDEZ DE LOS ELEMENTOS EN COORDENADAS LOCALES

Longitudes de los elementos

𝐿1 = 52 + 12 = 5.099𝑚 𝐿2 = 3.5𝑚 𝐿3 = 5𝑚

𝐿4 = 𝐿5 = 𝐿6 = 4.5𝑚 𝐿7 = 4.52 + 1.52 = 4.743𝑚

Cálculo de las matrices de rigidez

Matriz de rigidez del elemento 1

[k1]=k_EP_corte(A,f,Iz,Ix,L1,E,G)

Cabe mencionar que 𝑓 es el factor de forma, que en este caso

por tratarse de una sección rectangular 𝑓 = 1.2.
Análisis Matricial de Estructuras, VII Nivel, Carrera de Ingeniería Civil, ULEAM. Docente: Ing. Abel Zambrano Mejía, MSc. 42

MATRICES DE RIGIDEZ DE LOS ELEMENTOS EN COORDENADAS LOCALES

[k2]=k_EP_corte(A,f,Iz,Ix,L2,E,G) [k3]=k_EP_corte(A,f,Iz,Ix,L3,E,G)

[k4]=k_EP_corte(A,f,Iz,Ix,L4,E,G) [k7]=k_EP_corte(A,f,Iz,Ix,L7,E,G)
Análisis Matricial de Estructuras, VII Nivel, Carrera de Ingeniería Civil, ULEAM. Docente: Ing. Abel Zambrano Mejía, MSc. 43

MATRIZ DE TRANSFORMACIÓN DE LOS ELEMENTOS

Se calculan los ángulos α, 𝛽1 y 𝛽2 considerando el eje local de cada elemento mostrado

Elemento α β1 β2
1 -11.31° 0 0
2 0 0 0
3 0 0 0
4 -90° 0 0
5 -90° 0 0
6 -90° 0 0
7 -71.57° 0 0
Análisis Matricial de Estructuras, VII Nivel, Carrera de Ingeniería Civil, ULEAM. Docente: Ing. Abel Zambrano Mejía, MSc. 44

MATRIZ DE TRANSFORMACIÓN DE LOS ELEMENTOS

Cálculo de las matrices de transformación
Se utiliza la rutina “Th_EP” para calcular la matriz de transformación de cada elemento de la
siguiente manera:

%[T,h]=Th_EP(L,alfa,beta1,beta2)
[T1,h1]=Th_EP(L1,-11.31,0,0);
[T2,h2]=Th_EP(L2,0,0,0);
[T3,h3]=Th_EP(L3,0,0,0);
[T4,h4]=Th_EP(L4,-90,0,0);
T5=T4; h5=h4;
T6=T4; h6=h4;
[T7,h7]=Th_EP(L7,-71.57,0,0);
Análisis Matricial de Estructuras, VII Nivel, Carrera de Ingeniería Civil, ULEAM. Docente: Ing. Abel Zambrano Mejía, MSc. 45

CÁLCULO DE LAS MATRICES DE TRANSFORMACIÓN

La rutina “Th_EP” permite calcular dos variables de salida [T,h] para cada elemento. En total son 7 matrices
“T” y 7 matrices “h”, cuyos resultados se muestran a continuación
Análisis Matricial de Estructuras, VII Nivel, Carrera de Ingeniería Civil, ULEAM. Docente: Ing. Abel Zambrano Mejía, MSc. 46

CÁLCULO DE LAS MATRICES DE TRANSFORMACIÓN

Se recuerda que en este ejercicio T5=T6=T4 y h5=h6=h4. Debido a que los elementos 4, 5 y 6, poseen los mismos
datos en longitud y ángulos α, β1, β2.
 Análisis Matricial de Estructuras, VII Nivel, Carrera de Ingeniería Civil, ULEAM. Docente: Ing. Abel Zambrano Mejía, MSc. 47

MATRICES DE RIGIDEZ DE LOS ELEMENTOS EN COORDENADAS DE NUDO

Triple producto matricial

𝑲 = 𝑻 × 𝒌 × 𝑻𝑻

𝑲 = matriz de rigidez del elemento en coordenadas de nudo

𝒌= matriz de rigidez del elemento en coordenadas locales reducidas

𝑻 = matriz de transformación de un sistema local reducido a un sistema de nudos.

K1=T1*k1*T1'; K5=K4;
La matriz transpuesta en MATLAB se
K2=T2*k2*T2'; K6=K4;
coloca con un apóstrofe ( ‘ )
K3=T3*k3*T3'; K7=T7*k7*T7’;
K4=T4*k4*T4';
 Análisis Matricial de Estructuras, VII Nivel, Carrera de Ingeniería Civil, ULEAM. Docente: Ing. Abel Zambrano Mejía, MSc. 48

MATRICES DE RIGIDEZ DE LOS ELEMENTOS EN COORDENADAS DE NUDO

 Análisis Matricial de Estructuras, VII Nivel, Carrera de Ingeniería Civil, ULEAM. Docente: Ing. Abel Zambrano Mejía, MSc. 49

MATRICES DE RIGIDEZ DE LOS ELEMENTOS EN COORDENADAS DE NUDO

Análisis Matricial de Estructuras, VII Nivel, Carrera de Ingeniería Civil, ULEAM. Docente: Ing. Abel Zambrano Mejía, MSc. 50

VECTORES DE ENSAMBLAJE DE LOS ELEMENTOS

ve_1=[1 17 2 3 4 5];
ve_2=[3 4 5 6 7 8];
ve_3=[6 7 8 9 18 10];
Vectores de ensamblaje
ve_4=[11 12 13 3 4 5];
ve_5=[3 4 5 19 20 21];
ve_6=[14 15 16 6 7 8];
ve_7=[6 7 8 22 23 24];

Permite el ensamblaje de la matriz de rigidez de la estructura

Análisis Matricial de Estructuras, VII Nivel, Carrera de Ingeniería Civil, ULEAM. Docente: Ing. Abel Zambrano Mejía, MSc. 51

ENSAMBLAJE DE LA MATRIZ DE RIGIDEZ DE LA ESTRUCTURA

Procedimiento de ensamblaje directo utilizando la rutina “Ens”

Se crea primero la matriz de rigidez de la estructura “S” rellena de ceros, cuya dimensión es
de 24x24 considerando los grados de libertad y reacciones del sistema.

S=zeros(24,24) %[S]=Ens(ve,K,S,ngl)
[S]=Ens(ve_1,K1,S,24);
[S]=Ens(ve_2,K2,S,24);
[S]=Ens(ve_3,K3,S,24);
[S]=Ens(ve_4,K4,S,24);
[S]=Ens(ve_5,K5,S,24);
[S]=Ens(ve_6,K6,S,24);
[S]=Ens(ve_7,K7,S,24);
 Análisis Matricial de Estructuras, VII Nivel, Carrera de Ingeniería Civil, ULEAM. Docente: Ing. Abel Zambrano Mejía, MSc. 52

MATRICES DE RIGIDEZ DE LOS ELEMENTOS EN COORDENADAS DE NUDO

1 17 2 3 4 5

3 4 5 6 7 8
 Análisis Matricial de Estructuras, VII Nivel, Carrera de Ingeniería Civil, ULEAM. Docente: Ing. Abel Zambrano Mejía, MSc. 53

MATRICES DE RIGIDEZ DE LOS ELEMENTOS EN COORDENADAS DE NUDO

6 7 8 9 18 10

11 12 13 3 4 5

6 7 8 22 23 24
Análisis Matricial de Estructuras, VII Nivel, Carrera de Ingeniería Civil, ULEAM. Docente: Ing. Abel Zambrano Mejía, MSc. 54

MATRIZ DE RIGIDEZ DE LA ESTRUCTURA

A continuación, se muestra la matriz de rigidez de la estructura “S”:
Análisis Matricial de Estructuras, VII Nivel, Carrera de Ingeniería Civil, ULEAM. Docente: Ing. Abel Zambrano Mejía, MSc. 55

MATRIZ DE RIGIDEZ DE LA ESTRUCTURA

Fin de la presentación