Deflexiòn en viga simplemente apoyada

Programa para hallar la deflexión de una viga simplemente apoyada

Introduzca al matlab en un M-file los siguientes comandos:

--------------------------------------------------------------------------------

clc
clear all
%Programa para calcular la deflexión en la viga
%Datos de ingreso:
F=input('introduzca la carga que se aplicará a la viga[N]------------')
a=input('Ingrese la distancia[m] a la que se ubicará la fuerza vista desde el extremo izquierdo------------')
b=input('Ingrese la distancia[m] a la que se ubicará la fuerza vista desde el extremo derecho------------')
l=a+b;
d=input('Ingrese el diámetro[mm] de la sección------------')
E=input('Ingrese el módulo de Young[GPa]------------')
I=3.14159*(d/1000)^4/64;
E_=E*1000000000;
%Fórmulas de deflexión
x=0:0.01:a;
yAB=F.*b.*x.*(x.^2+b.^2-l.^2)/(6.*E_.*I.*l);
x1=a:0.01:l;
yBC=F.*a.*(l-x1).*(x1.^2+a.^2-2*l.*x1)/(6.*E_.*I.*l);
subplot(3,1,3),plot(x1,yBC,x,yAB)
xlabel('Tramo')
ylabel('Deflexión')
defl=F.*b.*a.*(a.^2+b.^2-l.^2)/(6.*E_.*I.*l);
disp('La deflexión en el punto donde se aplica la carga es de')
disp(defl)
%Fuerza cortante:
R1=F*b/(a+b);
R2=-F*a/(a+b);
R=R1+R2;
subplot(3,1,1),plot(x,R1,x1,R2)
xlabel('Tramo')
ylabel('Fuerza Cortante')
%Momento flector:
Ma=R1*x;
Mb=R1*x1-F*(x1-a);
subplot(3,1,2),plot(x,Ma,x1,Mb)
xlabel('Tramo')
ylabel('Momento Flector')
disp('Cortesía de Carlos León Chacón')

--------------------------------------------------------------------

Se obtiene algo como: