CAPITULO 8

MATRIZ DE RIGIDEZ Y DE FLEXIBILIDAD DE

UN ESTRUCTURA A PARTIR DEL CONCEPTO

RESUMEN
Se presenta el clculo de las matrices de rigidez K y de flexibilidad F de una estructura
usando el concepto. En los captulos posteriores a partir del concepto se obtendr en forma prctica
la matriz de rigidez.
Posteriormente se trabaja con la matriz de transformacin de coordenadas para encontrar las
matrices K y F . Finalmente se presenta un algoritmo orientado al uso del computador para
estructuras con elementos flexibles para encontrar las matrices indicadas en otro sistema de
coordenadas.

8. 1 MATRIZ DE RIGIDEZ DE UNA ESTRUCTURA K

8.1.1 Definicin
En el captulo 6 se indic que un trmino cualquiera de la matriz de rigidez de una estructura
k ij , es el valor de la carga generalizada Qi correspondiente a la deformada elemental q j = 1 y
dems nulas. Por consiguiente si se desea calcular, por ejemplo, los elementos de la primera
columna de la matriz de rigidez de una estructura se deber calcular el vector de cargas
generalizadas que corresponde al estado de desplazamiento elemental q1 = 1 y qi = 0, i 1 . De
igual forma se proceder con las dems columnas de K .

8.1.2

Procedimiento de clculo

El procedimiento de clculo para hallar la matriz de rigidez de una estructura K a partir del
concepto es el siguiente:

Roberto Aguiar Falcon

CEINCI-ESPE

234

1) Construir la deformada elemental cuya columna se desea calcular.

2) Encontrar las deformaciones p en cada uno de los elementos asociados a la deformada
elemental. Es un problema de geometra.
3) Transformar las deformaciones p de cada elemento en cargas internas P por medio de la
matriz de rigidez del elemento k . La ecuacin matricial que se utiliza es: P = k p .

4) Usando la esttica se realiza el equilibrio de cada uno de los elementos que conforman la
estructura.
5) Encontrar el equilibrio de cada una de las juntas de la estructura.
6) En el paso anterior se obtienen las cargas que actan sobre la estructura y el vector de
cargas generalizadas que son los elementos de la matriz de rigidez de la estructura.

8.1.3

Primera forma de clculo numrico

Por didctica nicamente se denomina primera forma de clculo de la matriz de rigidez de

una estructura a aquella en que se utiliza como sistema de coordenadas del elemento el indicado en
la figura 8.1

Figura 8.1 Sistema P p para la primera forma de clculo.

Se recuerda que las deformaciones p se obtienen con las siguientes ecuaciones:

p1 = 1

v 2 v1
L

p2 = 2

v 2 v1
L

p3 = u 2 u1

La matriz de rigidez de un elemento para el sistema de coordenadas indicado es:

4 EI
L
2 EI
k=
L
0

2 EI
L
4 EI
L
0

0
0
EA
L

EJEMPLO N.- 1

Encontrar la primera columna de la matriz de rigidez de la estructura indicada en la figura

8.2.1 si todos los elementos tienen la misma seccin transversal y la misma longitud. Encontrar
aplicando el concepto y no considerar el efecto de corte.

ANLISIS MATRICIAL DE ESTRUCTURAS

Figura 8.2.1

235

Figura 8.2.2

Figura 8.2.3

SOLUCIN

Por ser todos sus elementos flexibles se tiene seis grados de libertad. En consecuencia el
sistema de coordenadas de la estructura Q q es el indicado en la figura 8.2.2. Para una mejor
compresin en la figura 8.2.3 se indica la numeracin de los elementos.
En la figura 8.2.4 se indica el sistema de coordenadas de miembro P p de cada uno de
los elementos de la estructura.

Figura 8.2.4 Sistema P p de la estructura.

Para diferenciar las deformaciones y cargas internas se escribe entre parntesis y como
subndice el nmero del elemento al cual corresponde.
De acuerdo al procedimiento indicado en el apartado anterior para encontrar los elementos de
la primera columna de la matriz de rigidez de la estructura se procede de la siguiente manera:
1)

Deformada elemental

q1 = 1 y qi = 0, i 1 .

Figura 8.2.5 Diagrama elemental q1

Roberto Aguiar Falcon

CEINCI-ESPE

236
2) Clculo de las deformaciones p .

1
L
1
=
L
=0

p1(1) =

p1( 2 ) = 0

p1(3) = 0

p 2(1)

p 2( 2 ) = 0

p 2(3) = 0

p3( 2) = 1

p3(3) = 0

p3(1)

Para encontrar las deformaciones del elemento 1, se procedi de la siguiente manera:

u1 = 0

u2 = 0

v1 = 0

v 2 = 1

1 = 0

2 = 0

v 2 v1
1 0 1
= 0
=
L
L
L
v v
1 0 1
= 2 2 1 = 0
=
L
L
L
= u 2 u1 = 0 0 = 0

p1(1) = 1
p 2(1)
p3(1)

Se puede obtener stos mismos valores directamente sin necesidad de aplicar las
ecuaciones si se recuerda que las deformaciones p1 y p 2 son los ngulos comprendidos entre la
cuerda y la tangente de los nudos inicial y final respectivamente como se indic en el captulo 6. Para
el elemento 1 en la figura 8.2.6 se ha separado la deformada elemental correspondiente a la columna
izquierda en ella se ha unido la cuerda entre B y B, luego se han trazado las tangentes en el nudo
inicial y final se observa que el ngulo entre la cuerda y la tangente es antihorario luego las
deformaciones p1 y p 2 son positivas.

Figura 8.2.6 Regla para encontrar las deformaciones del elemento 1.

En la figura 8.2.6 se aprecia que el ngulo p1 es igual al ngulo p 2 por ser alternos e
internos y que p1 es igual al cateto opuesto sobre el cateto adyacente. Por lo tanto se tiene:

p1(1) =

1
L

p 2(1) =

1
L

3) Se obtienen las cargas internas P en cada uno de los elementos del prtico plano. Como
todos los elementos tienen la misma seccin transversal y la misma longitud, la matriz de
rigidez de cada uno de ellos es la misma y vale:

ANLISIS MATRICIAL DE ESTRUCTURAS

k (1) = k ( 2 ) = k (3)

237

4 EI
L
2 EI
L

2 EI
L
4 EI
L

0
0
EA
L

Las cargas internas del elemento 1 se obtienen del producto matricial P

Donde P

(1)

(1)

= k (1) p (1) .

es de la siguiente forma:

(1)
1

P2(1)
(1)
3

4 EI
L
2 EI
=
L

2 EI
L
4 EI
L

0
0
EA
L

1
L
1
L
0

De donde:

P1(1) =

6 EI
L2

P2(1) =

6 EI
L2

P3(1) = 0

Lo propio se realiza con los elementos 2 y 3 es decir se multiplica la matriz de rigidez del
respectivo elemento por su vector de deformaciones, los resultados que se obtienen, son:

P1( 2 ) = 0

P2( 2 ) = 0

P3( 2 ) =

P1(3) = 0

P2(3) = 0

P3(3) = 0

EA
L

4) Se encuentra el equilibrio de cada elemento.

En la figura 8.2.7 se indica con lnea continua las cargas internas P obtenidas en el paso
anterior y con lnea entrecortada las diferentes fuerzas que equilibran cada uno de los elementos.

Figura 8.2.7 Equilibrio de elementos.

El cortante en la columna se obtiene sumando los momentos y dividiendo para la longitud.

Roberto Aguiar Falcon

CEINCI-ESPE

238
5) Se realiza el equilibrio de juntas.

En las juntas de la estructura las fuerzas y momentos internos de cada elemento actan con
sentido contrario, stas se han representado con lnea continua y las cargas exteriores que equilibran
cada uno de los nudos se presentan con lnea entrecortada en la figura 8.2.8.

Figura 8.2.8 Equilibrio de juntas o nudos.

Ntese que para equilibrar la fuerza horizontal del nudo B se ha sumado el cortante
proveniente de la columna y la fuerza axial que viene de la viga.
6) Finalmente se determinan las cargas exteriores y el vector de cargas generalizadas.
En la figura 8.2.9 se presentan las fuerzas y momentos exteriores que se deben aplicar a la
estructura para que produzcan el diagrama elemental de la figura 8.2.5.

Figura 8.2.9 Valores de la primera columna de K

Por consiguiente el vector de cargas generalizadas Q resulta:

12 EI EA
+
L
L3
K 11
0
K 21
6 EI
K 31
=
Q=
L2
K 41
EA

K 51
L
K 61
0
0

ANLISIS MATRICIAL DE ESTRUCTURAS

239

Por definicin los elementos de Q son los trminos de la primera columna de la matriz de
rigidez de la estructura. Se deja al estudiante el obtener las dems columnas de K aplicando el
concepto. El resultado total es el siguiente:

12 EI EA
+
L
L3
0

K=

6 EI
L2
EA

0
0

8.1.4

12 EI EA
+
L
L3
6 EI
L2

6 EI
L2
6 EI
L2
8EI
L

12 EI
L3
6 EI
L2

6 EI
L2
2 EI
L

EA
L
0
0

12 EI EA
+
L
L3
0
6 EI
L2

0
12 EI
L3
6 EI
2
L

0
12 EI EA
+
L
L3
6 EI
2
L

0
6 EI
L2
2 EI
L
6 EI
L2
6 EI
2
L
8 EI
L

Segunda forma de clculo numrico

Se puede considerar otro sistema de coordenadas del elemento para encontrar la matriz de
rigidez de la estructura, por ejemplo el presentado en la figura 8.3

Figura 8.3 Otro sistema de coordenadas del elemento P p

EJEMPLO N.- 2

Calcular la primera columna de la matriz de rigidez de la estructura mostrada en la figura

8.2.1, considerando el sistema de coordenadas del elemento el indicado en la figura 8.3. En las
figuras 8.2.2 y 8.2.3 se indican el sistema de coordenadas de la estructura Q q y la numeracin de
los elementos respectivamente.

SOLUCIN

En la figura 8.4.1 se presenta las coordenadas de cada uno de los elementos que forman la
estructura.

Roberto Aguiar Falcon

CEINCI-ESPE

240

Figura 8.4.1 Sistema de coordenadas de los elementos del ejemplo 2.

Se resume a continuacin el procedimiento de clculo.
1) Deformada elemental q1

Figura 8.4.2 Deformada elemental q1

Las frmulas con las cuales se obtienen las deformaciones de un elemento se obtuvieron en
el captulo 6, en el apartado 6.3.4, stas son:

p1 = u 2 u1
p 2 = v 2 v1 L 1
p3 = 2 1
2) Deformaciones de los elementos
Las deformaciones en cada uno de los elementos son las siguientes:

