The intersection of three cylinders

Hi everyone. I want to plot and create the intersection of 3 perpendicular cylinders.

Your help would be greatly appreciated!!

[darova]

Think how to draw one simple segment

Then just copy it

[sylvain jack]

I still can't understand it..Would you be more specific about this? Much thanks..

[Maria Aizaga]

Hi, I have to graph the solid of the 3 cylinders, but I don't know how to do it. I need help, I was guiding myself with your code but the graph does not come out

Answer 1

[darova]

That result of intersection can be build with one simple element

Define boundaries of element:

phi = 0 .. 45 degree
% R - radius of cylinder
X = R*cos(phi);
Y = R*sin(phi);

Z = 0 .. Y;        % Z varies (because of 45 degree)

Then just just changing data with rotating system coordinates:

Then rotate around Z and Y

[Joseph Barreiro]

Could you pass me the code of your example?I am new to this and would like to know how to structure it.

Thanks

Answer 2

[darova]

Here is another example

• start with usual meshgrid to create mesh

• modify Y coordinate

• calculate Z coordinate

• copy object

[Joseph Barreiro]

Your answer helped me, I make the figures separately and in separate they look as if they would combine well, but when one of them simply does not come out the figure I would like .. I am new to matlab, sorry if I ask silly questions

this is the code:

R = 4;
xx = linspace(-1,1,11)*R/sqrt(2);   
[x,y] = meshgrid(xx);               
y = x.*y/sqrt(2)*R;             
z = sqrt(R^2-x.^2);             
h(1) = surf(x,y,z);
h(2) = surface(z,y,x);
hold on
R = 4;
xx = linspace(-1,1,11)*R/sqrt(2);   
[x,y] = meshgrid(xx);               
y = x.*y/sqrt(2)*R;             
z = sqrt(R^2-x.^2);            
u(1) = surf(x,y,z);
u(2) = surface(z,y,x);
direction = [0 0 1];
rotate(u,direction,90)
hold on
R = 4;
xx = linspace(-1,1,11)*R/sqrt(2);   
[x,y] = meshgrid(xx);               
y = x.*y/sqrt(2)*R;        
z = sqrt(R^2-x.^2);
b(1) = surf(x,y,z);
b(2) = surface(z,y,x);
direction = [0 0 1];
rotate(b,direction,180)
view(45,45)
hold on
R = 4;
xx = linspace(-1,1,11)*R/sqrt(2);   
[x,y] = meshgrid(xx);               
y = x.*y/sqrt(2)*R;             
z = sqrt(R^2-x.^2);             
a(1) = surf(x,y,z);
a(2) = surface(z,y,x);
direction = [0 0 1];
rotate(a,direction,270)
view(45,45)

[darova]

You don't need to create x,y,z data each time. Just use copyobj

clc,clear
cla
R = 4;
xx = linspace(-1,1,11)*R/sqrt(2);   
[x,y] = meshgrid(xx);               
y = x.*y*sqrt(2)/R;
z = sqrt(R^2-x.^2);           

    % top part
h1 = surf(x,y,z);
h2 = surface(y,x,z);

    % bottom part
h3 = surf(x,y,-z);
h4 = surt(y,x,-z);

    % create sides
for i = 0:3
    h11 = copyobj(h1);
    h22 = copyobj(h2);
    rotate(h11,[1 0 0],90,[0 0 0])
    rotate(h22,[1 0 0],90,[0 0 0])
    rotate(h11,[0 0 1],90*i,[0 0 0])
    rotate(h22,[0 0 1],90*i,[0 0 0])
end

view(45,45)

[Maria Aizaga]

Hi, I have to graph the solid of the 3 cylinders, but I don't know how to do it. I need help, I was guiding myself with your code but the graph does not come out .

[Joseph Barreiro]

@darova I almost forgot to thank you, you should share your knowledge on youtube, brilliant people like you should be on that kind of platform! Thanks for everything!

[darova]

You are welcome

Answer 3

[darova]

Another example

• create zigzag

• refine it using interp1

• convert it to cylindrical coordinates

• add bottom part

• use different combination (X,Y,Z), (Z,Y,X), (X,Z,Y)