Introducción a los Métodos
Computacionales
TRABAJO N° 3
Integrantes:
-Mayhua Paucar, Boris Zosimo
Docente:
-Cuba Samaniego, Lisbeth Mónica.
• TRABAJO 3
• ESTATICA MANUAL DEL PROBLEMA
• DIAGRAMAS DE FUERZA CORTANTE O
MOMENTO FLECTOR
• ECUACIONES DE MOMENTOS POR TRAMO
50 x3
f x = − + 265x, (0 ≤ X ≤ 3)
9
f x = −50x 2 + 415x − 150,(3 < X ≤ 6)
f x = −185x + 1650,(6 < X ≤ 10)
f x = 100x − 1200,(10 < X ≤ 12)
Remplazar la función cuando f(x)=0 , para calcular la posición en la cual el momento es cero.
Como evaluamos cada tramo y obtuvimos la gráfica del momento flector, podemos ubicar en
que función evaluar para determinar la posición donde el momento es cero:
𝐟 𝐱 = −𝟏𝟖𝟓𝐱 + 𝟏𝟔𝟓𝟎, (𝟔 < 𝐗 ≤ 𝟏𝟎)
𝟎 = −𝟏𝟖𝟓𝐱 + 𝟏𝟔𝟓𝟎
𝐱 = 𝟖. 𝟗𝟏𝟖𝟗 𝐟𝐭
• PROGRAMACION (MATLAB)
Graficando método 1 Graficando método 2
x=linspace(0,12,100); x1=0:0.1:3;
y1=((-50*x.^3)./9+265*x).*(0<=x & x2=3:0.1:6;
x<=3); x3=6:0.1:6;
y2=(-50*x.^2+415*x-150).*(3<x & x4=10:0.1:12;
x<=6); Mx1=(-50*x1.^3)./9+265*x1;
y3=(-185*x+1650).*(6<x & x<=10); Mx2=-50*x2.^2+415*x2-150;
y4=(100*x-1200).*(10<x & x<=12); Mx3=-185*x3+1650;
f=y1+y2+y3+y4; Mx4=100*x4-1200;
plot(x,f) x=[x1 x2 x3 x4];
grid on Mx=[Mx1 Mx2 Mx3 Mx4];
plot(x,Mx)
grid on
• ANALISIS DE LA POSICION EN EL
MOMENTO CERO
function [n,m] = bisec(a,b,er2)
Método Bisección er1=100;
m0=100;
for n=0:100
m=(a+b)/2;
%Fa=-185*a+1650
Fa=viga1(a);
%Fm=-185*m+1650
Fm=viga1(m);
if Fa*Fm<0
b=m;
else
a=m;
end
er1=(abs(m-m0)/m)*100;
m0=m;
if er1<=er2
break;
end
end
Método Regula Falsi
function [n,m] = RegulaFalsi(a,b,er2)
er1=100;
m0=100;
for n=0:100
%Fa=-185*a+1650;
Fa=viga1(a);
%Fb=-185*b+1650;
Fb=viga1(b);
m=(a*Fb-b*Fa)/(Fb-Fa);
%Fm=-185*m+1650;
Fm=viga1(m);
if Fa*Fm<0
b=m;
else
a=m;
end
er1= (abs(m-m0)/m)*100;
m0=m;
if er1<=er2
break;
end
end
Método Secante
function [n ,m ] =
secante(x0,x1,er2)
er1=100;
m0=100;
for n=0:100
Fx0=viga1(x0);
Fx1=viga1(x1);
m=(x0*Fx1-x1*Fx0)/(Fx1-Fx0);
x0=x1;
x1=m;
er1=abs((m-m0)/m)*100;
m0=m;
if er1<=er2
break
end
end
end
Método Newton Rapshon
function [n,m] = newton100(x0,er1)
er2=0;
m0=0;
for n=0:100
Fx0=-185*x0+1650;
dFx0=-185;
m=x0-Fx0/dFx0;
x0=m;
er2= abs((m-m0)/m)*100;
m0=m;
if er2<=er1
break
end
end
• INTERFAZ GUIDE (PROGRAMACION
MATLAB)
IMPORTAR GRAFICA
x=linspace(2,10,100)
axes(handles.axes1);
imshow(imread('[Link]'));
GRAFICA DE MOMENTO FLECTOR
x1=0:0.1:3;
x2=3:0.1:6;
x3=6:0.1:6;
x4=10:0.1:12;
Mx1=(-50*x1.^3)./9+265*x1;
Mx2=-50*x2.^2+415*x2-150;
Mx3=-185*x3+1650;
Mx4=100*x4-1200;
x=[x1 x2 x3 x4];
f=[Mx1 Mx2 Mx3 Mx4];
axes(handles.axes2);
plot(x,f)
grid on
title('DIAGRAMA DE MOMENTO FLECTOR')
RESULTADOS DE LOS METODOS CODIFICACIÓN
function popupmenu1_Callback(hObject, eventdata, case 'REGULA FALSI'
handles)
[n,m] = RegulaFalsi(a,b,er2);
global a; datos2=[0 0; n m0];
global b; set(handles.uitable1,'data',datos2);
a=eval(get(handles.edit1,'string')); case 'SECANTE'
b=eval(get(handles.edit2,'string')); x0=a;
er2=eval(get(handles.edit3,'string')); x1=b;
[n ,m0 ] = secante(x0,x1,er1);
contenido=get(hObject,'string');
a1=get(hObject,'value'); datos3=[0 0; 0 0; n m0];
tipo=contenido(a1); set(handles.uitable1,'data',datos3);
er1=100; case 'NEWTON RAPSHON'
m0=100; x0=a;
switch cell2mat(tipo) [n,m] = newton100(x0,er1)
datos4=[0 0; 0 0; n m0];
case 'BISECCION'
[n,m0] = bisec(a,b,er2); set(handles.uitable1,'data',datos4);
datos1=[n m0]; end
set(handles.uitable1,'data',datos1);
