Script para analizar el cambio de los esfuerzos normal y cortante en un plano oblicuo bajo carga axial

plano_oblicuo.matlab

%% Datos iniciales
%
% Script para analizar el cambio de los esfuerzos normal y cortante 
% en un plano oblicuo bajo carga axial

A_0 = 1; % Area de la seccion transversal
P = 1; % Fuerza externa

% La variacion en theta desde un corte transversal (0º - 0 rad) 
% a un corte longitudinal (90º - pi/2 rad)
theta = linspace(0,pi/2,100);
% theta = 30*pi/180;

% El area de la seccion oblicua
A_theta = A_0/cos(theta);

% Fuerza normal
F = P*cos(theta);

% Fuerza cortante
V = P*sin(theta);

% Esfuerzo normal
% El operador ./ significa division elemento a elemento
% entre los vectores F y A_theta
sigma = F./A_theta;

% Esfuerzo cortante
tau = V./A_theta;

% Graficar el esfuerzo normal vs theta
% theta se pasa de radianes a grados
figure,plot(theta*180/pi,sigma)
xlabel('theta')
ylabel('Esfuerzo normal')

figure,plot(theta*180/pi,tau)
xlabel('theta')
ylabel('Esfuerzo cortante')


figure,plot(theta*180/pi,sigma)
hold on
plot(theta*180/pi,tau,'r')
xlabel('theta')
ylabel('Esfuerzo')
legend('Normal','Cortante','Location','NE')
axis tight