How to combine two implicit functions

I want to give thickness to a Gyroid Function for Gyroid is

Gyroid=(cos(x) .* sin(y) + cos(y) .* sin(z) + cos(z) .* sin(x));
% Reverse Gyroid
R=sqrt( (-sin(x).*sin(y)+cos(x).*cos(z)).^2+(-sin(y).*sin(z)+cos(y).*cos(x)).^2+(-sin(z).*sin(x)+cos(z).*cos(y)).^2 );
ReverseGyroid=(Gyroid.*T./R);
isosurface(Gyroid)
isosurface(ReverseGyroid)

This is the output I expect but as one Isosurface def and not by using isosurface separately for both surfaces.

Thank you in advance

[Jim Joy]

Darpan,

Could you please clarify what you mean when you say that you expect this output with a single call to "isosurface"?

Do you expect the command below to give the output shown?

isosurface(Gyroid)

Best,

Jim

[Darpan Shidid]

Hi Jim,

Basically I want to add thickness to the Gyroid and close the gap between two gyroids using isocaps. So fill the gap between the blue and yellow surfaces. Hoping to find a simple solution.

[Jim Joy]

Hi Darpan,

Thank you for clarifying that you would like to fill the space between the two Gyroids.

To clarify further: are you constructing the 'reverse Gyroid' for visualization purposes, or is there a mathematical significance to the space between the Gyroids?

I am asking, because the "isocaps" function might be useful for you. Does the command below give a result similar to what you expect?

isocaps(x,y,z,Gyroid,0)

In general, filling the space between arbitrary surfaces is a bit trickier. The MATLAB Answers post here details one approach for doing so, which I believe could be adapted to this scenario. The basic idea is to create a 3D grid of data in which the value at a point (x,y,z) is the distance to the nearest surface. You can then use "isocaps" to connect the surfaces.

If you would like further suggestions about how to do this, please post your code for generating the 'reverse gyroid'. In particular, please post how to generate the variable 'T'. This will help users on this forum to experiment with ideas that may help provide a more complete solution.

Best Regards,

Jim

[Darpan Shidid]

Hi Jim, Yes I want to use ISOCAPS for filling the gap. The code is as supplied in the Question. The var T is just thickness that can be 0.5. Currently if I use the isocaps command it will fill the two isosurfaces separately. I want a thick Gyroid just like the image attached. How can I achieve that using Isocaps? Here is the code again.

Gyroid=(cos(x) .* sin(y) + cos(y) .* sin(z) + cos(z) .* sin(x));
% Reverse Gyroid
R=sqrt( (-sin(x).*sin(y)+cos(x).*cos(z)).^2+(-sin(y).*sin(z)+cos(y).*cos(x)).^2+(-sin(z).*sin(x)+cos(z).*cos(y)).^2 );
ReverseGyroid=(Gyroid.*T./R);
isosurface(Gyroid)
isosurface(ReverseGyroid)

Answer 1

[Paulo GOMES]

clear all

close all

clc

x_period = 4;

y_period = 4;

z_period = 4;

[x,y,z] = meshgrid(0:0.1:x_period*pi, 0:0.1:y_period*pi, 0:0.1:z_period*pi);

Gyroid=(cos(x) .* sin(y) + cos(y) .* sin(z) + cos(z) .* sin(x));

isosurface(x,y,z,Gyroid)

isocaps(x,y,z,Gyroid,0)

If you do something like this, I think you will get the solution you want!

[Yash Mistry]

hi, i am getting the error of too many input arguments for your code(isosurface) in your code. can you help me with it.

Answer 2

[Deison Preve]

hi Yash Mistry, did you manage to do so? I am trying to do same thing.. add some thickness on the TPMS isosurface.. If you did, could you share the lane?

Answer 3

[Fabian Günther]

The attached code generates a stl file of a Gyroid shell arrangement when you include the stlwrite function.

clc
clear
close all


n=200;
t=pi/n;

x_max=2;
y_max=2;
z_max=2;

xi = 0:t:x_max;
yi = 0:t:y_max;
zi = 0:t:z_max;

[x,y,z] = meshgrid(xi,yi,zi);

ri=-0.35;
ra=-ri;

F=cos(2.*pi.*x).*sin(2.*pi.*y)+cos(2.*pi.*y).*sin(2.*pi.*z)+cos(2.*pi.*z).*sin(2.*pi.*x); 
F=-(F+ri).*(F+ra);


[fs,v]=isosurface(x,y,z,F,0);
[fc,v2,c] = isocaps(x,y,z,F,0);
fn = [fs ; fc+length(v(:,1))];
vn = [v ; v2];
stlwrite('geometry.stl',fn,vn);

[NickKnack CAD]

hi @Fabian Günther your code just happens to be what ive been searching for.

Could you please provide a small explanation of whats in the last part of your code after the gyroid surface gets offset to just above the WRITESTL statement?

[Kamal Uddin Mohammad Miah]

thank you @Fabian Günther