How can I do a drawing of this shape؟
2 views (last 30 days)
Show older comments
%File of Result of stress by FEM v=0.35;% Poission ratio sv=66;% vertical stress sh=45;% min horizontal stress Pp=30;% pore pressure sH=75; % max horizontal stress a=0.08;
m=0:39; r=a+(5..*10.^-3+0.1.*m); n=linspace(1,120,40) theta= (3./2+3.*(n-1)) az=0; inc=0; b=((a.^2)./(r.^2)); N=40; % eqn 3.2.6 sx = sH.*(cos(az).* cos(r)).^2 + sh.*(sin(az).*cos(inc)).^2 + sv.*(-sin(inc).^2); sy = sH.*(-sin(az)).^2 + sh.*(cos(az)^2); sz = sH.*((cos(az).* sin(inc)).^2) + sh.*(sin(az).^2).*(sin(inc).^2)+ sv.*(cos(inc).^2); syz = (sH.*-sin(az).*sin(inc).*cos(az))+(sh.*cos(az).* sin(inc).*sin(az))+(sv.*cos(inc).*0); sxz = (sH.*cos(az).* cos(az).*sin(r).*sin(r))+(sh.*sin(az).*cos(inc).*sin(az).*sin(inc))+(sv.*- sin(r).*cos(inc)); sxy =(sH*cos(az)* cos(inc)*-sin(az))+ (sh*sin(az)*cos(inc)*cos(az))+(sv*-sin(inc)*0); sigma_tt = 0.5.*(sx + sy).*(1 +b) - 0.5.*(sx - sy).*(1 + 3.*b).*cos(2.*theta) - sxy.*(1 + 3.*b).*sin(2.*theta);
3 Comments
Image Analyst
on 17 Aug 2018
Not sure what the code does. I guess it does something but not what you want.
For what it's worth, I've attached my color wheel demo. Adapt as needed.
Accepted Answer
More Answers (0)
See Also
Categories
Find more on Stress and Strain in Help Center and File Exchange
Products
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!