How to develop MATLAB code for 3D Voronoi in cubical volume?

Question

Muhammad Zulkarnanin 2016-2-10

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/267415-how-to-develop-matlab-code-for-3d-voronoi-in-cubical-volume

评论： Oluwayemisi Ademola 2021-9-30

I just started using MATLAB for the research. My research is concentrated on composite material model. I want to develop fiber composite and porosity material model. From my literature study, voronoi technique is powerful to develop fiber and porosity material. In the MATLAB coding, in 2D or cell problem is already develope but I really need to understand in 3D or volume study case using MATLAB

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Mike Garrity 2016-2-10

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/267415-how-to-develop-matlab-code-for-3d-voronoi-in-cubical-volume#answer_209289

What have you tried so far? The voronoin function is the starting point.

3 个评论
显示 1更早的评论隐藏 1更早的评论

Muhammad Zulkarnanin 2016-2-11

在 MATLAB Online 中打开

Dear Mr Mike Garrity Thank you a lot to settle my problem. Now I can run the 3D voronoi, unfortunately when I come to visualization the voronoi perform in surface area. What I need the 3D voronoi look like construction network. Here I present my MATLAB coding below

clear all close all clc

r1=10; %distance between node [nm] rad=r1/2; % radius of distance [nm] numparticles = 200; %Number node launching top=50; %maximum distance of RVE box [nm] num=top/r1; %number of coordinate points along axis (x,y,z)

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % Data save ix_save=zeros(1,length(numparticles)); iy_save=zeros(1,length(numparticles)); iz_save=zeros(1,length(numparticles));

Ixstore=zeros(1,length(numparticles)); Iystore=zeros(1,length(numparticles)); Izstore=zeros(1,length(numparticles));

% Particle seed

    a = [rad:r1:top];
    b = [rad:r1:top]; 
    c = [rad:r1:top];
    xi=a(randint(1,1,[1,num]));
    yi=b(randint(1,1,[1,num]));
    zi=c(randint(1,1,[1,num]));
    Ixstore(1,1)=xi;
    Iystore(1,1)=yi;
    Izstore(1,1)=zi;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%   
%Launching Partikel terhadap partikel terdekat dengan menggunakan random

for i=2:numparticles; ix(i) = a(randint(1,1,[1,num])); iy(i) = b(randint(1,1,[1,num])); iz(i) = c(randint(1,1,[1,num]));

             Ixstore(1,i)=ix(i);
             Iystore(1,i)=iy(i);
             Izstore(1,i)=iz(i);

end X = [Ixstore(1,:)', Iystore(1,:)', Izstore(1,:)'];

[V,C] = voronoin(X); for i=1:length(C), disp(C{i}), end

for k=1:length(C) if all(C{k}~=1) VertCell = V(C{k},:); KVert = convhulln(VertCell); patch('Vertices',VertCell,Faces,KVert,'FaceColor','g','FaceAlpha',0.5)