3)

p1(1) = 0

p1( 2 ) = 1

p1(3) = 0

p 2(1) = 1

p 2( 2 ) = 0

p 2(3) = 0

p3(1) = 0

p3( 2 ) = 1

p3(3) = 0

Cargas Internas

La matriz de rigidez de elemento para el sistema de coordenadas P p de la figura 8.3 fue

deducida en el captulo anterior en el apartado 7.3.3 y es la siguiente:

ANLISIS MATRICIAL DE ESTRUCTURAS

EA
L
k=

0
0

241

12 EI
L3
6 EI
2
L

6 EI
L2
4 EI
L

Por los datos del ejemplo la matriz de rigidez es igual para todos los elementos. Al realizar el
producto matricial P = k p en cada uno de los elementos de la estructura se obtiene:

P1(1) = 0

4)

P1( 2 ) =

12 EI
L3
6 EI
= 2
L

EA
L

P1(3) = 0

P2(1) =

P2( 2 ) = 0

P2( 3) = 0

P3(1)

P3( 2 ) = 0

P3( 3) = 0

Equilibrio de elementos

En la figura 8.4.3 se presenta con lnea continua las cargas P obtenidas en cada uno de los
elementos y con lnea entrecortada las diferentes fuerzas que equilibran los elementos.

Figura 8.4.3 Equilibrio de elementos

Las figuras 8.4.3 y 8.2.7 son iguales. En consecuencia los siguientes pasos que faltan para
hallar los elementos de la primera columna de la matriz de rigidez de la estructura son los dados en el
ejemplo anterior. De idntica forma se obtendrn las dems columnas de la matriz de rigidez.
El haber obtenido la misma matriz de rigidez de la estructura, utilizando diferentes sistemas
de coordenadas de elemento, no es un hecho casual. Esto es debido a que la inmovilizacin de los
desplazamientos como cuerpo rgido, que se estudi en el captulo 6, utilizando diferentes tipos de
vnculos es un artificio. Si se utiliza otro sistema de coordenadas de miembro se obtendr la misma
matriz de rigidez K .
Finalmente, se puede resolver la estructura utilizando sistemas de coordenadas de elemento
P p , diferentes para cada uno de los elementos, por ejemplo el que se indica en la figura 8.5.

Roberto Aguiar Falcon

CEINCI-ESPE

242

Figura 8.5 Nuevo sistema de coordenadas de los elementos para la estructura de los ejemplos 1 y 2.

8. 2 MATRIZ DE FLEXIBILIDAD DE UNA ESTRUCTURA F

8.2.1 Definicin
Se conoce que q = F Q . Por lo tanto un elemento cualquiera de la matriz de flexibilidad de
una estructura Fij ser el valor del desplazamiento o giro q i correspondiente al estado de cargas

Q j = 1 y las dems nulas.

La matriz de flexibilidad F transforma las cargas generalizadas Q en coordenadas
generalizadas q .

8.2.2

pasos:

Procedimiento de clculo
Para encontrar la matriz de flexibilidad F a partir de su definicin, se realizan los siguientes

1) Para el estado de carga elemental Qi = 1 y las dems nulas, se debe hallar las cargas
internas P que actan en cada uno de los elementos. Esto es un procedimiento de esttica.
2) Conocido el vector P se encuentran las deformaciones p , por medio de la siguiente
ecuacin: p = f P .

3) A partir de las deformaciones internas en cada uno de los elementos de la estructura se

determina el vector de coordenadas generalizadas q . Este es un problema de geometra. Al
encontrar q se tienen ya los elementos de la columna de F que se est calculando.

EJEMPLO N.- 3

En la figura 8.6 se presenta un prtico plano con todos sus elementos completamente
flexibles, se pide calcular los elementos de la primera columna de la matriz de flexibilidad.

ANLISIS MATRICIAL DE ESTRUCTURAS

243

Figura 8.6 Estructura de ejemplo 3.

SOLUCIN

Por facilidad se ha considerado que los dos elementos tienen la misma seccin transversal y
la misma longitud. En consecuencia la matriz de flexibilidad de sus elementos es igual.
En la figura 8.7.1 se indican los 6 grados de libertad que tiene el prtico y en la figura 8.7.2 se
indica la numeracin de los elementos.

Figura 8.7.1 Sistema Q q

Figura 8.7.2

De acuerdo al procedimiento de clculo indicado en el apartado 8.2.2 se tiene:

1)

Q1 = 1 y Qi = 0 i 1

Figura 8.7.3 Momento unitario en A. Figura 8.7.4 Equilibrio de fuerzas externas

Roberto Aguiar Falcon

CEINCI-ESPE

244

Debido al momento unitario que acta en el nudo A de la estructura, se generan reacciones

en los vnculos A y C que valen
muestran stas reacciones.

1
con lo cual la estructura est en equilibrio. En la figura 8.7.4 se
L

Por efecto del sistema de cargas en cada uno de los elementos se tienen fuerzas y
momentos internos los mismos que se indican en la figura 8.7.5. Con lnea continua se indican las
acciones que vienen de las reacciones de los nudos y con lnea entrecortada las acciones que
conducen al equilibrio en los elementos.
Para el equilibrio de las juntas se colocan en primer lugar las fuerzas externas que actan en
la junta las mismas que estn indicadas en la figura 8.7.4. Luego se equilibra el nudo con fuerzas y
momentos. Estas fuerzas son las que pasan a los elementos con sentido contrario y son las que se
han indicado en la figura 8.7.5.

