Scribd
Cargar
67 vistas15 páginas

Ejemplos Finitos

Este documento presenta los pasos para resolver ecuaciones diferenciales ordinarias de segundo orden mediante el método de elementos finitos. Primero, se formula la ecuación débil y se aproxima la solución mediante funciones de forma lineales por tramos. Luego, se integran las ecuaciones ponderadas y se ensamblan en un sistema matricial que permite calcular los desplazamientos nodales. Finalmente, se resuelve el sistema para obtener la solución aproximada de la ecuación diferencial original.

Cargado por

Juan Carlos Fernandez
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd
67 vistas15 páginas

Ejemplos Finitos

Este documento presenta los pasos para resolver ecuaciones diferenciales ordinarias de segundo orden mediante el método de elementos finitos. Primero, se formula la ecuación débil y se aproxima la solución mediante funciones de forma lineales por tramos. Luego, se integran las ecuaciones ponderadas y se ensamblan en un sistema matricial que permite calcular los desplazamientos nodales. Finalmente, se resuelve el sistema para obtener la solución aproximada de la ecuación diferencial original.

Cargado por

Juan Carlos Fernandez
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd

http://claudiovz.github.

io/

Ejemplo 1

Resolver
d2 u
EA +q =0
dx2
u(0) = 0
P
u0 (L) =
EA

Solución exacta

qL + P q
u(x) = x− x2
EA 2EA

Solución aproximada generalizada


La forma débil de la ecuación diferencial es
Z L L
d2 û
Z 
R(x) W (x) dx = EA 2 + q W dx = 0
0 0 dx
multiplicando
L L
d2 û
Z Z
W EA dx + W q dx = 0
0 dx2 0
Usando el teorema de Gauss o integrando por partes
  L Z L Z L
dû dW dû
W EA − EA dx + W q dx = 0
dx 0 0 dx dx 0

Reordenando
Z L Z L   L
dW dû dû
EA dx = W q dx + W EA
0 dx dx 0 dx 0
Reemplazando F = EA dû
dx
Z L Z L
dW dû L
EA dx = W q dx + (W F )|0
0 dx dx 0

1
http://claudiovz.github.io/

Ejemplo 2

Resolver
d2 u
EA +q =0
dx2
u(0) = 0
P
u0 (L) =
EA

Solución exacta

qL + P q
u(x) = x− x2
EA 2EA

Solución aproximada lineal


La forma débil de la ecuación diferencial es
Z L L
d2 û
Z 
R(x) W (x) dx = EA 2 + q W dx = 0
0 0 dx
reduciendo el grado de las derivadas
Z L Z L
dW dû dû(L) dû(0)
EA dx = W q dx + W (L) EA − W (0) EA
0 dx dx 0 dx dx
usando bases lineales en coordenadas locales
   
x x
u(x) ≈ û(x) = u1 1 − + u2
L L
ûx es
dû u2 − u1
=
dx L
las funciones ponderadas son
dû x
W1 = =1−
du1 L
dû x
W2 = =
du2 L
formando el sistema de ecuaciones
Z L Z L
dW1 dû dû(L) dû(0)
EA dx = W1 q dx + W1 (L) EA − W1 (0) EA
0 dx dx 0 dx dx
Z L Z L
dW2 dû dû(L) dû(0)
EA dx = W2 q dx + W2 (L) EA − W2 (0) EA
0 dx dx 0 dx dx

1
http://claudiovz.github.io/

funciones ponderadas y sus derivadas


x dW1 1 x dW2 1
W1 = 1 − =− W2 = =
L dx L L dx L
valores de las funciones ponderadas en los nodos

W1 (L) = 0 W1 (0)= 1 W2 (L) = 1 W2 (0)= 0

fuerzas en los nodos


dû(L) dû(0)
EA = F2 EA = F1
dx dx
reemplazando
L    Z L 
u2 − u1
Z
1 x
− EA dx = 1− q dx + 0(F2 ) − 1(F1 )
0 L L 0 L
Z L    Z L 
1 u2 − u1 x
EA dx = q dx + 1(F2 ) − 0(F1 )
0 L L 0 L

reordenando
 Z L Z L
u1 − u2 q
EA dx = q − x dx − F1
L2 0 0 L
 Z L Z L
−u1 + u2 q
EA dx = x dx + F2
L2 0 0 L

integrando
 
u1 − u2 qL
EA = − F1
L 2
 
−u1 + u2 qL
EA = + F2
L 2

en forma matricial
      
qL
EA  1 −1 u1 2 −F1
  =   +  
L −1 1 u2 qL
F2
2

reemplazando fuerzas y desplazamientos


      
qL
EA  1 −1 0 −F 1
  =  2  +  
L −1 1 u2 qL
P
2

resolviendo
qL2 + 2P L
u2 =
2EA
F1 = qL + P

reemplazando en la solución aproximada


qL2 + 2P L x
 
qL + 2P
û(x) = = x
2EA L 2EA

2
http://claudiovz.github.io/

Ejemplo 3

Resolver
d2 u
EA +q =0
dx2
u(0) = 0
P
u0 (L) =
EA

Solución exacta

qL + P q
u(x) = x− x2
EA 2EA

Solución aproximada cuadrática


La forma débil de la ecuación diferencial es
Z L L
d2 û
Z 
R(x) W (x) dx = EA 2 + q W dx = 0
0 0 dx
reduciendo el grado de las derivadas
Z L Z L
dW dû dû(L) dû(0)
EA dx = W q dx + W (L) EA − W (0) EA
0 dx dx 0 dx dx
usando bases cuadráticas en coordenadas locales
     
3 2 4 4 1 2
u(x) ≈ û(x) = u1 1 − x + 2 x2 + u2 x − 2 x2 + u3 − x + 2 x2
L L L L L L
ûx es
dû 3u1 − 4u2 + u3 4u1 − 8u2 + 4u3
=− + x
dx L L2
las funciones ponderadas son
dû 3 2
W1 = = 1 − x + 2 x2
du1 L L
dû 4 4 2
W2 = = x − 2x
du2 L L
dû 1 2
W3 = = − x + 2 x2
du3 L L
formando el sistema de ecuaciones

1
http://claudiovz.github.io/

Z L Z L
dW1 dû dû(L) dû(0)
EA dx = W1 q dx + W1 (L) EA − W1 (0) EA
0 dx dx 0 dx dx
Z L Z L
dW2 dû dû(L) dû(0)
EA dx = W2 q dx + W2 (L) EA − W2 (0) EA
0 dx dx 0 dx dx
Z L Z L
dW3 dû dû(L) dû(0)
EA dx = W3 q dx + W3 (L) EA − W3 (0) EA
0 dx dx 0 dx dx

funciones ponderadas y sus derivadas


3 2 dW1 3 4 4 4 dW2 4 8
W1 = 1 − x + 2 x2 = − + 2x W2 = x − 2 x2 = − 2x
L L dx L L L L dx L L
1 2 dW3 1 4
W3 = − x + 2 x2 = − + 2x
L L dx L L
valores de las funciones ponderadas en los nodos

W1 (L) = 0 W1 (0)= 1 W2 (L)= 0 W2 (0) = 0


W3 (L) = 1 W3 (0)= 0

fuerzas en los nodos


dû(L) dû(0)
EA = F3 EA = F1
dx dx
reemplazando
Z L    Z L 
3 4 3u1 − 4u2 + u3 4u1 − 8u2 + 4u3 3 2 2
+ x EA − + x dx = 1 − x + 2 x q dx + 0(F3 ) − 1(F1 )
0 L L2 L L2 0 L L
Z L    Z L 
4 8 3u1 − 4u2 + u3 4u1 − 8u2 + 4u3 4 4 2
− 2 x EA − + x dx = x − x q dx + 0(F3 ) − 0(F1 )
0 L L L L2 0 L L2
Z L    Z L 
1 4 3u1 − 4u2 + u3 4u1 − 8u2 + 4u3 1 2 2
− + 2 x EA − + x dx = − x + x q dx + 1(F3 ) − 0(F1 )
0 L L L L2 0 L L2

reordenando
L   Z L 
3u1 − 4u2 + u3 4u1 − 8u2 + 4u3
Z
EA 4 3 2 2
3+ x −+ x dx = 1 − x + x q dx − F1
L 0 L L L2 0 L L2
EA L
Z    Z L 
8 3u1 − 4u2 + u3 4u1 − 8u2 + 4u3 4 4 2
4− x − + x dx = x − x q dx
L 0 L L L2 0 L L2
EA L
Z    Z L 
4 3u1 − 4u2 + u3 4u1 − 8u2 + 4u3 1 2 2
−1 + x − + x dx = − x + 2 x q dx + F3
L 0 L L L2 0 L L

integrando
 
EA 7 8 1 qL
u1 − u2 + u3 = − F1
L 3 3 3 6
 
EA 8 16 8 2qL
− u1 + u2 − u3 =
L 3 3 3 3
 
EA 1 8 7 qL
u1 − u2 + u3 = + F3
L 3 3 3 6

en forma matricial

2
http://claudiovz.github.io/

      
7 qL
3 − 83 1
3 u1 6 −F1
EA       
− 83  u2  =  2qL
 8 16
− +
     
0
L  3 3    3  
  

1 qL
3 − 83 7
3 u3 6 F3
reemplazando fuerzas y desplazamientos
      
7 qL
3 − 83 1
3 0 −F1
EA     6   
 8 16  2qL 
− − 83  u2  = + 0 
   
L  3 3
 3 
     
1 qL
3 − 83 7
3 u3 6 P
resolviendo
3qL2 + 4P L
u2 =
8EA
qL2 + 2P L
u3 =
2EA
F1 = qL + P

reemplazando en la solución aproximada


3qL2 + 4P L
    2  
4 4 2 qL + 2P L 1 2 2 qL + P q
û(x) = x − 2x + − x + 2x = x− x2
8EA L L 2EA L L EA 2EA

3
http://claudiovz.github.io/

Ejemplo 4

Resolver
d4 v
EI +q =0
dx4
v(0) = 0 EIv 00 (L) = 0
v 0 (0) = 0 EIv 000 (L)= P

Solución exacta

qL2 + 2P L 2 qL + P 3 q
v(x) = − x + x − x4
4EI 6EI 24EI

Solución aproximada generalizada


La forma débil de la ecuación diferencial es
Z L L
d4 v̂
Z 
R(x) W (x) dx = EI + q W dx = 0
0 0 dx4
multiplicando
L Z L Z L
d4 v̂ d4 v̂
Z 
EI + q W dx = W EI dx + W q dx = 0
0 dx4 0 dx4 0

usando el teorema de Gauss o integrando por partes


 L Z L 2 Z L
d3 v̂ d2 v̂ d2 v̂

dW d W
W EI 3 − EI 2 + EI dx + W q dx = 0
dx dx dx 0 0 dx2 dx2 0

reordenando
L L  L
d2 W d2 v̂ d3 v̂ d2 v̂
Z Z 
dW
EI dx = − W q dx − W EI 3 − EI 2
0 dx2 dx2 0 dx dx dx 0
cortante y momento
d3 v̂ d2 v̂
V = EI M = EI
dx3 dx2
reemplazando
L L  L
d2 W d2 v̂
Z Z 
dW
EI dx = − W q dx − WV − M
0 dx2 dx2 0 dx 0

1
http://claudiovz.github.io/

Ejemplo 5

Resolver
d4 v
EI +q =0
dx4
v(0) = 0 EIv 00 (L) = 0
v 0 (0) = 0 EIv 000 (L)= P

Solución exacta

qL2 + 2P L 2 qL + P 3 q
v(x) = − x + x − x4
4EI 6EI 24EI

Solución aproximada cúbica


La forma débil de la ecuación diferencial es
Z L L
d4 v̂
Z 
R(x) W (x) dx = EI + q W dx = 0
0 0 dx4
reduciendo el grado de las derivadas
L L
d2 W d2 v̂ d3 v̂(L) d3 v̂(0) dW (L) d2 v̂(L) dW (0) d2 v̂(0)
Z Z
EI dx = − W q dx − W (L) EI + W (0) EI + EI − EI
0 dx2 dx2 0 dx3 dx3 dx dx2 dx dx2
usando bases cúbicas en coordenadas locales
       
3 2 2 1 3 2 2 3 1 2 1 3
v(x) ≈ v̂(x) = v1 1 − 2 x2 + 3 x3 + θ1 x − x2 + 2 x3 + v2 x − x + θ 2 − x + x
L L L L L2 L3 L L2
v̂xx es
d2 v̂
       
6 12 4 6 6 12 2 6
= − 2 + 3 x v1 + − + x θ1 + − 3 x v2 + − + 2 x θ 2
dx2 L L L L2 L2 L L L
las funciones ponderadas son
dv̂ 3 2 2
W1 = =1− 2
x + 3 x3
dv1 L L
dv̂ 2 2 1 3
W2 = = x − x + 2x
dθ1 L L
dv̂ 3 2 2 3
W3 = = 2x − 3x
dv2 L L
dv̂ 1 2 1
W4 = = − x + 2 x3
dθ2 L L

1
http://claudiovz.github.io/

formando el sistema de ecuaciones


Z L 2 Z L
d W1 d2 v̂ d3 v̂(L) d3 v̂(0) dW1 (L) d2 v̂(L) dW1 (0) d2 v̂(0)
2
EI 2
dx = − W1 q dx − W1 (L) EI + W1 (0) EI + EI − EI
0 dx dx 0 dx3 dx 3 dx dx 2 dx dx2
Z L 2 L
d2 v̂ d3 v̂(L) d3 v̂(0) dW2 (L) d2 v̂(L) dW2 (0) d2 v̂(0)
Z
d W2
EI dx = − W 2 q dx − W 2 (L) EI + W2 (0) EI + EI − EI
0 dx2 dx2 0 dx3 dx3 dx dx2 dx dx2
Z L 2 2 Z L 3
d W3 d v̂ d v̂(L) d3 v̂(0) dW3 (L) d2 v̂(L) dW3 (0) d2 v̂(0)
2
EI 2 dx = − W3 q dx − W3 (L) EI + W3 (0) EI + EI − EI
0 dx dx 0 dx3 dx3 dx dx2 dx dx2
Z L 2 L
d2 v̂ d3 v̂(L) d3 v̂(0) dW4 (L) d2 v̂(L) dW4 (0) d2 v̂(0)
Z
d W4
2
EI 2 dx = − W4 q dx − W4 (L) EI + W4 (0) EI + EI − EI
0 dx dx 0 dx3 dx 3 dx dx 2 dx dx2
funciones ponderadas y sus derivadas
3 2 2 d 2 W1 6 12 2 2 1 d2 W2 4 6
W1 = 1 − x + 3 x3 = − 2 + 3x W2 = x − x + 2 x3 = − + 2x
L2 L dx2 L L L L dx2 L L
3 2 d 2 W3 6 12 1 1 d2 W4 2 6
W3 = 2 x2 − 3 x3 = 2 − 3x W 4 = − x2 + 2 x3 = − + 2x
L L dx2 L L L L dx2 L L
valores de las funciones ponderadas en los nodos
dW1 (L) dW1 (0)
W1 (L) = 0 W1 (0)= 1 =0 =0
dx dx
dW2 (L) dW2 (0)
W2 (L) = 0 W2 (0)= 0 =0 =1
dx dx
dW3 (L) dW3 (0)
W3 (L) = 1 W3 (0)= 0 =0 =0
dx dx
dW4 (L) dW4 (0)
W4 (L) = 0 W4 (0)= 0 =1 =0
dx dx
cortante y momento en los nodos
d3 v̂(L) d3 v̂(0)
EI = V2 EI = V1
dx3 dx3
d2 v̂(L) d2 v̂(0)
EI = M2 EI = M1
dx2 dx2
reemplazando
Z L          
6 12 6 12 4 6 6 12 2 6
− 2 + 3 EI − 2 + 3 x v1 + − + 2 x θ1 + − x v 2 + − + x θ2 dx
0 L L L L L L L2 L3 L L2
Z L 
3 2
=− 1 − 2 x2 + 3 x3 q dx − 0(V2 ) + 1(V1 ) + 0(M2 ) − 0(M1 )
0 L L
Z L          
4 6 6 12 4 6 6 12 2 6
− + 2 x EI − 2 + 3 x v1 + − + 2 x θ1 + − x v 2 + − + x θ2 dx
0 L L L L L L L2 L3 L L2
Z L 
2 1
=− x − x2 + 2 x3 q dx − 0(V2 ) + 0(V1 ) + 0(M2 ) − 1(M1 )
0 L L
Z L          
6 12 6 12 4 6 6 12 2 6
− x EI − + x v 1 + − + x θ 1 + − x v 2 + − + x θ2 dx
0 L2 L3 L2 L3 L L2 L2 L3 L L2
Z L 
3 2 2 3
=− x − x q dx − 1(V2 ) + 0(V1 ) + 0(M2 ) − 0(M1 )
0 L2 L3
Z L          
2 6 6 12 4 6 6 12 2 6
− + 2 x EI − 2 + 3 x v1 + − + 2 x θ1 + − 3 x v2 + − + 2 x θ2 dx
0 L L L L L L L2 L L L
Z L 
1 1
=− − x2 + 2 x3 q dx − 0(V2 ) + 0(V1 ) + 1(M2 ) − 0(M1 )
0 L L

2
http://claudiovz.github.io/

integrando
EI qL
(12v1 + 6Lθ1 − 12v2 + 6Lθ2 ) = − + V1
L3 2
EI 2 2 qL2
(6Lv 1 + 4L θ 1 − 6Lv2 + 2L θ2 ) = − − M1
L3 12
EI qL
3
(−12v1 − 6Lθ1 + 12v2 − 6Lθ2 ) = − − V2
L 2
EI qL2
3
(6Lv1 + 2L2 θ1 − 6Lv2 + 4L2 θ2 ) = + M2
L 12
en forma matricial
      
