ANÁLISIS MATRICIAL DE ESTRUCTURAS

EJERCICIO N.º 02 - PARRILLAS

01. Determinar la matriz de transformación de desplazamientos [A].

Solución:

1º Definimos las propiedades de "EI y EA" y asignamos numeración a la viga continua.

2º Determinamos los grados de libertad y definimos el sistema global {Q-D}.

3º Definimos el sistema local {q-d}

MAESTRÍA EN CIENCIA E INGENIERÍA INGENIERÍA ESTRUCTURAL

 ANÁLISIS MATRICIAL DE ESTRUCTURAS

4º Determinamos la matriz de transformación [A1], realizando D1=1, D2=D3=0

⎡0⎤
⎢1⎥
⎢ ⎥ 1
0⎥
A1 ≔ ⎢------------
⎢0⎥
⎢1⎥
⎢⎣ 0 ⎥⎦ 2

5º Determinamos la matriz de transformación [A2], realizando D2=1, D1=D3=0

⎡ 1 ⎤
⎢ 0 ⎥
⎢ ⎥ 1
⎢ _ _0_ _⎥_ _ _
A2 ≔ ⎢ -0.5 ⎥
⎢ 0 ⎥ 2
⎢ ⎥
- ‾‾
⎢ ―― 3⎥
⎢⎣ 2 ⎥⎦

6º Determinamos la matriz de transformación [A3], realizando D3=1, D1=D2=0

⎡ 0 ⎤
⎢ 0 ⎥
⎢ ⎥ 1
⎢ _ 1
_ _ ⎥_ _ _ _
A3 ≔ ⎢ ‾‾ 3 ⎥
⎢ ―― ⎥ 2
⎢ 2 ⎥
⎢ 0 ⎥
⎢⎣ -0.5 ⎥⎦

MAESTRÍA EN CIENCIA E INGENIERÍA INGENIERÍA ESTRUCTURAL

 ANÁLISIS MATRICIAL DE ESTRUCTURAS

10º Determinamos la matriz de transformación [A]

⎡0 1 0 ⎤
⎢1 0 0 ⎥ 1
⎢ ⎥
0 0 1
A ≔ augment ((A1 , A2 , A3)) A = ⎢ -------------------------------------
⎥
⎢ 0 -0.5 0.866 ⎥
⎢1 0 0 ⎥ 2
⎢⎣ 0 -0.866 -0.5 ⎥⎦

02. Determinación de la Matriz de Rigidez [K] por Ensamblaje

1º Propiedades de la Sección
Predimensionamiento de la sección

L1 ≔ 600 Dimensión elemento 01 (cm)

L2 ≔ 578 Dimensión elemento 02 (cm)

Lmax ≔ 600
Lmax
h ≔ ――― h = 75 Altura de la sección (cm)
8
h
b≔― Base de la sección (cm)
2
b = 37.5
Seleccionamos las dimensiones por criterio
bdis ≔ 40 Base de la sección de diseño
(cm)
hdis ≔ 80 Altura de la sección de diseño
⎛⎝hdis 3 ⎞⎠ (cm)
Ixx ≔ bdis 1 ――― Inercia de la sección dirección xx
12 (cm4)
Ixx = 1.707 ⋅ 10 6

MAESTRÍA EN CIENCIA E INGENIERÍA INGENIERÍA ESTRUCTURAL

 ANÁLISIS MATRICIAL DE ESTRUCTURAS

⎛⎝hdis 2 ⎞⎠
J ≔ bdis 1 ――― Modulo de sección (cm3)
6
J = 4.267 ⋅ 10 4

2º Propiedades de los Materiales

f´c ≔ 280 Resistencia a la compresión del
concreto (Kg/cm2)

E ≔ 15000 ⋅ ‾‾‾
f´c Modulo de Elasticidad del
Concreto

E = 2.51 ⋅ 10 5

u ≔ 0.25 Modulo de Poisson

