Creating Spherical Particle Agglomerates in Matlab in 2D or 3D

Question

Chris Dan 2022-3-27

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1681629-creating-spherical-particle-agglomerates-in-matlab-in-2d-or-3d

编辑： Chris Dan 2022-3-31

Hello,

I have a small question and I would like to know any tips or hints or a small code patch on how to proceed with it.

I have to create spherical particle agglomerates in matlab of different sizes all joined together. The number of particles are arbitary, given by me, it could be 5 ,10 ,50, 100 or 500 etc. The particles would be stationary. I have the diameters of the spherical particles and their <x,y> coordinates. At first I have to create particles in 2D and than in 3D. The output I would need after creating the particles would be their <x,y> coordinates again and their visualization so I can compare if they are created correctly or not.

Does anyone how how can I do it in matlab?

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

请先登录，再进行评论。

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

Answer 1

Chris Dan 2022-3-31

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1681629-creating-spherical-particle-agglomerates-in-matlab-in-2d-or-3d#answer_931374

编辑：Chris Dan 2022-3-31

在 MATLAB Online 中打开

close all; clear all;
%% Data for Agglomerates
N = 1000; % number of spheres
a = 1; % lowest diameter of sphere
b = 200 ; % highest diameter of sphere
Diam = a + (b-a).*rand(N,1);
Diam = round(Diam);
aaa= 0;  % minimum x and y coordinate limit for spheres
bbb= sum(Diam) ;% maximum x and y coordinat limit for sphere
M=3  ;
Axes= zeros(N,M);
Axes(:,1)=aaa+(bbb-aaa)*rand(N,1);
for k=2:M
    aaa=randperm(N);
    Axes(:,k)=Axes(aaa,1);
end
Axes_Label ={'X axis','Y axis','Z axis'};
Data_agglo = [Diam Axes];
Data_Label ={'Diameter','X axis','Y axis','Z axis'};
R = Diam ./2;
s = sum(Axes);
%% Plotting Agglomerates in 3D
f = figure('visible','on');
tic
for i =1:1:size(Data_agglo,1)
    
    p = linspace(-Diam(i),Diam(i),200);
    
    [X,Y,Z] = meshgrid(p,p,p); % box mesh
    active = X.^2 + Y.^2 + Z.^2 <= R(i)^2 ;
    plot3(Data_agglo(i,2)+X(active),Data_agglo(i,3)+Y(active),Data_agglo(i,4)+Z(active),'o');
    hold on
end
hold all
toc

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

请先登录，再进行评论。

Answer 2

Chris Dan 2022-3-31

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1681629-creating-spherical-particle-agglomerates-in-matlab-in-2d-or-3d#answer_931309

在 MATLAB Online 中打开

found the answer in 2D

close all; clear all;
%% Data for Agglomerates
N = 100; % number of spheres
a = 1; % lowest diameter of sphere
b = 100 ; % highest diameter of sphere
Diam = a + (b-a).*rand(N,1);
Diam = round(Diam);
aaa= 1;  % minimum x and y coordinate limit for spheres
bbb= sum(Diam) ;% maximum x and y coordinat limit for sphere
M=2  ;
Axes= zeros(N,M);
Axes(:,1)=aaa+(bbb-aaa)*rand(N,1);
for k=2:M
    aaa=randperm(N);
    Axes(:,k)=Axes(aaa,1);
end
Axes_Label ={'X axis','Y axis'};
Data_agglo = [Diam Axes];
Data_Label ={'Diameter','X axis','Y axis'};
R = Diam ./2;
s = sum(Axes);
%% Plotting Agglomerates in 2D
f = figure('visible','on');
for i =1:1:size(Data_agglo,1)
    
    p = linspace(-Diam(i),Diam(i),200);
    
    [X,Y] = meshgrid(p,p); % box mesh
    active = X.^2 + Y.^2  <= R(i)^2 ;
    plot(Data_agglo(i,2)+X(active),Data_agglo(i,3)+Y(active),'o');
    hold on
end