Scribd
Upload
90 views2 pages

1d Isotropic Linear Hardening

This document contains MATLAB code that is simulating plastic deformation behavior using the elastic-plastic constitutive model. The code defines material properties like Young's modulus and yield strength, then increments strain in small steps. It calculates the stress, plastic strain, and yield surface location at each step based on whether the material yields or remains elastic. The stress-strain response is plotted at the end.

Uploaded by

mohan
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as TXT, PDF, TXT or read online on Scribd
90 views2 pages

1d Isotropic Linear Hardening

This document contains MATLAB code that is simulating plastic deformation behavior using the elastic-plastic constitutive model. The code defines material properties like Young's modulus and yield strength, then increments strain in small steps. It calculates the stress, plastic strain, and yield surface location at each step based on whether the material yields or remains elastic. The stress-strain response is plotted at the end.

Uploaded by

mohan
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as TXT, PDF, TXT or read online on Scribd
You are on page 1/ 2

clc

clear all

% Input parameters
E = 200e6; % Young's modulus
sigma_yo = 200e5; % Yield strength
K = 25e5; % Hardening modulus
totalstrainforsimulation=1;

% strain increment
delstrain=0.001;
strain=0:delstrain:totalstrainforsimulation;

aa=size(strain);
bb=size(strain,2);

strainp=zeros(aa); %intital state plastic strain


alpha=zeros(aa); % intial yield srface location..
sigma_y=zeros(aa);
delgama=zeros(aa);
sigma=zeros(aa);

sigma_y(1,1)=sigma_yo;
strain(1,1)=0;
strainp(1,1)=0;
alpha(1,1)=strainp(1,1);
delgama(1)=0;
sigma(1,1)=0;

sigma_trail=zeros(aa);
strainp_trail=zeros(aa);
alpha_trail=zeros(aa);
F_trail=zeros(aa);

for i=2:bb

%trail state
sigma_trail(1,i)=sigma(i-1)+(E*delstrain);
%%sigma_trail(1,i)=E*(strain(1,i-1)+delstrain-strainp(1,i-1));
strainp_trail(1,i)=strainp(1,i-1);
alpha_trail(1,i)=alpha(1,i-1);
F_trail(1,i)=abs(sigma_trail(1,i))-(sigma_y(1,1)+(K*alpha(1,i-1))); %%%

if F_trail(1,i)<0
strainp(1,i)=strainp(1,i-1);
alpha(1,i)=alpha(1,i-1); %%%%%%%%%%
sigma(1,i)=sigma_trail(1,i);
sigma_y(1,i)=sigma_y(1,1)+(K*alpha(1,i));%%%%%%%%
%strain(1,i)=strain(1,i-1)+delstrain;

else
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
delgama(1,i)=F_trail(1,i)/(E+K);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
sigma(1,i)=sigma_trail(1,i)-(delgama(1,i)*E*sign(sigma_trail(1,i)));
strainp(1,i)=strainp(1,i-1)+(delgama(1,i)*sign(sigma_trail(1,i)));
alpha(1,i)=alpha(1,i-1)+delgama(1,i);%%%%
sigma_y(1,i)=sigma_y(1,1)+(K*alpha(1,i));%%%%%%%
%strain(1,i)=strain(1,i-1)+delstrain;
end

end
%plot(strain,delgama)

%plot(strain,alpha)

plot(strain, sigma)

You might also like