4D plot around a cylinder
1 visualización (últimos 30 días)
Mostrar comentarios más antiguos
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 comentarios
Respuesta aceptada
David Sanchez
el 1 de Ag. de 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.
Más respuestas (0)
Ver también
Categorías
Más información sobre Stress and Strain en Help Center y File Exchange.
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!