caLcuLo de estructurasdos

TEMA 07. INTRODUCCIÓN AL CÁLCULO MATRICIAL DE ESTRUCTURAS

1.CONTINUIDAD
1.1. Introducción
Cuando se trata de proyectar un elemento a flexión sobre varios apoyos consecutivos existe, en
general, una ventaja indudable en disponer una viga continua. Ello implica que la deformación tiene que
corresponder a una pieza enteriza, con curvatura continua, sin brusquedades. Si se disponen tramos
sueltos, cada uno apoyado en sus extremos, el régimen de cada tramo es de momentos POSITIVOS, y
la curvatura experimenta un cambio brusco, un quiebro, al pasar por el apoyo, ya que cada tramo posee
en ese punto una inclinación diferente. No existe la interacción de los diferentes tramos, limitándose
estos a compartir el apoyo común, que responde con la suma de reacciones que corresponderían a
cada tramo.

El comportamiento en continuidad permite el funcionamiento en régimen alternante de momentos

POSITIVOS y NEGATIVOS, como el que se produce en una viga con voladizos. Cada tramo interacciona
con los con los colindantes pactando una única inclinación sobre cada apoyo, sin giro relativo entre los
dos tramos, respondiendo el apoyo con la reacción que corresponde al régimen de momentos de cada
tramo que acomete a él.

Efectos de la continuidad

En la solución en continuidad, con un régimen de momentos positivos y negativos, los provenientes

de una nueva carga dada pueden ser menores que en la situación de tramos sueltos, resultando ser
válida una pieza de menor sección. En los casos en los que el diseño es por momentos o hay que
reducir la flecha, la viga continua la reduce drásticamente, siendo una solución ventajosa al construir en
acero u hormigón. En madera, donde el diseño es frecuentemente por cortante no suele haber diferencia
entre una u otra solución.

1
introduccion al calculo matricial de estructuras
 1.2.Hiperestatismo

En una viga continua no pueden obtenerse los

esfuerzos a partir de las fuerzas que actúan, ya que para
su cálculo necesitamos las reacciones de los apoyos,
que no pueden deducirse de las condiciones de equilibrio
(ecuaciones de la estática).
ΣV= 0
ΣM= 0

Sistema equivalente: si no existiesen los apoyos

intermedios, el sistema sería isostático, y las
incongruencias estarían en los puntos 2 y 3, que
experimentarían un descenso, luego hay una reacción
que se opone a ello. Tenemos dos nuevas condiciones
(ecuac. de compatibilidad) que nos permiten resolver la
viga.

ΣV= 0 δ =0
ΣM= 0 δ =0
2

2 Cálculo de vigas continuas. MÉTODO DE LA FLEXIBILIDAD

Para el cálculo de vigas continuas es corriente tomar

como incógnita los momentos en los apoyos, ya que
una vez conocidos, sabemos las leyes que tienen los
distintos tipos de cargas.

En este caso nos interesa conocer el valor de M2, de

donde deducimos el resto de los esfuerzos (cortante y
axil).

La continuidad obliga a la existencia de un momento

sobre el apoyo común, M, que produce un giro en cada
viga, siendo el ángulo ϕ a los lados igual y del mismo
sentido.

2
introduccion al calculo matricial de estructuras
 Para el análisis separamos la viga en sus dos tramos. M2 es la acción de un tramo sobre el otro:

Giro debido a la carga P Giro debido a la carga q

Pl 12 Pl23
ϕp = + ϕq = −
16 EI 24 EI
Giro debido al momento M2 Giro debido al momento M2
M 2l1 M 2 l2
ϕM = − ϕM = +
3EI 3EI

Sabemos los giros producidos por las diferentes solicitaciones, bien hallándolos o a través de las tablas.
La coherencia geométrica o compatibilidad en las deformaciones obliga a que el giro en un extremo y en
otro del apoyo sea el mismo, luego:

Pl12 M l ql23 M l M 2 (l1 + l2 ) Pl12 ql23

ϕ d = ϕi ⇒ − 21 = − + 22 ⇒ = +
16 EI 3EI 24 EI 3EI 3EI 16 EI 24 EI

Vemos como en el segundo término nos quedan la suma de los giros producidos por las cargas en el
apoyo 2. El primer término son los giros producidos por el momento.

2.2 Cálculo de vigas continuas. GENERALIZACIÓN DEL MÉTODO

En este caso sabemos que los momentos en 1 y 5 son cero, luego hay que hallarlos en los demás
puntos. Hacemos lo mismo que en el caso anterior: separamos la viga en sus tramos y dibujamos los
momentos Mi que actúan.
i
( NOTA: ϕc indica el giro sobre el nudo i producido por la carga correspondiente )

3
introduccion al calculo matricial de estructuras
 ϕ1c − ϕc2 + ϕc2 − ϕ3c ϕ3c − ϕc4

M 2l1 M 2 l2 M 3l 2 M 3l3 M 3l3 M 4 l4

− − −
3EI 3EI 3EI 3EI 6 EI 3EI
M 3l 2 M 2 l2 M 4l3 M 4l3
− −
6 EI 6 EI 6 EI 3EI

Tenemos que igualar los giros en los distintos puntos, de donde obtenemos las siguientes ecuaciones:

M 2l1 M l M l
2. ϕ1c − = −ϕc2 + 2 2 + 3 2
3EI 3EI 6 EI
M 2 l 2 M 3l2 Ml M l
3. ϕc2 − − = −ϕ3c + 3 3 + 4 3
6 EI 3EI 3EI 6 EI
M 3l3 M 4 l3 M l
4. ϕc3 − − = −ϕc4 + 4 4
6 EI 3EI 3EI

Ordenamos las ecuaciones, en un término los momentos y en el otro los giros:

M 2 (l1 + l2 ) M 3l 2
+ = ϕ1c+ 2
3EI 6 EI
M 2 l2 M 3 (l2 + l3 ) M 4 l3
+ + = ϕc2 + 3
6 EI 3EI 6 EI
M 3l3 M 4 (l3 + l4 )
+ = ϕ3c + 4
6 EI 3EI

- las incógnitas que queremos hallar son los distintos momentos Mi

- el término independiente es la suma de los giros debidos a la carga en los nudos

( NOTA: ∑ϕ i
: sumatorio de giros en el nudo i )

4
introduccion al calculo matricial de estructuras
 En el caso anterior, si suponemos EI = cte., la matriz del sistema queda:

l + l2 l2 
1 0   M 2   ∑ ϕ2 
 3 6     
   
1  l2 l2 + l3 l3  
M  = ϕ 
 3  ∑ 3 
 
EI  6 3 6 
    
l3 l3 + l4   M   ϕ 
 0
6 3 
  3  ∑ 4 


2.3Cálculo de vigas continuas. FORMACIÓN DE LA MATRIZ DE FLEXIBILIDAD

- Es siempre una matriz simétrica
- En la diagonal va la suma de las longitudes de los tramos que afectan al nudo en cuestión
dividido entre 3EI (Si EI = cte. sale fuera de la matriz)
- Los demás términos por encima de la diagonal son (recordamos que la matriz es simétrica)
- cero, si no influyen entre sí (M4 sobre el nudo 2)
- la longitud del tramo en la que ese momento actúa sobre el nudo (M3 sobre el nudo dos
→ L2) dividido entre 6EI .
- El término independiente es la suma de los giros en cada nudo.

2.4Cálculo de vigas continuas. OBSERVACIONES SOBRE EL MÉTODO

1. Los giros proporcionados por las cargas son siempre positivos con cargas verticales hacia abajo.
Tenemos que son siempre el sumatorio directo en valor absoluto.

ϕ1c + ϕc2 = ϕ1+ 2

Luego este es un método pensado para cargas VERTICALES. En el caso de una carga dirigida
hacia arriba:

En este caso el giro producido por esa carga se restaría.

2. La ventaja que tenemos con este método es que con la geometría de la viga ya podemos establecer
la matriz, y añadir al final el término independiente, que pueden ser varias hipótesis de carga.

L L L  M 2   ∑ 2 ∑ ∑ 
     
     
L L L  M  =  
   3  ∑ 3 ∑ ∑ 
     
   M   
L L L  4   ∑ 4 ∑ ∑ 

5
introduccion al calculo matricial de estructuras
2.5 Matriz de flexibilidad. EJEMPLO

La matriz de flexibilidad es inmediata:

4 + 5 5 
 0   M2 
 3 6   
 5  
5+8 8  M 
 
 6 3 6   3
  
 0 8 8 + 6   
6 3 
  M4 


Los términos independientes:

qL13 PL22 2 * 43 10 * 52 )
∑ ϕ2 = +
24 EI 16 EI
⇒
24 EI
+
16 EI
= 18,2916

10 * 5 2
1* 8 3
)
∑ ϕ3 = 16 EI + 24 EI = 36,9583
1 * 83 5 * 6 2 )
∑ ϕ4 = +
24 EI 16 EI
= 32,583

Podemos establecer el sistema de ecuaciones:

)
 3 0,833 0   M 2  18,2916 
     
     
 0,833 4,333 1,333   M  =  36,9583) 
   3  
     
 0 1,333 4,666   M   32,583) 
  4  

De donde sacamos los valores de los momentos directamente:

M2= 4,4129
M3= 6,0659
M4= 5,2501

Con estos datos ya podemos dibujar las gráficas de momentos, cortantes, axiles y hacer todos los
cálculos que sean necesarios.

6
introduccion al calculo matricial de estructuras
 3Cálculo de vigas continuas. MATRIZ DE RIGIDEZ
En el método anterior se hace un análisis de la viga continua a partir de otra equivalente que
difería de la original en que respetaba el equilibrio y no tenía más remedio que acusar incoherencia
geométrica, igualando los giros en los nudos. Ahora se intenta plantear el problema a la inversa, es
decir, suponer inicialmente coherencia total geométrica –giro nulo en los nudos- aunque se produzca un
desequilibrio de momentos que es preciso corregir. Este método nos permite una fácil generalización a
todo tipo de problemas, como es el de los pórticos planos.

En este caso se consideran las vigas como EMPOTRADAS.

Si uno de los apoyos gira es porque existe un momento que le ha

producido ese giro.

En el empotramiento el giro es CERO.

ML M 1 L M
− + =0 ⇒ M1 =
6 EI 3EI 2

Como tenemos un momento, en el otro apoyo tenemos un

momento que es la mitad.

Si hallamos el giro que se produce:

M
L
ML 2 ML 4 EI 2 EI
ϕ= − = ⇒M = ϕ ; M1 = ϕ
3EI 6 EI 4 EI L L

3.1Matriz de rigidez. EXPLICACIÓN DEL MÉTODO

Queremos resolver esta estructura.

La viga se deforma para soportar las cargas.

Para el cálculo descomponemos la viga en sus cargas y en sus

movimientos:

1. La estructura sin la deformación (impedimos el giro en el

extremo apoyado) y con la carga.

2. Consideramos únicamente la deformación, el giro en el

apoyo.

7
introduccion al calculo matricial de estructuras
 En 1, los momentos al no existir el giro son inmediatos: momentos de empotramiento
perfecto.

En 2 hemos calculado anteriormente el valor.

a) Calculamos el movimiento:
Por las condiciones de sustentación sabemos que el momento en el apoyo es cero:

qL2 4 EI L qL2 qL3

− + ϕ=0 ⇒ ϕ= * = EI
12 L 4 EI 12 48

Una vez calculado en giro, hallamos el momento:

qL2 2 EI qL3 qL2

MT = + * =
12 L 48 EI 8

Observamos que con este método el problema se resuelve de un modo indirecto; aunque
lo que se buscan son las solicitaciones, se plantean en función de las deformaciones, que son las
incógnitas del método; una vez resueltas éstas se vuelve tras los pasos anteriores para calcular las
solicitaciones. La ventaja está en que en general conocemos los momentos de empotramiento, con lo
que no hace falta calcularlos cada vez que se plantea el problema.

3.2 Matriz de rigidez. PLANTEAMIENTO

El planteamiento del método es como en el caso anterior;

-consideramos todos los nudos como empotrados

-consideramos las deformaciones en cada nudo
individualmente, es decir, los giros impedidos uno a uno.

No hemos tenido en cuenta el momento exterior aplicado en

2. Sólamente lo consideraremos al sumar todos los momentos,
que en ese apoyo no serán igual a cero, sino que tendrán que
ser igual a Mext.

La suma de momentos:

qL12 PL2 4 EI 4 EI 2 EI
2. − + + ϕ1 + ϕ1 + ϕ2 = M
12 8 L1 L2 L2