12 6L −12 v16L − qL V
    2   1 
2
2
−6L 2L  θ1  − qL
2  −M1 
     
EI  6L 4L
    =  12  +  
L3 
−12 −6L 12 −6L v2   − qL
    
−V 2

    22  
 
 
qL
6L 2L2 −6L 4L2 θ2 12 M 2

reemplazando fuerzas y desplazamientos


      
12 6L −12 6L 0 − qL V 1
    2   
2
6L 4L2 −6L 2L2   0  − qL −M
     
EI     =  12  + 
  1

L3 
−12 −6L 12 −6L v2   − qL
    
−P

    22  
  
 
2 2 qL
6L 2L −6L 4L θ2 12 0
resolviendo
3qL4 + 8P L3
v2 = −
24EI
qL3 + 3P L2
θ2 = −
6EI
V1 = qL + P
qL2 + 2P L
M1 = −
2
reemplazando en la solución aproximada

3qL4 + 8P L3 qL3 + 3P L2 5qL2 + 12P L 2 qL + 2P 3


     
3 2 2 3 1 2 1 3
v̂(x) = − x − x + − − x + x = − x + x
24EI L2 L3 6EI L L2 24EI 12EI

3
http://claudiovz.github.io/

Ejemplo 6

Resolver
d4 v
EI +q =0
dx4
v(0) = 0 EIv 00 (L) = 0
0
v (0) = 0 EIv 000 (L)= P

Solución exacta

qL2 + 2P L 2 qL + P 3 q
v(x) = − x + x − x4
4EI 6EI 24EI

Solución aproximada de quinto orden


La forma débil de la ecuación diferencial es
Z L L
d4 v̂
Z 
R(x) W (x) dx = EI + q W dx = 0
0 0 dx4
reduciendo el grado de las derivadas
L L
d2 W d2 v̂ d3 v̂(L) d3 v̂(0) dW (L) d2 v̂(L) dW (0) d2 v̂(0)
Z Z
2
EI 2 dx = − W q dx − W (L) EI 3
+ W (0) EI 3
+ EI 2
− EI
0 dx dx 0 dx dx dx dx dx dx2
usando bases de quinto orden en coordenadas locales
   
23 66 68 24 6 13 12 4
v(x) ≈ v̂(x) = v1 1 − 2 x2 + 3 x3 − 4 x4 + 5 x5 + θ1 x − x2 + 2 x3 − 3 x4 + 4 x5
L L L L L L L L
   
16 2 32 3 16 4 8 32 40 16
+ v2 x − 3 x + 4 x + θ 2 − x2 + 2 x3 − 3 x4 + 4 x5
L2 L L L L L L
   
7 2 34 3 52 4 24 5 1 2 5 3 8 4 4 5
+ v3 x − x + x − x + θ 3 − x + x − x + x
L2 L3 L4 L5 L L2 L3 L4
v̂xx es
d2 v̂
   
46 396 816 480 12 78 144 80
= − 2 + 3 x − 4 x2 + 5 x3 v1 + − + 2 x − 3 x2 + 4 x3 θ1
dx2 L L L L L L L L
   
32 192 192 16 192 480 320
+ − 3 x + 4 x2 v 2 + − + 2 x − 3 x2 + 4 x3 θ 2
L2 L L L L L L
   
14 204 624 2 480 3 2 30 96 2 80 3
+ − x + x − x v 3 + − + x − x + x θ3
L2 L3 L4 L5 L L2 L3 L4

1
http://claudiovz.github.io/

las funciones ponderadas son