E
G ≔ ――― Modulo de Corte del concreto
2 (1 + u))
(
(Kg/cm2)

G = 1.004 ⋅ 10 5

3º Matriz de Rigidez Local

⎡ J ⎤
⎢ G ⋅ ―― 0 0 ⎥
⎢ L1 ⎥
K1 ≔ ⎢ 0
Ixx
12 ⋅ E ⋅ ――
Ixx ⎥
6 ⋅ E ⋅ ――
⎢ L1 3 L1 2 ⎥
⎢ ⎥
⎢ 0 Ixx Ixx ⎥
6 ⋅ E ⋅ ―― 4 ⋅ E ⋅ ――
⎢⎣ L1 2 L1 ⎥⎦

⎡ 7.139 ⋅ 10 6 0 0 ⎤
⎢ 4 ⎥
K1 = ⎢ 0 2.38 ⋅ 10 7.139 ⋅ 10 6 ⎥
⎢⎣ 0 7.139 ⋅ 10 6 2.856 ⋅ 10 9 ⎥⎦

⎡ J ⎤
⎢ G ⋅ ―― 0 0 ⎥
⎢ L2 ⎥
K2 ≔ ⎢ 0
Ixx
12 ⋅ E ⋅ ――
Ixx ⎥
6 ⋅ E ⋅ ――
⎢ L2 3 L2 2 ⎥
⎢ ⎥
⎢ 0 Ixx Ixx ⎥
6 ⋅ E ⋅ ―― 4 ⋅ E ⋅ ――
⎢⎣ L2 2 L2 ⎥⎦

MAESTRÍA EN CIENCIA E INGENIERÍA INGENIERÍA ESTRUCTURAL

 ANÁLISIS MATRICIAL DE ESTRUCTURAS

⎡ 7.411 ⋅ 10 6 0 0 ⎤
⎢ ⎥
K2 = ⎢ 0 2.662 ⋅ 10 4 7.693 ⋅ 10 6 ⎥
⎢⎣ 0 7.693 ⋅ 10 6 2.964 ⋅ 10 9 ⎥⎦

4º Ensamblaje de la Matriz de Rigidez

⎡A A A ⎤
⎢ 1,1 1,2 1,3⎥ ⎡0 1 0⎤
A1 ≔ ⎢ A A A ⎥ A1 = ⎢ 1 0 0 ⎥
⎢ 2,1 2,2 2,3⎥ ⎢ ⎥
⎣0 0 1⎦
⎢A A A ⎥
⎣ 3,1 3,2 3,3⎦

⎡A A A ⎤
⎢ 4,1 4,2 4,3⎥ ⎡ 0 -0.5 0.866 ⎤
A2 ≔ ⎢ A A A ⎥ ⎢
A2 = 1 0 0 ⎥
⎢ 5,1 5,2 5,3⎥ ⎢ ⎥
⎣ 0 -0.866 -0.5 ⎦
⎢A A A ⎥
⎣ 6,1 6,2 6,3⎦

K ≔ ((A1)) T ⋅ ((K1)) ⋅ A1 + A2 T ⋅ K2 ⋅ A2
⎡ 5.042 ⋅ 10 4 -6.663 ⋅ 10 6 3.293 ⋅ 10 6 ⎤
Matriz de Rigidez de la
⎢ ⎥
K = ⎢ -6.663 ⋅ 10 6 2.232 ⋅ 10 9 1.28 ⋅ 10 9 ⎥ Estructura
6
⎢⎣ 3.293 ⋅ 10 1.28 ⋅ 10 9 3.602 ⋅ 10 9 ⎥⎦

5º Solución por superposición de cargas

5.1. Solución del problema primario

MAESTRÍA EN CIENCIA E INGENIERÍA INGENIERÍA ESTRUCTURAL

 ANÁLISIS MATRICIAL DE ESTRUCTURAS

A) Cargas Nodales

⎡0⎤
F≔⎢0⎥
⎢ ⎥
⎣0⎦

⎡0⎤
R1 ≔ -F R1 = ⎢ 0 ⎥
⎢ ⎥
⎣0⎦

B) Cargas En las Barras

P ≔ 12000 Cargas Puntuales (Kg)
w ≔ 180 Cargas Distribuidas (Kg/cm)

⎡ 0 ⎤
⎢ -P ⎥
⎢ ⎥ ⎡ 0 ⎤
L1 ⎥
⎢ -2 ⋅ P ⋅ ―― ⎢ -1.2 ⋅ 10 4 ⎥ ⎡ qemp ⎤
⎢ 9 ⎥ ⎢ 6
⎥ ⎢ 1,1
⎥
qemp ≔ ⎢ 0 ⎥ qemp = ⎢ -1.6 ⋅ 10 ⎥ qemp1 ≔ ⎢ qemp ⎥
⎢ ⎢ ⎥
L2 ⎥ 0 2,1
⎢ ⎥
⎢ -w ⋅ ―― ⎥ ⎢ -5.202 ⋅ 10 4 ⎥ ⎢ qemp ⎥
⎢ 2 ⎥ ⎢ 6
⎥ ⎣ 3,1⎦
⎢ L2 2 ⎥ ⎣ -5.011 ⋅ 10 ⎦
⎢ -w ⋅ ―― ⎥ ⎡ qemp ⎤
⎣ 12 ⎦ ⎢ 4,1
⎥
qemp2 ≔ ⎢ qemp ⎥
5,1
⎢ ⎥
⎢ qemp ⎥
⎣ 6,1⎦

MAESTRÍA EN CIENCIA E INGENIERÍA INGENIERÍA ESTRUCTURAL

 ANÁLISIS MATRICIAL DE ESTRUCTURAS

R2 ≔ A T ⋅ qemp

⎡ -6.402 ⋅ 10 4 ⎤
⎢ ⎥
R2 = ⎢ 4.34 ⋅ 10 6 ⎥
⎢⎣ 9.056 ⋅ 10 5 ⎥⎦
⎡ -6.402 ⋅ 10 4 ⎤
⎢ ⎥
R ≔ R1 + R2 R = ⎢ 4.34 ⋅ 10 6 ⎥
⎢⎣ 9.056 ⋅ 10 5 ⎥⎦

5.2. Solución del problema Complementario

A) Vector de Cargas Nodales

⎡ 6.402 ⋅ 10 4 ⎤
⎢ ⎥
Q ≔ -R Q = ⎢ -4.34 ⋅ 10 6 ⎥
⎢⎣ -9.056 ⋅ 10 5 ⎥⎦

B) Desplazamientos Nodales
D ≔ K -1 ⋅ Q

⎡ 3.658 ⎤
D = ⎢ 0.014 ⎥
⎢ ⎥
⎣ -0.009 ⎦

C) Deformación en los Extremos de los elementos

⎡ 0.014 ⎤
d≔A⋅D ⎢ 3.658 ⎥ ⎡d ⎤ ⎡d ⎤
⎢ ⎥ ⎢ 1,1⎥ ⎢ 4,1⎥
⎢ -0.009 ⎥
d= d1 ≔ ⎢ d ⎥ d2 ≔ ⎢ d ⎥
⎢ -0.014 ⎥ ⎢ 2,1⎥ ⎢ 5,1⎥
⎢ 3.658 ⎥ ⎢d ⎥ ⎢d ⎥
⎢⎣ -0.008 ⎥⎦ ⎣ 3,1⎦ ⎣ 6,1⎦

⎡ 0.014 ⎤ ⎡ -0.014 ⎤
d1 = 3.658 ⎥
⎢ d2 = ⎢ 3.658 ⎥
⎢ ⎥ ⎢ ⎥
⎣ -0.009 ⎦ ⎣ -0.008 ⎦

D) Solución del Problema Complementario (Fuerzas en los elementos en la

dirección de las coordenadas locales)

qcomp1 ≔ K1 ⋅ d1 ⎡ 9.897 ⋅ 10 4 ⎤
⎢ ⎥
qcomp1 = ⎢ 2.621 ⋅ 10 4 ⎥
⎢⎣ 1.78 ⋅ 10 6 ⎥⎦

qcomp2 ≔ K2 ⋅ d2 ⎡ -1.061 ⋅ 10 5 ⎤
⎢ ⎥
qcomp2 = ⎢ 3.781 ⋅ 10 4 ⎥
⎢⎣ 5.187 ⋅ 10 6 ⎥⎦

MAESTRÍA EN CIENCIA E INGENIERÍA INGENIERÍA ESTRUCTURAL

 ANÁLISIS MATRICIAL DE ESTRUCTURAS

5.3. Respuesta Estructural Unidades (kg, cm y rad)

⎡ 9.897 ⋅ 10 4 ⎤ Kg-cm/rad
⎢ ⎥
q1 ≔ qemp1 + qcomp1 q1 = ⎢ 1.421 ⋅ 10 4 ⎥ Kg
⎢⎣ 1.796 ⋅ 10 5 ⎥⎦ Kg-cm

⎡ -1.061 ⋅ 10 5 ⎤
⎢ ⎥
q2 ≔ qemp2 + qcomp2 q2 = ⎢ -1.421 ⋅ 10 4 ⎥
⎢⎣ 1.755 ⋅ 10 5 ⎥⎦

MAESTRÍA EN CIENCIA E INGENIERÍA INGENIERÍA ESTRUCTURAL