ANALISIS MATRICIAL DE ESTRUCTURAS

Msc. Ing. Gerald Astorga

Temario

1. Conceptos básicos
2. Método de la rigidez
3. Estructuras articuladas
4. Transformación de coordenadas
5. Ensamblaje de la matriz de rigidez de la estructura
6. Apoyos inclinados y apoyos elásticos
7. Apoyos con desplazamiento inicial
8. Análisis tridimensional de barras articuladas
9. Estructuras tipo marco rígido
10. Matriz de rigidez de elementos tipo marco rígido en 3D
11. Rigidez lateral de pórticos

Bibliografía

Bazán Enrique, Meli Roberto. (2002). Diseño sísmico de edificios. Limusa. México.

Gonzalo Cuevas Oxar. (2002). Análisis estructural. Limusa. México.

Hannon Bruce, Ruth Matthias. (2001). Dynamic Modeling. Springer. USA.

Kassimali Aslam. (2001). Análisis Estructural. Thompson. Australia.

Mc cormac Jack, Nelson Jammek. (2002). Análisis de estructuras; Método elástico y matricial.
Alfa Omega. México.

Blanco Conera Suarez. (2015). Análisis matricial de estructuras. CIMNE. España

Roberto Aguilar. (2014). Análisis matricial de estructuras. Universidad Fuerzas Armadas.

Ecuador.

Hayrettin Kardentuncen. (1975). Introducción al análisis estructural con matrices. Mc Graw Hill.
USA.
TEMA 1

CONCEPTOS BÁSICOS

Mecánica: Estudia el comportamiento de las partículas debido a las acciones.

Cuerpo rígido: INDEFORMABLES, clásica newtoniana.

Cuerpo deformable: Teoría Relatividad EINSTEN

Fluidos: Comprensibles-gases.

Incompresibles.

Acciones: Materiales

Diseño H° A°

Metálicas

Maderas

La resistencia de un material está definida por un previo análisis estructural modelando barras
rectas en 2 o 3 dimensiones, básicamente se debe realizar el procedimiento para determinar
los esfuerzos internos y externos lo que implicar cuantificar en sus distintas unidades.

Básicamente los parámetros fundamentales son:

ε = Deformación unitaria del medio continuo.

σ = Esfuerzo asociado a las deformaciones sufridas por el medio continuo.

δ = Desplazamiento de los diferentes puntos del medio continuo que forma la estructura.

Medio continuo ↔ Cuerpo deformable

Proceso de análisis

1.-Que las acciones no superen a los materiales

2.-Que la estructura no se deforme más de lo normal.

 El método de la flexibilidad → reacciones

 El método de la rigidez → deformaciones

{ F }=[ K ] {δ }
{ Vector fuerza }=[ Matriz de rigidez ] { Vector desplazamiento }
Esfuerzo: Fuerza que se aplica en una determinada área

F
σ=
A
Deformación:

δ
ε=
Lo

Mat: → Dúctiles

→ Frágil

Ley de Hooke

Límite de proporcionalidad donde el esfuerzo versus la deformación responden a una función

lineal. El límite de elasticidad es el punto máximo de deformidad y vuelve pero no a su estado
real.

El punto de fluencia es el punto donde sin aumentar carga se sigue deformando.

En la región elástica

σ
tg ϴ= =E
ε
Región plástica

P
σ A PL
E= = = =E
ε δ δA
L

P=Fuerza axial
AE
P= δ
L
P=K δ

Por cualquier método resolver:

 TEMA 2

METODO DE LA RIGIDEZ

Dentro de los métodos clásicos (Cross, Cani) es necesario comprender el comportamiento de

los distintos tipos de estructuras que se tienen.

El método de rigidez o matricial a diferencia de los métodos elásticos considera las

deformaciones como las reacciones a las acciones externas producidas sobre una estructura.

El método de la rigidez o matricial es la base para el desarrollo de software en la actualidad

(SAP, Robots, Lira, Adina).

El método de la rigidez permite la resolución de todo tipo de estructuras y se basa en la

construcción y operación de matrices de rigidez de cada elemento y global de la estructura,
dada la simplicidad de la metodología y lo estructural de los algoritmos de solución y es que
podemos resolver cualquier estructura independiente de su magnitud.

2.1 Origen (1945- 1955)

Análisis de matrices por la aeronáutica.

(1956) Turner, Clough, Martin y Topp)

Análisis

El método consiste en reemplazar la estructura original de elementos estructurales finitos

cuyas propiedades pueden expresarse de forma matricial.
2.2 Consideraciones y conceptos básicos.

Se supone módulo de elasticidad mc constante en toda la longitud de los elementos, es

