Создание объемной 3D-сетки в форме арахиса в GMSH ⇐ Python

[Anonymous]

Я пытаюсь создать трехмерную объемную сетку в форме арахиса в GMSH (но я мог бы использовать другую платформу) для моделирования методом конечных элементов (FEM).
В Python я использовал gmsh.mode.occ, чтобы создать две сферы и объединить их, а затем расширить их вдоль одной оси:

Код: Выделить всё

s1 = gmsh.model.occ.addSphere(0, 0, r1-off, r1)
s2 = gmsh.model.occ.addSphere(0, 0, -r2 + off, r2)
peanut, _ = gmsh.model.occ.fuse([(3, s1)], [(3, s2)])
gmsh.model.occ.dilate(peanut, 0,0,0,1,1,-dilat)
gmsh.model.occ.synchronize()

Это работает, но приводит к излому поверхности, что проблематично для последующего FEM.
Я решил использовать функцию установки уровня:
f(x,y,z) = ((y-a)² + x² + z²)((y+a)² + x² + y²) - b + c * y
Это дает очень гладкую форму арахиса (см. рисунок).
Используя gmsh из Python, я не смог найти способ создать трехмерную объемную сетку из этой функции набора уровней. Есть какие-нибудь советы?



Подробнее здесь: https://stackoverflow.com/questions/798 ... sh-in-gmsh