RESOLUCIÓN DE LOSAS 1

Resolución de Losas
de Hormigón Armado mediante el
Método de las Diferencias Finitas

Ing. Sebastian Stefanini

U.T.N. – F.R.B.A.

Estructuras de hormigón

Agosto 11 de 2023

Sebastian Stefanini
RESOLUCIÓN DE LOSAS 2

Índice
Método de las diferencias finitas .......................................................................................... 3

Resolución de la ecuación diferencial de la elástica de una viga por el método de las

diferencias finitas ............................................................................................................................. 5

Planteo del método de las diferencias finitas con funciones de dos variables ..................... 7

Resolución de la ecuación diferencial de la elástica de una losa en un punto de su dominio

........................................................................................................................................................ 10

Tablas resumen de valores significativos ........................................................................... 16

Bibliografía......................................................................................................................... 16

Sebastian Stefanini
RESOLUCIÓN DE LOSAS 3

Método de las Diferencias Finitas

Dada una función y  y ( x ) suficientemente diferenciable y

dados tres puntos de coordenadas x  a  x , x  a , x  a  x pertenecientes al dominio

de la función.

De acuerdo al teorema de Taylor y siempre que x  0

1
y ( a  x )  y ( a )  y( a )  x  y( a )  x 2
2

1
y ( a  x )  y ( a )  y( a )  x  y( a )  x 2
2

Restando miembro a miembro

y ( a  x )  y ( a  x )  2 y ( a )  x

y ( a  x )  y ( a  x )
y ( a )  (1)
2 x

Sumando miembro a miembro

y (a  x)  y (a  x)  2 y (a)  y(a)  x 2

y ( a  x )  y ( a   x )  2 y ( a )
y (a )  (2)
x 2

Aplicando la expresión (1)

y ( a  x )  y (a  x)

y (a ) 
2x

Reemplazando por la expresión (2)

y ( a  2  x )  y ( a )  2 y ( a   x ) y ( a )  y ( a  2 x )  2 y ( a   x )

y (a ) 
  x 2
x 2
2 x

y ( a  2 x )  2  y ( a  x )  y ( a   x )   y ( a  2 x )
y ( a )  (3)
2 x 3

Aplicando la expresión (2)

y ( a  x )  y (a  x)  2 y (a )

y ( IV ) ( a ) 
x 2

y ( a  2 x )  4 y ( a   x )  6 y ( a )  4 y ( a  x )  y ( a  2 x )
y ( IV ) ( a )  (4)
x 4

Sebastian Stefanini
RESOLUCIÓN DE LOSAS 4

Se pueden identificar los valores de la función y en las proximidades del punto x  a con la
siguiente notación:

y (a  2x)  yi  2

y (a  x)  yi 1

y (a )  yi

y (a  x)  yi 1

y (a  2x)  yi  2

Resultando

yi 1  yi 1
y ( a ) 
2 x

yi 1  yi 1  2 yi
y (a ) 
x 2

yi  2  2  yi 1  yi 1   yi  2
y ( a ) 
2x 3

yi  2  4 yi 1  6 yi  4 yi 1  yi  2
y ( IV ) ( a ) 
x 4

Esta notación se puede utilizar para resolver en forma aproximada algunas ecuaciones
diferenciales como la de la elástica de una viga

d4y q

dx 4 EI

Sebastian Stefanini
RESOLUCIÓN DE LOSAS 5

Resolución de la Ecuación Diferencial de la Elástica de una Viga por el Método de las

Diferencias Finitas

Particionando la longitud de una viga en n partes iguales, quedan determinados los puntos de
coordenadas  x0 , x1 ,..., xn1 , xn  tales que x  xi  xi 1 para i  1,..., n

Para una viga apoyada en sus extremos se verifica que y0  yn  0

Si la viga se encuentra simplemente apoyada en el borde izquierdo se verifica que M 0  0 o

sea

y0 1  y01  2 y0
y (0)  0  y0 1   y0 1
x

Análogamente en el borde derecho.

Si la viga se encuentra empotrada en el borde izquierdo se verifica que

y01  y0 1
y (0)   0 o sea y01  y0 1
2 x

Planteando la ecuación diferencial en cada uno de los puntos  x1 ,..., xn1 resulta en un

sistema de n  1 ecuaciones lineales con n  1 incógnitas.

Resolviendo dicho sistema se obtienen las deformaciones de la viga en los puntos

 x1,..., xn1
Por ejemplo, dada la ecuación diferencial de la elástica de una viga

d4y q

dx 4 EI

Planteando la ecuación en un punto de coordenada i

yi  2  4 yi 1  6 yi  4 yi 1  yi  2 q

x 4
EI

Despejando la incógnita yi

q
4( yi 1  yi 1 )  yi  2  yi  2  x 4
yi  EI
6

En una hoja de cálculo se resuelve fácilmente el sistema de n  1 ecuaciones lineales con

n  1 incógnitas.

Una vez calculados los valores de la elástica y dado que

Sebastian Stefanini
RESOLUCIÓN DE LOSAS 6

d2y d3y
M   EI y V   EI
dx 2 dx 3

Entonces

yi 1  yi 1  2 yi yi  2  2  yi 1  yi 1   yi  2
M i   EI y Vi   EI
x 2 2 x 3

Ejemplo 1

Determinar la elástica, los esfuerzos de flexión y corte en una viga de rigidez unitaria EI  1
, carga distribuida unitaria q  1 y longitud unitaria lc  1 por el método de las diferencias finitas
suponiendo:

a) Viga simplemente apoyada

b) Viga apoyada – empotrada

c) Viga empotrada – empotrada

Adoptar una partición de la luz de la viga igual al 10% de la luz

d) Comparar las flechas obtenidas con los valores exactos.

e) Comparar los esfuerzos de flexión máximos obtenidos con los exactos.

Sebastian Stefanini
RESOLUCIÓN DE LOSAS 7

Planteo del Método de las Diferencias Finitas con Funciones de Dos Variables

Dada una función z  w( x, y ) suficientemente diferenciable y

dados cinco puntos de coordenadas ( x, y)  (a  x, b) , ( x, y)  (a, b) , ( x, y)  (a  x, b) ,

( x, y)  (a, b  y) , ( x, y)  (a, b  y) pertenecientes al dominio de la función.

De acuerdo al teorema de Taylor y siempre que x  0

1
w( a  x, b)  w( a, b)  wx ( a, b)  x   wxx ( a, b)  x 2 
2

1
w( a  x, b)  w( a, b)  wx ( a, b)  x   wxx ( a, b)  x 2 
2

Restando miembro a miembro

w(a  x, b)  w( a  x, b)  2 wx (a, b)  x

w( a  x, b)  w( a  x, b)
wx (a, b)  (1)
2 x

Sumando miembro a miembro

w( a  x, b)  w( a  x, b)  2 w(a, b)  wxx ( a, b)  x 2

w( a  x, b)  w( a  x, b)  2 w( a, b)
wxx (a, b )  (2)
x 2

Análogamente

w(a, b  y )  w(a, b  y )
wy (a, b)  (3)
2y

w( a, b  y )  w(a, b  y )  2 w(a, b)
wyy (a, b)  (4)
y 2

Sebastian Stefanini
RESOLUCIÓN DE LOSAS 8

Además

wx ( a, b  y )  wx (a, b  y )

wxy (a, b) 
2y

Reemplazando (1)

w(a  x, b  y )  w( a  x, b  y ) w(a  x, b  y )  w( a  x, b  y )


wxy (a, b)  2 x 2x
2 y

w(a  x, b  y )  w(a  x, b  y )  w(a  x, b  y )  w(a  x, b  y )

wxy (a, b)  (5)
4xy

De la ecuación (1)

wxx (a  x, b)  wxx

 (a  x, b)
xxx ( a , b ) 
w
2 x

w( a  2x, b)  w( a, b)  2 w( a  x, b) w(a, b)  w(a  2x, b)  2 w(a  x, b)


w ( a , b )   x 2
x 2
xxx
2x

w( a  2x, b)  2 w( a  x, b)  w( a  2x, b)  2 w( a  x, b)

xxx ( a , b ) 
w (6)
2 x 3

Análogamente