dv̂ 23 2 66 3 68 4 24 5
W1 = =1− x + 3x − 4x + 5x
dv1 L2 L L L
dv̂ 6 2 13 3 12 4 4 5
W2 = = x − x + 2x − 3x + 4x
dθ1 L L L L
dv̂ 16 2 32 3 16 4
W3 = = 2x − 3x + 4x
dv2 L L L
dv̂ 8 2 32 3 40 4 16 5
W4 = = − x + 2x − 3x + 4x
dθ2 L L L L
dv̂ 7 2 34 3 52 4 24 5
W5 = = 2x − 3x + 4x − 5x
dv3 L L L L
dv̂ 1 2 5 3 8 4 4
W6 = = − x + 2 x − 3 x + 4 x5
dθ3 L L L L
formando el sistema de ecuaciones
Z L 2 Z L
d W1 d2 v̂ d3 v̂(L) d3 v̂(0) dW1 (L) d2 v̂(L) dW1 (0) d2 v̂(0)
EI dx = − W 1 q dx − W 1 (L) EI + W1 (0) EI + EI − EI
0 dx2 dx2 0 dx3 dx3 dx dx2 dx dx2
Z L 2 L
d2 v̂ d3 v̂(L) d3 v̂(0) dW2 (L) d2 v̂(L) dW2 (0) d2 v̂(0)
Z
d W2
2
EI 2 dx = − W2 q dx − W2 (L) EI + W2 (0) EI + EI − EI
0 dx dx 0 dx3 dx 3 dx dx 2 dx dx2
Z L 2 L
d2 v̂ d3 v̂(L) d3 v̂(0) dW3 (L) d2 v̂(L) dW3 (0) d2 v̂(0)
Z
d W3
2
EI 2 dx = − W3 q dx − W3 (L) EI + W3 (0) EI + EI − EI
0 dx dx 0 dx3 dx 3 dx dx 2 dx dx2
Z L 2 L
d2 v̂ d3 v̂(L) d3 v̂(0) dW4 (L) d2 v̂(L) dW4 (0) d2 v̂(0)
Z
d W4
EI dx = − W 4 q dx − W 4 (L) EI + W4 (0) EI + EI − EI
0 dx2 dx2 0 dx3 dx3 dx dx2 dx dx2
Z L 2 2 Z L 3
d W5 d v̂ d v̂(L) d3 v̂(0) dW5 (L) d2 v̂(L) dW5 (0) d2 v̂(0)
2
EI 2 dx = − W5 q dx − W5 (L) EI + W5 (0) EI + EI − EI
0 dx dx 0 dx3 dx3 dx dx2 dx dx2
Z L 2 L
d2 v̂ d3 v̂(L) d3 v̂(0) dW6 (L) d2 v̂(L) dW6 (0) d2 v̂(0)
Z
d W6
2
EI 2 dx = − W6 q dx − W6 (L) EI + W6 (0) EI + EI − EI
0 dx dx 0 dx3 dx 3 dx dx 2 dx dx2
funciones ponderadas y sus derivadas
23 2 66 3 68 4 24 5 d 2 W1 46 396 816 480
W1 = 1 − x + 3x − 4x + 5x = − 2 + 3 x − 4 x2 + 5 x3
L2 L L L dx2 L L L L
6 2 13 3 12 4 4 5 d2 W2 12 78 144 2 80 3
W2 = x − x + 2x − 3x + 4x = − + 2x − 3 x + 4x
L L L L dx2 L L L L
16 2 32 3 16 4 d2 W3 32 192 192 2
W3 = 2x − 3x + 4x = 2 − 3 x+ 4 x
L L L dx2 L L L
2
8 2 32 3 40 4 16 5 d W4 16 192 480 320
W4 = − x + 2x − 3x + 4x 2
= − + 2 x − 3 x2 + 4 x3
L L L L dx L L L L
7 2 34 3 52 4 24 5 d2 W5 14 204 624 2 480 3
W5 = 2x − 3x + 4x − 5x = 2 − 3 x+ 4 x − 5 x
L L L L dx2 L L L L
1 2 5 3 8 4 4 d2 W6 2 30 96 2 80 3
W6 = − x + 2 x − 3 x + 4 x5 = − + 2x − 3x + 4x
L L L L dx2 L L L L
valores de las funciones ponderadas en los nodos
dW1 (L) dW1 (0)
W1 (L) = 0 W1 (0)= 1 =0 =0
dx dx
dW2 (L) dW2 (0)
W2 (L) = 0 W2 (0)= 0 =0 =1
dx dx
dW3 (L) dW3 (0)
W3 (L) = 0 W3 (0)= 0 =0 =0
dx dx

2
http://claudiovz.github.io/

dW4 (L) dW4 (0)


W4 (L) = 0 W4 (0)= 0 =0 =0
dx dx
dW5 (L) dW5 (0)
W5 (L) = 1 W5 (0)= 0 =0 =0
dx dx
dW6 (L) dW6 (0)
W6 (L) = 0 W6 (0)= 0 =1 =0
dx dx
cortante y momento en los nodos
d3 v̂(L) d3 v̂(0)
EI = V3 EI = V1
dx3 dx3
d2 v̂(L) d2 v̂(0)
EI = M3 EI = M1
dx2 dx2
reemplazando
Z L   
46 396 816 2 480 3 46 396 816 2 480 3
− 2 + 3 x − 4 x + 5 x EI − 2 + 3 x − 4 x + 5 x v1
0 L L L L L L L L
     