Figura 8.7.5 Equilibrio de los elementos.

Se deja al lector la verificacin del equilibrio de juntas y de elementos. Por otra parte se
trabaja con el sistema de coordinas de miembro de la figura 8.1, en la cual se tiene que P1 es el
momento en el nudo inicial y es positivo si es antihorario, P2 es el momento en el nudo final y es
positivo si es antihorario y P3 es la fuerza axial en el nudo final y es positiva si produce traccin en el
elemento. Con stas indicaciones las cargas internas en los elementos que se obtienen de la figura
8.7.5 son las siguientes:

P1(1) = 1
(1)
2

P1( 2 ) = 1

= 1

P3(1) =

1
L

( 2)
2

=0

(1)

1
= 1

P3( 2 ) = 0

1
L

p= f P

f (1) = f ( 2 )

L
6 EI
L
3EI
0

1
= 0
0

2) Clculo de las deformaciones de los elementos.

L
3EI
L
=
6 EI

( 2)

0
0
L
EA

ANLISIS MATRICIAL DE ESTRUCTURAS

p (1) = f (1) P (1)

245

L
3EI
L
=
6 EI

p ( 2) = f ( 2) P ( 2)

L
3EI
L
=
6 EI
0

L
6 EI
L
3EI

0
0
L
EA

L
6 EI
L
3EI

0
0
L
EA

L
2 EI
1
L
1 =
2 EI
1
1

L
EA

L
3EI
1
L
0 =
6 EI
0
0

En consecuencia se tiene:

L
2 EI
L
=
2 EI
1
=
EA

L
3EI
L
=
6 EI

p1(1) =

p1( 2) =

p 2(1)

p 2( 2)

p3(1)

p3( 2 ) = 0

3) Clculo de los desplazamientos y giros q

Esta es la parte ms difcil del clculo de la matriz de flexibilidad de una estructura. Se quiere
hallar los desplazamientos y giros que experimentan las juntas de tal forma que los elementos tengan
las deformaciones p indicadas. Para facilitar la explicacin, se denomina:

L
2 EI

1
EA

Con sta nomenclatura las deformaciones en los elementos valen:

p1(1) =

2
p1( 2 ) =
3

p 2(1) =

p 2( 2 ) =

(1)
3

( 2)
3

=0

En la figura 8.7.6 se indican stas deformaciones y en la misma se aprecia que se est

violando dos reglas bsicas en el nudo B que son:
i)

El nudo B es discontinuo ; se aprecia que el nudo B por parte de la

columna ha bajado.

Roberto Aguiar Falcon

CEINCI-ESPE

246

ii)

El nudo B antes de deformarse meda 90 grado, ahora no.

Figura 8.7.6 Deformaciones de los elementos obtenidos.

Para solucionar stos dos problemas se deben dar ciertos desplazamientos y giros a la
estructura. En primer lugar para corregir la anomala de la discontinuidad en el nudo B, al elemento
BC se desplaza verticalmente una magnitud igual a como se indica en la figura 8.7.7.

Figura 8.7.7 Solucin de la discontinuidad en el nudo B pero no en el ngulo.

Si bien se ha solucionado un problema, ahora se ha creado otro que es la posicin del nudo
C ya que debido al tipo de vnculo existente en esa junta el nudo C no debe bajar. Esto se va a
resolver posteriormente, se concentra la atencin en resolver el problema de que la junta B antes de
deformarse meda 90 grados y despus de deformarse debe medir 90 grados.
En la figura 8.7.7 se aprecia que el ngulo B vale
B debe valer 90 grados luego el elemento BC se rota
en la figura 8.7.8.

2
3

5
3

+ 90 + = 90 + . Pero el ngulo

5
haciendo centro en B esto se presenta
3

ANLISIS MATRICIAL DE ESTRUCTURAS

247

Luego de la rotacin el punto C se ha desplazado verticalmente

5
L hasta llegar a C.
3

Figura 8.7.8 Solucin de la rotacin del nudo B.

La posicin final del nudo C no puede ser C como est indicado en la figura 8.7.8 debe estar
en cualquier parte de la recta BC, para lograr ste objetivo se debe rotar el nudo A un ngulo .

5
3
L

+ L

5
+
L 3

Esta rotacin se indica en la figura 8.7.9, se observa que la posicin final del nudo C es
C. El nudo C se ha desplazado horizontalmente L igual corrimiento experimenta el nudo B como
se aprecia en la figura mencionada.

Figura 8.7.9 Posicin final de la estructura.

En resumen, los desplazamientos y giros de los nudos son los siguientes:

Roberto Aguiar Falcon

CEINCI-ESPE

248

q1 = +
q 2 = L
q3 =
q 4 = +
q 5 = L
5
q6 = +
3 3

Rotacin del nudo A.

Desplazamiento horizontal del nudo B.
Desplazamiento vertical del nudo B.
Rotacin del nudo B.
Desplazamiento horizontal del nudo C.
Rotacin del nudo C.

Al reemplazar , , , se encuentra:

5
8
4L
1
+ = + =
+
L 3
L 3EI EAL
3
5
5
1
5 L2
= + L = L =

L 3
EA 6 EI
3
1
= =
EA
5
2
L
1
= + + = + =
+
L 3
L 3
EAL 3EI
2
1
5L
= L =

EA 6 EI
5

L
1
= 2 + + = =

