Comment assurer des points nodaux dans gmsh?

Question

J'essaye de construire un maillage assez simple. J'ai une boîte:

box_size = 50; 
lb = 10.; 
Point(1) = {-box_size/2, -box_size/2, -box_size/2, lb}; 
Point(2) = {box_size/2, -box_size/2, -box_size/2, lb}; 
Point(3) = {box_size/2, box_size/2, -box_size/2, lb}; 
Point(4) = {-box_size/2, box_size/2, -box_size/2, lb}; 
Line(1) = {1, 2}; 
Line(2) = {2, 3}; 
Line(3) = {3, 4}; 
Line(4) = {4, 1}; 
Line Loop(5) = {1, 2, 3, 4}; 
Plane Surface(6) = {5}; 
Extrude {0, 0, box_size} { 
    Surface{6}; 
} 

Cela fonctionne très bien et gmsh est content de le mailler (3D). Le problème est que je voudrais m'assurer que certains points à l'intérieur de la boîte sont des points nodaux. Donc, ma question est, comment puis-je faire en sorte que les points, comme ceux-ci

lc = 10; 
Point(5) = {7.150548, 1.000000, -6.990684, lc}; 
Point(6) = {-4.438894, 1.000000, -8.960816, lc}; 
Point(7) = {-9.893936, 1.000000, 1.452595, lc}; 
Point(8) = {-1.675894, 1.000000, 9.858569, lc}; 
Point(9) = {8.858176, 1.000000, 4.640336, lc}; 
Point(10) = {1.675894, 4.750000, -9.858569, lc}; 
Point(11) = {-8.858176, 4.750000, -4.640336, lc}; 
Point(12) = {-7.150548, 4.750000, 6.990684, lc}; 
Point(13) = {4.438894, 4.750000, 8.960816, lc}; 
Point(14) = {9.893936, 4.750000, -1.452595, lc}; 
Point(15) = {7.150548, 8.500000, -6.990684, lc}; 

font partie du maillage? La raison pour laquelle j'ai besoin de ceci est que j'ai besoin d'imposer des conditions aux limites à ces points spécifiques.

Si cela est plus facile dans un autre logiciel, je suis également heureux de changer. J'espère que quelqu'un peut aider.

Merci d'avance.

Answer 1

Je crois que la seule option est de diviser la structure de telle sorte qu'ils comportent vos points et maillent la structure après cela. Vous pouvez maintenant appliquer vos charges et conditions sur les points physiques ou les lignes physiques. Exemple: si vous avez un cube que vous voulez mailler. Et la condition aux limites est sur le plan au centre, puis divise le cube sur ce plan. Faites de l'avion une entité physique, par exemple, surface physique (14) = {numéro de plan intermédiaire} .mesh tout et vous êtes prêt à partir!

Answer 2

Une réponse assez tardive mais qui pourrait aider de toute façon. Si l'indice du point est p et celle du volume ou de la surface est q, alors:

Point{p} In Volume {q}; 

Ou si elle est sur une surface:

Point{p} In Surface {q};