decir elementos prismáticos.

Las relaciones fundamentales de equilibrio, compatibilidad, fuerza – desplazamiento se

mantienen vigentes.

Los grados de libertad son los desplazamientos independientes en traslación y rotación de los
nodos que son necesarios para especificar la forma deformada de una estructura cuando se
somete a una carga unitaria.

2.3 Comprensiones en el método matricial.

Se considera lo siguiente:
 a) Los nodos se cuentan con un número dentro de un círculo empezando de un nodo
libre.
b) Los elementos se cuentan con un número dentro de un rectángulo o cuadrado
empezando del nodo de menor numeración.
c) Los grados de libertad se representan con flechas en caso de traslación y flechas curvas
en caso de rotación.
Las estructuras se trabajan en sistemas cartesianos ortogonales bajo un eje local y un
eje global en caso de estructuras espaciales.

TEMA 3

ESTRUCTURAS ARTICULADAS

También conocidas como celosías. Las tipologías más usadas por triángulos son:

3.1 Matriz de rigidez de elementos articulados.

3.1.1 Barra simple

{ F }=[ K ] {δ }
 F1 X =K (U 1 −U 2) → F1= K U 1 −K U 2

F2 X =K (U 2−U 1) → F2 =K U 2−K U 1

AE
K=
L
Fx1 - U1
= AE AE
L L
Fx2 - U2
AE AE
L L

F1x 0 - 0 U1x
AE AE
= L L
F1y 0 0 0 0 U1y

Fx2 - 0 0 U2x
AE AE
L L
F2y 0 0 0 0 U2y

3.1.2 Caso barra compuesta (2 tramos).

Barra 1

F1 −K 1U1
= K1
F´1 - −K 1U2
K1
Barra 2

F2´´ −K 2U2
= K2
F3 - −K 2U3
K2

Entonces:

F1 −K 1 0 U1
K1
F2 = - K 1+ K 2 −K 2 U2
K1
F3 0 −K 2 K2 U3

Caso espacial

Matriz de rigidez coplanar para coplanares

F1 K1 −K 1 0 . . . 0

F2 - K 1+ K 2 −K 2 . . . 0
K1
F3 0 −K 2 K 2+ K 3 −K 3 . . 0
=
. . . −K 3 K 3+ K 4 −K 4 . 0
. . . . −K 4 .
. . . . . .
Fn 0 0 0 0 K n+1 + K n

Ejercicio
Barra 1

- -
K1
- EA
EA 3 c m2 2
K1 3cm
L 50 cm
L 50 cm
= =E
- - 2
- EA 3cm
K1
K1 EA 3 c m2
L 50 cm
L 50 cm
Barra 2

- -
- EA 2 c m2
K2
K2
EA 2 c m2
L 70 cm
L 70 cm Barra 3
= =E
- -
- EA 2 c m2
K3
K - K2 = EA - =E 21 c m 22 -
2 L 50 cm2
K3 L EA 70 cm 1cm
L 70 cm
- -
-
K3 EA
EA 2 1 c m2
K3 1cm
L 70 cm
L 70 cm

-
F4 0 0 U1
K1 K1
-
K1+ K2 - K2 0 U2
F3 =
K1
-
F2 0 K2+K3 - K3 U3
K2
F1 0 0 - K3 K3+ K4 U4
Así:
 -
0 0,0886 0 U2
0,028
- -
5000 = E 0,042 U3
0,028 0,014
-
2000 0 0,014 U4
0,014

- - - -
0,042 0,042
0,014 0,028 0,014 0,028
-
- -
0,014 0 0,014 0
0,014 0,014

-
0,028 0 0,0886 0 0,0886
0,028
- -
- -
0,014 0 0,014 0
0,014 0,014

- -
0 0,0886 0 0,0886
0,028 0,028
-
- - - -
0,042 0,042
0,014 0,028 0,014 0,028

0,000392 0,000392 0,000392

0,000392 0,00124 0,00124

0,000392 0,00124 0,002937

Invirtiendo

U2 16,6384 16,6384 16,6384 0 0,0582

=
U3 1 16,6384 51,6407 51,6407 5000 0,1807
L
U4 16,6384 51,6407 121,6197 2000 0,2507

TAREA 3

Para un edificio de 3 niveles de altura total de 8 mts., determinar los desplazamientos

horizontales por nivel para un efecto sísmico replicable estáticamente con fuerzas de acuerdo
al diagrama, hallar la rotación del vínculo.
6
E=2.1∗10 kg/cm2
Hº 0,007

Aº 0,010
Albañilería 0,005

Madera 0,010

D/Hpiso (cm)