L 3
L 3 EAL 6 EI

q1 = +
q2
q3
q4
q5
q6

Por definicin se tiene:

4L
1
+
3EI EAL
1
5 L2
= F21 =

EA 6 EI
1
= F31 =
EA
1
L
= F41 =
+
EAL 3EI
1
5 L2
= F51 =

EA 6 EI
1
L
= F61 =

EAL 6 EI

q1 = F11 =
q2
q3
q4
q5
q6

Las restantes columnas de la matriz de flexibilidad F se obtienen en forma similar. Con el

propsito de que el estudiante calcule cualquiera de las columnas de la matriz de flexibilidad a partir
de su definicin, se indica el resultado completo pero nicamente se presenta la matriz triangular
superior debido a que la matriz es simtrica.

ANLISIS MATRICIAL DE ESTRUCTURAS

4L
1
+
3EI EAL

1 5 L2

EA EI
2 L3
L
+
3EI EA

EA
L
EA
L
EA

F=

8.2.3

249

1
L
1
5 L2
1
L
+

EAL 3EI
EA 6 EI
EAL 6 EI
2
3
L
1
2L
L
L2
1

3EI EA
3EI EA
6 EI EA
1
L
1

EA
EA
EA
2
L
1
L
1
L
1
+

+
3EI LEA
3EI EA
6 EI LEA
3
2L
2L
L2
1
+

3EI EA
6 EI EA
L
1
+
3EI LEA

Principio de superposicin

Como se indic la parte ms difcil del clculo de la matriz de flexibilidad de una estructura a
partir de su definicin es calcular los desplazamientos y giros q luego que se han obtenido las
deformaciones p . Si se emplea el principio de superposicin se puede hacer esto de una forma
sencilla como se ilustra con el ejemplo 4.

EJEMPLO N.- 4

Encontrar la primera columna de la matriz de flexibilidad del ejemplo 3, aplicando el principio

de superposicin para calcular el vector de coordenadas q .

SOLUCIN
Para el ejemplo que se est analizando se tiene:

p1(1) =

2
p1( 2 ) =
3

p 2(1) =

p 2( 2 ) =

(1)
3

( 2)
3

=0

En el clculo del vector de coordenadas generalizadas q se consideran dos etapas, a saber:

Etapa 1.- Actan slo las deformaciones del elemento uno, por lo tanto se considera al elemento dos
como totalmente rgido.
Etapa 2.- El elemento uno es totalmente rgido y slo hay deformaciones en el miembro dos.

Roberto Aguiar Falcon

CEINCI-ESPE

250

A los desplazamientos y giros de la etapa uno se los denomina q

giros de la etapa dos q

( 2)

(1)

y a los desplazamientos y

. En consecuencia:

q = q (1) + q ( 2 )
Etapa 1.- En la figura 8.8.1 se presenta sta fase del clculo que consiste en tener
deformaciones nicamente en el elemento vertical.

Figura 8.8.1 nicamente se deforma elemento AB

Al igual que en el ejemplo anterior, en la figura 8.8.1 hay que rectificar dos errores que se
tienen y son los siguientes:
i)

El nudo B es discontinuo.

ii)

El ngulo del nudo B antes de la deformacin es diferente al obtenido despus de la

deformacin.

Para corregir el primer error se desplaza el elemento BC una cantidad como lo indica la
figura 8.8.2 y para solucionar el segundo error haciendo centro en B se rota el elemento BC una
magnitud igual a como lo ilustra la figura 8.8.3. Ntese que el punto C se ha desplazado hasta C
una cantidad L .

Figura 8.8.2

Figura 8.8.3

Pero la posicin final del nudo C tiene que estar a lo largo de BC ya que no puede
desplazarse verticalmente para lograr esto haciendo centro en el nudo A se rota la estructura un
ngulo 1 .

ANLISIS MATRICIAL DE ESTRUCTURAS

1 =

251

+ L
L

Por lo tanto la ubicacin final del nudo C al terminar sta etapa es C como se indica en la
figura 8.8.4.

Figura 8.8.4 Fin de etapa uno

Por lo tanto los desplazamientos y giros encontrados en la etapa uno son:

q1(1) = + 1 = 2 +
q

(1)
2

L
= 1 L = L

q3(1) =
q 4(1) = + 1 =
q

(1)
5

L
= 1 L = L

q 6(1) = + 1 =

Etapa 2.- En sta etapa actan las deformaciones en el elemento 2 ya que el miembro 1 se
considera como elemento rgido, estas deformaciones iniciales se indican en la
figura 8.9.1

Figura 8.9.1 Slo se deforma elemento BC.

Roberto Aguiar Falcon

CEINCI-ESPE

252

El ngulo final del nudo B en la figura 8.9.1 se aprecia que es 90 +

o

2
. Esto no puede ser.
3

Se deja al estudiante que justifique las figuras 8.9.2 y 8.9.3 con las cuales se corrige sta anomala.

Figura 8.9.2

Figura 8.9.3

Luego de la geometra realizada se tiene:

2
q1( 2 ) = 2 =
3
2
q 2( 2 ) = 2 L = L
3
( 2)
q3 = 0
2
2
2
q 4( 2 ) = + 2 =
3
3
3
2
q5( 2 ) = 2 L = L
3
2

q 6( 2 ) = + 2 =
3 3
3
Al sumar los desplazamientos y giros obtenidos en las dos etapas se tiene:

q = q (1) + q ( 2 )
2
8

+ = +
L 3
L
3
2
5
q 2 = L L = L
3
3
q 3 = + 0 =
q1 = 2 +

q4 =

2
+
L 3

2
5
q 5 = L L = L
3
3
q6 =

ANLISIS MATRICIAL DE ESTRUCTURAS

Al sustituir

253

L
1
y =
se obtienen los valores de la primera columna de la matriz
2 EI
EA

de flexibilidad anotados en el ejemplo 3.

8. 3 TRANSFORMACIN DE COORDENADAS DE UNA ESTRUCTURA

8.3.1

Clculo de la matriz de rigidez y de flexibilidad

Conocida la matriz de rigidez K para un determinado sistema de coordenadas Q q de una

estructura, se desea ahora calcular la matriz de rigidez K

en otro sistema de coordenadas

Q q . Este clculo se lo va a realizar utilizando la matriz de transformacin de coordenadas.

EJEMPLO N.- 5

El prtico plano de la figura 8.10.1 est compuesto por una columna de altura H y una viga de
longitud L, los dos elementos se consideran axialmente rgidos y el sistema de coordenadas Q q
es el indicado en la figura 8.10.2. La matriz de rigidez K asociado a ste sistema de coordenadas
es la siguiente:

12 EI o
H3
6 EI o
K=
H2
0

6 EI o
H2
4 EI o 4 EI 1
+
H
L
2 EI 1
L

Se desea encontrar la matriz de rigidez K

0
2 EI 1
L
4 EI 1
L

para el sistema de coordenadas Q q de la

figura 8.10.3 por medio de la matriz de transformacin de coordenadas T .

Figura 8.10.1

Figura 8.10.2 Sistema Q q

Figura 8.10.3 Sistema Q q

Roberto Aguiar Falcon

CEINCI-ESPE

254

SOLUCIN
Se establece una relacin entre los dos sistemas de coordenadas de la siguiente manera:

q = T q
La relacin as definida fue estudiada en el captulo 5 por consiguiente el estudiante deber
dibujar cada una de las deformadas elementales y encontrar la siguiente matriz T .

T= 1
0

0
1

0
0

En el apartado 7.3.2 del captulo anterior se demostr que la matriz de rigidez de un elemento
en otro sistema de coordenadas se obtiene con la siguiente ecuacin:

k = T t k T

( 8.1 )

La ecuacin ( 8.1 ) que fue deducida para un elemento es aplicable a una estructura. Por lo
tanto para hallar la matriz de rigidez de la estructura K
ecuacin.

(mayscula) se utilizar la siguiente

K =Tt K T
Para el ejemplo se tiene:

K =

0
1

0
0

1
0

12 EI o
H3
6 EI o
H2

6 EI o
H2
4 EI o 4 EI 1
+
H
L
2 EI 1
L

0
Por lo tanto la matriz de rigidez K

4 EI o 4 EI 1
+
H
L
2 EI 1
K =
L
6 EI o
H2

0
0
1
0

2 EI 1
L
4 EI 1
L

0
0
1

1
0
0

para el sistema de coordenadas de la figura 8.10.3 es:

2 EI 1
L
4 EI 1
L
0

6 EI o
H2
0
12 EI o
H3

Si se desea calcular la matriz de flexibilidad F para el sistema de coordenadas de la figura

8.10.3, conocida la matriz de flexibilidad F para el sistema Q q de la figura 8.10.2, hay que
hacerlo por medio de la matriz T1 estudiada en el captulo 5.

Q = T1 Q

ANLISIS MATRICIAL DE ESTRUCTURAS

255

La matriz de flexibilidad buscada se obtiene del siguiente triple producto matricial.

F = T1t F T1
Se deja al lector la demostracin de las formulas con las cuales se obtiene K
a la teora presentada en el captulo 5.

8.3.2

y F en base

Regla prctica

Se puede obtener directamente la matriz de rigidez y de flexibilidad de una estructura cuyos

elementos son totalmente flexibles, cuando se cambia el sistema de coordenadas generalizadas, por
consiguiente no es necesario calcular las matrices T y T1 respectivamente.
En efecto, para hallar K

para el sistema de coordenadas Q q

de la figura 8.10.3 a

partir de la matriz K calculado para las coordenadas de la figura 8.10.2 se han de intercambiar las
filas y columnas de sta matriz de acuerdo al cambio de numeracin del nuevo sistema de
coordenadas.
saber:

La forma de pasar las coordenadas de la figura 8.10.2 a la 8.10.3 se presenta en dos fases a

1) Se cambian la numeracin de los dgitos uno y tres en el sistema Q q esto se presenta en

la figura 8.10.4. Por lo tanto ahora q1 es la rotacin en el rodillo C y q 3 es el desplazamiento
horizontal de la junta B.

Figura 8.10.4 Sistema de coordenadas generalizadas de fase uno.

La matriz de rigidez para este nuevo sistema de coordenadas se obtiene intercambiando las
columnas y las filas uno y tres de la matriz de rigidez, con lo que se halla:

4 EI 1
L
2 EI 1
K de la fase uno =
L
0

2 EI 1
L
4 EI o 4 EI 1
+
H
L
6 EI o
H2

0
6 EI o
H2
12 EI o
H3

Roberto Aguiar Falcon

CEINCI-ESPE

256

