FINAL INDEPENDENT PROJECT

This code is written to calculate the elastic deformation when a disc is subjected to centrifugal force, i have to perform this calulation to find out the deformation at 5mm each for the disc at 900 angular velocity. after getting values i have compared it by FEA (using ansys).

clear all
close all
clc

p= 2710;
w=900;
w1=w^2;
mu= 0.33;
E= 68000000
%hoop stress
for r= linspace(0.005,0.1,20)
X = ((p*w1)/8).*[(3+mu)*(0.1*0.1) - (1+ 3*mu).*(r*r)]

%radial stress

Y = ((p*w1)/8).*[(3 + mu)*((0.1*0.1)- (r*r))]
%elastic strain due to radial
Z = (1/E)+(mu/E)*(X - Y)
Z1 = (1/E)+(mu/E)*(Y - X)
total_strain = Z+Z1;

end