12 78 144 2 80 3 32 192 192 2 16 192 480 2 320 3
+ − + 2 x − 3 x + 4 x θ1 + − 3 x + 4 x v2 + − + 2 x − 3 x + 4 x θ 2
L L L L L2 L L L L L L
    
14 204 624 480 2 30 96 80
+ − 3 x + 4 x2 − 5 x3 v3 + − + 2 x − 3 x2 + 4 x3 θ3 dx =
L2 L L L L L L L
Z L 
23 66 68 24
− 1 − 2 x2 + 3 x3 − 4 x4 + 5 x5 q dx − 0(V3 ) + 1(V1 ) + 0(M3 ) − 0(M1 )
0 L L L L
Z L   
12 78 144 80 46 396 816 480
− + 2 x − 3 x2 + 4 x3 EI − 2 + 3 x − 4 x2 + 5 x3 v1
0 L L L L L L L L
     
12 78 144 2 80 3 32 192 192 2 16 192 480 2 320 3
+ − + 2 x − 3 x + 4 x θ1 + − x + x v 2 + − + x − x + x θ2
L L L L L2 L3 L4 L L2 L3 L4
    
14 204 624 480 2 30 96 80
+ 2
− 3 x + 4 x2 − 5 x3 v3 + − + 2 x − 3 x2 + 4 x3 θ3 dx =
L L L L L L L L
Z L 
6 2 13 3 12 4 4 5
− x − x + 2 x − 3 x + 4 x q dx − 0(V3 ) + 0(V1 ) + 0(M3 ) − 1(M1 )
0 L L L L
Z L     
32 192 192 2 46 396 816 2 480 3 12 78 144 2 80 3
− x + x EI − + x − x + x v 1 + − + x − x + x θ1
0 L2 L3 L4 L2 L3 L4 L5 L L2 L3 L4
     
32 192 192 16 192 480 320 14 204 624 480
+ − 3 x + 4 x2 v2 + − + 2 x − 3 x2 + 4 x3 θ2 + − 3 x + 4 x 2 − 5 x 3 v3
L2 L L L L L L L2 L L L
   Z L 
2 30 96 80 16 2 32 3 16 4
+ − + 2 x − 3 x2 + 4 x3 θ3 dx = − x − 3 x + 4 x q dx − 0(V3 ) + 0(V1 ) + 0(M3 ) − 0(M1 )
L L L L 0 L2 L L
Z L      
16 192 480 2 320 3 46 396 816 2 480 3 12 78 144 2 80 3
− + 2 x − 3 x + 4 x EI − 2 + 3 x − 4 x + 5 x v1 + − + 2 x − 3 x + 4 x θ1
0 L L L L L L L L L L L L
     
32 192 192 2 16 192 480 2 320 3 14 204 624 2 480 3
+ − 3 x + 4 x v2 + − + 2 x − 3 x + 4 x θ2 + − 3 x + 4 x − 5 x v3
L2 L L L L L L L2 L L L
   Z L 
2 30 96 80 8 32 40 16
+ − + 2 x− 3 x2 + 4 x3 θ3 dx = − − x2 + 2 x3 − 3 x4 + 4 x5 q dx−0(V3 )+0(V1 )+0(M3 )−0(M1 )
L L L L 0 L L L L
Z L     
14 204 624 2 480 3 46 396 816 2 480 3 12 78 144 2 80 3
− 3 x + 4 x − 5 x EI − 2 + 3 x − 4 x + 5 x v1 + − + 2 x − 3 x + 4 x θ1
0 L2 L L L L L L L L L L L
     
32 192 192 2 16 192 480 2 320 3 14 204 624 2 480 3
+ − 3 x + 4 x v2 + − + 2 x − 3 x + 4 x θ2 + − 3 x + 4 x − 5 x v3
L2 L L L L L L L2 L L L
   Z L 
2 30 96 80 7 2 34 3 52 4 24 5
+ − + 2 x− 3 x2 + 4 x3 θ3 dx = − 2
x − 3 x + 4 x − 5 x q dx−1(V3 )+0(V1 )+0(M3 )−0(M1 )
L L L L 0 L L L L

3
http://claudiovz.github.io/

Z L     
2 30 96 2 80 3 46 396 816 2 480 3 12 78 144 2 80 3
− + 2 x − 3 x + 4 x EI − 2 + 3 x − 4 x + 5 x v1 + − + 2 x − 3 x + 4 x θ1
0 L L L L L L L L L L L L
     
32 192 192 16 192 480 320 14 204 624 480
+ − 3 x + 4 x2 v2 + − + 2 x − 3 x2 + 4 x3 θ2 + − 3 x + 4 x 2 − 5 x 3 v3
L2 L L L L L L L2 L L L
   Z L 
2 30 96 80 1 5 8 4
+ − + 2 x− 3 x2 + 4 x3 θ3 dx = − − x2 + 2 x3 − 3 x4 + 4 x5 q dx−0(V3 )+0(V1 )+1(M3 )−0(M1 )
L L L L 0 L L L L
integrando
 
EI 5092 1138L 512 384L 1508 242L 7qL
3
v1 + θ1 − v2 + θ2 − v3 + θ3 = − + V1
L 35 35 5 7 35 35 30
332L2 64L2 38L2 qL2
 
EI 1138L 128L 242L
v 1 + θ 1 − v 2 + θ 2 − v 3 + θ 3 = − − M1
L3 35 35 5 7 35 35 60
 
EI 512 128L 1024 512 128L 8qL
3
− v1 − θ1 + v2 − v3 + θ3 = −
L 5 5 5 5 5 15
64L2 256L2 64L2
 
EI 384L 384L
v1 + θ1 + θ2 − v3 + θ3 = 0
L3 7 7 7 7 7
 
EI 1508 242L 512 384L 5092 1138L 7qL
3
− v1 − θ1 − v2 − θ2 + v3 − θ3 = − − V3
L 35 35 5 7 35 35 30
38L2 64L2 332L2 qL2
 
EI 242L 128L 1138L
v 1 + θ 1 + v 2 + θ 2 − v 3 + θ 3 = + M3
L3 35 35 5 7 35 35 60
en forma matricial
      
5092
35
1138L
35 − 512
5
384L
7 − 1508
35
242L
35 v1 − 7qL V1
     30   
 1138L 332L2 64L2 38L2     qL2  
− 128L − 242L − −M

θ
  1   60   1
 35 35 5 7 35 35
     
 
 512    8qL  
− − 128L 1024
− 512 128L 
−

EI  0 v 2 0
 5 5 5 5 5    =  15  + 
      

3
L  384L 64L2 256L2 64L2 
− 384L
      
0 θ 0 0
 7 7 7 7 7   2 
      
  
 1508 7qL
− 35 − 242L − 512 − 384L 5092
− 1138L − −V
      
35 5 7 35 v
35    3 3
 30
      
   
242L 38L2 128L 64L2 332L2 qL2
35 35 5 7 − 1138L
35 35 θ 3 60 M 3

reemplazando fuerzas y desplazamientos


      
5092 1138L
− 512 384L
− 1508 242L
0 − 7qL V
 35 35 5 7 35 35     30   1 
2 2
 1138L 332L 128L 64L2 38L2   0  − qL
− 5 − 242L  −M1 
    
 35 35 7 35 35
     60   
 512 128L 1024    8qL  
− − 5 − 512 128L 
−

EI  0 v 0
 5 5 5 5   2  =  15  + 
      

L3  384L
 7 64L 2
0 256L2
− 384L 64L 2  
θ 2
 
0
 
0

7 7 7 7
      
      
 1508 242L    7qL  
− 35 − 35 − 512 − 384L 5092 1138L 
− 35  v3  − 30   −P 

 5 7 35    2   
242L 38L2 128L 64L2 332L2 qL
35 35 5 7 − 1138L
35 35 θ 3 60 0
resolviendo
17qL4 + 40P L3
v2 = −
384EI
7qL + 18P L2
3
θ2 = −
48EI
3qL4 + 8P L3
v3 = −
24EI
qL3 + 3P L2
θ3 = −
6EI

4
http://claudiovz.github.io/

V1 = qL + P
qL2 + 2P L
M1 = −
2
reemplazando en la solución aproximada
17qL4 + 40P L3 7qL3 + 18P L2
     
16 2 32 3 16 4 8 2 32 3 40 4 16 5
v̂(x) = − x − x + x + − − x + x − x + x
384EI L2 L3 L4 48EI L L2 L3 L4
3qL4 + 8P L3 qL3 + 3P L2
     
7 2 34 3 52 4 24 5 1 2 5 3 8 4 4 5
+ − x − x + x − x + − − x + x − x + x
24EI L2 L3 L4 L5 6EI L L2 L3 L4
qL2 + 2P L 2 qL + P 3 q
=− x + x − x4
4EI 6EI 24EI

También podría gustarte