FEM MeshGmshFromShape/en: Difference between revisions

From FreeCAD Documentation
(Updating to match new version of source page)
(Updating to match new version of source page)
Line 20: Line 20:
==Description==
==Description==


For a finite elements analysis the geometry needs to be discretized into a [[FEM_Mesh|FEM Mesh]]. This command uses gmsh (which needs to be installed on the system) for calculating the mesh.
For a finite elements analysis the geometry needs to be discretized into a [[FEM_Mesh|FEM Mesh]]. This command uses the program [https://en.wikipedia.org/wiki/Gmsh Gmsh] (which needs to be installed on the system) for calculating the mesh.

Gmsh is bundled with the FreeCAD installation binaries. Alternatively you can install it separately from FreeCAD and then use the menu {{MenuCommand|Edit → Preferences → FEM → Gmsh}} to set the path to the ''gmsh.exe''.


==Usage==
==Usage==


# Select the shape you want to analyze. For volume FEM this needs to be a solid or compsolid. A compsolid is necessary if your part is made from multiple materials. (A compsolid can be created with the [[Part_BooleanFragments|BooleanFragments]] command.) -- For shell and beam FEM somebody needs to fill in the details here.
# Select the shape you want to analyze. For volume FEM this needs to be a solid or compsolid. A compsolid is necessary if your part is made from multiple materials. (A compsolid can be created with the [[Part_BooleanFragments|BooleanFragments]] command.)
# There are several ways to invoke the command:
#* Press the {{Button|[[Image:FEM_MeshGmshFromShape.svg|16px]] [[FEM_MeshGmshFromShape|FEM MeshGmshFromShape]]}} button.
#* Press the {{Button|[[Image:FEM_MeshGmshFromShape.svg|16px]] [[FEM_MeshGmshFromShape|FEM MeshGmshFromShape]]}} button.
#* Select the {{MenuCommand|Mesh → [[Image:FEM_MeshGmshFromShape.svg|16px]] FEM mesh from shape by Gmsh}} option from the menu.
#* Select the {{MenuCommand|Mesh → [[Image:FEM_MeshGmshFromShape.svg|16px]] FEM mesh from shape by Gmsh}} option from the menu.
Line 32: Line 33:
# Close the task. You now should see a new FEMMeshGMSH object in your active analysis container.
# Close the task. You now should see a new FEMMeshGMSH object in your active analysis container.


After the mesh has been crated you can change its properties using the [[Property_editor|property editor]]. After you changed a property, you must reopen the Gmsh dialog again and click the {{Button|Apply}} button. (You can leave the dialog open while changing properties.)

== Properties ==

* {{PropertyData|Algorithm2D}}: The algorithm to create 2D meshes. The different algorithms are [https://gmsh.info/doc/texinfo/gmsh.html#Choosing-the-right-unstructured-algorithm explained here]. For Delaunay, see [https://en.wikipedia.org/wiki/Delaunay_triangulation Delaunay triangulation].
* {{PropertyData|Algorithm3D}}: The algorithm to create 3D meshes. The different algorithms are [https://gmsh.info/doc/texinfo/gmsh.html#Choosing-the-right-unstructured-algorithm explained here].
* {{PropertyData|Characteristic Length Max}}: The maximal size of the mesh elements. If set to ''0.0'', the size will be set automatically. This property can also be changed in the Gmsh dialog in the field '''Max element size'''.
* {{PropertyData|Characteristic Length Min}}: The minimal size of the mesh elements. If set to ''0.0'', the size will be set automatically. This property can also be changed in the Gmsh dialog in the field '''Min element size'''.
* {{PropertyData|Coherence Mesh}}:
** true (default); duplicate mesh nodes will be removed
** false
* {{PropertyData|Element Dimension}}: The dimension of the mesh elements. This property can also be changed in the Gmsh dialog in the field '''Mesh element dimension'''.
** From Shape (default); the dimension will be determined from the dimension of the object that is meshed
** 1D
** 2D
** 3D
* {{PropertyData|Element Order}}: The [https://www.comsol.de/support/knowledgebase/1270 mesh element order]. This property can also be changed in the Gmsh dialog in the field '''Mesh order'''. {{Version|0.20}}
** 1st
** 2nd (default)
* {{PropertyData|Mesh Size From Curvature}} {{Version|0.20}}: The number of mesh elements per <math>2\pi</math> times the radius of the curvature. To get a finer mesh at small corners or holes, this value can be increased for better results
[[Image:FEM_Gmsh-MeshSizeFromCurvature.png|450px]]
{{Caption|Effect of ''Mesh Size From Curvature'''; left: set to 12, right: deactivated}}
* {{PropertyData|Second Order Linear}}: Option if second order nodes (if {{PropertyData|Element Order}} set to ''2nd'') and/or mesh refinement points are created by linear interpolation.
** true; linear interpolation is used
** false (default); curvilinear interpolation is used

==Notes==

=== Nonpositive Jacobians ===

When you get a meshing erro about nonpositive Jacobians, you can try out the following strategies:

* Set {{PropertyData|Second Order Linear}} to ''true'' but keep {{PropertyData|Element Order}} at ''2nd''.
* Set {{PropertyData|Element Order}} to ''1st''.
* Use a smaller element size by reducing the {{PropertyData|Characteristic Length Max}}.

===Mesh Growth===

At edges and small geometric entities the mesh has to be smaller than in areas without edges. So the mesh element size grows away from edges. The growing strategy of Gmsh is to grow between edges with different sizes. So the growing fails when an area has the same sized edges like for example this tube:

[[Image:FEM_Gmsh-MeshGrowth-failing.png|400px]]
{{Caption|Failing mesh growing because the cylindrical area is surrounded be the same edges}}

To enable a sensible mesh growing, you must in this case add an edge to the area. In the example this would be a circle in the middle of the cylinder. The circle is added as part of a [[Part_BooleanFragments|BooleanFragments]] compound (to form a CompSolid), see [https://forum.freecadweb.org/download/file.php?id=146255 the project file] of the example.

[[Image:FEM_Gmsh-MeshGrowth-success.png|400px]]
{{Caption|Sensible mesh growing due to the additional edge in the middle of the cylindrical aread}}


{{Docnav
{{Docnav

Revision as of 20:08, 1 April 2021

This documentation is not finished. Please help and contribute documentation.

GuiCommand model explains how commands should be documented. Browse Category:UnfinishedDocu to see more incomplete pages like this one. See Category:Command Reference for all commands.

See WikiPages to learn about editing the wiki pages, and go to Help FreeCAD to learn about other ways in which you can contribute.

FEM MeshGmshFromShape

Menu location
Mesh → FEM mesh from shape by Gmsh
Workbenches
FEM
Default shortcut
None
Introduced in version
-
See also
FEM tutorial

Description

For a finite elements analysis the geometry needs to be discretized into a FEM Mesh. This command uses the program Gmsh (which needs to be installed on the system) for calculating the mesh.

Gmsh is bundled with the FreeCAD installation binaries. Alternatively you can install it separately from FreeCAD and then use the menu Edit → Preferences → FEM → Gmsh to set the path to the gmsh.exe.

Usage

  1. Select the shape you want to analyze. For volume FEM this needs to be a solid or compsolid. A compsolid is necessary if your part is made from multiple materials. (A compsolid can be created with the BooleanFragments command.)
    • Press the FEM MeshGmshFromShape button.
    • Select the Mesh → FEM mesh from shape by Gmsh option from the menu.
  2. Optionally edit the minimal and maximal element size. (Autodetection works fine unless you apply complicated boundary conditions.)
  3. Click the Apply button and wait for the computation of the mesh to complete
  4. Close the task. You now should see a new FEMMeshGMSH object in your active analysis container.

After the mesh has been crated you can change its properties using the property editor. After you changed a property, you must reopen the Gmsh dialog again and click the Apply button. (You can leave the dialog open while changing properties.)

Properties

  • DataAlgorithm2D: The algorithm to create 2D meshes. The different algorithms are explained here. For Delaunay, see Delaunay triangulation.
  • DataAlgorithm3D: The algorithm to create 3D meshes. The different algorithms are explained here.
  • DataCharacteristic Length Max: The maximal size of the mesh elements. If set to 0.0, the size will be set automatically. This property can also be changed in the Gmsh dialog in the field Max element size.
  • DataCharacteristic Length Min: The minimal size of the mesh elements. If set to 0.0, the size will be set automatically. This property can also be changed in the Gmsh dialog in the field Min element size.
  • DataCoherence Mesh:
    • true (default); duplicate mesh nodes will be removed
    • false
  • DataElement Dimension: The dimension of the mesh elements. This property can also be changed in the Gmsh dialog in the field Mesh element dimension.
    • From Shape (default); the dimension will be determined from the dimension of the object that is meshed
    • 1D
    • 2D
    • 3D
  • DataElement Order: The mesh element order. This property can also be changed in the Gmsh dialog in the field Mesh order. introduced in version 0.20
    • 1st
    • 2nd (default)
  • DataMesh Size From Curvature introduced in version 0.20: The number of mesh elements per times the radius of the curvature. To get a finer mesh at small corners or holes, this value can be increased for better results

Effect of Mesh Size From Curvature'; left: set to 12, right: deactivated

  • DataSecond Order Linear: Option if second order nodes (if DataElement Order set to 2nd) and/or mesh refinement points are created by linear interpolation.
    • true; linear interpolation is used
    • false (default); curvilinear interpolation is used

Notes

Nonpositive Jacobians

When you get a meshing erro about nonpositive Jacobians, you can try out the following strategies:

  • Set DataSecond Order Linear to true but keep DataElement Order at 2nd.
  • Set DataElement Order to 1st.
  • Use a smaller element size by reducing the DataCharacteristic Length Max.

Mesh Growth

At edges and small geometric entities the mesh has to be smaller than in areas without edges. So the mesh element size grows away from edges. The growing strategy of Gmsh is to grow between edges with different sizes. So the growing fails when an area has the same sized edges like for example this tube:

Failing mesh growing because the cylindrical area is surrounded be the same edges

To enable a sensible mesh growing, you must in this case add an edge to the area. In the example this would be a circle in the middle of the cylinder. The circle is added as part of a BooleanFragments compound (to form a CompSolid), see the project file of the example.

Sensible mesh growing due to the additional edge in the middle of the cylindrical aread