Ecuación Diferencial de la Viga

Diego Alejandro Varelas

Luis Miguel Hernández E.

Ecuaciones Diferenciales Parciales

Vibraciones Libres

Haremos el estudio de las vibraciones en vigas elásticas, definidas por su longitud L y su área transversal A,
y constituidas por un material de módulo de elasticidad E y densidad ρ en la teorı́a elemental. Supongamos
que:

i) La barra es homogénea.

ii) Todas las tensiones se consideran cero, excepto en la tensión normal σ, que se considera que varı́a
linealmente a lo largo de la sección transversal, con la coordenada y.

La tensión normal σ se puede escribir como:

σ = Ky, con K 6= 0. (1)

Donde K-constante y la coordenada y = 0 contiene al eje neutro a lo largo del cual la tensión normal es
nula.
La suposición de que el resto de tensiones son nulas requiere de que la resultante de todos los esfuerzos
internos a lo largo de toda la superficie transversal A de la viga sea cero, y que los momentos generados por
dichos esfuerzos en tanto al eje neutro igualen al momento de flexión interno M , esto es:
Z Z
σ dA = 0, K y dA = 0 (2)
A A

1
 Z Z
yσ dA = −M, K y 2 dA = −M (3)
A A

Como K 6= 0, la ecuación (2) implica que el eje neutro y el eje geométrico de la sección transversal coinciden.
La ecuación (3) nos sirve para definir el valor de K como:

M
K=−
Iz

Donde Iz es el momento de inercia de la sección transversal respecto al eje z.

Z
Iz = y 2 dA
A

Por lo tanto:
My My
yK = − , esto es σ = −
Iz Iz
De la ley de Hooke sabemos que:
δ F
ε= = (∗)
L AE
δ: Alargamiento.
L: Longitud original.
E: Módulo Young (Elasticidad longitudinal).
A: Sección transversal de la pieza estirada.
F : Fuerza aplicada sobre la viga.
En este caso F = −M y, A = Iz usando ley de Hooke en (∗) obtenemos

σ My
ε= =−
E Iz E

Por la ecuación de la elástica, esta una ecuación diferencial que permite encontrar la forma concreta de la
curva elástica para una viga de eje recto.

d2 v(t, x) ε M
=− = , entonces
dx2 y EIz

M = EIv 00 (4)

donde v(t, x) es el desplazamiento transversal de la viga.

A (4) se le conoce como ley de Euler-Bernoulli.
Ecuación del movimiento

Se considera un volumen infinitecimal de dx a una distancia x del final de la viga. Donde V y M son
respectivamente el esfuerzo cortante y el momento flector y F (t, x) la carga por unidad de longitud de la
viga.

Aplicando el equilibrio de momentos en el extremo izquierdo del elemento diferencial se tiene

∂M ∂V 2 d2 x
M+ dx − M − V dx − d x − F (x, t) =0
∂x ∂x 2
∂M ∂V 2 d2 x
dx − V dx − d x − F (x, t) =0
∂x ∂x 2
∂M ∂V dx
dx( −V − dx − F (x, t) ) = 0, como dx > 0
∂x ∂x 2
∂M ∂V dx
−V − dx − F (x, t) = 0.
∂x ∂x 2

Ası́ cuando dx → 0
∂M
=V (5)
∂x
La condición de equilibrio dinámico para las vibraciones transversales de la viga, se obtienen aplicando la
2da Ley de Newton

∂2v ∂V
ρAdx 2
=V −V − dx − F (t, x)dx
∂t ∂x
∂2v ∂V
ρA 2 = − − F (t, x)
∂t ∂x

sustituyendo (5) y (4) se tiene

∂2v ∂EIz v 00
 
∂
ρA 2 = − − F (t, x)
∂t ∂x ∂x
por tanto
∂2v ∂4v
ρA 2
= −EIz 4 − F (t, x) (6)
∂t ∂x
La función v(t, x) se considera lo suficientemente suave para que las derivadas apropiadas existan.
 Vibraciones Libres: Se considera en (6) F (t, x) = 0, se tiene

∂2v ∂4v
ρA 2
= −EIz 4
∂t ∂x

Luego
∂2v EIz ∂ 4 v
= − .
∂t2 ρA ∂x4
r
EIz
Sea C = , entonces
ρA
vtt + C 2 vxxxx = 0 (7)

Ejemplo 0.1. Supongamos v(t, x) = X(x)T (t) solución de (7), con X, T no triviales, demostrar que

X (4) T 00
= − 2 = λ4 = const
X c T
X(x) = A cos λx + B sen λx + C cosh λx + D senh λx

T (t) = a cos cλ2 t + b sen cλ2 t

Demostración. Resolvamos por variables separables, donde v(t, x) = X(x)T (t), entonces obtenemos lo si-
guiente. vtt (x, t) = X(x)T 00 (t) y vxxxx = X (4) (x)T (t) donde

X(x)T 00 (t) + c2 X (4) (x)T (t) = 0

X(x)T 00 (t) = −c2 X (4) (x)T (t)

T 00 (t) X (4) (x)
− = = λ4 , con λ ∈ R
c2 T (t) X(x)

Entonces

 X (4) (x) − λ4 X(x) = 0 (1)
 T 00 (t) + c2 λ4 T (t) = 0 (2)
Hallemos las soluciones de (1)
m4 − λ4 = 0 ⇔ (m2 − λ2 )(m2 + λ2 ) = 0, de donde

i) m2 − λ2 = 0 ⇔ m2 = λ2 ⇔ m = ±λ, de forma que X1 (x) = C 0 eλx + D0 e−λx

ii) m2 + λ2 = 0 ⇔ m = ±λi, entonces X2 (x) = A cos λx + B sen λx y por tanto

X(x) = A cos λx + B sen λx + C 0 eλx + D0 e−λx
C +D C −D Ceλx Deλx Ce−λx De−λx
Y haciendo C 0 = y D0 = entonces (1’) C 0 eλx = + y (2’)D0 e−λx = −
2 2 2 2 2 2
sumando (1’) y (2’)
(eλx + e−λx ) (eλx − e−λx )
C 0 eλx + D0 e−λx = C +D = C cosh λx + D senh λx
2 2
Entonces
X(x) = A cos λx + B sen λx + C cosh λx + D senh λx
Ahora hallemos las soluciones de (2)
m2 + c2 λ4 = 0 ⇔ m2 = −c2 λ4 ⇔ m = ±cλ2 i entonces
T (t) = a cos(cλ2 t) + b sen(cλ2 t)
 Ejemplo 0.2. Encontrar las soluciones vn = Xn (x)Tn (t) correspondientes a la velocidad inicial cero y que
satisfagan las condiciones de la frontera
v(t, 0) = 0, v(t, L) = 0 (extremos simplemente apoyados para todos los instantes t)
vxx (t, 0) = 0, vxx (t, L) = 0 (momento cero, por consiguiente curvatura cero, en los extremos)

Demostración. Supongamos que la solución es la forma v(t, x) = X(x)T (t), con X y T no triviales, entonces
igual
 que el ejemplo anterior, obtendremos que
 X (4) (x) − λ4 X(x) = 0 (1)
 T 00 (t) + c2 λ4 T (t) = 0 (2)
Y por las condiciones que nos dan, sabemos que
0 = vxx (t, 0) = X 00 (0)T (t) ⇒ X 00 (0) = 0
0 = vxx (t, L) = X 00 (L)T (t) ⇒ X 00 (L) = 0
Luego

 X (4) (x) − λ4 X(x) = 0
 X 00 (0) = X 00 (L) = 0 = X(0) = X(L)
Del ejemplo anterior sabemos que las solución de X(x) es

X(x) = A cos λx + B sen λx + C cosh λx + D senh λx

X 0 (x) = −Aλ sen λx + Bλ cos λx + Cλ senh λx + Dλ cosh λx

X 00 (x) = −Aλ2 cos λx − Bλ2 sen λx + Cλ2 cosh λx + Dλ2 senh λx

Utilizando las continuación condiciones de frontera obtendremos lo siguiente

0 = X 00 (0) = −λ2 A + Cλ2 ⇔ −A + C = 0 ⇔ A = C

x(0) = 0 = A + C Reemplazando lo anterior tenemos

A + A = 0 ⇔ A = 0 ⇔ C = 0, Entonces
X(x) = B sen λx + D senh λx
X 00 (x) = −Bλ2 sen λx + Dλ2 senh λx

i) X(L) = 0 = B sen λL + D senh λL

ii) X 00 (L) = 0 = −Bλ2 sen λL + Dλ2 senh λL

De (i)-(ii) tenemos
2B sen λL = 0 ⇔ sen λL = 0 ⇔ λL = nπ, n ∈ Z
nπ
λn = ,n ∈ Z
L  nπx   nπx 
Xn (x) = Bn sen + Dn senh
L L
Ahora por el ejemplo anterior tenemos
T (t) = a cos(cλ2 t) + b sen(cλ 2
t)
h nπ i2  h i 
nπ 2
Tn (t) = an cos ct + bn sen ct
L L
Por tanto
vn (t, x) = Xn (x)Tn (t)