Dcm=Ucm* R *0.75

U 1 = 0,0007

D1 = 0,0007 * 0,75*8
Cm m
3
D1 = 4,2 * 10 −¿ /300cm
¿
3
D1 = 4,2*10 −¿ /300cm = 0,000014
¿
0,000014 ≤ 0,0007 → (NORMA)

Ejemplo:

F= K*U
 KN
20 ∗250 m2
6OKN = mm 2 * U
150 mm

U =1,8mm
Elemento 1

1 - 1
K1 = EA
- 1 1
L

Elemento 2

1 - 1
K2 = EA
- 1 1
L

Matriz General

- 0
1
1
- -
K = EA 2
1 1
L
- 1
0
1 F = K U

-
F1 1 0 U1
1
- - AE
F2 = 2 = U2
1 1 L
-
F3 0 1 U3
1

-
F2 2 U2
1 AE
= =
- L
F3 1 1,2 mm
1
 -
60 2 U2
EA 1
= =
L -
0 1 1,2 mm
1

100 100 ( 1 )∗6

60= (2 ) U 2−
3 3 5
200 120
60= U 2−
3 3

U 2=1.5 mm

F1 1 -1 0 0
EA
60 = -1 2 1 1.5
L
F3 0 -1 1 1.2

100
F 1= (1∗0+ (−1 )∗1.5+0∗1.2)
3
F 1=−50 KN
100
F 3= (0∗0+ (−1 )∗1.5+1∗1.2)
3
F 3=−10 KN

Ejercicio

Barra 1
 K1 -K1 1 -1
EA 25∗300 2 -2
= =
L 150
-K1 K1 -1 1 -2 2

Barra 2

300 −300 3 −3
K2 -K2
200 200 2 2
=E =E
−300 300 −3 3
-K2 K2
200 200 2 2

Barra 3

100 −100
K3 -K3 1 -1
100 100
=E =E
−100 100
-K3 K3 -1 1
100 100

F1 2 -2 0 0

F2 -2 3,5 -3/2 0
=E
F3 0 -3/2 5/2 -1

F4 0 0 -1 1

10 3,5 -1,5 0 U2

0 =E -1,5 2,5 -1 U3

-17 0 -1 1 U4
Hallando inversa

3,5 -1 -1,5 -1 -15 3,5

-
-1 1 0 1 0 -1
-
-1,5 0 3,5 0 3,5 1,5
Adj= -
-1 1 0 1 0 -1
-
-1,5 0 3,5 0 3,5 1,5

3,5 -1 -1,5 -1 -1,5 3,5

 1,5 1,5 1,5
Adj t 1,5 3,5 3,5
=
-1,5 3,5 6,5

Det (A) = 3

1
A−1= Adj t
IAI

0,5 0,5 0,5

−1
A 0,5 1,17 1,17
=
-0,5 1,17 2,17

U2 -0,14

U3 = -0,5956 cm

U4 -1,6756

TEMA 4

TRANSFORMACION DE COORDENADAS
 X, Y COORDENADAS GLOBALES

X, Y COORDENADAS LOCALES

Globales para saber desplazamientos globales

x=P cos ( ɵ+ α ) =Pcos ɵ∗cos α−Psen ɵ∗Psen α

y=P sen ( ɵ+α )=Psen ɵ∗cos α + Pcos ɵ∗Psen α

x=P cos α
y=Psen α

Reemplzando en “x” y “y”

x=x cos θ− y sen θ

y=x sen θ+ y cos θ

{xy}=[ cos θ
sen θ ]{ }
−sen θ x
cos θ y

[ T ] → [ T ] =1

[ T ]−1=[ T ]t

{xy}=[ cos θ
sen θ ]{ }
−sen θ x
cos θ y
→ T Globales o locales
 −1

{}[
x = cos θ −sen θ
y sen θ cos θ ]{} x
y

{xy}=[−sen
cos θ
θ
sen θ x
cos θ y ]{ } → J Locales o globales

{ U }= [ T ] {U }
F=[ K ] {U }
F=[ T ] { δ }

4.2 Matriz de rigidez global de una barra articulada

{ U }= [ T ] {U } → { U }= [ J ] { U } ①

{ F }= [ T ] { F } → { F }=[ J ] { F } ②

{ F }= [ K ] {U } Local ③

① y ② en ③

[ J ] { F }= [ K ] [ J ] { U }
−1
{ F }=[ J ] [ K ][ J ] {U }
{F } −1
=[ J ] [ K ][ J ]
{U }
[ J ] Traspuesta de la matriz de coordenadas
Rigidez global

[ K ] =[ J ]−1 [ K ] [ J ]
Para axial

