Universidad Privada Domingo Savio

Facultad de Ingeniería
Carrera de Ingeniería Civil
Análisis Matricial de Estructuras

ELEMENTOS SOMETIDOS A FUERZAS

AXIALES POR EL MÉTODO DE LA RIGIDEZ

1 INTRODUCCIÓN

El desarrollo de los computadores ha estimulado enormemente la investigación en muchas

ramas de la ciencia permitiendo desarrollar procedimientos numéricos apropiados para el uso
de los mismos. En el campo del análisis de estructuras, el ordenador ha conducido al desarrollo
de métodos que utilizan las ideas del álgebra matricial.

La teoría matricial del análisis de estructuras aparece en la literatura técnica en la década de

los 50. Tras una confusión inicial en el mundo de la ingeniería estructural práctica que
oscureció en una primera etapa la relación existente entre el nuevo procedimiento y los
métodos estructurales clásicos, el desarrollo del nuevo método sufrió un impulso tal, que al
principio de la década de los 60 ya estaba perfectamente establecido.

Este impulso se debió a la confluencia de unas necesidades de cálculo, muchas veces tan
complejas, que los métodos clásicos resultaban claramente insuficientes con el desarrollo y
operatividad del ordenador.

El empleo de la notación matricial presenta dos ventajas en el cálculo de estructuras:

a) Permite desde el punto de vista teórico, utilizar métodos de cálculo de una forma más
compacta, precisa y al mismo tiempo completamente general. Los principios
fundamentales no se ven oscurecidos por las operaciones de cálculo o diferencias
geométricas en las tipologías estructurales analizadas.
b) Proporciona en la práctica, un sistema adecuado de análisis y determina las bases
idóneas para el desarrollo de programas de ordenador.

Introducción Página 1 de 11
 Universidad Privada Domingo Savio
Facultad de Ingeniería
Carrera de Ingeniería Civil
Análisis Matricial de Estructuras

2 CONCEPTOS

2.1 RIGIDEZ

𝐹
𝐾=
𝐷

Fuerza requerida en un elemento estructural para producir una unidad de desplazamiento o

deformación.

2.2 FLEXIBILIDAD

Es la inversa de la rigidez, es decir, es el desplazamiento requerido por unidad de fuerza.

1 𝐷
𝐹= =
𝐾 𝐹

2.3 SUPOSICIONES DE CÁLCULO

Para el cálculo de elementos sometidos a fuerzas axiales, se toman en cuenta las siguientes
suposiciones:

✓ Linealidad geométrica (desplazamientos pequeños).

✓ Linealidad material.
✓ Superposición.
✓ Material isotrópico.
✓ Sección prismática.

Introducción Página 2 de 11
 Universidad Privada Domingo Savio
Facultad de Ingeniería
Carrera de Ingeniería Civil
Análisis Matricial de Estructuras

Los resultados deberán cumplir:

✓ Condiciones de compatibilidad.
✓ Condiciones de equilibrio.
✓ Contornos.

2.4 MIEMBROS Y NODOS

Uno de los primeros pasos para aplicar el método de la rigidez consiste en identificar los
elementos o miembros de la estructura y sus nodos. Cada elemento se especifica por un
numero encerrado en un cuadrado, y para identificar los nodos se usa un numero dentro de
un círculo. También se identificarán los extremos inicial y final de cada elemento mediante
una flecha indicada a lo largo del elemento, con la punta de la fecha dirigida hacía el extremo
final.

2.5 GRADOS DE LIBERTAD

Los grados de libertad (GDL) Son el número de desplazamientos y rotaciones independientes

de un sistema; por ejemplo, para la cercha mostrada, el sistema cuenta con 8 GDL, 2 GDL
restringidos y 6 GDL no restringidos. En este caso, representan las incógnitas principales en el
método de la rigidez.

2.6 COORDENADAS GLOBALES Y LOCALES

Dado que las cargas y los desplazamientos son cantidades vectoriales, es necesario establecer
un sistema de coordenadas a fin de precisar el sentido correcto de la dirección. Se emplean
dos tipos diferentes de sistemas coordenados.

Introducción Página 3 de 11
 Universidad Privada Domingo Savio
Facultad de Ingeniería
Carrera de Ingeniería Civil
Análisis Matricial de Estructuras

2.6.1 Coordenadas globales

Con notación (x, y), el cual es único, y sirve para especificar el sentido de cada uno de los
componentes de la fuerza externa y el desplazamiento en los nodos.

2.6.2 Coordenadas locales

Para especificar el sentido de dirección de sus desplazamientos y las cagas internas en el

