Draw partial spheroid include a spheroid

7 次查看（过去 30 天）

KOU DU 2019-7-29

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/473950-draw-partial-spheroid-include-a-spheroid

评论： Kim 2019-10-17

I want to draw 1/8 spheroid include a small spheroid and output the geometry for mesh. But my current coding always have discontinue in the cutting plan.

Can anyone help provide a idea of cutting the spheroid in 1/8 not for showing but get the data.

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

请先登录，再进行评论。

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

采纳的回答

Bruno Luong 2019-8-21

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/473950-draw-partial-spheroid-include-a-spheroid#answer_388516

编辑：Bruno Luong 2019-8-21

在 MATLAB Online 中打开

The code bellow us this FEX to generate mesh points.

% radius of the inner/outer spherical parts
r1 = 1;
r2 = 2;
n = 20; % discretization parameters of spherical parts
W = allVL1(3, n); % FEX file
% Connectivity
XY = W*[0 1 0;
        0 0 1].';
F=delaunay(XY);
% Points in S2
Tri3 = eye(3);
W = W/n;
XYZ = W*Tri3';
XYZ = XYZ ./ sqrt(sum(XYZ.^2,2));
% inner/outer sphericals
XYZ1 = XYZ*r1;
XYZ2 = XYZ*r2;
% TRUE for points on the boundary
ibdr = W==0; 
close all
patcharg = {'FaceColor', 'g', 'FaceAlpha', 0.5};
patch('Faces', F, 'Vertices', XYZ1, patcharg{:});
patch('Faces', F, 'Vertices', XYZ2, patcharg{:});
for k=1:3
    XYZk = [XYZ1(ibdr(:,k),:); flipud(XYZ2(ibdr(:,k),:))];
    % You are free to mesh XYZk, I leave it as polygonal shape (quarter of a rings)
    patch('XData', XYZk(:,1), 'YData', XYZk(:,2), 'ZData', XYZk(:,3), patcharg{:});
end
view(3)
axis equal

13 个评论
显示 11更早的评论隐藏 11更早的评论

Bruno Luong 2019-8-22

编辑：Bruno Luong 2019-8-22

在 MATLAB Online 中打开

Well, I'll be kind for once and ended up doing the whole thing for you

 % radius of the outer spherical
r1 = 1;
% parameter of inner part
a=0.2; b=0.2; c=0.2; p=0.7;
n = 20; % discretization parameters of spherical parts
W = allVL1(3, n); % FEX file
% Connectivity
XY = W*[0 1 0;
        0 0 1].';
F = delaunay(XY);
% Points in S2
XYZ = W/n;
% inner/outer sphericals
XYZ1 = XYZ .* (r1./ sqrt(sum(XYZ.^2,2)));
q = 2*p;
qnorm = sum((XYZ ./ [a,b,c]).^q,2).^(1/q); % add abs if components are negative
XYZ2 = XYZ ./ qnorm;
% TRUE for points on the boundary
ibdr = W==0;
% Group vertices
XYZ = [XYZ1; XYZ2];
m = size(XYZ1,1);
Fout = 0 + F;
Fin = m + fliplr(F);
Fwall = cell(1,3);
for k=1:3
    bk = find(ibdr(:,k));
    Fk = [[bk(1:end-1) bk(2:end) bk(1:end-1)]+[0 0 m];
          [bk(2:end) bk(1:end-1) bk(2:end)]+[m m 0]];
    % Re-orienting triangular faces so that they are consistently oriented
    T = reshape(XYZ(Fk,:),[],3,3);        
    N = cross(T(:,2,:)-T(:,1,:),T(:,3,:)-T(:,1,:),3);
    reverse = N(:,:,k) > 0;
    Fk(reverse,:) = fliplr(Fk(reverse,:));
    Fwall{k} = Fk;
end
Faces = cat(1,Fin,Fout,Fwall{:});
close all
patcharg = {'FaceColor', 'g', 'FaceAlpha', 0.7};
patch('Faces', Faces, 'Vertices', XYZ, patcharg{:});
view(3)
axis equal
stlwrite('watermeloon.stl', Faces, XYZ)

KOU DU 2019-8-22

Thank you!

Kim 2019-10-17

Since I have a similar problem, I tried to compile this programm but as in my case, I always get the error "Input argument must be a triangulation object."

Could anybody point out, what mistake I've been making?

请先登录，再进行评论。

更多回答（2 个）

darova 2019-8-19

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/473950-draw-partial-spheroid-include-a-spheroid#answer_388144

在 MATLAB Online 中打开

Use patch() to generate planes

clc,clear
R = 10;
r = 3;
t = linspace(0,pi/2,20);
x = [r*cos(t) fliplr(R*cos(t))];
y = [r*sin(t) fliplr(R*sin(t))];
patch(x,y,x*0,'b')
hold on
patch(x,x*0,y,'b')
patch(x*0,x,y,'b')
hold off
alpha(0.5)
view(3)

13 个评论
显示 11更早的评论隐藏 11更早的评论

KOU DU 2019-8-22

Thanks darova. Sorry I don't see the question of shapes. The inner shape is not only a ellispoid.

KOU DU 2019-8-22

And as you said. I try to introduce p but not success. Whatever, thank you very much.

请先登录，再进行评论。

Abhisek Pradhan 2019-8-7

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/473950-draw-partial-spheroid-include-a-spheroid#answer_386540

在 MATLAB Online 中打开

Following code may be used as an alternative to draw a sphere. Theta and Phi can be varied to get the desired result.

R=10;
Phi=linspace(-pi,pi);
Theta=linspace(0,2*pi);
[Phi,Theta]=meshgrid(Phi,Theta);
Z=R*sin(Phi);
X=R*cos(Phi).*cos(Theta);
Y=R*cos(Phi).*sin(Theta);
hSurface  = surf(X,Y,Z);
set(hSurface,'FaceColor',[0 0 1], 'FaceAlpha',0.5,'FaceLighting','gouraud','EdgeColor','none');