Création d'une sphère ouverte des 2 cotés

Version imprimable

Salut,

Afin de créer une sphère centrée en (x,y,z) de rayon r=1.5mm alors j'ai fait ce qui suit:
Code:

1 2 3 4 5 6 7 8 9 10 11 12 13 [xo,yo,zo] = sphere(30); r=1.5; phi=linspace(0,pi,30); theta=linspace(0,2*pi,40); [phi,theta]=meshgrid(phi,theta); xo=r*sin(phi).*cos(theta); yo=r*sin(phi).*sin(theta); zo=r*cos(phi); surf(xo+x, yo+y, zo+z)
Maintenat je voudrais faire des ouvertures des deux coté de la sphère. comment pourrais-je le faire?
merci:fleur:

Tu ne nous dis rien sur la forme de ces ouvertures...

Essaie ceci :

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
r = 1.5;

phi = linspace(pi*.1,pi*.9,30);
theta = linspace(0,2*pi,40);

[phi,theta] = meshgrid(phi,theta);

x = r*sin(phi).*cos(theta);
y = r*sin(phi).*sin(theta);
z = r*cos(phi);

xo = rand;
yo = rand;
zo = rand;

figure
surf(xo+x, yo+y, zo+z)
axis equal vis3d

Merci beaucoup Dut :ccool:

Voila, grâce à vous j'ai obtenu le ressort que je voulais avoir:
Code:

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 t = 0:pi/30:3*pi; xs=sin(t); ys=cos(t); t1(:,1)=xs; t1(:,2)=ys; t1(:,3)=t; for i=1: length(t1) r = 1; phi = linspace(pi/2,pi/2,30); theta = linspace(0,2*pi,40); [phi,theta] = meshgrid(phi,theta); xo = r*sin(phi).*cos(theta); yo = r*sin(phi).*sin(theta); zo = r*cos(phi); %xo = rand; %yo = rand; %zo = rand; mesh(xo+t1(i,1), yo+t1(i,2), zo+t1(i,3)) end
Maintenant j'essaie de transformer ce mesh en isosurface mais je n'y arrive pas voila ce que j'ai ajouté mais ça ne donne rien:

Code:

patch(isosurface(xo+t1(i,1), yo+t1(i,2), zo+t1(i,3),-3))

Pouvez-vous m'aidez s'il vous plait?

Voici ton code optimisé :

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
figure
hold on
 
t = 0:pi/30:3*pi;
xs = sin(t);
ys = cos(t);
 
r = 1;
phi = linspace(pi/2,pi/2,30);
theta = linspace(0,2*pi,40);
 
[phi,theta] = meshgrid(phi,theta);
 
xo = r*sin(phi).*cos(theta);
yo = r*sin(phi).*sin(theta);
zo = r*cos(phi);
 
n = numel(xs);
 
X = zeros(40,n);
Y = zeros(40,n);
Z = zeros(40,n);
 
for i = 1:n    
    X(:,i) = xo(:,1)+xs(i);
    Y(:,i) = yo(:,1)+ys(i);
    Z(:,i) = zo(:,1)+t(i);
end
 
mesh(X,Y,Z)
 
axis equal vis3d

Que cherches-tu à faire ensuite ?

20/10/2012, 14h51
Décembre

Merci Dut pour l'optimisation :applo:

En fait je cherche à exporter ce ressort en format *.STL, pour le lire dans un logiciel de modélisation 3D tel que blender par exemple.

Pour cela il faudrait que je crée l'isosurface pour avoir la liste des vertex et des faces, j'ai lu sur le help de matlab que c'est faisable avec la fonction isosurface mais l'exemple montre le cas d'une matrice 3D or dans mon cas je n'ai qu'un mesh plotté :aie:
20/10/2012, 16h51
Jerome Briot

Citation:

Envoyé par Décembre

Pour cela il faudrait que je crée l'isosurface pour avoir la liste des vertex et des faces,

Utilise la fonction SURF2PATCH :

Code:

f = surf2patch(X,Y,Z,'triangles');
20/10/2012, 17h21
Décembre

Je ne sais pas comment te remercier Dut toutes les formules de remerciement du monde ne peuvent pas exprimer ma gratitude :rose2:

Merci :fleur:
22/10/2012, 20h23
Décembre
Bonsoir,

Je n'arrive pas à arrêter de poser des questions :help:

Grâce à l'aide de dut j'ai réussi à modéliser un ressort à partir d'une sphère, en le comparant avec les ressorts que j'obtiens avec les logiciels de modélisation 3D j'ai constater que le modéliser avec des sphères était une erreur alors j'ai penser à refaire la même chose avec un cylindre:
Code:

1 2 3 4 5 6 7 8 9 10 11 t = 0:0.01:pi; xs = sin(pi*t); ys = cos(pi*t); [x,y,z] = cylinder(0.5) for i=1:length(t) hold on surf(xs(i)+x, ys(i)+y, t(i)+z) end
Si vous exécuter le programme vous allez comparer que les cylindres sont mis dans le mauvais sens, alors j'aimerais bien vous demander s'il y a moyen de les inverser :question: