4D plot around a cylinder
1 vue (au cours des 30 derniers jours)
Afficher commentaires plus anciens
Hi, I have a formula for stress distribution around a borehole, sigtet(r,T)=(((sig1+sig2)*(1+R^2/a^2))/2)-(((sig1-sig2)*(1+(3*R^4)/a^4)*cos(T*2))/2)-sig12*(1+3*(R/a)^4)*sin(2*T)-delp*(R/a)^2
for a special depth I can enter sig1, sig2 and delp values and get the stress distribution around the wellbore by this code: R=10.5, sig1=100, sig2=35, sig12=0, delp=10,
for ii=1:500
for jj=1:201
a=R+(ii+1)/100;
T=((jj-1)/10)*pi;
sigtet(ii,jj)=(((sig1+sig2)*(1+R^2/a^2))/2)-(((sig1-sig2)*(1+(3*R^4)/a^4)*cos(T*2))/2)-sig12*(1+3*(R/a)^4)*sin(2*T)-delp*(R/a)^2;
x(ii,jj)=(a)*cos(T);
y(ii,jj)=(a)*sin(T);
end
end
pcolor(x,y,sigtet)
shading interp
but sig1 and sig2 and delp values are functions of depth (z), as u see for 2D plot we have stress values in different radius and angles, and the color shows stress intensity for a known depth, now I want to graph it in 3D plot which shows the stress distribution not only for different radius and angles but for different depths. so the new function will be:
sigtet(r,T,z)=(((sig1(z)+sig2(z))*(1+R^2/a^2))/2)-(((sig1(z)-sig2(z))*(1+(3*R^4)/a^4)*cos(T*2))/2)-sig12*(1+3*(R/a)^4)*sin(2*T)-delp(z)*(R/a)^2
0 commentaires
Réponse acceptée
David Sanchez
le 1 Août 2013
I am sorry to disappoint you, but the 4D plot is not yet available. You should choose 3 of your variables at a time and plot them together. Visit (type it)
doc plot3
doc mesh
doc surf
I am sure you will find a plotting function that suits your needs.
Plus de réponses (0)
Voir également
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!