PRINCIPIO DE LOS TRABAJOS VIRTUALES

f(x)
 dσ 1 d 2σ 2 
 σ + dx + 2
dx ....A
σA  dx 2 dx 

Ecuación diferencial

dσ du d  du 
Equilibrio f ( x )dx + dxA = 0 σ=E  E  A = −f ( x )
dx dx dx  dx 
dσ
A = −f ( x ) Para E=cte a lo largo de la barra
dx
Compatibilidad du EAu ' '+f ( x ) = 0
ε= d 2u
dx EA 2 = −f ( x )
dx
Ecuación constitutiva σ = Eε Ecuación diferencial barra con carga axial

Ahora haremos lo siguiente:

1) Multiplicar ambos lados de la ecuación por una cantidad arbitraria q(x)

[EAu' '+f ( x )]q(x ) = 0

2) Integrar sobre el dominio del problema, de 0 a L
L L L

∫ [EAu ' '+f ( x )]q( x )dx = 0

0
⇒ ∫ EAu ' ' q( x )dx + ∫ f ( x)q( x )dx = 0
0 0

RSM/2009 1
3) Integrar por partes: ∫ rds = − ∫ sdr +rs
ds = u ' ' dx → s = u '
r = q → dr = q' dx
L L
− ∫ EAu ' q' ( x )dx + EAu ' q ( x ) 0 + ∫ f ( x )q ( x )dx = 0
L

0 0
Término asociado a condiciones de borde.
L L Arbitrariamente se define como nulo.
− ∫ EAu ' q ' ( x )dx + ∫ f ( x )q ( x )dx = 0
0 0

EAu ' q ( x ) 0 = EA[u ' (L)q ( L) − u ' (0)q(0)] = 0

L

d
Definamos q=δu, una función de desplazamiento arbitraria. Por lo tanto, q' = δu = δε
dx
Notar que du
u' = =ε y Eu’=σ L L

∫ σAδεdx = ∫ f (x )δudx + Fδu 0

L
dx
0 0

Principio de los trabajos virtuales: Si a una estructura deformable que está en equilibrio
bajo la acción de un grupo de cargas (f(x)) se le da una pequeña deformación virtual (δu),
entonces el trabajo hecho por las cargas externas actuando sobre los desplazamientos
L
virtuales ( ∫0 f ( x)δudxvirtual hecho porL las fuerzas internas actuando sobre
) es igual al trabajo
las deformaciones virtuales ( ). ∫ σAδεdx
RSM/2009 0 2
Este principio es independiente de la relación constitutiva del material. Lo podemos re-
derivar de la ecuación de equilibrio sin usar la ley de Hooke.
L
 dσ  ∫ rds = −∫ sdr +rs L L
∫0 
 f ( x ) + A q( x )dx = 0
dx  dσ − ∫ σAq' ( x )dx + q ( x )σA 0 + ∫ f ( x )q( x )dx = 0
L

ds = dx → s = σ 0 0
dx
L L
r = q → dr = q' dx
∫ σAδεdx = ∫ f ( x )δudx + Fδu 0
L

0 0

L
Nos queda resolver el problema con el término Fδu 0 = Aσδu L − Aσδu 0

Supongamos las siguientes condiciones de borde: f(x) F

Si F=0, (Aσδu) L =0 ya que Aσ=0. Luego, ; (Aσδu)0 =0 también debe ser 0

→ δu =0 donde los desplazamientos reales son conocidos.
Ahora si F≠0, (Aσδu)L – (Fδu)L =0 → (Aσ – F)δuL=0, pero Aσ – F=0 por equilibrio.
Luego, δu es arbitrario donde el esfuerzo es conocido.
Por lo tanto, si se especifica un desplazamiento, no se puede prescribir una
fuerza, y si especifica una fuerza, el desplazamiento es incógnita.

RSM/2009 3
 ENERGÍA POTENCIAL DE UNA BARRA CARGADA AXIALMENTE

Energía de deformación interna

Trabajo externo

Energía potencial total

La energía potencial de una barra cargada axialmente es un funcional que solo depende
de u(x), que es la variable primaria

1
Π = U −W =
2 ∫ EAu '2 dx − ∫ f ( x )udx 4
 CONCEPTO DE VARIACION CINEMATICAMENTE ADMISIBLE

