# Unintended plot while using the cylinder function

조회 수: 2 (최근 30일)
Kevin Hanekom . 2021년 10월 17일
댓글: Kevin Hanekom . 2021년 10월 18일
Hello!
I have been attempting to create a 3d object by rotating a 2d plot utilizing the cylinder function. But for some reason my plot comes out with this funnel attached to the bottom. Could you possibly bring some insight on how I may be able to create the shape located above the funnel? To be clear the shape above is exactly what we would expect.
Heres the code.
Yieldstr = 1000; % KPa, Sets the str for all following 2d Plots
minValue = -4*Yieldstr;
maxValue = 4*Yieldstr;
fidelity = Yieldstr/125; %able to use much higher fidelity in 2d
sigx = minValue:fidelity:maxValue;
sigy = sigx.';
VM = (1/(sqrt(2)) .* sqrt((sigx-sigy).^2 + (sigy).^2 +(-sigx).^2));
[x,y,~] = find(VM<Yieldstr);
PStress = [sigx(y).' sigy(x)];
A = boundary(PStress,0);
plot(PStress(A,1),PStress(A,2),'LineWidth',1.75);
view(0,90);
colormap jet
[X,Y,Z] = cylinder(PStress);
mesh(X,Y,Z)
view([-32.636 3.808])
Please let me know if you can think of a better way to achieve this 3d shape.
Thank you,
Kevin.

댓글을 달려면 로그인하십시오.

### 채택된 답변

Robert U 2021년 10월 18일
Hi Kevin Hanekom,
it seems your vector PStress has two colums which are both interpreted as radius by the function cylinder and merged accordingly. If just one column is of interest, you would have to call it explicitely.
Yieldstr = 1000; % KPa, Sets the str for all following 2d Plots
minValue = -4*Yieldstr;
maxValue = 4*Yieldstr;
fidelity = Yieldstr/125; %able to use much higher fidelity in 2d
sigx = minValue:fidelity:maxValue;
sigy = sigx.';
VM = (1/(sqrt(2)) .* sqrt((sigx-sigy).^2 + (sigy).^2 +(-sigx).^2));
[x,y,~] = find(VM<Yieldstr);
PStress = [sigx(y).' sigy(x)];
A = boundary(PStress,0);
plot(PStress(A,1),PStress(A,2),'LineWidth',1.75);
view(0,90);
colormap jet
[X,Y,Z] = cylinder(PStress(:,2));
mesh(X,Y,Z)
view([-32.636 3.808])
Kind regards,
Robert
##### 댓글 수: 1없음 표시없음 숨기기
Kevin Hanekom 2021년 10월 18일
Amazing, Thank you!!

댓글을 달려면 로그인하십시오.

### 카테고리

Help CenterFile Exchange에서 Orange에 대해 자세히 알아보기

R2021a

### Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by