The document contains code for modeling the stress and displacement distributions in the plastic and elastic regions around a cylindrical hole. Variables are defined for material properties and geometry. Stresses and displacements are calculated in loops moving radially out from the center of the hole. Results are plotted showing the transition from plastic to elastic behavior with increasing distance from the center.
The document contains code for modeling the stress and displacement distributions in the plastic and elastic regions around a cylindrical hole. Variables are defined for material properties and geometry. Stresses and displacements are calculated in loops moving radially out from the center of the hole. Results are plotted showing the transition from plastic to elastic behavior with increasing distance from the center.
The document contains code for modeling the stress and displacement distributions in the plastic and elastic regions around a cylindrical hole. Variables are defined for material properties and geometry. Stresses and displacements are calculated in loops moving radially out from the center of the hole. Results are plotted showing the transition from plastic to elastic behavior with increasing distance from the center.
The document contains code for modeling the stress and displacement distributions in the plastic and elastic regions around a cylindrical hole. Variables are defined for material properties and geometry. Stresses and displacements are calculated in loops moving radially out from the center of the hole. Results are plotted showing the transition from plastic to elastic behavior with increasing distance from the center.
Download as DOCX, PDF, TXT or read online from Scribd
Download as docx, pdf, or txt
You are on page 1/ 10
a=1;E=6777.9312;vp=.2103448;cohe=3.
45; phi=30;kphi=0;P0=30;Pi=0;reg=5;nx1=100;nx2=100; pai=pi/180; G=E/2/(1+vp); kp=(1+sin(phi*pai))/(1-sin(phi*pai)); q=2*cohe*tan((45+.5*phi)*pai); r00=1/(kp-1); r01=q/(kp-1); r0=a*((2/(kp+1))*(P0+r01)/(Pi+r01))^r00; sre=(2*P0-q)/(kp+1); kps=(1+sin(kphi*pai))/(1-sin(kphi*pai)); c10=(r0/a)^(kp-1); c13=(1-vp)*(kp*kp-1)*(Pi+r01)/(kp+kps); c14=((1-vp)*(kp*kps+1)/(kp+kps)-vp)*(Pi+r01); c15=(2*vp-1)*(P0+r01); delta=(r0-a)/nx1; %stresses in the plastic zone for i=0:nx1 r=a+delta*i; ri(i+1)=r; c11=(r0/r)^(kps+1); c12=(r/a)^(kp-1); pur(i+1)=(r/2/G)*(c15+c13*c10*c11+c14*c12); psxx(i+1)=-r01+(Pi+r01)*c12; psyy(i+1)=-r01+kp*(Pi+r01)*c12; end delta=(reg-r0)/nx2; r1=ri; %stresses in the elastic zone for i=0:nx2 r=r0+delta*i; ri(i+1)=r; eur(i+1)=(P0-(2*P0-q)/(kp+1))*(r0/2/G)*(r0/r); esxx(i+1)=P0-(P0-sre)*(r0/r)^2; esyy(i+1)=P0+(P0-sre)*(r0/r)^2; end rr=[r1,ri];
plot(rr,[pur eur]) hold on plot(ri,eur,'r') line([ri(1) ri(1)],[0 max(pur)+.002],'LineStyle','--','Color','m') text(ri(2),max(pur),'\rightarrow Elastic Region','VerticalAlignment','top') text(r1(end-2),max(pur),'Plastic Region \leftarrow','VerticalAlignment','top','HorizontalAlignment','right') xlabel('Radial distance from center (m)') ylabel('Radial displacement (m)') title('Radial displacement u_r for the cylindrical hole in plastic and elastic regions') figure
plot(rr,[psxx esxx],'k') hold on plot(rr,[psyy esyy]) hold on plot(ri,esxx,'r') hold on plot(ri,esyy,'c') line([ri(1) ri(1)],[0 50],'LineStyle','--','Color','m') text(ri(2),48,'\rightarrow Elastic Region','VerticalAlignment','top') text(r1(end-2),48,'Plastic Region \leftarrow','VerticalAlignment','top','HorizontalAlignment','right') xlabel('Radial distance from center (m)') ylabel('Stress(MPa)') legend('\sigma3','\sigma1','\sigma3','\sigma1',1) title('\sigma1 ,\sigma3 for the cylindrical hole in plastic and elastic regions') figure
z = [psxx esxx]; theta = 0:pi/20:2*pi; xx = bsxfun(@times,rr',cos(theta)); yy = bsxfun(@times,rr',sin(theta)); zz = repmat(z',1,length(theta)); surf(xx,yy,zz,'LineStyle','none') colorbar
view(2) for i=2:9 d=find(.1*fix(10*rr)>.5*i); dd(i-1)=d(1); end title('Minor Principal strss \sigma3 distribution')