Supongamos que u(x) cambia a u(x)+δu(x). El funcional cambia de Π a Π+δΠ. La función δu(x) y el
escalar δΠ se denominan variaciones de u(x) y Π, respectivamente. La variación δu(x) no debe
confundirse con la derivada ordinaria du(x)=u´(x)dx ya que al tomar la variación la variable
independiente queda fija.
Una variación de desplazamiento δu(x) se dice admisible cuando ambos u(x) y u(x)+δu(x) son
cinemáticamente admisible en el sentido del principio de los trabajos virtuales. Un desplazamiento
axial cinemáticamente admisible obedece las siguientes dos condiciones:
1) Es continuo sobre el largo de la barra, esto es, u(x) es de clase C0 para 0≤x≤L.
2) Satisface exactamente las condiciones de borde de desplazamiento, en la figura por ejemplo,
u(0)=0.

5
Cálculo variacional
Considerar la función F(x,u,u’). Para un valor fijo arbitrario de la variable independiente x, F depende de u
y u’. El cambio εv en u, donde ε es constante y v es una función se llama la variación de u y se denota
por δu: δu≡ εv. El operador δ se llama operador variacional. La variación δu de una función u representa
un cambio admisible en la función u(x) a un valor fijo de la variable independiente x. Si u se especifica en
un punto (usualmente en la frontera), la variación de u es cero porque el valor especificado no puede
cambar. Por lo tanto, la variación de una función u debe satisfacer las condiciones de contorno. La
variación δu representa un cambio virtual pero admisible en u. Asociado a este cambio en u (es decir,
cuando u cambia en u+εv), hay un cambio en F,
∆F=F(x, u+ εv, u’+ εv’)-F(x, u, u’)

Primer variacional de F

Notar la analogía entre el primer variacional y el diferencial total.

Ya que x no varía durante la variación de u a u+δu, dx=0 y por lo tanto δ actúa como un operador diferencial
con respecto a las variables dependientes.
RSM/2009 6
 PRINCIPIO DE ENERGIA POTENCIAL
Para materiales hiperelásticos, existe una función de densidad de energía de deformación U0
tal que:
∂U 0 U0 es una función única de las deformaciones y se asume como definida positiva.
σ ij =
∂ε ij
El PTV se puede expresar en términos de U0:
∂U 0
∫Ω ∂ε ij δε ij dv − [∫ ρf δu dv + ∫ T δu ds] = 0
Ω
i i
S
i i

La primera integral es igual a ∫

Ω
σ ij δε ij dv = ∫ δU 0 dv = δU
Ω

 εij  1
U = ∫ U 0 dv = ∫  ∫ σ ij dε ij dv = ∫ σ ij ε ij dv U: energía de deformación
Ω Ω  2 Ω
0 
Válido para material lineal elástico
La expresión en corchetes se puede identificar como –δW E. Se supone que existe un
potencial V=W E cuya primera variación es: δU = δW = − ρf δu dv + T δu ds
E [∫Ω i i ∫
S i i ]
El PTV toma la forma: δU + δV = δ(U + V ) = δΠ = 0
En que Π se conoce como energía potencial total. El principio de mínima energía potencial
establece que de todos los desplazamientos posibles, aquellos que satisfacen las ecuaciones
de equlibrio minimizan la energia potencial.

7
Energía Potencial de una barra cargada axialmente: Π = 1 ∫ EAu ' 2 dx − ∫ f ( x )udx
2
1
Π=
2 ∫ EAu ' 2 dx − ∫ f ( x )udx Minimizar Π es equivalente a satisfacer el PTV en
estructuras elásticas lineales.
EA L L

2 ∫0 ∫0 f ( x)δudx
δΠ = 2 u ' δu ' dx −
Como usamos esto para resolver problemas reales?
L L
δΠ = ∫ σAδεdx − ∫ f ( x )δudx = 0
0 0
barra 1 barra 2
Fi L1, (EA)1, ε1 L2, (EA)2, ε2
2
1 
Π = ∑  ∫ EAε 2 dx  − ∑ Fi u i Fk
i =1  2  nodo i Fj
nodo j nodo k
 Fi  ui 
FFFF

;;;;

uuuu

=  Fj  =  u j 

uuuu
BBBB
Bu

uuuu
FFFF
1 

TTTT
TTTT

TTTT
2 Li

Π = ∑  EA ∫