[ ]
AE −AE
K= L L
−AE AE
L L
 [ ]
AE −AE
0 0
L L
0 0 0 0
K=
−AE AE
0 0
L L
0 0 0 0

[ ]
1 0 −1 0
AE 0 0 0 0
K=
L −1 0 1 0
0 0 0 0

−1
[ K ] =[ J ] [ K ] [ J ]
K11 K12
K=
K21 K22

C2 Cs -C2 -Cs
Cs S2 -Cs -S2
K = AE/L
-C2 -Cs C2 Cs
-Cs -S2 Cs S2

Matriz de transformación

cosθ=c
sen θ=s

Ejercicio
barra de x1 y1 a x2 y2 Long (cm)
Area(cm^2)
l m lm l^2 m^2 AE/L

1,0 1,0 0,0 0,0 2,0 1,5 1,0 180,3 25,0 0,8 0,6 0,5 0,7 0,3 277350,1

2,0 1,0 0,0 0,0 3,0 3,0 0,0 300,0 25,0 1,0 0,0 0,0 1,0 0,0 166666,7

3,0 2,0 1,5 1,0 3,0 3,0 0,0 180,3 25,0 0,8 -0,6 -0,5 0,7 0,3 277350,1

4,0 2,0 1,5 1,0 4,0 3,0 2,0 180,3 25,0 0,8 0,6 0,5 0,7 0,3 277350,1

5,0 3,0 3,0 0,0 4,0 3,0 2,0 200,0 25,0 0,0 1,0 0,0 0,0 1,0 250000,0

6,0 3,0 3,0 0,0 5,0 4,5 1,0 180,3 25,0 0,8 0,6 0,5 0,7 0,3 277350,1

7,0 4,0 3,0 2,0 5,0 4,5 1,0 180,3 25,0 0,8 -0,6 -0,5 0,7 0,3 277350,1

8,0 3,0 3,0 0,0 6,0 6,0 0,0 300,0 25,0 1,0 0,0 0,0 1,0 0,0 166666,7

9,0 5,0 4,5 1,0 6,0 6,0 0,0 180,3 25,0 0,8 -0,6 -0,5 0,7 0,3 277350,1
 l^2 lm (-l^2) (-lm)
lm m^2 (-lm) (-m^2)
K = AE/L
(-l^2) (-lm) l^2 lm
(-lm) (-m^2) lm m^2

0,69230769 0,46153846 0,692307692 0,461538462

0,46153846 0,30769231 0,461538462 0,307692308
K1 = 277350,098
0,69230769 0,46153846 0,692307692 0,461538462
0,46153846 0,30769231 0,461538462 0,307692308

1 0 -1 0
0 0 0 0
K2 = 166666,667
-1 0 1 0
0 0 0 0

0,692 -0,462 -0,692 0,462

-0,462 0,308 0,462 -0,308
K3 = 277350,098
-0,692 0,462 0,692 -0,462
0,462 -0,308 -0,462 0,308

0,6923 0,4615 -0,6923 -0,4615

0,4615 0,3077 -0,4615 -0,3077
K4 = 277350,098
-0,6923 -0,4615 0,6923 0,4615
-0,4615 -0,3077 0,4615 0,3077

0 0 0 0
0 1 0 -1
K5 = 250000
0 0 0 0
0 -1 0 1

0,6923 0,4615 -0,6923 -0,4615

0,4615 0,3077 -0,4615 -0,3077
K6 = 277350,098
-0,6923 -0,4615 0,6923 0,4615
-0,4615 -0,3077 0,4615 0,3077

0,6923 -0,4615 -0,6923 0,4615

-0,4615 0,3077 0,4615 -0,3077
K7 = 277350,098
-0,6923 0,4615 0,6923 -0,4615
0,4615 -0,3077 -0,4615 0,3077

K8 1 0 -1 0
= 166666,667 0 0 0 0
 -1 0 1 0
0 0 0 0

0,6923 -0,4615 -0,6923 0,4615

-0,4615 0,3077 0,4615 -0,3077
K9 = 277350,098
-0,6923 0,4615 0,6923 -0,4615
0,4615 -0,3077 -0,4615 0,3077

MATRIZ GENERAL

MATRIZ INVERSA (9X9)

MULTIPLICANDO POR LAS FUERZAS ACTUANTES DE LA ARMADURA

 DESPLAZAMIENTOS

desplazamientos
1 0,049504382 x
2 -0,07620958 y
3 0,0375 x
4 -0,080545104 y
5 0,047124056 x
6 -0,076545104 y
7 0,035119674 x
8 -0,075021608 y
9 0,063 x