PL2 2 EI 4 EI
3. + ϕ1 + ϕ2 = 0
8 L2 L2

8
introduccion al calculo matricial de estructuras
 Agrupamos términos:

 4 EI 4 EI  2 EI  qL2 PL2 
 + ϕ1 + ϕ2 = M −  − 1 + 
 L1 L2  L2  12 8 
2 EI 4 EI  PL2 
ϕ1 + ϕ2 = 0 − − 
L2 L2  8 

Si suponemos EI=cte

4 4 2
 + 
 ϕ1   ∑ M E − ∑ Mˆ 1 
1
L L2 L2 
EI  1   =
 2 4  ϕ2   ∑ M E − ∑ Mˆ 2 
2

 L2 
 L2

En este caso el término independiente es la suma de momentos exteriores menos el sumatorio

de los momentos de empotramiento.
Al resolver el sistema obtendemos los giros de los nudos. Para conocer el valor de los MOMENTOS,
hay que sumar todos los momentos que tenemos hallados en función de los GIROS.

En una viga cualquiera, la suma de momentos es:

4 EI 2 EI
Mˆ + ϕ1 * + ϕ2 *
L1 L1

O sea,

Mempotramiento+ giro del nudo*4EI/L1+ giro del nudo de enfrente*2EI/L1

9
introduccion al calculo matricial de estructuras
3.3Matriz de rigidez. EJEMPLO

La matriz de RIGIDEZ está formada por los nudos que se mueven, luego podemos sustituir el vuelo
por su acción.

Formamos la matriz de rigidez. Es una matriz simétrica:

ϕ1 ϕ2 ϕ3

4 EI 2 EI
ϕ1 0
L1 L1

2 EI 4 EI 4 EI 2 EI
ϕ2 +
L1 L1 L2 L2

ϕ3 2 EI 4 EI 4 EI
0 +
L2 L1 L2

Los momento de empotramiento son:

El sistema queda:

 0 , 66 0 , 33 0   ϕ1  4 − 6 
      Soluciones:
     
EI  0 , 33 1 , 66 0 ,5   ϕ  =  3 ,3  ϕ1= -5/EI
      ϕ2= 3.93/EI
2

     
 0
 0 ,5 1 , 66  ϕ 
 3 
 − 3 ,3


 ϕ3= -3.17/EI

10
introduccion al calculo matricial de estructuras
Una vez calculados los giros, los momentos son:

4 EI 2 EI
Mˆ + ϕ1 * + ϕ2 *
L1 L1

nudo2

3.93 4 EI  3.17 2 EI 
2 . 2 .7 + * + − *  = 5.04
EI 4  EI 4 

nudos 3 y 4

 3.17 4 EI 
3. 6 +  − *  = 3 .9
 EI 6 

 3.17 2 EI 
4. − 6 + 0 +  − *  = −7.06
 EI 6 
En este caso el nudo 4 no gira por ser un empotramiento.

1. En el nudo 1 sabemos que M=4 (momento debido al

voladizo).

Con estos datos obtenemos la gráfica de momentos, de donde sacamos los cortantes, axiles...

11
introduccion al calculo matricial de estructuras
4. análisis de pórticos. MATRIZ DE RIGIDEZ
El método de la rigidez se basaba en plantear un giro nulo en los nudos -suponiéndolos todos
empotrados- para luego equilibrar los momentos. La existencia de nudos rígidos en los dinteles de los
pórticos nos permite plantear las condiciones de deformación ya que el giro del soporte va a depender
del de la viga o dintel.
Para el análisis de los pórticos vamos a considerar inicialmente dos casos:
- pórticos intraslacionales: los nudos sólo pueden girar
- pórticos traslacionales: además pueden experimentar desplazamientos.

4.1matriz de rigidez. PÓRTICOS INTRASLACIONALES

En estos los nudos no se desplazan, únicamente pueden girar.

1. Si tomamos como hipótesis de partida el comportamiento de cada pieza con los extremos fijos
y sin girar, de hecho estamos ignorando la estructura como tal; el problema es sólo el de un conjunto de
barras por separado. El dintel está sometido en ambos extremos al que hemos denominado momento
de empotramiento perfecto, que para sección constante y carga distribuída era de valor M0=q*L2/12.

2. Estudiamos el comportamiento geométrico de la estructura, es decir, los giros ϕ de cada nudo.

