Discussion :

[conanran]

Bonjour,
En effet, afin de générer une trajectoire spirale d’outil d’usinage j’ai utilisé les solutions d’une PDE, pour avoir les contours. Et à partir ces contours je veux créer un spirale (comme montre l’image).
J’ai utilisé ce code et j’ai obtenu les points des contours (ContourMatrix) mais il apparait qu’il faut créer des point équi-angulaire sur chaque contour pour créer un spirale.
Merci pour votre aide

ci joint l'image et le code


le code

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
%% Création du model
g0 = [3,4,-.6,.6,.6,-.6,-.6,-.6,.6,.6]'; 
sf = 'C1';
g = decsg(g0,sf,sf'); % decomposed geometry matrix
problem = allzerobc(g); % zero Dirichlet conditions
[p,e,t] = initmesh(g,'Hmax',0.02);
c = 1;
a = 0;
f = 1;
%% solve the PDE
u = assempde(problem,p,e,t,c,a,f);
x = linspace(-.6,.6);
y = linspace(-.6,.6);
uxy=tri2grid(p,t,u,x,y);
[C,h]=contour(x,y,uxy,6);
set(h,'TextStep',get(h,'LevelStep')*0.1)
clabel(C)