( )dx  −
Fk   u k  i =1  2 0 

uuuu
BBBB
Bu TTTT

uuuu
BBBB
BBBB
uuuu

uuuu
FFFF
2  

TTTT
TTTT

TTTT
TTTT

TTTT
Li
ui  δΠ = ∑ EA ∫ δ
 1
( )
δ dx  − δ
BBBB 1111
uuuu

+
 
ε1 =
1
[− 1 1 0] u j  = i =1  2 0 
uuuu
BBBB
BBBB
uuuu
BBBB
uuuu
uuuu
BBBB

uuuu
BBBB
Bu
TTTT
TTTT

TTTT

TTTT
TTTT
L1
 u k  δ =( δ )T ( ) T
=δ
uuuu

BBBB
BBBB

uuuu
FFFF
 2  Li  
TTTT

TTTT
ui  ∑  EA ∫ ( )
BBBB
uuuu

δΠ = δ dx  −  = 0 Pero δu es arbitrario
1
[0 1 − 1] u j  = 
2222

ε2 =  i =1  0  
L2   Para una barra
 u k 

BBBB
2 
BBBB
uuuu iiii
BBBB
uuuu iiii
uuuu
BBBB
BBBB
uuuu iiii

BBBB
BBBB

uuuu
FFFF
TTTT

  EA ( 
TTTT

Li

)  − =0 = 1 [− 1 1]
(ε1 )2 = ε T ε = ( )T ( )= ∑ ∫0
T
iiii

dx L
 i =1   
 

BBBB
BBBB
KKKK

 14442444 

TTTT
3 1 − 1  1 − 1
  = 2  [− 1 1] =  
L 1 − 1 1 

kkkk
eeee
EA  1 − 1
RSM/2009 = 8
L − 1 1 
Viga de Euler-Bernoulli

d4u Ecuación de la elástica

u i  EI ⋅ =q
v  dx 4
 i

uuuu
 θi  d2v
=  φ= Curvatura
dx 2
u j 
v j  1
2∫
  Πb = σεdv Energía potencial interna
 θ j 
d2v
σ = Eε ε = −φy = − 2 y
dx
1
Πb =
2 ∫ E ( v' ' ) 2 y 2 dAdx
123 E, v’’ son independientes del área
dv
Π c arg as = ∫ q( x)vdx Por el momento no
consideramos el término
1 L de energía axial pues está
Π = Π b − Π c arg as =
2 ∫0
EI ⋅ ( )
v '' 2
⋅ dx −
L
∫0 q( x) ⋅ v ⋅ dx desacoplado de la flexión

Si EI es constante,
1
Π=
L
( )
2 L
EI ∫ v '' ⋅ dx − ∫ q ( x ) ⋅ v ⋅ dx
2 0 0

L
( ) L
δΠ = EI ∫ v '' δv' '⋅dx − ∫ q ( x ) ⋅ δv ⋅ dx
0 0

Necesitamos escribir los desplazamientos de la forma N(x)=f(vi,θi, vj,θj). Supongamos: v(x)=ax3+bx2+cx+d

(mismo orden que solución de la ecuación diferencial de la elástica). v puede diferenciarse dos veces asi
que funcionará en δπ=0. Tenemos 4 desplazamientos nodales y cuatro coeficientes.
RSM/2009 9
 Viga de Euler-Bernoulli
v(0)=vi → d= vi; v’(0)=θi → c= θi; v(L)=vj ; v’(L)=θj
Resolviendo para el resto de las constantes y reordenando se obtiene:
 x
3
x 
2
  x 
2
  x 2 x 
3
 x 2 x 
v( x ) = 1 + 2  − 3   v i + x 1 −   θi + 3  − 2   v j + x   −  θ j
 L  L     L     L   L    L  L 
v( x ) = N1 ( x ) vi + N 2 ( x )θi + N 3 ( x ) v j + N 4 ( x )θ j
ui 
v  Las funciones N se denominan Polinomios de Hermite

NNNN
 i

uuuu

uuuu
Nu
 θi  dV d
v( x ) = [0 N1 N2 0 N3 N 4 ]  = v' = = = [0 N '1 N '2 0 N'3 N '4 ]
u j 

uuuu
dx dx
v j  v' ' = [0 N ' '1 N' '2 0 N ' '3 N ' '4 ]
 
 θ j 
N' N' uuuu
uuuu '''' N'
''''
N'
'''' N' '
uuuu '''' ''''

NNNN
uuuu
L L
δΠ = EI ∫ ( )( δ ) ⋅ dx − ∫0 q( x ) ⋅ ( δ ) ⋅ dx = 0
0
uuuu ''''

uuuu uuuu

NNNN
uuuu

L L
δΠ = EI ∫ ( δ )T ( ) ⋅ dx − ∫0 ( δ )T q( x ) ⋅ dx = 0
0
N

uuuu
NNNN
TTTT

TTTT

TTTT

L L
δΠ = EI ∫ δ ⋅ dx − ∫ δ T
q ( x ) ⋅ dx = 0
0 0
uuuu

N'
''''
N'
''''
uuuu
----
uuuu

NNNN
TTTT

TTTT

TTTT

EI L dx ( ) L
∫ ∫ 42
T
δ δ q ( x )dx = 0
 0 
3
KKKK

0
1 44 42444 1 43
f
Ku
----
ffff
0000

=
RSM/2009 10
POLINOMIOS DE HERMITE

 x
3
x 
2
  x 2    x 2 x 
3
 x 2 x 
v( x ) = 1 + 2  − 3   vi +  x 1 −   θi + 3  − 2   v j + x   −  θ j
 L  L     L     L   L    L  L 
v( x ) = N1 ( x ) v i + N 2 ( x )θi + N 3 ( x ) v j + N 4 ( x )θ j

RSM/2009 11
  6x
2
6x     12 x 6 
   −    12  x  6
LL LL    −   −
  L2  L  L2 L3 L2 
   
x    6x − 4 
2
x  2x   2 x  2x 2  x
 1 −  − 1 −  

− 1 −  + 2 − 1 − 
  L2 L 
  L L  L  L L L L L
N' = 2 ; N' ' =   = 
     
 −  3 − 2 
6x 6x 6 12 x 12 x 6
  −     − 2 
   L L 
2
 LL LL   L L L

  x   2  x  1   2  x  1   2    6x − 2 
x
2
 2  x  1     −  +    −  + x  2    L2 L 
  L  − L  + x  L  L  − L    L  L  L   L  L  L   L 
        
L 2
 12x 6  12EI
Primer término de la matriz de rigidez: K11 = EI ∫  3 − 2  dx = ... = 3
0
L L  L
Tarea: Continuar con la operatoria y verificar que se obtiene la matriz de rigidez conocida del
análisis de estructuras:

NNNN
L
∫
T
Verificar además que de la integral q(x )dx
0

se obtiene el vector de fuerzas nodales equivalentes visto

en análisis de estructuras, para el caso de carga
uniformemente distribuida.
Como se puede ocupar las funciones de forma para
obtener el vector de fuerzas nodales equivalentes para el
caso de una carga puntual?

RSM/2009 12
RSM/2009 13
RSM/2009 14
 EFECTOS DE LA TEMPERATURA

∆L=α ∆TL ∆T

∆T: gradiente uniforme de temperatura εi= α ∆T

εi= α ∆T= ∆L/L; σ=0 Deformación real=0

σ=E εi =Eα ∆T
L

∫ σAδεdx = δu P + δu P
0
1 1 2 2 Trabajo virtual interno
L L L

σ = E(ε − εi ) εi= deformación inicial ∫ EA(ε − ε )δεdx = ∫ EAεδεdx − ∫ EAε δεdx

0
i
0 0
i
Bu

u 2 − u1 Para gradiente uniforme de temperatura,

ε= =
L
BBBB
uuuu

uuuu
BBBB
TTTT

δε = ( δ ) = δ
T T

[− 1 1]T dx = 
L
− EAα∆T 
uuuu

BBBB
BBBB

uuuu
----
uuuu

BBBB

L  L T  1
Fi = EAα∆T ∫
TTTT

TTTT

EA ∫ dx  EA ∫ ε dx  
iiii

T
Wi = δ δ
0  4 L  EAα ∆ T 
14 0
2443 0
Fi
uuuu
KuFFFF iiii

FFFF iiii

uuuu
PPPP
TTTT

TTTT

TTTT

Wi = δ [ − ]= δ
Ku

PPPP

= +
RSM/2009 15
 GRADIENTE TERMICO ATRAVES DEL ESPESOR EN VIGAS
εt=α ∆T/2 Gradiente uniforme de temperatura ∆T entre fibra
superior en inferior
φi
d

εb=-α
∆T/2

Wint = EI ∫ (v '' )δv' '⋅dx

L

εt − εb α ∆T
φi = − =−
L d
  α ∆T    α ∆T 
M = EI(v' '−φi ) = EI v' '− −   = EI v' '+ 
  d   uuuu d 

N'
''''
N'
''''
uuuu

N'
''''
TTTT

TTTT

TTTT
L α ∆T EIα∆T L
δv' '⋅dx = δ EI ∫ ( )dx 
L L
Wint = EI ∫ v' ' δv' '⋅dx + EI ∫ + ∫ ⋅dx
0 0 d  0 d 0
144 424443
Ku

FFFF
----
FFFF iiii

 0  Vector de fuerzas nodales equivalentes

=  
0
 EIα∆T 
− 
Fi =  d 
 0 
 0 
 EIα∆T 
 
 d 
RSM/2009 16
 EJEMPLO DE OBTENCION DE MATRIZ DE RIGIDEZ ELEMENTO DE VIGA
POR CONDENSACION ESTATICA

 12 EI 6 EI 12 EI 6 EI 
−  qL   6 2 6 4  2
 L3 L2 L3 L2  v EI  2 v1 + θ 1 − 2 v2 + θ 2  − qL =0
 
2 EI   1   qL2
 6 EI 2 L L L L  12
4 EI 6 EI 
 − 2  θ  qL2 6 2 6
 L
2
L L L  1  =  12  − 2 v1 − θ 1 + 2 v 2
6 EI    qL  12 EI L
− 12 EI 6 EI
− 2
12 EI
− 2   v 2   θ2 = L L
 L3 L L3 L θ 2   2  4
 6 EI 2 EI 6 EI 4 EI  − qL 
2
L
 − 2   12
 L2 L L L 

12 EI 6 EI 12 EI 6 EI L  qL2 6 2 6 
 L3 1 v + θ − v +  − 2 v1 − θ1 + 2 v2  + qL = V1
 L 2 1
L 3 2
L 4
2 12 EI L L L  2
 12 EI 9 EI   6 EI 3EI   12 EI 9 EI 
v1  3 − 3  + θ1  2 − 2  + v 2  − 3 + 3  +  qL + qL  = V1
 L L   L L   L L   8 2
3EI 3EI 3EI
v1 + 2 θ1 − 2 v 2 + 5qL = V1
L 3
L L 8
144444 42444444 3
Primera fila matriz de rigidez barra empotrada - apoyada

17
MATRIZ DE RIGIDEZ ELEMENTO DE ENREJADO A PARTIR DE MATRIZ
ELEMENTO DE VIGA BI-EMPOTRADA

RSM/2009 18
MATRIZ DE RIGIDEZ ELEMENTO DE ENREJADO A PARTIR DE MATRIZ
ELEMENTO DE VIGA BI-EMPOTRADA

RSM/2009 19
FORMULA GENERAL DE MATRIZ DE RIGIDEZ Y VECTOR DE FUERZAS
NODALES EQUIVALENTES CONDENSADOS PARA ELEMENTO DE VIGA

Kc kip c

f: grados de libertad asociados a fuerzas nodales diferente a cero

l: : grados de libertad asociados a fuerzas nodales nulas

RSM/2009 20
 DEFORMACIONES POR CORTE. ELEMENTO DE VIGA DE TIMOSHENKO.

En esta teoría se relaja la hipótesis de que las

secciones planas permanezcan perpendiculares al eje
de la viga. Esto representa una mayor aproximación a la
deformación real de la sección transversal en vigas de
gran altura o bajo módulo de corte transversal. Se
considera un giro medio de la sección de manera que
para efectos prácticos puede seguir considerándose
plana.

θ=
dv
+φ φ : Giro adicional de la sección
cuando se considera el efecto de
dx
las deformaciones de corte dv
θ≠
dx

du dθ Esto es mejor que la teoría de Euler-Bernoulli, pero aún es constante a través

εx = = −y del espesor (en teoría es cuadrática) Debe introducirse un factor de
dx dx
corrección, k.
dθ
σ x = − Ey
dx dθ
M = EI
dv du dv dx
γ xy = + = − θ = −φ
dx dy dx  dv 
V = ∫ τ xy dA = − GA − θ  = −GAγ xy
 dv   dx 
τ xy = G  − θ  A

 dx  dv V ( x)
− θ ( x) = −
dx kGA
RSM/2009 21
 MATRIZ DE RIGIDEZ ELEMENTO DE VIGA DE TIMOSHENKO.

Definiendo

Se puede comprobar que la matriz de rigidez del elemento de viga de Timoshenko está dada por:

RSM/2009 22
 MATRIZ DE RIGIDEZ ELEMENTO CON CACHOS RIGIDOS

El elemento de viga posee extremos rígidos de longitud ∆x1 y ∆x2. Al ser la curvatura nula, no hay
contribución en estos tramos a la energía potencial. Por lo tanto, es necesario definir el desplazamiento
vertical solo en el tramo deformable de longitud l =L- ∆x1 - ∆x2.
Determinemos el polinomio de Hermite N2 asociado a un giro unitario del extremo 1 (el resto de los
desplazamientos es nulo).

θ θ∆x1

∆x1 0≤x≤l
∆x2
l

 1 2∆x 
3 2 3 2 a =  2 + 3 1 θ
v(x) = ax + bx + cx + d v(x = 0) = d = θ∆x1 v(x = l) = al + bl + θ l + θ∆x = 0 l l 
v' (x) = 3ax 2 + 2bx + c v' (x = 0) = c = θ v' (x = l) = 3al2 + 2bl + θ = 0  2 3∆x 
b = − + 2 1  θ
l l 

 1 2∆x   2 3∆x 
v(x) = N 2 = ax 3 + bx 2 + cx + d =  2 + 3 1  x 3 −  + 2 1  x 2 + x + ∆x1
l l  l l 

Comprobación: haciendo ∆x1=0, l=L, se obtiene N2(x) determinado previamente

RSM/2009 23
  1 2 ∆x   2 3∆x 
N 2 ' = 3 2 + 3 1  x 2 − 2 + 2 1  x + 1
l l  l l 
 1 2 ∆x   2 3∆x 
N 2 ' ' = 6 2 + 3 1  x − 2 + 2 1 
l l  l l 

2
  1 2∆x   2 3∆x 
L
K 22 = EI ∫ 6 2 + 3 1  x − 2 + 2 1  dx = .
0 
l l  l l 
2
L
 
2 1 4∆x1 4∆x1 
2
 2 3∆x1 4∆x1 6∆x12   4 12∆x1 9∆x12 
= EI ∫ 36 x  4 + 5 + 6  − 24 x 3 + 4 + 4 + 5  + 4 2 + 3 + 4  dx
0 l l l  l l l l  l l l 
 3  1 4∆x1 4∆x12  
2 2 7∆x1 6∆x1 
2
 4 12∆x1 9∆x12 
= EI 12l  4 + 5 + 6  − 12l  3 + 4 + 5  + 4l  2 + 3 + 4 
 l l l  l l l  l l l 
4 ∆x1 ∆x1 
2
= EI  + 12 2 + 12 3 
l l l 

RSM/2009 24
 Determinemos ahora el polinomio de Hermite N1 asociado a un desplazamiento vertical unitario del
extremo 1 (el resto de los desplazamientos es nulo).

∆x1
1

0≤x≤l
∆x2
l
2
3 2
v(x) = ax + bx + cx + d v(x = 0) = d = 1
3 2
v(x = l) = al + bl + 1 = 0 a=
l3
v' (x) = 3ax 2 + 2bx + c v' (x = 0) = c = 0 v' (x = l) = 3al2 + 2bl = 0 3
b=−
l2
2 3 3 2
v(x) = N1 = ax 3 + bx 2 + cx + d = x − 2 x +1
l3 l
Ya que este N1 es similar al caso de un elemento de viga de sin extremos rígidos, debe obtenerse la
misma expresión para K11, reemplazando sólo L por l
L 2
 12 x 6  12 EI
K11 = EI ∫  3 − 2  dx = ... = 3
0
l l  l

RSM/2009 25
 MATRIZ DE RIGIDEZ ELEMENTO CON CACHOS RIGIDOS Y
DEFORMACIONES POR CORTE

r=L/GAc+L^3/(12*EI);

KL=[EA/L 0 0 -EA/L 0 0
0 1/r dx1/r+L/(2*r) 0 -1/r dx2/r+L/(2*r)
0 dx1/r+L/(2*r) (dx1^2+L*dx1)/r+EI/L+L^2/(4*r) 0 -dx1/r-L/(2*r)
dx1*dx2/r+L*(dx1+dx2)/(2*r)+L^2/(4*r)-EI/L
-EA/L 0 0 EA/L 0 0
0 -1/r -dx1/r-L/(2*r) 0 1/r -dx2/r-L/(2*r)
0 dx2/r+L/(2*r) dx1*dx2/r+L*(dx1+dx2)/(2*r)+L^2/(4*r)-EI/L 0 -dx2/r-L/(2*r) (dx2^2+L*dx2)/r+EI/L+L^2/(4*r)];
KE=AA'*KL*AA;

RSM/2009 26
 VIGA SOBRE FUNDACION ELASTICA

Consideremos una viga de largo infinito

apoyada sobre todo su largo sobre una
fundación elástica. Se aplica una carga
puntual P en el origen de los ejes x,y,z. La
carga provoca que la viga se deflecte, lo
que a su vez desplaza la fundación
elástica. Como resultado, una fuerza
distribuida se desarrolla entre la viga y la
fundación. Por lo tanto, relativa a la viga, la
rigidez de la fundación produce una fuerza
uniformemente distribuida q en la viga.
Para esta deducción, la fuerza q se toma
positiva cuando actúa hacia arriba. Cuando
la deflexión de la viga es hacia abajo
(positiva), el suelo es comprimido y empuja
hacia arriba la viga (q positivo). Cuando la
deflexión es hacia arriba (negativa), se
produce tensión en el suelo.

RSM/2009 27
Consideremos el diagrama de cuerpo libre de un elemento de viga delimitado por dos secciones
transversales distanciadas en ∆z. Para el convenio de signos indicado, y pequeños desplazamientos
se obtienen las siguientes relaciones diferenciales:

Para una fundación elástica lineal, la carga distribuida q es proporcional a la deflexión de la viga: q =k y
Donde el coeficiente de resorte k puede escribirse como: k= b k0
En que b es el ancho uniforme de la viga y k0 es el módulo de la fundación (Westergaard's modulus)
cuyas dimensiones son [F/L3]. Cuando la viga descansa sobre una capa de suelo, el módulo de la
fundación puede ser altamente variable. A menudo en diseño de fundaciones de edificios, las
estimaciones iniciales del módulo k0 se basan en descripciones cualitativas del suelo o algún
conocimiento de la resistencia a compresión no confinada qu.

Valores de módulo de suelo para arenas y arcillas

Tipo de suelo Rango de k0 [N/mm3]
Arena suelta 0.005 – 0.016
Arena media 0.010 – 0.080
Arena densa 0.063 – 0.126
Arena arcillosa (media) 0.031 – 0.080
Arena limosa (media) 0.024 – 0.048
Arcilla, qu < 0.2 N/mm2 0.012 – 0.024
Arcilla, 0.2 N/mm2<qu < 0.4 N/mm2 0.024 – 0.048
28
Arcilla, qu > 0.4 N/mm2 > 0.048
 La ecuación diferencial de la viga sobre fundación elástica está dada por: (2)

Definiendo (3)

La solución general puede expresarse como:

(4)

Esta ecuación representa la solución general para la respuesta de una viga de largo infinito sobre una
fundación elástica sometida a una carga puntual. Las constantes se determinan por las condiciones de
borde. Como las deflexión tiende a cero para valores grandes de z, las constantes C1 y C2 son cero, y
la ecuación anterior se reduce a:
(5)

Para determinar las dos constantes en el caso de una viga con carga puntual, usamos las siguientes
condiciones: a) la pendiente de la viga es cero bajo la carga debido a la simetría; b) la mitad de la carga P
debe ser soportada por la fundación elástica bajo la mitad de la longitud de la viga, especificada por
valores positivos de z. Se obtienen las relaciones:

en y (6)

La condición de pendiente cero en z=0 en ecuación (5) implica C3=C4=C, por lo que la ec. (5) queda
como:
(7)
Sustituyendo (7) en (6) se obtiene: (8)

La ecuación para la deflexión de la viga queda como:

(9)

Esta ecuación es válida para valores positivos de z. para valores negativos, por simetría, y(-z)=y(z). La
pendiente, momento y corte se obtienen reemplazando ec. (9) en (1).
 Como podemos aplicar lo que ya sabemos de métodos energéticos y elementos finitos en el análisis
de este tipo de problemas?
Para ello derivaremos, a través del principio de Energía Potencial Mínima una matriz de rigidez para
una viga apoyada sobre un medio elástico, la cual será luego implementada en una rutina de análisis
de vigas y marcos, y estudiaremos la convergencia a la solución teórica.
y, v

E, I, A
mi , θi
mj , θj x, u
Vi
Vj

k [N/mm2]
El funcional de energía potencial se puede escribir como:

1L 1 L 2 L

π = ∫ EI (v' ' ) dx + axial + k ∫ v dx − ∫ fvdx

2

20 2 0 0

Notar que aparece un término extra, el cual viene del trabajo hecho por la fundación elástica sobre la
viga deflectada. Si k es el coeficiente de rigidez del suelo o fundación, la fuerza restituyente del suelo
sobre un elemento diferencial de viga de largo dx es: F = kdx ⋅ v = kvdx
El trabajo hecho por esta fuerza es: W = 1 kvdx ⋅ v = 1 kv 2 dx
2 2
Asumiendo que EI es constante a través del largo del elemento, e ignorando los términos axiales, ya
que en un análisis lineal su contribución queda desacoplada de los términos del flexión, se tiene:
31
 L L L
δπ = EI ∫ V' ' δV' ' dx + k ∫ VδVdx − ∫ fδVdx
0 0 0

 U1 
V  ⇒ V ' ' = N' ' U
 1
θ  δV = N δU
V ( x ) = [0 N 1 N2 0 N3 N 4 ] 1 
U 2 
 V2  δV' ' = N' ' δU
 
 θ 2 

L L L
⇒ δπ = EI ∫ ( N' ' U)( N' ' δU)dx + k ∫ ( NU)( NδU)dx − ∫ fNδUdx
0 0 0
; L L L
δπ = EI ∫ δU N' ' N' ' Udx + k ∫ δU N NUdx − ∫ fNδUdx
T T T T

0 0 0

L T  L T  L

δπ = EIδU  ∫ N ' ' N ' ' dx U + kδU T

T
 ∫0 N Ndx  U − ∫0 fNδUdx
10
44244 3 14243
K debido a flexion K debido a fundacion elastica

K = K debido a flexion + K debido a fundacion elastica

RSM/2009 32
Utilizando los polinomios de interpolación de Hermite, el primer término de esta ecuación corresponde
a la matriz de rigidez regular:

 EA EA 
 0 0 − 0 0 
L L
 12 EI 6 EI 12 EI 6 EI 
 0 0 − 3 
 L3 L2 L L2 
 6 EI 4 EI 6 EI 2 EI 
0 0 − 2
 L2 L L L 
K debido a flexion = 
EA EA
− 0 0 0 0 
 L L 
 0 12 EI 6 EI 12 EI 6 EI 
− 3 − 2 0 − 2
 L L L3 L 
 6 EI 2 EI 6 EI 4 EI 
 0 0 − 2
L2 L L L 
El segundo término, luego de desarrollar las integrales, queda como:

 0 0 0 0 0 0 
 13L 11L 2
9L 13L2 
 0 35 0 −
210 70 420 
 11L2 L3 13L2 L3 
 0 0 − 
K debido a fundacion elastica = k 210 105 420 140 
0 0 0 0 0 0 
 9L 13L2 13L 11L2 
0 0 −
70 420 35 210 
 13L 2
L 3
11L 2
L3 
0 − − 0 − 
 420 140 210 105 

RSM/2009 33
 Se puede implementar esta matriz en una rutina de análisis de marcos planos. Por
ejemplo, para los siguientes valores numéricos:

P=100000; k=5000; E=29000000; I=100; L=480 in (largo total)

=0.02562288

En el centro, para x=0, la solución teórica para deflexión y momento dan:

v= 0.2562288 in
M=975690.34 lb-in

Dividiendo la viga en dos elementos se obtiene:

Momento en el centro:
M= 1378948 lb-in

El error en la estimación de la deflexión máxima es de un 70.1%, y de 29.2 % en el

cálculo del momento.???? Muy mala solucion. Es necesario discretizar más el modelo
de elementos finitos.

RSM/2009 34
35
RSM/2009 36