elemento. Este sistema se identificará con ejes (x',y') con el origen en el nodo inicial y
extendiéndose hacia el nodo final.

3 APLICACIÓN DEL MÉTODO DE RIGIDEZ

La aplicación de este método requiere subdividir a la estructura en una serie de elementos

finitos e identificar sus puntos extremos como nodos. Para el análisis de armaduras, los
elementos finitos se representan por cada uno de los miembros que forman la armadura y los
nodos representan los nudos de unión de los elementos. La información a obtener
corresponderá a cada punto de discretización.

Los pasos necesarios para resolver una estructura mediante los métodos matriciales son:

1) Definir la geometría de la estructura y las acciones, así como las condiciones de apoyo.
2) Identificar el número de movimientos incógnitas que determinan la deformación de la
estructura, teniendo en cuenta las correspondientes condiciones de compatibilidad en
los nudos.
3) Resolver las piezas individuales, en función de los movimientos de sus extremos,
satisfaciendo las condiciones de equilibrio y compatibilidad de las piezas.
4) Imponer las necesarias condiciones de equilibrio en los nudos.
5) Imponer las condiciones de apoyo de la estructura.
6) Determinar los movimientos incógnita resolviendo el sistema de ecuaciones resultante.
7) Determinar los esfuerzos y las reacciones en la estructura.

Introducción Página 4 de 11
 Universidad Privada Domingo Savio
Facultad de Ingeniería
Carrera de Ingeniería Civil
Análisis Matricial de Estructuras

4 MATRIZ DE RIGIDEZ DE UN ELEMENTO EN COORDENADAS LOCALES

Relaciona los desplazamientos incógnita de una estructura con las fuerzas exteriores
conocidas, lo cual permite encontrar las reacciones, esfuerzos internos y tensiones en
cualquier punto de la estructura.

La matriz de rigidez es una propiedad intrínseca del sistema estructural, no cambia en función
del estado de cargas o de condiciones de contorno a que se someta la estructura.

Sea la barra 1–2 sometida a las fuerzas 𝑃𝑥1 y 𝑃𝑥2 en los extremos, se tendrán desplazamientos
en estos nudos 𝛿𝑥1 y 𝛿𝑥2 , donde se cumplirá que:

Introducción Página 5 de 11
 Universidad Privada Domingo Savio
Facultad de Ingeniería
Carrera de Ingeniería Civil
Análisis Matricial de Estructuras

𝛿𝑡 = 𝛿𝑥1 − 𝛿𝑥2

Además, para que la barra esté en equilibrio:

𝑃𝑥1 = −𝑃𝑥2

En el caso de elementos sometidos a fuerzas axiales, solamente:

𝑃𝑦1 = 𝑃𝑦2 = 0 y 𝜃1 = 𝜃2 = 0

De la ley de Hooke se conoce:

𝑃
𝜎 = 𝜀𝐸 𝜎 = , esfuerzo
𝐴
𝛿
𝜀 = 𝐿 , deformación unitaria
𝑃 𝛿𝑇
= 𝐸 𝐸, método de elasticidad
𝐴 𝐿

𝐸𝐴
𝑃= 𝛿
𝐿 𝑇

Reemplazando 𝛿𝑡 en la ecuación anterior:

𝐸𝐴 𝐸𝐴
𝑃= 𝛿𝑥1 − 𝛿
𝐿 𝐿 𝑥2

Relacionando la ecuación anterior con 𝑃𝑥1 = −𝑃𝑥2 :

Introducción Página 6 de 11
 Universidad Privada Domingo Savio
Facultad de Ingeniería
Carrera de Ingeniería Civil
Análisis Matricial de Estructuras

𝐸𝐴 𝐸𝐴
𝑃𝑥1 = 𝛿𝑥1 − 𝛿 𝑃𝑦1 = 0
𝐿 𝐿 𝑥2

𝐸𝐴 𝐸𝐴 𝑃𝑦2 = 0
𝑃𝑥2 = − 𝛿𝑥1 + 𝛿
𝐿 𝐿 𝑥2

Las expresiones anteriores se pueden escribir de forma matricial, denotadas con un apostrofe
para indicar que son términos en coordenadas locales:

′ ′
𝑃𝑥1 𝐸𝐴⁄𝐿 0 − 𝐸𝐴⁄𝐿 0 𝛿𝑥1
′ ′
𝑃𝑥2 0] 𝛿𝑥2
=[ 0 0 0
′ ′
𝑃𝑦1 − 𝐸𝐴⁄𝐿 0 𝐸𝐴⁄𝐿 0 𝛿𝑦1
′ 0 0 0 ′
0 {𝛿𝑦2
{𝑃𝑦2 } }

Esta ecuación es más común denotarla como:

′ ′
𝐹𝑥1 𝐸𝐴⁄𝐿 0 − 𝐸𝐴⁄𝐿 0 𝛿𝑥1
′
𝐹𝑥2 0] 𝛿′𝑥2
=[ 0 0 0
′ ′
𝐹𝑦1 − 𝐸𝐴⁄𝐿 0 𝐸𝐴⁄𝐿 0 𝛿𝑦1
′ 0 0 0 ′
0 { 𝛿𝑦2
{𝐹𝑦2 } }

Vector de Matriz de rigidez Vector de

fuerza en para elementos desplazamiento
coordenadas biarticulados en en coordenadas
locales {𝐹 ′ } coordenadas locales locales {𝐷′ }
[𝐾 ′ ]

{𝐹 ′ } = [𝐾 ′ ]{𝐷 ′}

5 TRANSFORMACIÓN DEL DESPLAZAMIENTO Y FUERZA

Ahora se muestra la forma para transformar las fuerzas del elemento y los desplazamientos
definidos en coordenados locales a coordenadas globales mediante la matriz de
transformación T.

5.1 MATRIZ DE TRANSFORMACIÓN

Relaciona a los datos de rigidez en coordenadas locales a matrices de coordenadas globales.

Introducción Página 7 de 11
 Universidad Privada Domingo Savio
Facultad de Ingeniería
Carrera de Ingeniería Civil
Análisis Matricial de Estructuras

cos 𝛽 0 0 0
sin 𝛽 0 0 0
[𝑇 ] = [ ]
0 0 cos 𝛽 0
0 0 sin 𝛽 0

Una de las principales características de las matrices de transformaciones, es que estas son
ortogonales.

[𝑇]−1 [𝑇] = [𝐼 ] [𝑇]−1 , inversa de [𝑇]

[𝑇]𝑇 , transpuesta de [𝑇]
[𝑇]−1 = [𝑇]𝑇 [𝐼 ], matriz identidad

[𝑇]𝑇 [𝑇] = [𝐼 ]

Por tanto:

cos 𝛽 sin 𝛽 0 0
[𝑇]−1 0 0 0 0
= [𝑇 ] 𝑇 = [ ]
0 0 cos 𝛽 sin 𝛽
0 0 0 0

5.2 TRANSFORMACIÓN DEL DESPLAZAMIENTO

A partir del vector de los desplazamientos en coordenadas locales, multiplicando por la matriz
de transformación se tienen los desplazamientos en coordenadas locales:

′
𝛿𝑥1 cos 𝛽 0 0 0 𝛿𝑥1
′
𝛿𝑥2 sin 𝛽 0 0 0 𝛿𝑥2
𝛿𝑦1 = [ 0 0 cos 𝛽 0
] ′
𝛿𝑦1
′
{𝛿𝑦2 } 0 0 sin 𝛽 0
{𝛿𝑦2 }

Vector de Matriz de rotación o Vector de

desplazamiento transformación de desplazamiento
en coordenadas coordenadas locales a en coordenadas
globales {𝐷} coordenadas globales [𝑇] locales {𝐷′ }
{𝐷 } = [𝑇]{𝐷 ′}

Despejando el vector de desplazamientos en coordenadas locales {𝐷 ′}:

Introducción Página 8 de 11
 Universidad Privada Domingo Savio
Facultad de Ingeniería
Carrera de Ingeniería Civil
Análisis Matricial de Estructuras

[𝑇]−1 {𝐷 } = [𝑇]−1 [𝑇]{𝐷 ′} = [𝐼 ]{𝐷 ′}

{𝐷 ′} = [𝑇]−1 {𝐷 }

5.3 TRANSFORMACIÓN DE LA FUERZA

A partir del vector de las fuerzas en coordenadas locales, multiplicando por la matriz de
transformación se tienen las fuerzas en coordenadas locales:

′
𝐹𝑥1 cos 𝛽 0 0 0 𝐹𝑥1
′
𝐹𝑥2 sin 𝛽 0 0 0 𝐹𝑥2
𝐹𝑦1 = [ ] ′
0 0 cos 𝛽 0 𝐹𝑦1
{𝐹𝑦2 } 0 0 sin 𝛽 0 ′
{𝐹𝑦2 }

Vector de Matriz de rotación o Vector de

fuerzas en transformación de fuerzas en
coordenadas coordenadas locales a coordenadas
globales {𝐹 } coordenadas globales [𝑇] locales {𝐹 ′}

{𝐹 } = [𝑇]{𝐹 ′ }

Despejando el vector de fuerzas en coordenadas locales {𝐹 ′}:

[𝑇]−1 {𝐹 } = [𝑇]−1 [𝑇]{𝐹 ′ } = [𝐼 ]{𝐹 ′}

{𝐹 ′} = [𝑇]−1 {𝐹 }

6 MATRIZ DE RIGIDEZ DE UN ELEMENTO EN COORDENADAS

GLOBALES

Reemplazando [𝑇]−1 = [𝑇]𝑇 en las ecuaciones anteriores:

{𝐷 ′} = [𝑇]𝑇 {𝐷 }

{𝐹 ′} = [𝑇]𝑇 {𝐹 }

Reemplazando las anteriores matrices en la ecuación fundamental:

{𝐹 ′ } = [𝐾 ′ ]{𝐷′}

Introducción Página 9 de 11
 Universidad Privada Domingo Savio
Facultad de Ingeniería
Carrera de Ingeniería Civil
Análisis Matricial de Estructuras

[𝑇]𝑇 {𝐹 } = [𝐾 ′ ][𝑇]𝑇 {𝐷 }

[𝑇][𝑇]𝑇 {𝐹 } = [𝑇][𝐾 ′ ][𝑇]𝑇 {𝐷 }

[𝐼 ]{𝐹 } = [𝑇][𝐾 ′ ][𝑇]𝑇 {𝐷 }

{𝐹 } = [𝑇][𝐾 ′ ][𝑇]𝑇 {𝐷 }

Se deduce que la matriz de rigidez de un elemento en coordenadas globales es:

[𝐾 ] = [𝑇][𝐾 ′ ][𝑇]𝑇

Desarrollando la anterior ecuación, la matriz de rigidez en coordenadas globales es:

𝑐2 𝑐𝑠 −𝑐 2 −𝑐𝑠
𝐴𝐸 𝑐𝑠 𝑠2 −𝑐𝑠 −𝑠 2 ]
[𝐾 ] = [
𝐿 −𝑐 2 −𝑐𝑠 𝑐2 𝑐𝑠
−𝑐𝑠 −𝑠 2 𝑐𝑠 𝑠2

Donde: 𝑐 = cos 𝛽

𝑠 = sin 𝛽

7 MATRIZ DE RIGIDEZ ENSAMBLADA DE LA ESTRUCTURA

Una vez obtenidas las matrices de rigidez de cada elemento del sistema estructural, la matriz
ensamblada corresponde a la unión de todos los términos semejantes en cada nodo,
estableciendo de esta manera la compatibilidad de los elementos en sus extremos.

7.1 MATRIZ DE RIGIDEZ REDUCIDA

Al tener la matriz de rigidez de la estructura [𝐾 ], la matriz de rigidez reducida [𝐾 ]𝑟𝑒𝑑

representa a los datos después de la eliminación de las condiciones de desplazamiento nula,
es decir, la matriz de rigidez reducida se obtiene a partir de la eliminación de las condiciones
de apoyo que tiene la estructura.

8 VECTOR DE DESPLAZAMIENTO NODAL

Introducción Página 10 de 11
 Universidad Privada Domingo Savio
Facultad de Ingeniería
Carrera de Ingeniería Civil
Análisis Matricial de Estructuras

El vector de desplazamientos nodales se lo determina a partir de las matrices de fuerzas

externas y rigidez reducidas.

Obtenida la matriz de rigidez, operando en la ecuación:

{𝐹 } = [𝐾 ]{𝐷 }

[𝐾 ]−1 {𝐹 } = [𝐾 ]−1 [𝐾 ]{𝐷 }

[𝐾 ]−1 {𝐹 } = [𝐼 ]{𝐷 }

{𝐷 } = [𝐾 ]−1 {𝐹 }

Como se operan con matrices reducidas:

{𝐷 }𝑟𝑒𝑑 = [𝐾 ]−1
𝑟𝑒𝑑 {𝐹 }𝑟𝑒𝑑

9 EJEMPLO

Determinar:

✓ Vector de fuerzas externas ✓ Desplazamientos nodales.

aplicadas. ✓ Reacciones en los apoyos.
✓ Matriz de rigidez global de cada ✓ Esfuerzos en cada barra.
barra. ✓ Comprobar los resultados en
✓ Matriz de rigidez ensamblada. SAP2000.

Introducción Página 11 de 11