Es en este paso en el que realmente se considera a la estructura como a un todo, en el que unas partes
interaccionan con las demás.
Para que se produzca un giro ϕ en el soporte, debe actuar en su extremo un momento,

4 EI p
M p = ϕ*
L
y en la viga:

2 EI v
Mv = ϕ*
L
La ecuación de equilibrio es:

M p + Mv = M0
o sea;

4 EI p
2 EI v qL2
ϕ* + ϕ* =
L L 12
Ecuación con una única incógnita ϕ, que puede resolverse fácilmente. En el caso de la existencia
de más nudos, existirán más giros ϕn, que nos darán tantas incógnitas como ecuaciones.

12
introduccion al calculo matricial de estructuras
4.2 matriz de rigidez. APLICACIÓN PRÁCTICA EN EL CASO DE PÓRTICOS PLANOS
Este es un pórtico intraslacional que puede únicamente
girar en los nudos 2 y 3.

Estudiamos qué ocurre cuando gira el nudo 2. Estos son

los momentos debidos al giro estudiados anteriormente:

ϕ2 ϕ3
1+ 2 2
ϕ2  4 EI   2 EI 
   
 L   L 

ϕ3

Giramos el otro nudo:

ϕ2 ϕ3
1+ 2 2
ϕ2  4 EI   2 EI 
   
 L   L 
2 2+3
 2 EI   4EI 
ϕ3    
 L   L 

Los momentos a los que hay que igualar los momentos producidos por los giros son los exteriores
más los producidos por las cargas.

∑M E − ∑ Mˆ i
En este caso:

 qL2 PL 
∑ M E −  12 − 8 
 
 qL2 
∑ M E −  − 12 
 

13
introduccion al calculo matricial de estructuras
 Luego las ecuaciones reales son:

1+ 2 2
 4 EI   2 EI   qL22 PL1 
  * ϕ2 +   * ϕ3 = 0 −  − 
 L   L   12 8 

2 2+3
 2 EI   4 EI   qL22 
  * ϕ2 +   * ϕ3 = 0 −  − 
 L   L   12 

El sistema con la matriz es:

  4 EI 1+ 2  2 EI  
2
  qL22 PL1  
     0 −  − 
 L   L    ϕ2    12 8  
  = 
  2 EI  2  4 EI  
2 +3
 ϕ3   qL22  
       0 −  −  
  L   L     12  

NOTA: En el caso del mismo pórtico con un apoyo articulado, habría que considerar el giro en el
mismo al no estar impedido, luego en el nuevo sistema habría que incluirlo.

ϕ2 ϕ3 ϕ4
1+ 2 2 ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○
ϕ2  4 EI   2 EI 
   
 L   L  matriz
2 2+3
 2 EI   4EI  simétrica
ϕ3    
 L   L 
○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○

3 3
ϕ4  2 EI   4 EI 
0    
 L   L 

14
introduccion al calculo matricial de estructuras
4.3 matriz de rigidez. PÓRTICOS TRASLACIONALES
desplazamiento de un apoyo.MOMENTO DE EMPOTRAMIENTO LOCAL.

Si desplazamos uno de los apoyos de la viga una

cantidad, esta sufrirá una deformación.

Debido a ello aparecerán unos momentos que

justiquen la existencia de esos giros, luego estarán
provocando en los apoyos unos esfuerzos debidos
al desplazamiento además de los que puedan existir
debido a la propia carga.

- El giro en el nudo 1 será la suma de los producidos

por cada momento:
ML ML ML
ϕtotal = − =
3EI 6 EI 6 EI

Si despejamos el momento:
6 EI
M= *ϕ
L

Como estamos en la hipótesis de deformaciones

pequeñas,
∂
tgϕ = ϕ ⇒ ϕ =
L

Sustituimos en el momento:

6 EI ∂ 6 EI
M= * = 2
L L L
Debido al desequilibrio de momentos tiene que
existir un cortante:

M + M 12 EI
T= = 3 *∂
L L

15
introduccion al calculo matricial de estructuras
4.4Planteamiento de la matriz de rigidez en PORTICOS TRASLACIONALES

En este caso, los nudos 2 y 3 pueden girar,

como en el ejemplo anterior, y además se
DESPLAZAN.

Posibilidades de movimiento:

Giro del nudo 2 Giro del nudo 3 Desplazamiento del dintel

Formación de la MATRIZ DE RIGIDEZ:

○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○

ϕ1 ϕ2 ∂
1+ 2 2 1
 4 EI   2 EI   6 EI 
ϕ1      2 
 L   L   L 

2 2+ 3 3
ϕ2  2 EI   4 EI   6EI 
     2 
 L   L   L 
○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○

matriz  12 EI 
1+ 3

∂  3 
simétrica  L 

A continuación vemos como se han obtenido los valores de la matriz debidos al desplazamiento:

16
introduccion al calculo matricial de estructuras
Cortante producido en el dintel (barra 2) debido al giro del nudo 2

Aparece un CORTANTE debido al momento producido

por el giro, pero no se considera ya que no actúa en la
dirección del desplazamiento δ

Cortante producido en la barra 1 debido al giro del nudo 2

Para equilibrar los momentos en la barra, tiene

que existir un cortante, que será el que actúe
desplazando el dintel. En este caso la dirección
es en la del movimiento supuesto (hacia la
derecha), luego lo tomamos como positivo.
NOTA: sería como en el caso de un Cross
traslacional, donde en el estado 0 hemos
colocado las bielas pertinentes para impedir los
desplazamientos. En este caso actuaríamos
comprimiendo esa biela ficticia.

4 EI 2 EI
+
M1 + M 2 L1 L1 6 EI
T= = = 2
L L1 L1

Cortante producido en la barra 3 debido al giro del nudo 3

Es el mismo razonamiento:

4 EI 2 EI
+
M1 + M 2 L3 L3 6 EI
T= = = 2
L L3 L3

17
introduccion al calculo matricial de estructuras
 Cortantes producidos por el desplazamiento del dintel (barra 2)

Se producen los momentos de

empotramiento local, ya calculados, debidos
al desplazamiento de los nudos.

El valor del cortante es también conocido:

12 EI
T= *∂
L3

Para el cortante, en la matriz colocamos únicamente los coeficientes, ya que δ es una incógnita
que despejamos. En este caso, actúan dos cortantes en el desplazamiento, luego el término que irá en
la matriz será la suma de ambos:
1+ 3
12 EI 12 EI  12 EI 
T = 3 + 3 = 3 
L1 L3  L 

Es decir, que en este término actuará el cortante de todas las barras que se desplazan.
(NOTA: mucho cuidado con los signos).

Finalmente, teniendo en cuenta que la matriz de rigidez es simétrica:

 1+ 2 2
 6 EI  
1
  4 EI   2 EI 
   
 L   L   L  
  2 EI 2  4 EI 
2 +3
 6 EI  
3

     
 L   L   L  
  6 EI 1  6 EI 
3
 12 EI  
1+ 3

      3  
  L   L   L  

18
introduccion al calculo matricial de estructuras
4.5 matriz de rigidez. APLICACIÓN PRÁCTICA EN EL CASO DE PÓRTICOS TRASLACIONALES

En primer lugar hemos numerado todos los nudos y las barras.

Estamos en el caso de un pórtico traslacional, donde existe un posible desplazamiento δ del
dintel. Veamos los movimientos de la estructura:
- Los nudos 1, 2 y 3 son empotramientos, luego por definición tienen impedidos todos los
movimientos.
- Los nudos 4, 5 y 6 pueden girar.
- Como además se ha dicho, hay un desplazamiento del dintel.

NOTA: para hallar los términos de la matriz de rigidez es conveniente analizar los movimientos
uno a uno como en el caso anterior, y viendo como actúan, ya rellenar la matriz.

En este caso queda:

ϕ4 ϕ5 ϕ6 ∂
○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○

1+ 4 4 1
 4 EI   2 EI   6 EI 
ϕ4     0  2 
 L   L   L 

4 2+ 3 5 2
ϕ5  2 EI   4 EI   2 EI   6EI 
       2 
 L   L   L   L 

5 3+5 3
ϕ6  2 EI   4 EI   6EI 
0      2 
 L   L   L 
○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○
1+ 2 +3
 12 EI 
∂ este término es SIMÉTRICO  3 
 L 

19
introduccion al calculo matricial de estructuras
 Si consideramos la misma inercia y el mismo material (mismo E):

 1 .8 0 .4 0 0.375 
 
 0 .4 2.47 0.33 0.375 
 0 0.33 1.67 0.375 
 
 0.375 0.375 0.375 0.5625 

