Discussion :

[stage75]

Bonjour à tous,

Je débute avec Matlab et je sollicite votre aide, je veux créer via un code numérique Matlab une géométrie carrée (4 cm/4cm) et de faire le maillage de ce carrée (un maillage de 40 elements/40 elements, et j'en ai aucune aidée sur la méthode à utiliser !

Est ce quelqu'un pourrait m'aider svp

[thecrazydonut]

Bonjour,

Ça peut être ?

Code :

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

1
2
[x,y]=meshgrid(1:(4-1)/40:4,1:(4-1)/40:4);
plot(x,y,x',y');

[Jerome Briot]

Il faudrait être plus précis sur ce que tu souhaites obtenir au final.

As-tu uniquement besoin d'afficher le maillage ou bien cherches-tu à obtenir la table de connectivité des éléments ?

[stage75]

Merci bcp pour vos réponses,

Effectivement je veux afficher aussi la connectivité

[Kangourou]

salut,

pour la connexité entre les noeuds de ton maillage, c'est un peu plus dur. Une solution est de calculer une table des indices des noeuds, puis d'extraire les indices adjacents dans la table. 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
17
18
19
20
21
22
23
24
inds = [...
 1  5  9  13 17 ; ...    
 2  6 10 14 18 ; ...
 3  7 11 15 19 ; ...
 4  8 12 16 20 ];
 
% tableaux des indices des coins de chaque face
v1 = inds(1:end-1, 1:end-1);
v2 = inds(1:end-1, 2:end);
v3 = inds(2:end, 1:end-1);
v4 = inds(2:end, 2:end);
 
% adjacences horizontales
vx1 = inds(:, 1:end-1);
vx2 = inds(:, 2:end);
adjX = [vx1(:) vx2(:)];
 
% adjacences verticales
vy1 = inds(1:end-1, :)';
vy2 = inds(2:end, :)';
adjY = [vy1(:) vy2(:)];
 
% toutes les adjacences
adjList = [adjH ; adjV];

à tester...

[Jerome Briot]

Une solution :

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
18
19
20
21
22
23
24
25
nElmtsOnSide = 40;
squareSize = 4; %cm
seeds = linspace(0, squareSize, nElmtsOnSide+1);
[coordX,coordY] = meshgrid(seeds);
 
nodes = [coordX(:) coordY(:)];
 
indices = reshape(1:(nElmtsOnSide+1)*(nElmtsOnSide+1), nElmtsOnSide+1, []);
 
elmts = zeros(nElmtsOnSide*nElmtsOnSide,4);
 
for i = 1:nElmtsOnSide
    for j = 1:nElmtsOnSide
        idx = sub2ind(size(indices), [i i i+1 i+1], [j j+1 j+1 j]);
        elmts((i-1)*nElmtsOnSide+j,:) = indices(idx);
    end
end
 
figure
 
patch('vertices', nodes, 'faces', elmts, 'facecolor', 'g')
hold on
plot(coordX, coordY, 'r*');
 
axis equal tight