PDE meshing (generateMesh) failing after solid operations

Question

Pablo J. Rosado Junquera 2020-6-28

1
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/556270-pde-meshing-generatemesh-failing-after-solid-operations

评论： Juan Peña 2021-9-11

Hey there,

I'm trying to develop an algorithm to define, edit and solve finite elements problems using constructive solid geometry. The libs I'm using work well, but when it comes to generate the mesh, some complex geometries return the "Meshing failed for Hmax..." error whatever the Hmax, Hmin values are.

After some tests, It seems it is the very solid operation which makes the meshing fail, even if the operation is somehow "simplifying" the model, or if it leaves the model exactly the same way. An example:

Original model: Successful mesh generation

Modified (intersection with a cube to obtain a half): failed generation

How could I avoid this? For me, since the non-modified models are of the same complexity, it does not seem a "too complex geometry" problem.

Any workaround appreciated.

Thank you

10 个评论
显示 8更早的评论隐藏 8更早的评论

Pablo J. Rosado Junquera 2020-6-29

在 MATLAB Online 中打开

Sure!

The model I am working with is loaded from a STL file through "stlread". Then, I use "mesh_boolean" from the gptoolbox to perform the solid operation (union, intersect, difference), which outputs a new set of faces (connectivity list) and points. Using this script requires a previous MEX compilation, including CGAL, gmp, libigl and Eigen.

After that, I use the script "triangulation" to generate a triangulation struct from faces and points, which is the required variable type to feed the "stlwrite" function. When the new stl is saved, it is instantly re-loaded using "importGeometry", from the PDE Toolbox. Finally, "generateMesh" is used, also from the PDE Toolbox, to generate the pde model mesh.

Thank you

% Let V2, F2, V1, F1 be the points and faces of the two selected STL's. 
% Intersect solids
[semiV,semiF] = mesh_boolean(V2,F2,V1,F1, 'intersect');
% Write the altered STL model
triangulation_fish = triangulation(semiF,semiV); 
stlwrite(triangulation_fish,'fish.stl');
% Create PDE Model and load the STL
mymodel = createpde;
importGeometry(mymodel,'fish.stl');
% pdegplot(mymodel,'CellLabels','on')
generateMesh(mymodel,'Hmax',0.1);
pdeplot3D(mymodel);

Pablo J. Rosado Junquera 2021-9-5

编辑：Pablo J. Rosado Junquera 2021-9-5

Sure!

I used raw C++ functions compiled by MATLAB using the mex command. The function is called "mesh_boolean".

You will need:

Alecjacobson's gptoolbox .- It contains the MATLAB mesh_boolean function call you want to use. Also you can find useful documentation on how to compile the souce files.
libigl .- It contains the source files needed for you compilation (the actual mesh_boolean fuction writen in C++).
Eigen .- The classic linear algebra C++ lib. It is used by almost any other math - related library.

Theoretically, libigl just needs Eigen, although I used CGAL and Boost too (and their dependencies) because I think this function needs them specifically. Check the documentation carefully to be aware of the libs you need to use.

The last step would be to compile the libraries source files, and then compile the mesh_boolean function directly from MATLAB. In case you are not familiar with C++, to compile your libs:

Download Cmake .- It is a software which reads a txt input file called "Cmakelists" to arrange some compilation options and then generates a solution. (.sln). Let us say it saves us a headache due to the many libs we would need to include and configure to compile by ourselves.
Download Visual Studio .- It contains an easy-to-use compiler. Just open your previously generated .sln and BUILD your libs.

And you are done! Anyway, as I have already said, the gptoolbox documentation is pretty helpful and it is a good starting point.

Cheers!

Juan Peña 2021-9-11

Thanks!!! : )

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Pablo J. Rosado Junquera 2020-9-7

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/556270-pde-meshing-generatemesh-failing-after-solid-operations#answer_491137

I figured it out a coupe weeks ago. There is a script called tetgen.m, from the gptoolbox, which uses an external program, also called Tetgen, to generate the meshes.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

C 2021-3-29

Hi there, Could you please let me know how you used gptoolbox? Did you directly input the model geometry object? I can only see that tetgen() requires a triangle mesh and then triangle() requires a .poly or .node file.

请先登录，再进行评论。

PDE meshing (generateMesh) failing after solid operations

10 个评论
显示 8更早的评论隐藏 8更早的评论

采纳的回答

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

更多回答（0 个）

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

PDE meshing (generateMesh) failing after solid operations

10 个评论 显示 8更早的评论隐藏 8更早的评论

采纳的回答

1 个评论 显示 -1更早的评论隐藏 -1更早的评论

更多回答（0 个）

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

10 个评论
显示 8更早的评论隐藏 8更早的评论

1 个评论
显示 -1更早的评论隐藏 -1更早的评论