2) Por ltimo para pasar de las coordenadas de la figura 8.10.4 a las coordenadas de la figura
8.10.3 se intercambian los dgitos uno y dos. Por lo tanto en la matriz de rigidez de la fase
uno se intercambian los elementos de la fila uno a la fila dos y luego los elementos de la
columna uno a la columna dos,
presentada en el ejemplo anterior.

obteniendo de sta manera la matriz K

que ya fue

Un procedimiento similar se sigue para calcular la matriz de flexibilidad de una estructura en

un nuevo sistema de coordenadas. El procedimiento de clculo presentado es muy fcil programarlo.

8. 4 EJERCICIOS RESUELTOS
En los ejercicios que se van a resolver en ste captulo al igual que en el prximo, el lector
deber justificar cada uno de los pasos dados, toda vez que ya se han indicado la teora respectiva.

EJEMPLO N.- 6

Obtener directamente a partir de su definicin los elementos de la segunda columna de la

matriz de rigidez de la estructura mostrada en la figura 8.11.1, las vigas son totalmente rgidas y las
columnas son axialmente rgidas. Todas las columnas tienen la misma seccin transversal y longitud.

Figura 8.11.1 Estructura de ejemplo 6.

Figura 8.11.2 Sistema Q q

Figura 8.11.3 Sistema P p

ANLISIS MATRICIAL DE ESTRUCTURAS

257

SOLUCIN

La estructura tiene dos grados de libertad que son los desplazamientos horizontales de cada
uno de los pisos. Estos grados de libertad se indican en la figura 8.11.2. Por ser los elementos
verticales axialmente rgidos no existe deformacin axial luego hay dos coordenadas en el sistema
P p de cada elemento. En la figura 8.11.3 se muestra el sistema de coordenadas de los
elementos de la estructura analizada y la numeracin de los elementos. Ntese que no existe sistema
P p en los elementos horizontales esto se debe a que son elementos completamente rgidos.
De acuerdo al procedimiento de clculo indicado en el apartado 8.1.2 se construye la
deformada elemental q 2 ya que se va a encontrar los elementos de la segunda columna de la matriz
de rigidez.

q 2 = 1 y qi = 0 i 2

Figura 8.11.4 Deformada elemental q 2

1
L
1
=
L

1
L
1
=
L

p1(1) = 0

p1( 2 ) =

p1(3) = 0

p1( 4 ) =

p 2(1) = 0

p 2( 2 )

p 2(3) = 0

p 2( 4 )

Cargas Internas P = k p . nicamente los elementos dos y cuatro tienen deformaciones.

P ( 2) = P ( 4)

4 EI o
L
=
2 EI o
L

2 EI o
L
4 EI o
L

6 EI 0
1
L
L2
=
1
6 EI 0
L
L2

Por lo tanto:

6 EI o
L2
6 EI
= 2o
L

6 EI o
L2
6 EI
= 2o
L

P1(1) = 0

P1( 2 ) =

P1( 3) = 0

P1( 4 ) =

P2(1) = 0

P2( 2 )

P2( 3) = 0

P2( 4 )

Roberto Aguiar Falcon

CEINCI-ESPE

258
Equilibrio de elementos

Figura 8.11.5 Equilibrio de los elementos

Equilibrio de Juntas

Figura 8.11.6 Junta B

o

Figura 8.11.7 Junta C

Junta B

12 EI o
N5
L3
u + u'
FY = 0 = N 1 N 2
L
6 EI
M = 0 = 2o +u
L
FX = 0 =

Junta C

12 EI o
L3
u '' + u '''
FY = 0 = N 2
L
6
EI
M = 0 = u '' + 2 o
L
FX = 0 = N 6 +

ANLISIS MATRICIAL DE ESTRUCTURAS

Figura 8.11.8 Junta D

o

259

Figura 8.11.9 Junta E

Junta D

FX = 0 = K 22 + N 6

12 EI o
L3

u '' + u '''
FY = 0 = N 4 +
L
6 EI
M = 0 = 2 o + u ' ''
L
o

Junta E

12 EI o
+ N 5 + K 12
L3
u + u'
FY = 0 = N 3 N 4 +
L
6 EI
M = 0 = 2 o + u'
L
FX = 0 =

Al resolver las 12 ecuaciones con 12 incgnitas que se han presentado, se encuentra:

24 EI o
L3
12 EI o
N5 =
L3
6 EI
u '' = 2 o
L
N1 =

12 EI o
L3
12 EI
N6 = 3 o
L
6 EI
u ''' = 2 o
L
N2 =

24 EI o
L3
6 EI
u = 2o
L
24 EI o
K 12 =
L3
N3 =

Fuerzas exteriores

Figura 8.11.10 Elementos de la segunda columna de K

12 EI o
L3
6 EI
u' = 2 o
L
24 EI o
K 22 =
L3
N4 =

Roberto Aguiar Falcon

CEINCI-ESPE

260

De sta manera se ha obtenido los elementos de la segunda columna de la matriz de rigidez

a partir de su definicin, se deja al estudiante el obtener los elementos de la primera columna de K .
El resultado completo es:

48 EI o
L3
K=
24 EI o

L3

24 EI o
L3
24 EI o
L3

EJEMPLO N.- 7

Encontrar a partir de su definicin los elementos de la segunda columna de la matriz de

flexibilidad de la estructura presentada en la figura 8.11.1

SOLUCIN
De acuerdo al procedimiento de clculo indicado en el apartado 8.2.2 se tiene:

Q2 = 1 y Qi = 0 i 2

Figura 8.12.1 Fuerza unitaria Figura 8.12.2 Equilibrio de estructura.

Equilibrio de elementos

Figura 8.12.3 Fuerzas internas

ANLISIS MATRICIAL DE ESTRUCTURAS

261

Por consiguiente:

L
4
L
=
4

L
4
L
=
4

L
4
L
=
4

L
4
L
=
4

P1(1) =

P1( 2 ) =

P1(3) =

P1( 4) =

P2(1)

P2( 2 )

P2(3)

P2( 4)

Clculo de las deformaciones

f (1) = f ( 2 ) = f ( 3) = f ( 4 ) =

L 2
6 EI o 1

1
2

p= f P
p1(1) =
p

(1)
2

L2
24 EI o

p1( 2 ) =

L2
=
24 EI o

( 2)
2

L2
24 EI o

L2
=
24 EI o

p1(3) =
p

( 3)
2

L2
24 EI o

L2
=
24 EI o

p1( 4 ) =
p

( 4)
2

L2
24 EI o

L2
=
24 EI o

Clculo del vector q

Se denomina
elementos valen

L2
. Por lo tanto las deformaciones a flexin en cada uno de los
24 EI o

toda vez que tienen ese valor.

Figura 8.12.4 Deformaciones obtenidas

La estructura de la figura 8.12.4 no cumple con la geometra de deformacin, para esto se
procede de la siguiente forma:
i)

El nudo D se desplaza verticalmente una magnitud igual a L como lo indica la figura

8.12.5 de esa manera las vigas giran un ngulo y ya cumplen con el principio de Williot.
Pero el nudo F se ha desplazado a F. En consecuencia se solucion un problema pero se
cre otro y para solucionarlo se da el siguiente paso.

Roberto Aguiar Falcon

CEINCI-ESPE

262
ii)

Con centro en el nudo A, se rota a la estructura en sentido horario un ngulo

ilustra la figura 8.12.6.

Figura 8.12.5

como lo

Figura 8.12.6

De la figura 8.12.6 se tiene:

q1 = L =

L3
24 EI o

q2 = 2 L =

L3
12 EI o

Pero

L3
q1 = F12 =
24 EI o

L3
q 2 = F22 =
12 EI o

EJEMPLO N.- 8
Si en conexin con el problema 6, la matriz T define una transformacin de coordenadas de

la forma q = T q . Interpretar las nuevas coordenadas Q q y encontrar la matriz de rigidez

usando la ley de transformacin de coordenadas de ste nuevo sistema.

T=

1
1

0
1

SOLUCIN

q = T q
q1
q2

1
1

0
1

q1
q 2

q1 = 1 q1 + 0 q 2 = q1
q 2 = 1 q1 + 1 q 2 = q1 + q 2
De donde:

q1 = q1

ANLISIS MATRICIAL DE ESTRUCTURAS

263

q 2 = q 2 q1

Por lo tanto las coordenadas q miden desplazamientos relativos.

Figura 8.13.1 Sistema Q q

Clculo de K

Figura 8.13.2 Sistema Q q

K =Tt K T
1
K =
0

1
1

48 EI o
L3
24 EI o

L3

24 EI o
L3
K =
0

24 EI o
L3
24 EI o
L3

1
1

0
1

0
24 EI o
L3

EJEMPLO N.- 9
La matriz de flexibilidad F para la estructura del problema 7 es:

F=

L3
24 EI o

L3
24 EI o

L3
24 EI o

L3
12 EI o

En este caso el sistema de coordenadas Q q es el presentado en la figura 8.11.2. Se pide

calcular la matriz de flexibilidad F
regla prctica.

para el sistema de coordenadas de la figura 8.14, utilizando la

Roberto Aguiar Falcon

CEINCI-ESPE

264

Figura 8.14 Coordenadas generalizadas del ejemplo 14

SOLUCIN
Se cambia la fila uno a la fila dos.

L3
24 EI o

L3
12 EI o

L3
24 EI o

L3
24 EI o

Finalmente se cambia la columna uno por la columna dos.

L3
12 EI o

L3
24 EI o

L3
24 EI o

L3
24 EI o

8. 5 EJERCICIOS PROPUESTOS
Ejercicio N.- 1
Encontrar la matriz de rigidez utilizando el concepto, de la estructura del ejemplo 8,

empleando las coordenadas Q q .

Ejercicio N.- 2
Generar directamente de su definicin los trminos de la primera columna de la matriz de
flexibilidad de la estructura mostrada en la figura 8.11.1
Ejercicio N.- 3
Calcular la matriz de rigidez para el prtico de la figura 8.2.1, empleando la transformacin de
coordenadas. Si el nuevo sistema de coordenadas es:

ANLISIS MATRICIAL DE ESTRUCTURAS

265

Ejercicio N.- 4
Para la estructura de la figura 8.6 obtener la matriz de flexibilidad si el sistema de
coordenadas generalizadas es el mostrado a continuacin. Utilice la regla prctica.

Ejercicio N.- 5
Hallar la matriz de rigidez usando el concepto del siguiente prtico plano.

Emplear como sistemas de coordenadas del elemento:

a.

Roberto Aguiar Falcon

CEINCI-ESPE

266

b.

c. Sistema a. para elemento vertical y sistema b. para elemento horizontal.

Ejercicio N.- 6
Obtener la matriz de rigidez para la armadura presentada.

Usar como sistema de coordenadas del elemento el siguiente:

Sistema P-p