Programa para eje de múltiples elementos

PRÁCTICA CALIFICADA N°4

Alarcón Ortiz, Joel Efraín

20163004

DIAGRAMA DE FLUJO

[pic 1]

CÓDIGO DE MATLAB

clear;

clc;

%Introducción de datos

fprintf('ING215 Resistencia de Materiales 1\n');

fprintf('Joel Alarcón 20163004\n');

fprintf('Eje de múltiples elementos\n');

fprintf('momentos torsores, esfuerzos cortantes, ángulos de giro\n');

fprintf('Ingrese los datos\n');

%

n=input('Número de elementos=');

%Datos de las barras

for i=1:n     %para cada elemento de izquierda a derecha

fprintf('Elemento\n',i);

L(i)=input('Longitud del elemento (mm)=');

Ri(i)=input('Diámetro interno del elemento (mm)=');

Re(i)=input('Diámetro externo del elemento (mm)=');

G(i)=input('Módulo de corte del elemento (GPa)=');

T(i)=input('Momento torsor en el extremo izquierdo del elemento (N.m)=');

endfor

%Cálculos

i=1

J(i)=(Re(i)^4-Ri(i)^4)*pi/32        %Momento polar de cada elemento

K(i)=L(i)/J(i)/G(i)                 %Factor para el ángulo de giro Ang=K*T

for i=1:n-1

J(i+1)=(Re(i+1)^4-Ri(i+1)^4)*pi/32

K(i+1)=L(i+1)/J(i+1)/G(i+1)+K(i)    %Suma de todos los factores

Endfor

Kt=K(n)

i=1

J(i)=(Re(i)^4-Ri(i)^4)*pi/32

K(i)=L(i)/J(i)/G(i)

S(i)=0

for i=1:n-1

K(i+1)=L(i+1)/J(i+1)/G(i+1)

T(i+1)=T(i+1)+T(i)

S(i+1)=T(i+1)*K(i+1)+S(i)       %Suma del producto de los factores y el

                                   %total de momentos torsores aplicados

endfor

St=S(n)

Ma=St/Kt        %Simplificación de la ecuación de suma de giros = 0

Mb=T(n)-Ma      %La suma de las reacciones es igual a la suma de los torsores

for i=1:n

Ang(i)=K(i)*(Ma-T(i))*180/pi     %Se calcula el ángulo de giro de cada elemento

Esf(i)=(Ma-T(i))*Re(i)/2/J(i)    %Se calcula el esfuerzo de cada elemento

endfor

EsfM=1000*max(Esf)      %Esfuerzo cortante máximo

%Impresión de datos

fprintf('RESULTADOS DEL EJE\n');

fprintf('Momento torsor en A =%6.2f(N.m)\n', Ma);

fprintf('Momento torsor en B =%6.2f(N.m)\n', Mb);

fprintf('Esfuerzo cortante máximo =%6.2f(MPa)\n', EsfM);

fprintf('Elemento     Ángulo de giro (°)\n');

for i=1:n

fprintf('%5d          %10.4f\n', i, Ang(i));

Endfor

RESULTADOS DEL EJERCICIO 3.155

ING215 Resistencia de Materiales 1

Joel Alarcón 20163004

Eje de múltiples elementos

momentos torsores, esfuerzos cortantes, ángulos de giro

Ingrese los datos

Número de elementos=2

Elemento

Longitud del elemento (mm)=200

Diámetro interno del elemento (mm)=0

Diámetro externo del elemento (mm)=50

Módulo de corte del elemento (GPa)=77.2

Momento torsor en el extremo izquierdo del elemento (N.m)=0

Elemento

Longitud del elemento (mm)=250

Diámetro interno del elemento (mm)=0

Diámetro externo del elemento (mm)=38

Módulo de corte del elemento (GPa)=77.2

Momento torsor en el extremo izquierdo del elemento (N.m)=1400

...