Taller de torcion matlab

PRIMER TALLER TORCION.
clc
clear all
close all
% variables de entrada del sistema
%longitud de el 1 cilindro
a=250;
%longitud de el 2 cilindro
b=150;
%longitud de el 3 cilindro
c=200;%longitud total
long=600;
%diámetro de el cilindro 1
Diametro1=60;
%diámetro de el cilindro 2
Diametro2=30;
%diámetro de el cilindro 3
Diametro3=15;
%momento
m=65E3;
%módulo de rigidezG=125E6;
%número de elementos
n=3;
% variables de estudio
%ángulo
teta=zeros([1,n]);
%longitudes
Long=zeros([1,n]);
%diámetro
Dia=zeros([1,n]);
% K global
K_G=zeros([1,n]);
%Declarar matricesde rigidez elementales para cada elemento
K1_g_e=zeros(n+1);
K2_g_e=zeros(n+1);
K3_g_e=zeros(n+1);
%vector de longitudes
Long(1)=a;
Long(2)=b;
Long(3)=c;
%vector de diámetrosDia(1)=Diametro1;
Dia(2)=Diametro2;
Dia(3)=Diametro3;
%para la operación del momento polar de inercia
for i=1:n
J(i)=(1/4)*pi*((Dia(i)/2)^4);
end
%para calcular la rigidez de cada elemento
for i=1:nK(i)=(G)*J(i)/Long(i);
end



%matrices de rigidez
%matriz de rigidez de primer elemento
K1_g_e(1,1)=K(1);
K1_g_e(2,2)=K(1);
K1_g_e(1,2)=-K(1);
K1_g_e(2,1)=-K(1);
%matriz de rigidez desegundo elemento
K2_g_e(2,2)=K(2);
K2_g_e(3,3)=K(2);
K2_g_e(2,3)=-K(2);
K2_g_e(3,2)=-K(2);
%matriz de rigidez de tercer elemento
K3_g_e(3,3)=K(3);
K3_g_e(4,4)=K(3);
K3_g_e(3,4)=-K(3);K3_g_e(4,3)=-K(3);
%matriz global
K_G=K1_g_e+K2_g_e+K3_g_e;
%vector de cargas
P=[0;0;m;0];
%para aplicar condiciones de frontera
for i=2:3
for j=2:3
K2G(i-1,j-1)=K_G(i,j);
end
end%para reducir el vector de carga
P2=P(2:3);
teta2=K2G\P2;
teta(1)=0;
teta(4)=0;
%para hallar las torsiones
for i=1:n-1
teta(i+1)=teta2(i);
end
%para hallar reacción en 1
R1=K(1)*teta(2);%para hallar reacción en 3
R3=-K(3)*teta(3);
%para la reacción total
RT=R1-R3;
%variables y resultados a mostrar
J
K
P
K1_g_e
K2_g_e
K3_g_e
K_G
K2G
P2
teta
R1
R3
RT


J =...