w( a, b  2y )  2 w( a, b  y )  w( a, b  2y )  2 w(a, b  y )

yyy ( a , b ) 
w (7)
2 y 3

De la ecuación (3)

wxx (a, b  y )  wxx (a, b  y )

xxy ( a , b ) 
w
2y

w( a   x , b   y )  w( a   x , b   y )  2 w( a , b   y )
x 2
xxy ( a , b ) 
w 
2 y

w(a  x, b  y )  w(a  x, b  y )  2 w( a, b  y )

 x 2
2 y

w(a  x, b  y )  w(a  x, b  y )  2 w(a, b  y )

xxy ( a , b ) 
w 
2x 2 y

 w(a  x, b  y )  w(a  x, b  y )  2 w(a, b  y )

 (8)
2x 2 y

Sebastian Stefanini
RESOLUCIÓN DE LOSAS 9

Análogamente

w(a  x, b  y )  w(a  x, b  y )  2 w(a  x, b)

yyx ( a , b ) 
w 
2xy 2

 w(a  x, b  y )  w(a  x, b  y )  2 w(a  x, b)

 (9)
2xy 2

De la ecuación (6)

wxx (a  x, b)  wxx (a  x, b)  2 wxx ( a, b)

( IV )
wxxxx (a, b) 
x 2

w(a  2x, b)  w(a, b)  2 w(a  x, b) w(a, b)  w(a  2x, b)  2 w(a  x, b)

( IV )
wxxxx ( a, b )  x 2 x 2 
x 2

w(a  x, b)  w(a  x, b)  2 w(a, b)

2
 x 2
x 2

w( a  2x, b)  4 w( a  x, b)  6 w( a, b)  4 w( a  x, b)  w(a  2x, b)

( IV )
wxxxx ( a, b)  (10)
x 4

Análogamente

w(a, b  2y )  4 w(a, b  y )  6 w(a, b)  4 w(a, b  y )  w(a, b  2y )

w(yyyy
IV )
( a, b)  (11)
y 4

De la ecuación (6)

wxx (a, b  y )  wxx

 (a, b  y )  2 wxx (a, b)
( IV )
wxxyy ( a, b) 
y 2

w( a   x , b   y )  w( a   x , b   y )  2 w( a , b   y )

wxxyy ( a, b) 
( IV )  x 2

y 2

w( a   x , b   y )  w( a   x , b   y )  2 w( a , b   y ) w( a   x , b )  w( a   x , b )  2 w( a , b )
 2
x 2
x 2
y 2
w(a  x, b  y )  w(a  x, b  y )  w(a  x, b  y )  w(a  x, b  y )
( IV )
wxxyy ( a, b)  
x 2 y 2

 2  w( a   x , b )  w( a   x , b )  w( a , b   y )  w( a , b   y )   4 w( a , b )
 (12)
x 2 y 2

Sebastian Stefanini
RESOLUCIÓN DE LOSAS 10

Resolución de la Ecuación Diferencial de la Elástica de una Losa Delgada en Un Punto de su

Dominio

La ecuación diferencial de la elástica de una losa es:

p
( IV )
wxxxx (a, b)  2 wxxyy
( IV )
( a, b)  wyyyy
( IV )
(a, b) 
B

De las ecuaciones (10), (11) y (12)

w( a  2x, b)  4 w( a  x, b)  6 w(a, b)  4 w( a  x, b)  w( a  2x, b)


x 4

w(a  x, b  y )  w(a  x, b  y )  w(a  x, b  y )  w(a  x, b  y )

2 
x 2 y 2

2  w(a  x, b)  w(a  x, b)  w( a, b  y )  w(a, b  y )   4 w(a, b)

2 
 x 2 y 2

w(a, b  2y )  4 w(a, b  y )  6 w(a, b)  4 w(a, b  y )  w(a, b  2y ) p

 
y 4 B

Despejando z (a, b)

6 w( a, b) w(a  2x, b)  4 w(a  x, b)  4 w(a  x, b)  w( a  2x, b)

 
x 4 x 4

w(a  x, b  y )  w(a  x, b  y )  w(a  x, b  y )  w(a  x, b  y )

2 
x 2 y 2

8w( a , b ) 2  w(a  x, b)  w(a  x, b)  w( a, b  y )  w(a, b  y ) 

 2 
x y
2 2
x 2  y 2

6 w( a, b) w(a, b  2y )  4 w(a, b  y )  4 w( a, b  y )  w( a, b  2y ) p

  
x 4 x 4 B

6w(a, b)y 4   w(a  2x, b)  4w(a  x, b)  4w(a  x, b)  w(a  2x, b) y 4 

2  w(a  x, b  y )  w(a  x, b  y )  w(a  x, b  y )  w(a  x, b  y )  x 2 y 2 

 8 w ( a , b )  x 2  y 2  4  w ( a   x , b )  w( a   x , b )  w ( a , b   y )  w ( a , b   y )   x 2  y 2 

Sebastian Stefanini
RESOLUCIÓN DE LOSAS 11

p 4 4
6 w( a, b)x 4   w( a, b  2y )  4 w( a, b  y )  4 w( a, b  y )  w( a, b  2y ) x 4  x y
B

w(a, b) 6 y 4  8x 2 y 2  6x 4  

  w(a  2x, b)  4w(a  x, b)  4w(a  x, b)  w(a  2x, b) y 4 

2  w(a  x, b  y )  w(a  x, b  y )  w(a  x, b  y )  w(a  x, b  y )  x 2 y 2 

 4  w ( a   x , b )  w( a   x , b )  w ( a , b   y )  w ( a , b   y )   x 2  y 2 

p 4 4
  w( a, b  2y )  4 w( a, b  y )  4 w( a, b  y )  w( a, b  2y ) x 4  x y
B

w( a , b ) 

  w(a  2x, b)  4w(a  x, b)  4w(a  x, b)  w(a  2x, b) y 4



2  w(a  x, b  y )  w(a  x, b  y )  w(a  x, b  y )  w(a  x, b  y )  x 2 y 2 

 4  w ( a   x , b )  w( a   x , b )  w ( a , b   y )  w ( a , b   y )   x 2  y 2 

p 4 4
  w(a, b  2y )  4 w(a, b  y)  4w(a, b  y)  w(a, b  2y)  x 4  x y  
B 

1

 6y  8x 2 y 2  6x 4 
4

Si x  y

w( a , b ) 

  w(a  2x, b)  4w(a  x, b)  4w(a  x, b)  w(a  2x, b) x 4



2  w(a  x, b  y )  w(a  x, b  y )  w(a  x, b  y )  w(a  x, b  y ) x 4 

 4  w ( a   x , b )  w( a   x , b )  w ( a , b   y )  w ( a , b   y )   x 4 

p 8 1
  w( a, b  2y )  4 w(a, b  y )  4 w(a, b  y )  w( a, b  2y ) x 4  x  
 20x
4
B

Sebastian Stefanini
RESOLUCIÓN DE LOSAS 12

w( a , b ) 

w(a  2x, b)  8w(a  x, b)  8w(a  x, b)  w(a  2x, b) 

 2 w( a   x , b   y )  2 w( a   x , b   y )  2 w( a   x , b   y )  2 w( a   x , b   y ) 

p 4 1
 w( a, b  2y )  8w( a, b  y )  8w( a, b  y )  w( a, b  2y )  x  
B  20

La última función corresponde a la deformada de un punto de la losa en función de la

deformada de los puntos próximos, de la carga en dicho punto y de la rigidez de la losa.

Sebastian Stefanini
RESOLUCIÓN DE LOSAS 13

Para hallar la elástica en cada punto de la losa teniendo en cuenta las condiciones de borde
se procede del siguiente modo:

a) La losa se particiona en cuadrados de dimensión x

b) En el perímetro de la losa, supuesto que se encuentra apoyada en forma rígida se asigna
el valor nulo o sea
w( a, b)  0
c) Si la losa se encuentra simplemente apoyada o sea
w( a  x, b)  w( a  x, b)  2 w( a, b)
wxx ( a, b)  0  0 
x 2
w( a   x , b )   w( a   x , b )
Análogamente
w ( a , b   y )   w ( a , b  y )
d) Si la losa se encuentra empotrada o sea
w( a  x, b)  w( a  x, b)
wx (a, b)  0   0  w( a   x , b )  w( a   x , b )
2 x
Análogamente
w( a , b   y )  w( a , b   y )
e) En los puntos interiores de la losa y en cada vértice de dichos cuadrados se expresa el
valor de la elástica en función de los valores contiguos
f) Se resuelve el sistema de ecuaciones y se obtiene el valor de la elástica en el vértice de
cada cuadrado de la partición.

Seguido se proceden a calcular los momentos flexores en las distintas direcciones

 2w 2w 
mx   B   2   2 
 x y 

 w(a  x, b)  w(a  x, b)  2 w(a, b) w(a, b  y )  w(a, b  y )  2 w(a, b) 

mx ( a, b)   B    
 x 2
y 2 

Análogamente

 w(a, b  y )  w(a, b  y )  2 w(a, b) w(a  x, b)  w(a  x, b)  2 w(a, b) 

m y ( a, b)   B    
 y 2
x 2 

Sebastian Stefanini
RESOLUCIÓN DE LOSAS 14

Los momentos torsores

2w
mxy  (1   ) B 
xy

w(a  x, b  y )  w(a  x, b  y )  w(a  x, b  y )  w(a  x, b  y )

mxy  (1   ) B 
4xy

Los esfuerzos de corte

 3w 3w 
tx   B   3  2 
 x y x 

 w( a  2x, b)  2 w( a  x, b)  w( a  2x, b)  2 w(a  x, b) 

tx   B    ... 
 2 x 3


 w ( a   x , b   y )  w ( a   x , b  y )  2 w ( a   x , b ) 
t x   B   ...   ... 
 2 xy 2


  w( a   x , b   y )  w( a   x , b   y )  2 w( a   x , b ) 
t x   B   ...  
 2 x y 2 

Si x  y

 w( a  2x, b)  4 w( a  x, b)  w( a  2x, b)  4 w(a  x, b)  ... 

tx   B   
 2 x 3 

 ...  w( a  x, b  y )  w( a  x, b  y )  w( a  x, b  y )  w( a  x, b  y ) 

tx   B   
 2 x 3 

Análogamente

 w( a, b  2 y )  4 w( a, b  y )  w( a, b  2y )  4 w( a, b  y )  ... 
ty  B  
 2 y 3 

 ...  w(a  x, b  y )  w( a  x, b  y )  w(a  x, b  y )  w(a  x, b  y ) 

ty  B  
 2 y 3 

Sebastian Stefanini
RESOLUCIÓN DE LOSAS 15

Ejemplo 2

Determinar la elástica, los esfuerzos de flexión y corte en una losa rectangular de rigidez
unitaria B  1 , carga distribuida unitaria p  1 y cuyos lados miden 1 y 1.5 respectivamente por el
método de las diferencias finitas suponiendo:

a) Lados simplemente apoyados

b) Lados empotrados

Adoptar una partición de la luz de la viga igual al 10% de la luz menor

Sebastian Stefanini
RESOLUCIÓN DE LOSAS 16

Tablas Resumen de Valores para Diseñar

En las siguientes tablas se resumen coeficientes que permiten calcular fácilmente las flechas,
momentos máximos y reacciones en los apoyos en losas rectangulares con sus aristas simplemente
apoyadas o empotradas. Se adoptó una rigidez y una carga uniformemente distribuida unitaria. El
lado menor de la losa es unitario. El coeficiente de Poisson se adoptó igual a 0.16

A la longitud menor de una losa se la denomina a y a la longitud mayor se la denomina b

p  a4
Para calcular la flecha de una losa se multiplicará el coeficiente W por siendo p una
B
Ec h 3
carga uniformemente distribuida y B  su rigidez
12(1   2 )

Para calcular los momentos flexores multiplicará el coeficiente correspondiente por p  a 2

Para calcular las reacciones en los apoyos se multiplicará el coeficiente correspondiente por
p  a b

Bibliografía

Belluzzi, O. (1970). Ciencia de la construcción. Madrid: Aguilar S.A.

Kalmanok, A. S. (1961). Manual para cálculo de placas. Montevideo: Editora Inter Ciencia.

Sebastian Stefanini