Metafor has a few integrated meshers, quite simple but usually capable of meshing parts that are not too complicated. However, 2D meshers require that the
wire delimiting the domain has the right orientation. The wire must be defined with its “area to the left”, which means that when the wire is followed along its orientation (defined by the succession of its curves), the matter of the part is on the left-hand
The auto-detection is working if the
Sides is made of 4 edges, each of which is only made of one
| ||= number of the
In the example above, the auto-detection is doable.
A similar function for higher degrees is:
HighDegreeTransfiniteMesher2D(sideset(number), degree).execute(type, sideNode, tri)
| ||= number of the side
| ||= mesh degree (2 or 3)|
For 9-nodes mesh elements, the central node is linked to the
For 16-nodes mesh elements, central nodes are linked to the
mat = ((mat11, mat12, ..., mat1nbmax), (mat21, mat22, ..., mat2nbmax), (mat31, mat32, ..., mat3nbmax), (mat41, mat42, ..., mat4nbmax) ) TransfiniteMesher2D(sideset(number)).execute2(mat, type)
mat is a python tuple with 4 components. Each component is a tuple which contains the number of the curves which constitute each edge. Each of these four edges can be made of a different number of curves.
This face can be meshed by the command:
TransfiniteMesher2D(sideset(1)).execute2( (1,2,(3,4),5) )
The order of the lines does not matter in 2D (in opposition to the 3D case). Therefore, the following command also works:
TransfiniteMesher2D(sideset(1)).execute2( ((4,3),1,2,5) )
The mesh can be projected on the surface
surfNo (whether with or without auto-detection) if this surface is associated to the
Side before the meshing operation:
This is done with the commands:
MesherTFI2D =TransfiniteMesher2D(sideset(number)) MesherTFI2D.setEnableDistribution() MesherTFI2D.execute(type, tri) or MesherTFI2D.execute2(mat, type)
when one if the edges is already discretized using a mesh elements distribution.
Note : By default, the parameter of the function
setEnableDistribution is set to
Quadrangles (offset + frontal):
d is the average length of the mesh elements edges to generate.
This quadrangular mesher is based on Sarrate algorithm. This method consists in diving the domain in a recursive way until only one quadrangle remains. It is very robust, and enables efficient meshing of complex domains.
Gen4 can be used to mesh:
These lines are used to mesh the face #1 of the domain, whose vertices (points #1, #2, #3 and #4) densities are 0.1, 0.1, 0.2 and 0.05.
from toolbox.meshers import Gen4Mesher defaultDensity = 0.1 mesher = Gen4Mesher(sideset(1), domain, defaultDensity) mesher.setPointD(3,0.2) mesher.setPointD(4,0.05) mesher.execute()
A default density, used unless explicitly stated differently, can be defined. This way, a density is only defined on specific mesh points. If no value is assigned in the data set, $0.1$ is taken as default density value.
CAREFUL : The face #1 must be defined in the XY plane.
Create triangles on a planar face, as a function of