Cálculo del TÉRMINO INDEPENDIENTE:

En el caso de los términos correspondientes a los giros ϕ 4 , ϕ 5 , ϕ 6 son como ya sabemos la

suma de los momentos exteriores menos los de aquellos producidos por las cargas.

∑M E − ∑ Mˆ i
Entonces:
) )
nudo 4 ⇒
q1 L21 q2 L24
∑ M E − ∑ Mˆ i = 0 −− (
= 0 − − 4 + 4.16 = −0,16 )
12 12
) )
(
nudo 5 ⇒ ∑ M E − ∑ Mˆ i = 0 − 6 − 4.16 = −1,83 )
nudo 6 ⇒ ∑M E − ∑ Mˆ i = 0 − (− 6 ) = +6

NOTA: para ver mejor estos cálculos, conviene dibujar la estructura con los MEMP, lo que permite
una visión más rápida:

Los MEMP son los debidos a la carga distribuída:

qL2
ME =
12
La carga puntual no produce momentos en las
barras al estar aplicada en un nudo.

En el caso del término independiente debido al desplazamiento, corresponderá a la suma de los

cortantes que influyen en el movimiento de la estructura:

∑T = ∑T − ∑TE i

Es decir, el sumatorio de las fuerzas exteriores actuando en la dirección del movimiento menos el
sumatorio de los cortantes debidos a las cargas exteriores en la dirección del movimiento.
Para ver mejor como actúa, los dibujamos:

qL 3 * 4
T1 = = =6
2 2

El signo negativo es porque hemos supuesto

el movimiento δ hacia la derecha.

20
introduccion al calculo matricial de estructuras
Finalmente, nuestro sistema de ecuaciones quedará como:

 1 .8 0 .4 0 0.375   ϕ4   0 − 0,17 
     
 0 .4 2.47 0.33 0.375   ϕ5   0 − 1,8 
EI 
0 0.33 1.67 0.375  ϕ  =  0 + 6  =
   6   
 0.375 0.375 0.375 0.5625   ∂   3 − (− 6 )

Solución:

3,74 0,476
ϕ4 = − ϕ6 = −
EI EI
3,24 20,975
ϕ5 = − ∂=
EI EI

Con estos datos ya podemos hallar los momentos:

)  4 EI  4  2 EI 
4

M 45 = M +   * ϕ4 +   * ϕ5
 L   L 
 4 EI 3,74   2 EI 3.24 
M 45 = 4,17 +  − *  + − * 
 5 EI   5 EI 

M 45 = 4,17 − 2,992 − 1,296 = −0,118mT

 2 EI 3,74   4 EI 3.24 
M 54 = −4,17 +  − *  + − *  = − 8,258mT
 5 EI   5 EI 

 4 EI 3,24   2 EI 0,476 
M 56 = 6 +  − *  + − *  = 3,681mT
 6 EI   6 EI 

 2 EI 3,24   4 EI 0,476 
M 65 = −6 +  − *  + − *  = −7,397mT
 6 EI   6 EI 

Para el momento M52 podemos hacerlo de dos formas:

- equilibrando el nudo:

-8,258 + 3,681 + M52= 0 ; M52= 4,577mT

- o como los casos anteriores:

 4 EI 3,24   6 EI 
M 52 = 0 +  − *  +  − 2 * ∂  = 4 ,625 mT
 L EI   L 

21
introduccion al calculo matricial de estructuras
 El momento en los APOYOS será: (tenemos en cuenta que ϕ 1 , ϕ 2 , ϕ 3 son cero, por definición
de empotramiento)

2 EI 4 EI 6 EI
M 1 = Mˆ + * ϕ4 + * ϕ1 + *∂ =
L1 L1 L2
2 EI 3,74 4 EI 6 EI 20,975
= 4− * + *0 + * = 9,995 ≈ 10mT
4 EI L1 16 EI

2 EI 3,24 6 EI 20,975
M2 = 0 − * + * = 6,245mT
4 EI 16 EI

2 EI 0,476 6 EI 20,975
M3 = 0 − * + * = 7,627mT
4 EI 16 EI

La gráfica de momentos será:

22
introduccion al calculo matricial de estructuras