Discussion :

[yoannthevenoux]

Bonjour à tous,

Je cherche actuellement à représenter la distribution de pression dans un cylindre avec MATLAB.
Une image sera plus utile qu'un long discours, je voudrais avoir un résultat à peu près comme ça :


J'ai une expression de la pression, en coordonnées cylindriques. Et après avoir essayer plusieurs solutions trouvées sur internet je n'ai pas trouvé mon bonheur.

Dans l'absolu, je n'ai pas besoin de représenter tout mon volume. Une solution à laquelle j'ai pensé est de créer un cylindre comme sur cette image
auquel on ajouterait les surfaces à chaque extrémité. Mais je n'ai pas réussi à réaliser cela.

Quelqu'un a-t-il une solution?

PS:Mon niveau en MATLAB est plutôt débutant

[Jerome Briot]

Quelle est la forme de tes données ?
Sont-elles réparties sur une grille régulière le long du cylindre ?

[yoannthevenoux]

J'ai simplement une equation de la forme (j'ai mis les variables en gras) :

p(r,ϕ,z)= cos⁡(lπz/L).cos⁡(nϕ).J(α.r/R)
avec J une fonction de bessel du premier ordre.

donc pour la répartition des valeurs c'est moi qui décide.
Je peut donc créer les vecteurs R, ϕ et z avec un espacement régulier.

[Jerome Briot]

Un exemple :

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
[theta,z] = meshgrid(0:pi/50:2*pi,1:30);
rho = ones(size(theta));
p = z.*cos(theta);
 
[x,y,z] = pol2cart(theta,rho,z);
 
figure
 
surf(x,y,z,p,'edgecolor','none')
 
hold on
 
patch(x(1,:),y(1,:),z(1,:),p(1,:)); % Face inferieure
patch(x(end,:),y(end,:),z(end,:),p(end,:)); % Face superieure
 
axis vis3d

[yoannthevenoux]

Merci beaucoup pour ta réponse et ta rapidité,
Je testerais ça lundi matin! (je n'ai pas matlab chez moi!)

[yoannthevenoux]

J'ai aujourd'hui essayer le code proposé, mais malheuresement un problème persiste. En effet, sur les deux surfaces patch, la répartition des couleurs ne se fait pas suivant mon équation. Sur l'image suivante, il y a à gauche le résultat que j'ai actuellement. Or je voudrais que la surface "patch" soit comme celle de droite (image qui a été créé avec surf, donc c'est une surface 3D, qui est ici vu du dessus)


Est-il possible de modifier la fonction patch pour répondre à mon besoin?
(Je n'ai pas encore vraiment compris comment fonctionne la fonction patch)

Edit: Dans le code au dessus il y a :

Code :

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

rho = ones(size(theta));

Dans ce cas rho a toujours la meme valeur, pour mes surfaces inférieurs et supérieur, dois-je modifié cela?