7/23/2019 Exemplo de Utilizao do Matlab em Vigas e Fluidos

1/3

Aula 3

%Programa 1 - Fluidoclear allclose all

g = 9.8;deltah = 0.1;d = 0.2;deltad = 0.1;

forj = 1:1:5;h = 0.5;

fori = 1:1:5;v2 = sqrt(2*g*h);a2 = (pi*d^2)/4;q2 = v2*a2;hv(i,j) = h;v2v(i,j) = v2;dv(i,j) = d;q2v(i,j) = q2;h = h+deltah;

endd = d+deltad;end

figure(1)plot(hv,q2v);grid;xlabel('altura(m)')ylabel('vazo(m^3/s)')figure(2)plot3(hv,dv,q2v);grid;xlabel('altura(m)')ylabel('dimetro(m)')zlabel('vazo(m^3/s)')

%Programa 2 - Vigaclear allclose all

l = 1.0; %em mme = 200e9; %mdulo de elasticidademi = 0.5; %momento de inrcia

w = 300.0; %carga concentradadeltax = 0.1;deltaw = 100;

fori = 1:1:5;x = 0.0;

forj = 1:1:11;y = (((-w*x^4)/24)+((w*l^3*x)/6)-((w*l^4)/8))/(me*mi);xv(j,i) = x;yv(j,i) = y;wv(j,i) = w;x = x+deltax;

end

w = w+deltaw;end

7/23/2019 Exemplo de Utilizao do Matlab em Vigas e Fluidos

2/3

figure(1)plot(xv,yv);grid;xlabel('comprimento(m)')ylabel('deflexo(m)')figure(2)

plot(wv,yv);grid;xlabel('carga(N)')ylabel('deflexo(m)')figure(3)plot3(xv,wv,yv);grid;xlabel('comprimento(m)')ylabel('carga(N)')zlabel('deflexo(m)')

%Programa 3 - Vigaclear allclose all

l = 1.0; %em mme = 200e9; %mdulo de elasticidademi = 0.5; %momento de inrciaw = 300.0; %carga concentradadeltax = 0.1;deltami = -0.1;

fori = 1:1:3;x = 0.0;

forj = 1:1:11;y = (((-w*x^4)/24)+((w*l^3*x)/6)-((w*l^4)/8))/(me*mi);xv(j,i) = x;yv(j,i) = y;miv(j,i) = mi;x = x+deltax;

endmi = mi+deltami;end

%Programa 4 - Fluidoclear allclose all

g = 9.8;h = 0.5;

deltah = 0.1;d = 0.2;

fori = 1:1:5;v2 = sqrt(2*g*h);a2 = (pi*d^2)/4;q2 = v2*a2;hv(i) = h;v2v(i) = v2;a2v(i) = a2;q2v(i) = q2;

h = h+deltah;end

7/23/2019 Exemplo de Utilizao do Matlab em Vigas e Fluidos

3/3

figure(1)plot(hv,q2v);grid;xlabel('altura(m)')ylabel('vazo(m^3/s)')

%Programa 5 - Fluidoclear allclose all

g = 9.8;deltah = 0.1;d = 0.2;deltad = 0.1;

forj = 1:1:5;h = 0.5;

fori = 1:1:5;v2 = sqrt(2*g*h);

a2 = (pi*d^2)/4;q2 = v2*a2;hv(i,j) = h;v2v(i,j) = v2;dv(i,j) = d;q2v(i,j) = q2;h = h+deltah;

endd = d+deltad;end

figure(1)plot(hv,q2v);grid;xlabel('altura(m)')