Análisis Matricial

Tema 1
Relaciones Matriciales fuerza – desplazamiento en una Barra

©Copyright 2015 Francisco Pellicer Martínez. Todos los derechos reservados. 1

Índice: Relaciones Matriciales fuerza–desplazamiento en una Barra

1- Respuesta como sólido rígido. Matriz de Equilibrio (H).

2- Vector Estado (e) y Matriz de Transferencia (G).
3- Relación de transferencia de una barra cargada.
4- Matriz de Flexibilidad (S).
5- Relación de rigidez de una barra cargada (K).
6- Anexo. La viga de Timoshenko.
7- Apéndice.
8- Bibliografía

Objetivos

• Conocer la expresión matricial convencional de las ecuaciones de equilibrio de una barra.

• Determinar la matriz de equilibrio de una barra en cualquier sistema de referencia.
• Identificar los componentes del vector estado en el extremo de una barra.
• Conocer la relación de transferencia de una barra cargada.
• Ser capaz de obtener la matriz de transferencia de la barra y el vector que representa el
efecto de las fuerzas exteriores aplicadas en puntos intermedios.
• Conocer las variables que relaciona la matriz de flexibilidad y las condiciones empleadas
para establecer esta relación.
• Determinar la matriz de flexibilidad de una barra.
• Calcular la matriz de rigidez de una barra a partir de la de flexibilidad, a partir de la matriz de
transferencia y directamente (resolviendo en flexibilidad el problema que proporciona cada
columna).
• Comprender porque las fuerzas de empotramiento perfecto se llaman así, y saber
calcularlas resolviendo en flexibilidad el problema adecuado y también a partir de la relación
de transferencia.
• Determinar los términos asociados a las fuerzas exteriores en las relaciones de rigidez y
transferencia a partir de las correspondientes a una fuerza puntual, recurriendo al teorema
de superposición.
• Conocer los rasgos más característicos de la teoría de vigas de Timoshenko en
contraposición de la teoría de vigas de Navier-Bernoulli.
• Determinar los desplazamientos de la directriz de una viga recta de sección constante según
la teoría de Timoshenko.

©Copyright 2015 Francisco Pellicer Martínez. Todos los derechos reservados. 2

1. Respuesta como sólido rígido. Matriz de Equilibrio

Dada una barra cualquiera bajo unos ejes cualesquiera OXYZ, las relaciones entre las fuerzas del
nudo dorsal f1 y f2 frontal se pueden expresar con las siguientes expresiones.

m2Z

f2Z
m1Z 2 ∆Z

f1Z m1X f2Y m2Y

f1Y f2X
Y
m1Y
m2X
1
f1X
∆X

∆Y
X

=0 → + =0 , =0 → + − ∆ + ∆ =0

=0 → + =0 , =0 → + + ∆ − ∆ =0

=0 → + =0 , =0 → + − ∆ + ∆ =0

Las expresiones anteriores se pueden expresar con notación matricial de la siguiente forma,

+ · =

Donde:

• y son los vectores columna que agrupan las fuerzas que actúan en los nudos dorsal y
frontal, respectivamente, y

• H es, por definición, la matriz de equilibrio de la barra, ya que establece el equilibrio de

fuerzas entre los dos nudos.

Si la barra responde como un sólido rígido se cumple

= · =

©Copyright 2015 Francisco Pellicer Martínez. Todos los derechos reservados. 3

Donde

• y son los vectores columna que agrupan los desplazamientos de los nudos dorsal y
frontal, respectivamente, - el subíndice R recuerda que esta relación solamente es válida si
la barra se comporta como un sólido rígido - y

• H es la matriz de equilibrio de la barra.

La relación existente entre las ecuaciones de equilibrio y las que definen el movimiento como sólido
rígido se llama contragradiencia.

Las ecuaciones de equilibrio del sólido rígido 1-2 son

+ =0
+ =0
+ =0
+ − ∆ + ∆ =0
+ + ∆ − ∆ =0
+ − ∆ + ∆ =0
Y matricialmente se expresan de la siguiente forma:
+ · =
1 0 0 0 0 0 0
# #
0 1 0 0 0 0# 0#
" " " "
"+ 0 0 1 0 0 0" "= 0"
" 0 −∆ ∆ 1 0 0" " 0"
" ∆ 0 −∆ 0 1 0" " 0"
! −∆ ∆ 0 0 0 1! ! 0!
∆ =% −%
∆ =& −&
∆ =' −'
Matriz de Equilibrio que es diferente para las barras de los diferentes tipos de estructuras

1 0 0 0 0 0
0 1 0 0 0 0#
"
0 0 1 0 0 0"
Barra de estructura espacial → = 0 −∆ ∆ 1 0 0"
∆ 0 −∆ 0 1 0"
−∆ ∆ 0 0 0 1!
1 0 0
Barra de pórtico plano → =7 0 1 08
−∆ ∆ 1
1 0 0
Barra de emparrillado plano → =7 ∆ 1 08
−∆ 0 1

©Copyright 2015 Francisco Pellicer Martínez. Todos los derechos reservados. 4

Relaciones entre los movimientos de los puntos 1 y 2 del mismo sólido rígido

2 ∆Z

∆X

∆Y
X

Relación entre los giros –infinitesimales–

: =:

Relación entre los desplazamientos –infinitesimales–

; =; +: ×= − >

Donde ; = =? @ A > y ; = =? @ A >.

C D E
: × = − > = BF F F B = =∆ F − ∆ F >C + =∆ F − ∆ F >D + =∆ F − ∆ F >E
∆ ∆ ∆

; = ; + =∆ F − ∆ F >C + =∆ F − ∆ F >D + =∆ F − ∆ F >E

? 1 0 0 0 ∆ −∆ ?
@ # −∆ 0 ∆ # @ #
0 1 0
A " " A "
"= 0 0 1 ∆ −∆ 0 " "
F " 0 0 0 1 0 0 " F "
F " 0 0 0 0 1 0 " F "
F ! 0 0 0 0 0 1 ! F !

= · =

©Copyright 2015 Francisco Pellicer Martínez. Todos los derechos reservados. 5

2. Vector Estado y Matriz de Transferencia

Se llama vector estado en el extremo i de la barra el que agrupa los desplazamientos di de este
extremo y las fuerzas fi que actúan en él. Lo representaremos por ei.

G =H I ; G =H I

Se llama matriz de transferencia de la barra, y la representaremos por G, la que relaciona los

vectores estado en los extremos, cuando la barra solamente esté cargada en ellos, según la
expresión.

G =K·G

La composición del vector estado para diferentes tipos de estructuras es.

GL = H L I
L

Barra de estructura espacial

? L
@# #
A"
L
"
L"
L = F " ; L =
" L"
F " L"
F ! L!

Barra de pórtico plano

? L
L =7@8 ; L =7 L 8
F L

Barra emparrillado plano

A L
L = 7F 8 ; L =7 L8
F L

©Copyright 2015 Francisco Pellicer Martínez. Todos los derechos reservados. 6

Matriz de transferencia

La relación de transferencia de una barra solicitada exclusivamente en los extremos se puede

desarrollar como

KMM KMN
H I=H I·H I
KNM KNN

Comparando la segunda ecuación matricial con la ecuación de equilibrio

+ · = → =− O
·
Resulta evidente que

KNM = P ; KNN = − O
; KMM = Q

donde Ω es una matriz de ceros de la dimensión adecuada.

Si las fuerzas fueran nulas no habría ningún esfuerzo ni ninguna deformación y, entonces, la
relación de movimientos de los extremos habría de ser la de un sólido rígido.

Por lo que la matriz

KMN

Es una matriz propia de cada tipo de viga, que se puede determinar a partir de los teoremas de
Mohr, las fórmulas de Navier-Bresse o las fórmulas de cálculo de los giros de torsión.

Para una barra espacial recta en la que X es la dirección de la barra e Y y Z son ejes principales de
inercia, a matriz KMN es la siguiente:

R
− 0 0 0 0 0 #
ST R "
RU 0 0 − "
0 0 2SW
6SW "
R "
RU 0 0
0 0 2SW "
6SW "
KMN =
R "
0 0 0 − 0 0
R YZ "
0 0 − R "
2SW 0 − 0 "
SW "
R
0 0 R "
2SW 0 0 −
SW !

©Copyright 2015 Francisco Pellicer Martínez. Todos los derechos reservados. 7

Para una barra de un pórtico plano la matriz de transferencia es

R
− 0 0 #
ST "
RU R "
KMN = 0 − "
6SW 2SW
"
R R "
0 −
−2SW SW !
Para una barra de un emparrillado plano la matriz de transferencia es

