Ex_5_8_1 - Example 5.8-1 Clear the work space clear clear...

Info iconThis preview shows page 1. Sign up to view the full content.

View Full Document Right Arrow Icon
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: Example 5.8-1 Clear the work space clear clear all; close all; clc; Inner radius of the cylinder (inch) a=2/2 a=2/2 a= 1 Outer radius of the cylinder (inch) b=6/2 b=6/2 b= 3 Inside pressure (psi) ip=5000 ip=5000 ip = 5000 Outside pressure (psi) op=1000 op=1000 op = 1000 Range for plotting the stress r=[a:0.001:b]'; r=[a:0.001:b]'; Stress distribution (subjected to either/both inside and outside pressure) ri=a; ri=a; ro=b; ip=5000; op=0; Radial Stress Distribution sr1=( (ip*ri*ri)((ri*ro./r).^2)* (op)/(ro^2sr1=( (ip*ri*ri)- (op*ro*ro) + ( ((ri*ro./r).^2)* (op-ip) ) )/(ro^2-ri^2); Tangential Stress Distribution st1=( (ip*ri*ri)(op)/(ro^2st1=( (ip*ri*ri)- (op*ro*ro) - ( ((ri*ro./r).^2)* (op-ip) ) )/(ro^2-ri^2); Longitudinal Stress Distribution sz1=( (ip*ri*ri))/(ro^2sz1=( (ip*ri*ri)- (op*ro*ro) )/(ro^2-ri^2); Von-Mises Stress Distribution sVM1=(0.5*( (st1(sz1(sr1sVM1=(0.5*( (st1-sz1).^2 + (sz1-sr1).^2 + (sr1-st1).^2 ) ).^0.5; figure; plot(r,st1/1000,'--',r,sr1/1000,'--',r,sr1/1000,' plot(r,st1/1000,'--',r,sr1/1000,'-.') legend('Tangential','Radial') title('Stress distribution in thick cylinder') xlabel('Radius') ylabel('Stress, (in kpsi)') grid on figure; plot(r,sVM1/1000,' t(r,sVM1/1000,'plot(r,sVM1/1000,'-') legend('Vonlegend('Von-Mises Stress') title('Vontitle('Von-Mises Stress distribution in thick cylinder') xlabel('Radius') ylabel('Stress, (in kpsi)') grid on ...
View Full Document

Ask a homework question - tutors are online