Application Center - Maplesoft

App Preview:

Resolucion de la ecuación de la elástica para una viga simplemente apoyada con carga distribuida q(x), por el método aproximado de las series de Fourier

You can switch back to the summary page by clicking here.

Learn about Maple
Download Application


 

Image 

 

 

Resolucion de la ecuaci?n de la el?stica para una viga simplemente apoyada con carga distribuida q(x), por el m?todo aproximado de las series de Fourier 

 

 

> restart; 1; `*`(La, `*`(ecuacion, `*`(de, `*`(la, `*`(elastica, `*`(es)))))); 1; diff(diff(diff(diff(y(x), x), x), x), x) = `/`(`*`(q(x)), `*`(E, `*`(J))); 1; `*`(Tomamos, `*`(los, `*`(par?metros, `*`...
restart; 1; `*`(La, `*`(ecuacion, `*`(de, `*`(la, `*`(elastica, `*`(es)))))); 1; diff(diff(diff(diff(y(x), x), x), x), x) = `/`(`*`(q(x)), `*`(E, `*`(J))); 1; `*`(Tomamos, `*`(los, `*`(par?metros, `*`...
 

 

 

 

 

 

 

`*`(La, `*`(ecuacion, `*`(de, `*`(la, `*`(elastica, `*`(es))))))
diff(diff(diff(diff(y(x), x), x), x), x) = `/`(`*`(q(x)), `*`(E, `*`(J)))
`*`(Tomamos, `*`(los, `*`(par?metros, `*`(en, `*`(forma, `*`(gen?rica))))))
N
E
J
l (1.1)
 

Donde : 

 

N = N?mero de t?rminos de la serie a desarrollar 

E = M?dulo de elasticidad longitudinal 

J = Momento de inercia respecto del eje sobre el que act?a el momento flector 

l = Longitud de la viga 

 

Las condiciones de borde para el caso de la viga simplemente apoyada son: 

 

y(0) = y(l) = 0 el descenso en los apoyos es nulo 

 

M(0) = -EJ y''(0) = M(l) = -EJ y''(l) = 0 el momento en los apoyos es nulo 

 

Desarrollando en series de Fourier la funci?n y(x), la cual satisface automaticamente las condiciones de borde, tenemos: 

> assume(n, integer); 1; `:=`(y, proc (x, N) options operator, arrow; sum(`*`(w[n], `*`(sin(`/`(`*`(n, `*`(Pi, `*`(x))), `*`(l))))), n = 1 .. N) end proc); -1; y(x, N); 1
 

sum(`*`(w[n], `*`(sin(`/`(`*`(n, `*`(Pi, `*`(x))), `*`(l))))), n = 1 .. N) (1.2)
 

 

Desarrollando tambi?n  en series de Fourier la funci?n q(x), tenemos: 

> `:=`(qf, proc (x, N) options operator, arrow; sum(`*`(q[n], `*`(sin(`/`(`*`(n, `*`(Pi, `*`(x))), `*`(l))))), n = 1 .. N) end proc); -1; qf(x, N); 1
 

sum(`*`(q[n], `*`(sin(`/`(`*`(n, `*`(Pi, `*`(x))), `*`(l))))), n = 1 .. N) (1.3)
 

Donde hemos llamado qf(x,N) a la funcion q(x) expresada en series de Fourier 

 

Teniendo en cuenta la ecuacion de la el?stica nos queda: 

 

> diff(diff(diff(diff(y(x, N), x), x), x), x) = `/`(`*`(qf(x, N)), `*`(E, `*`(J))); 1
 

sum(`/`(`*`(w[n], `*`(sin(`/`(`*`(n, `*`(Pi, `*`(x))), `*`(l))), `*`(`^`(n, 4), `*`(`^`(Pi, 4))))), `*`(`^`(l, 4))), n = 1 .. N) = `/`(`*`(sum(`*`(q[n], `*`(sin(`/`(`*`(n, `*`(Pi, `*`(x))), `*`(l)))))... (1.4)
 

 

De la expesi?n anterior se observa que existe la siguiente igualdad: 

 

> `:=`(igualdad, `/`(`*`(w[n], `*`(sin(`/`(`*`(n, `*`(Pi, `*`(x))), `*`(l))), `*`(`^`(n, 4), `*`(`^`(Pi, 4))))), `*`(`^`(l, 4))) = `/`(`*`(q[n], `*`(sin(`/`(`*`(n, `*`(Pi, `*`(x))), `*`(l))))), `*`(E, `...
 

`/`(`*`(w[n], `*`(sin(`/`(`*`(n, `*`(Pi, `*`(x))), `*`(l))), `*`(`^`(n, 4), `*`(`^`(Pi, 4))))), `*`(`^`(l, 4))) = `/`(`*`(q[n], `*`(sin(`/`(`*`(n, `*`(Pi, `*`(x))), `*`(l))))), `*`(E, `*`(J))) (1.5)
 

 

De donde podemos obtener el valor de wn en funcion de qn: 

> solve({igualdad}, {w[n]}); 1
 

{w[n] = `/`(`*`(q[n], `*`(`^`(l, 4))), `*`(`^`(n, 4), `*`(`^`(Pi, 4), `*`(E, `*`(J)))))} (1.6)
 

 

La funci?n q(x) tiene la forma: 

> `:=`(q, proc (x) options operator, arrow; qo end proc); 1
 

proc (x) options operator, arrow; qo end proc (1.7)
 

 

Determinamos el valor del coeficiente qn para el desarrollo en series de Fourier de la funci?n q(x): 

> `:=`(p, q(x)); -1; `:=`(q[n], proc (x) options operator, arrow; `+`(`/`(`*`(2, `*`(int(`*`(p, `*`(sin(`/`(`*`(n, `*`(Pi, `*`(x))), `*`(l))))), x = 0 .. l))), `*`(l))) end proc); 1; `:=`(`#msub(mi(
`:=`(p, q(x)); -1; `:=`(q[n], proc (x) options operator, arrow; `+`(`/`(`*`(2, `*`(int(`*`(p, `*`(sin(`/`(`*`(n, `*`(Pi, `*`(x))), `*`(l))))), x = 0 .. l))), `*`(l))) end proc); 1; `:=`(`#msub(mi(
 

 

proc (x) options operator, arrow; `+`(`/`(`*`(2, `*`(int(`*`(p, `*`(sin(`/`(`*`(n, `*`(Pi, `*`(x))), `*`(l))))), x = 0 .. l))), `*`(l))) end proc
`+`(`-`(`/`(`*`(2, `*`(qo, `*`(`+`(`-`(1), `^`(-1, n))))), `*`(n, `*`(Pi))))) (1.8)
 

Por lo que wn valdr?: 

> `:=`(k, `#msub(mi(
 

`+`(`-`(`/`(`*`(2, `*`(qo, `*`(`+`(`-`(1), `^`(-1, n)), `*`(`^`(l, 4))))), `*`(`^`(n, 5), `*`(`^`(Pi, 5), `*`(E, `*`(J))))))) (1.9)
 

> `:=`(ps, proc (x, N) options operator, arrow; sum(`*`(k, `*`(sin(`/`(`*`(n, `*`(Pi, `*`(x))), `*`(l))))), n = 1 .. N) end proc); -1
 

Donde debimos sustituir a qf(x,N) por ps(x,N) para poder intoducir en ?l, el valor determinado del coeficiente qn al que llamamos k 

 

Entonces nuestro descenso y(x) ser?: 

> y(x, N); 1
 

sum(`+`(`-`(`/`(`*`(2, `*`(qo, `*`(`+`(`-`(1), `^`(-1, n)), `*`(`^`(l, 4), `*`(sin(`/`(`*`(n, `*`(Pi, `*`(x))), `*`(l)))))))), `*`(`^`(n, 5), `*`(`^`(Pi, 5), `*`(E, `*`(J))))))), n = 1 .. N) (1.10)
 

 

El momento flector ser?: 

> `:=`(M, proc (x, N) options operator, arrow; `+`(`-`(`*`(E, `*`(J, `*`(diff(diff(y(x, N), x), x)))))) end proc); 1; M(x, N); 1
 

 

proc (x, N) options operator, arrow; `+`(`-`(`*`(E, `*`(J, `*`(diff(diff(y(x, N), x), x)))))) end proc
`+`(`-`(`*`(E, `*`(J, `*`(sum(`+`(`/`(`*`(2, `*`(qo, `*`(`+`(`-`(1), `^`(-1, n)), `*`(`^`(l, 2), `*`(sin(`/`(`*`(n, `*`(Pi, `*`(x))), `*`(l)))))))), `*`(`^`(n, 3), `*`(`^`(Pi, 3), `*`(E, `*`(J)))))), ... (1.11)
 

La carga distrbuda, denomina ahora ps(x) ser?: 

> ps(x, N); 1
 

sum(`+`(`-`(`/`(`*`(2, `*`(qo, `*`(`+`(`-`(1), `^`(-1, n)), `*`(sin(`/`(`*`(n, `*`(Pi, `*`(x))), `*`(l))))))), `*`(n, `*`(Pi))))), n = 1 .. N) (1.12)
 

Asignando valores particulares a los par?metros tendremos: 

> `:=`(E, 1); 1; `:=`(J, 1); 1; `:=`(l, 1); 1; `:=`(qo, 1); 1
 

 

 

 

1
1
1
1 (1.13)
 

> y(x, N); 1
 

sum(`+`(`-`(`/`(`*`(2, `*`(`+`(`-`(1), `^`(-1, n)), `*`(sin(`*`(Pi, `*`(x, `*`(n))))))), `*`(`^`(n, 5), `*`(`^`(Pi, 5)))))), n = 1 .. N) (1.14)
 

> M(x, N); 1
 

`+`(`-`(sum(`+`(`/`(`*`(2, `*`(`+`(`-`(1), `^`(-1, n)), `*`(sin(`*`(Pi, `*`(x, `*`(n))))))), `*`(`^`(n, 3), `*`(`^`(Pi, 3))))), n = 1 .. N))) (1.15)
 

> ps(x, N); 1
 

sum(`+`(`-`(`/`(`*`(2, `*`(`+`(`-`(1), `^`(-1, n)), `*`(sin(`*`(Pi, `*`(x, `*`(n))))))), `*`(n, `*`(Pi))))), n = 1 .. N) (1.16)
 

> `:=`(Ty, [seq(y(x, N), N = 1 .. 7)]); -1; plot(Ty, x = 0 .. l, color = [red, green, yellow, blue, pink, magenta, black], title = `L?nea El?stica`); 1
`:=`(Ty, [seq(y(x, N), N = 1 .. 7)]); -1; plot(Ty, x = 0 .. l, color = [red, green, yellow, blue, pink, magenta, black], title = `L?nea El?stica`); 1
`:=`(Ty, [seq(y(x, N), N = 1 .. 7)]); -1; plot(Ty, x = 0 .. l, color = [red, green, yellow, blue, pink, magenta, black], title = `L?nea El?stica`); 1
 

Plot_2d
 

 

La gr?fica aproximada del momento flector ser?: 

> `:=`(TM, [seq(M(x, N), N = 1 .. 7)]); -1; plot(TM, x = 0 .. l, color = [red, green, yellow, blue, pink, magenta, black], title = `Momento Flexor`); 1
`:=`(TM, [seq(M(x, N), N = 1 .. 7)]); -1; plot(TM, x = 0 .. l, color = [red, green, yellow, blue, pink, magenta, black], title = `Momento Flexor`); 1
 

Plot_2d
 

> `:=`(Tps, [seq(ps(x, N), N = 1 .. 40)]); -1; plot(Tps, x = 0 .. l, color = [red, black], title = `Carga Distribuida`); 1
 

Plot_2d
 

Dando valores particulares a "x" y a "N" se tiene: 

 

> `:=`(x, `+`(`*`(.5, `*`(l)))); 1; `:=`(N, 101); 1
 

 

.5
101 (1.17)
 

El descenso en el punto "x" valdr?: 

> y(x, N); 1
 

0.1302083332e-1 (1.18)
 

El valor del Momento Flector ser?: 

> `:=`(M, `+`(`-`(sum(`/`(`*`(`*`(200, `+`(`-`(1), `^`(-1, n))), `*`(sin(`+`(`*`(`/`(1, 10), `*`(n, `*`(Pi, `*`(x)))))))), `*`(`^`(n, 3), `*`(`^`(Pi, 3)))), n = 1 .. N)))); 1
 

2.375037479 (1.19)
 

El valor de la carga ser?: 

> ps(x, 5); 1
 

1.103474272 (1.20)
 

>
 

>
 

 

Legal Notice: ? Maplesoft, a division of Waterloo Maple Inc. 2009. Maplesoft and Maple are trademarks of Waterloo Maple Inc. Neither Maplesoft nor the authors are responsible for any errors contained within and are not liable for any damages resulting from the use of this material.  This application is intended for non-commercial, non-profit use only. Contact the authors for permission if you wish to use this application in for-profit activities.