RU R
0 #
6SW 2SW "
R "
KMN = 0 − 0 "
YZ
"
R R "
− 0 −
2SW SW !

3. Relación de transferencia de una barra cargada

En el punto anterior se ha establecido la relación entre los vectores estado en los extremos de la
barra solicitada únicamente en estos puntos mediante la expresión

G =K·G

Donde G es la matriz de transferencia de la barra.

Ahora se muestra que esta relación, cuando actúan un sistema de fuerzas en puntos interiores de la
barra, se puede escribir como

G =K·G −[

Donde G sigue siendo la matriz de transferencia de la barra y Q es un vector que depende de las
características de la barra y de las fuerzas exteriores, que no recibe ningún nombre específico.

Determinación del vector Q

Para justificar la posibilidad de escribir la relación de transferencia de la forma indicada y mostrar la

manera de determinar el vector Q se presentará el cálculo de éste en un caso particular: el de la
barra de pórtico plano sometida a una fuerza puntual en un punto interior.

©Copyright 2015 Francisco Pellicer Martínez. Todos los derechos reservados. 8

 P

1
2
t

Para hacerlo se procede como en el caso anterior, es decir:

- Estableciendo las condiciones de equilibrio se determina la segunda fila matricial (f2).

- Calculando los desplazamientos del extremo frontal a partir de los del dorsal se obtendrá la
primera fila matricial (d2).

Condiciones de equilibrio

Y
P
MZ1 MZ2
fY1 fY2

fX1 1 t 2 X
fX2

=0 → + =0

=0 → + +\ =0

, =0 → + − R − \=R − ]> = 0
1 0 0 0
7 8 = 70 1 08 · 7 8−7 \ 8 ↔ =− O
− [N
0 −R 1 −\=R − ]>
La ecuación de equilibrio la proporciona la inversa de la matriz de equilibrio de una barra recta de
pórtico plano.

©Copyright 2015 Francisco Pellicer Martínez. Todos los derechos reservados. 9

Cálculo de los desplazamientos

P
MZ1 MZ2
fY1 fY2

fX1 1 t 2 X
fX2

fX1 -
+

fY1 + + P
+

MZ1 MZ1 - fY1·X

- +
-
P(L-t)

Determinación del vector Q

Cálculo de los desplazamientos

R
? =? −
ST

R RU =R − ]>U
@ =@ +F R− + + \
2SW 6SW 6SW

R R =R − ]>
F =F − + + \
SW 2SW 2SW

Se puede desarrollar en forma matricial

R
− 0 0 # 0
? ? ST " =R − ]>U #
1 0 0 RU R " − \"
7 8 = 70 1 R 8 · 7 @ 8 +
@ 0 − ·7 8− 6SW "
F 0 0 1 F 6SW 2SW " =R − ]> "
R R " − \!
0 2SW
2SW SW !

©Copyright 2015 Francisco Pellicer Martínez. Todos los derechos reservados. 10

Y expresada de forma matricial queda

= · + K_` · − [M

Hasta ahora se ha obtenido

= · + K_` · − [M

=− O
− [N

Que se puede escribir

KMN [M
H I=H I·H I − H[ I ↔ G =K·G −[
P − O N

Con esto queda comprobar que la estructura de la matriz de transferencia de la barra cargada es la
indicada al comienzo y que el vector Q se determina fácilmente a partir de los teoremas de Mohr,
fórmulas de Navier-Bresse, fórmulas de torsión,…etc.

Esta conclusión sigue siendo válida para otros tipos de barras y de cargas.

Forma de proceder con los otros tipos de cargas:

- Fuerza puntual o momento puntual → como el caso anterior.

- Sistema de fuerzas (horizontales o verticales) y de momentos puntuales → por
superposición.
- Fuerza repartida horizontal o vertical, o momento repartido → por superposición a partir
del resultado correspondiente a la fuerza puntual. En este caso, por sumar los efectos de
infinitas fuerzas (o momentos) puntual infinitesimales hay que recurrir a las integrales.

©Copyright 2015 Francisco Pellicer Martínez. Todos los derechos reservados. 11

4. Matriz de Flexibilidad

Definición

Se llama matriz de flexibilidad de una barra la que proporción los desplazamientos del nudo frontal
de ésta, considerándola empotrada por el nudo dorsal, a partir de las fuerzas que actúan en el nudo
frontal.
La matriz de flexibilidad se representa por S. Los desplazamientos del nudo frontal se representan
para recordar que se han obtenido considerando la viga empotrada en el nudo dorsal.
f2 2
d2 = δ

a =b·

La obtención es inmediata. Solamente es necesario recurrir a los teoremas de Mohr, las fórmulas
de Navier-Bresse y/o las fórmulas de la torsión.

©Copyright 2015 Francisco Pellicer Martínez. Todos los derechos reservados. 12

Matriz de flexibilidad

• Barra de estructura espacial

R
0 0 0 0 0 #
ST R "
RU 0 0
0 0 2SW "
3SW "
R
RU 0 − 0 "
0 0 2SW "
3SW "
b=
R
0 0 0 0 0 "
R YZ "
0 0 − R "
2SW 0 0 "
SW "
R
0 0 R "
2SW 0 0
SW !
• Barra de pórtico plano

R
0 0 #
ST "
RU R "
b= 0
3SW 2SW "
"
R R "
0
2SW SW !
• Barra de emparrillado plano

RU R
0 − #
3SW 2SW "
R "
b= 0 0 "
YZ
"
R R "
− 0
2SW SW !

©Copyright 2015 Francisco Pellicer Martínez. Todos los derechos reservados. 13

5. Relación de rigidez de una barra cargada

Se llama relación de rigidez de una barra cargada la que proporciona las fuerzas que actúan en los
extremos de la barra en función de los desplazamientos de estos puntos y de las fuerzas aplicadas
en puntos interiores de la viga. Se escribe como

d d f
H I=H I·H I + e fg → = d· +
d d

En esta relación:

K es la matriz de rigidez de la viga, y solamente depende de las características geométricas y

mecánicas de ésta, y

f0 son las fuerzas de empotramiento perfecto, que dependen de la viga y de las cargas exteriores.

Relación entre la matriz de flexibilidad y la de rigidez

Estado I Estado II

f2 d2 = δ
2
d2R
2

Deformación pura
a =b·
Movimiento de sólido rígido
1 = ·
1
+ · =
d1R

Superposición de los estados I y II

=
= +h
= · +b·

Despejando .

= bO = − > = bO − bO

©Copyright 2015 Francisco Pellicer Martínez. Todos los derechos reservados. 14

Y sustituyendo el resultado en las ecuaciones de equilibrio

=− · = − · bO + bO
Se llega finalmente a

H I = i bO − bO j · H
O
I
−b bO
d d
I = i bO − bO j
O
H
d d −b bO

Inclusión del efecto de las fuerzas exteriores (SUPERPOSICIÓN)

2
d2R
2
f2 d2R 2
f2 f2 0

= +
1
1 1
f1
d1R f1
d1R
f1 0

Inclusión del efecto de las fuerzas exteriores

• Basándonos en el teorema de superposición, se ha comprobado cómo afectan las fuerzas

que actúan sobre puntos interiores de la viga a las que lo hacen en los extremos.

• Esta justificación explica el nombre de fuerzas de empotramiento perfecto.

• Las fuerzas de empotramiento perfecto se pueden determinar fácilmente planteando el

problema en flexibilidad. También se puede recurrir a soluciones tabuladas.

©Copyright 2015 Francisco Pellicer Martínez. Todos los derechos reservados. 15

Obtención de la relación de rigidez a partir de la matriz de transferencia

• En un apartado anterior se ha obtenido

KMM KMN [M = KMM · + KMN · − [M

H I=H I·H I − H[ I → k
KNM KNN N = KNM · + KNN · − [N

• Para determinar la relación de rigidez (1) se despeja f1 de la primera ecuación, y (2) se

sustituye en la segunda

=1> = KMN O = − KMM · + [M > = KMN O − KMN O · KMM · + KMN O · [M

=2> = KNM · + KNN · lKMN O − KMN O · KMM · + KMN O · [M m − [N

=2> = lKNM − KNN · KMN O KMM m · + KNN · KMN O + lKNN · KMN O [M − [N m

y, finalmente, escribirlo en forma matricial

−KMN O · KMM KMN O KMN O · [M

H I=e g·H I+e g
KNM − KNN · KMN O · KMM KNN · KMN O KNN · KMN O · [M − [N

Es decir

d d f
H I=H I·H I + e fg
d d

Donde:

d = −KMN O · KMM

d = KMN O

d = KNM − KNN · KMN O · KMM

d = KNN · KMN O

No solo se ha establecido otra manera de determinar la matriz de rigidez sino que también de
obtener las fuerzas de empotramiento perfecto.

©Copyright 2015 Francisco Pellicer Martínez. Todos los derechos reservados. 16

También se puede obtener la matriz de rigidez de forma directa

• Si las fuerzas que actúan en puntos interiores de la viga son nulas, la relación de rigidez se
reduce a

d d f
d d
H I=H I·H I + e fg = H I·H I
d d d d

• Se supone la viga empotrada en los dos extremos e imponemos un desplazamiento unidad en

uno de los grados de libertad. Entonces:

- El vector desplazamiento se reducirá a un vector de ceros con un único elemento igual a 1.

- El producto de la matriz de rigidez por este vector será una columna de esta matriz.

- Esta columna tendrá que ser igual al vector que agrupa las reacciones del problema
planteado.

- Esto proporciona un nuevo método para determinar la matriz de rigidez.

Matriz de rigidez de una barra recta de sección constante (teoría de Navier-Bernoulli)

ST ST
0 0 − 0 0 #
R R "
12SW 6SW 12SW 6SW "
0 0 − U ?
# RU R R R "
6SW 4SW 6SW 2SW " @ #
" 0 0 − "
"= R R R R " F
"
ST ST " ? "
" − 0 0 0 0 "
" R R @ "
" F !
! 12SW 6SW 12SW 6SW "
0 − U − 0 −
R R RU R "
6SW 2SW 6SW 4SW "
0 0 −
R R R R !

Pórtico plano

©Copyright 2015 Francisco Pellicer Martínez. Todos los derechos reservados. 17

Matriz de rigidez de una barra recta de sección constante (teoría de Navier-Bernoulli)

12SWo 6SWo 12SWo 6SWo

0 − − 0 − #
RU R RU R "
YZ YZ
0 0 0 − 0 "
R R " A
# #
6SWo 4SWo 6SWo 2SWo " F
" − 0 0 "
o "= R R R R " F
"
" 12SW o 6SWo 12SWo 6SWo " A "
− U 0 0 " F "
" R R RU R "
! YZ YZ F !
0 − 0 0 0 "
R R "
6SWo 2SWo 6SWo 4SWo "
− 0 0 !
R R R R

Emparrillado plano

©Copyright 2015 Francisco Pellicer Martínez. Todos los derechos reservados. 18

6. Anexo. La viga de Timoshenko

Introducción

• En Mecánica de Estructuras y Teoría de Estructuras del Grado en Ingeniería Civil se estudia

la Teoría de Vigas de Navier-Bernoulli, que es aplicable a piezas suficientemente alargadas
(h/L ≤1/10*).

• La piezas que, siendo alargadas, no los son lo suficiente para poder aplicar la regla (1/5 ≤ h/L
≤ 1/10*) aún se pueden estudiar como vigas, recurriendo a la Teoría de Timoshenko, cuya
principal característica es que tiene en cuenta la deformación por cortante.

• La mayoría de los programas comerciales de cálculo de estructuras permiten considerar

barras que responden según esta segunda teoría.

• Este apartado tiene como finalidad mostrar los rasgos fundamentales de la viga de
Timoshenko con la profundidad necesaria para permitiros esta posibilidad.



(*) Estos límites son aproximados.

©Copyright 2015 Francisco Pellicer Martínez. Todos los derechos reservados. 19

A1 - La viga de Navier-Bernoulli
Hipótesis cinemática

Las secciones transversales de la viga permanecen planas, indeformables y normales a la directriz

durante el proceso de deformación.

Y
Z
Z
Y
X

Sección plana y
perpendicular a la directriz

Directriz inicial X Directriz deformada

Esfuerzos

• Las esfuerzos coherentes con la hipótesis de Navier-Bernoulli son:

- el axil y

- los momentos flectores

• Eso quiere decir que existen unas ecuaciones constitutivas que determinan estos esfuerzos
en función de parámetros que miden la deformación (las relaciones axil-alargamiento
unitario y momento-curvatura).

• Así mismo, implica que todas las ecuaciones que describen la respuesta de la viga se pueden
escribir recurriendo, exclusivamente, a estos tres esfuerzos.

• Los cortantes no son esfuerzos coherentes con esta teoría –por eso no hay ninguna
ecuación constitutiva del cortante –. Todo esto, suele emplearse como variables auxiliares –
representan derivadas de los momentos – para escribir algunas ecuaciones, como ahora las
de equilibrio.

©Copyright 2015 Francisco Pellicer Martínez. Todos los derechos reservados. 20

Deformación originada por el axil

• En un prisma elemental (una rebanada

de viga de espesor dX), el axil origina
un desplazamiento de la cara frontal
respecto de la dorsal en la dirección
de la directriz.

• El alargamiento unitario de la directriz

es proporcional al axil y queda
determinado por p = q · S · T.

• El axil no origina ningún

desplazamiento perpendicular a la
directriz ni ningún giro.

Deformación originada por el momento flector

• En un prisma elemental, el flector flector. No origina ningún

origina un giro de la cara frontal desplazamiento paralelo a la directriz.
respecto de la dorsal en la dirección
de la directriz.

• El giro por unidad de longitud es

proporcional al flector y queda
determinado por =r·S·W (en
ejes principales de inercia).

• Coincide con la curvatura de la

directriz.

• La figura muestra el desplazamiento

perpendicular a la directriz debido al

©Copyright 2015 Francisco Pellicer Martínez. Todos los derechos reservados. 21

Deformación originada por el cortante

• La teoría de Navier-Bernoulli desprecia la deformación por cortante, es decir, supone que el

cortante, por grande que sea, no origina ninguna deformación.

• Daros cuenta que esto explica que no haya ninguna ecuación constitutiva del cortante: ¿cómo
puede calcularse el cortante a partir de una deformación inexistente?

• Naturalmente, la realidad es que las tensiones tangenciales si que producen deformaciones,

pero los efectos de estos en los movimientos de la viga son despreciables si esta es muy
alargada (recordar, h/L ≤ 1/10) y, entonces, la teoría de Navier-Bernoulli proporciona resultados
razonables.

A2 - La viga de Timoshenko
Hipótesis cinemática

• Las secciones transversales de la viga permanecen planas e indeformables durante el proceso

de deformación.

Y
Z
Z
Y
X

Sección plana que NO es

perpendicular a la directriz

X Directriz deformada
Directriz inicial

• Daros cuenta: respecto a la viga de Navier-Bernoulli, desaparece la condición que exige que las
secciones permanezcan perpendiculares a la directriz.

©Copyright 2015 Francisco Pellicer Martínez. Todos los derechos reservados. 22

Esfuerzos

• Los esfuerzos coherentes con la hipótesis de Timoshenko son

- el axil

- los momentos flectores, y

- los cortantes

• Esto quiere decir que existen unas ecuaciones constitutivas que determinan estos esfuerzos en
función de parámetros que miden la deformación (las relaciones axil-alargamiento unitario,
momento-curvatura y cortante-distorsión angular).

• Así mismo, implica que se necesitan todos los cinco esfuerzos para poder escribir todas las
ecuaciones que describen la respuesta de la viga.

Deformaciones originadas por el axil y por el flector

• Todo lo que se ha dicho respecto a las deformaciones y movimientos originados por el axil de la
viga de Navier-Bernoulli sigue siendo válido en la de Timoshenko.

• Todo lo que se ha dicho respecto a las deformaciones y movimientos por el flector en la viga de
Navier-Bernoulli sigue siendo válido en la de Timoshenko, excepto que el giro por unidad de
longitud χ pueda identificarse con la curvatura de la directriz.

• En un problema plano, en teoría lineal y para la viga de Timoshenko, la curvatura resulta

s @ st
= +r
s% s%

donde t es la deformación originada por el cortante (ejes principales de inercia)

• En la teoría de Timoshenko, el giro es la derivada de la flecha menos la deformación por

cortante.

©Copyright 2015 Francisco Pellicer Martínez. Todos los derechos reservados. 23

Deformación originada por el cortante

• En un prisma elemental, el cortante

origina un desplazamiento de la carga
frontal respecto de la dorsal
perpendicular a la directriz (en la dirección
del cortante).

• La distorsión angular es proporcional al

cortante y queda determinada por
u = t · Y · Tv

• Tv es una característica de la sección que

sólo depende de la geometría.

• El cortante no origina desplazamiento

paralelo a la directriz ni giro.

Determinación del área de cortante wx

• Sección genérica.

©Copyright 2015 Francisco Pellicer Martínez. Todos los derechos reservados. 24

Fórmulas de Navier-Bresse generalizadas Problema plano de la viga recta

(Para calcular los desplazamientos y giros de la directriz, en la teoría de Timoshenko)

| p
?y = ?z − F z · =&y − &z > + { · s%
}
ST

| | u
~y = ~z − F · =%y − %z > + { • · =%y − %z >€ s% + { · s%
z
}
SW }
YTv

‚
Fy = Fz + { · s•
f SW

Ecuaciones referidas a los ejes principales de inercia

Una interpretación que nos permitiría generalizar el segundo teorema de Mohr y enunciarlo del
siguiente modo: la flecha en A es igual a la flecha en B más el giro de A por la distancia desde A
hasta la vertical que pasa por B, más el momento estático del área de flectores respecto de esta
recta dividido por EIz, más el área de cortantes dividida por GAQY.

©Copyright 2015 Francisco Pellicer Martínez. Todos los derechos reservados. 25

Energía de deformación de la viga -incluyendo la correspondiente al cortante-

(Para calcular los desplazamientos y giros de la directriz, en teoría de Timoshenko, mediante el

segundo teorema de Castigliano).

• Viga recta, problema plano

‚
p u
ƒ ={ •
∗
+ + € · s%
f ST SW YTv

Viga recta, caso general

‚
p u u …
ƒ∗ = { • + + + + ++ € · s%
f ST SW SW YTv YTv YZ

Ecuaciones referidas a los ejes principales de inercia

El factor de cortante α

• En la formulación matricial del problema de la viga de Timoshenko, la flexibilidad frente al

cortante suele representarse mediante el factor de cortante, que se define como

12SW 1 R
†= → =†
YTv R YTv 12SW

• Es adimensional

• El factor de cortante mide la sensibilidad de la pieza frente a la deformación por cortante:

cuanto más grande es, mayor resulta el efecto de esta deformación.

• Suele ser del orden de las centésimas o, como mucho, de las décimas.

• Haciendo α=0 en las expresiones determinadas según la teoría de Timoshenko se obtienen las
correspondientes a la de Navier-Bernoulli.

• Influencia de la deformación por cortante en los movimientos (Monleón, 1999).

©Copyright 2015 Francisco Pellicer Martínez. Todos los derechos reservados. 26

• EL factor de cortante para diferentes secciones y relaciones h/L

Secciones rectangular y circular, α=0,2

Secciones metálicas, α=0,3

©Copyright 2015 Francisco Pellicer Martínez. Todos los derechos reservados. 27

7. Apéndice. Matrices de rigidez y Fuerzas de Empotramiento Perfecto.

Matriz de rigidez de una barra recta de sección constante (teoría de Navier-Bernoulli)

Estructura espacial

©Copyright 2015 Francisco Pellicer Martínez. Todos los derechos reservados. 28

Matriz de rigidez de una barra recta de sección constante (teoría de Navier-Bernoulli)

ST ST
0 0 − 0 0 #
R R "
12SW 6SW 12SW 6SW "
0 0 − U
RU R R R " ?
# #
6SW 4SW 6SW 2SW " @
" 0 0 − "
"= R R R R " F
"
ST ST " ? "
" − 0 0 0 0 "
" R R @ "
" F !
! 12SW 6SW 12SW 6SW "
0 − U − 0 −
R R RU R "
6SW 2SW 6SW 4SW "
0 0 −
R R R R !

Pórtico plano

Matriz de rigidez de una barra recta de sección constante (teoría de Navier-Bernoulli)

12SWo 6SWo 12SWo 6SWo

0 − − 0 #
−
RU R RU R "
YZ YZ
0 0 0 − 0 " A
R R "
# 6SWo 4SWo 6SWo 2SWo " F #
" − 0 0 "
"= R R R R " F
"
" 12SW o 6SWo 12SWo 6SWo " A "
− U 0 0 " F "
" R R RU R " F !
! YZ YZ "
0 − 0 0 0
R R "
6SWo 2SWo 6SWo 4SWo "
− 0 0 !
R R R R

Emparrillado plano

©Copyright 2015 Francisco Pellicer Martínez. Todos los derechos reservados. 29

Fuerzas de empotramiento perfecto de una barra recta de sección constante sometida a una
fuerza puntual vertical. (teoría de Navier-Bernoulli)

Fuerzas de empotramiento perfecto de una barra recta de sección constante sometida a una
fuerza puntual horizontal. (teoría de Navier-Bernoulli o Timoshenko)

©Copyright 2015 Francisco Pellicer Martínez. Todos los derechos reservados. 30

Apéndice

Matriz de rigidez de una barra recta de sección constante

(pórtico plano, teoría de Timoshenko)

ST ST
0 0 − 0 0 #
R R "
1 12SW 1 6SW 1 12SW 1 6SW
0 · U · 0 − · U · "
1+† R 1+† R 1+† R 1+† R " ?
# 1 6SW 4 + † 4SW 1 6SW 4 + † 2SW " #
@
" 0 · · 0 − · "
1+† R 1+† R 1+† R 1+† R " F
"= "
ST ST " ? "
" − 0 0 0 0
" " @ "
R R
! 1 12SW 1 6SW 1 12SW 1 6SW " F !
0 − · U − · 0 · − · "
1+† R 1+† R 1+† RU 1+† R "
1 6SW 4 + † 2SW 1 6SW 4 + † 4SW "
0 · 0 − ·
1+† R 1+† R 1+† R 1+† R !

Pórtico plano

Matriz de rigidez de una barra recta de sección constante

(emparrillado plano, teoría de Timoshenko)

1 12SW& 1 6SW& 1 12SW& 1 6SW&

· 0 − · − · 0 − · #
1 + †' R 3
1 + †' R 2
1 + †' R 3
1 + †' R2 "
YZ YZ "
0 0 0 − 0
R R " A
# 1 6SW& 4 + †' 4SW& 1 6SW& 4 + †' 2SW& " #
F
" − · 2 0 · · 0 · "
"= 1 + †' R 1 + †' R 1 + †' R2 1 + †' R " F
o
"
" 1 12SW & 1 6SW & 1 12SW& 1 6SW& " A "
− · 0 · 2 · 0 · 2 " F
" 1 + †' R 3
1 + † R 1 + †' R3 1 + †' R "
'
" F !
! YZ YZ
0 − 0 0 0 "
R R "
1 6SW& 4 + †' 2SW& 1 6SW& 4 + †' 4SW& "
− · 0 · · 0 ·
1 + †' R 2
1 + †' R 1 + †' R2 1 + †' R !

©Copyright 2015 Francisco Pellicer Martínez. Todos los derechos reservados. 31

Matriz de rigidez de una barra recta de sección constante

(estructura espacial, teoría de Timoshenko)

Fuerzas de empotramiento perfecto de una barra recta de sección constante sometida a una
fuerza puntual vertical (teoría de Timoshenko).

©Copyright 2015 Francisco Pellicer Martínez. Todos los derechos reservados. 32

Fuerzas de empotramiento perfecto de una barra recta de sección constante sometida a un
momento puntual. (teoría de Timoshenko)

8. Bibliografía

• Albiges, M., Coín. Y Journet, H. (1970), Estudio de las estructuras por métodos matriciales,
Editores Técnicos Asociados, Barcelona.

• Casanova-Colón, J. 1999. Apuntes de Cálculo de Estructuras III. Ingeniería de Caminos, Canales

y Puertos. Universidad Politécnica de Valencia.

• Casanova-Colón, J. 2007. Apuntes de Cálculo de Estructuras III. Ingeniería de Caminos, Canales

y Puertos. Universidad Politécnica de Valencia.

• González, J.R. y Samartín A. (1999), Cálculo de Estructuras, Colegio de Ingenieros de Caminos,

Canales y Puertos, Madrid.

• Livesley, R.K. 1964. Métodos matriciales para el cálculo de estructuras, Blume, Madrid, 1970.

• Monleón, S. 1999. Análisis de Vigas, Arcos, Placas y Láminas: una Presentación Unificada,
Editorial de la Universidad Politécnica de Valencia, Valencia.

• Sáez-Benito, J.M. 1975. Cálculo Matricial de Estructuras Formadas por Piezas Prismáticas,
Fondo Editorial de Ingeniería Naval, Madrid.

©Copyright 2015 Francisco Pellicer Martínez. Todos los derechos reservados. 33