Body Fitting 3D meshing algorithm

Body Fitting algorithm generates hexahedrons of a Cartesian grid in the internal part of geometry and polyhedrons and other types of elements at the intersection of Cartesian cells with the geometrical boundary. The algorithm supports construction of viscous boundary layers (use Viscous Layers hypothesis to define them).

_images/cartesian3D_sphere.png

A sphere meshed by Body Fitting algorithm

The meshing algorithm is as follows.

  1. Lines of a Cartesian structured grid defined by Body Fitting Parameters hypothesis are intersected with the geometry boundary, thus nodes lying on the boundary are found. This step also allows finding out for each node of the Cartesian grid if it is inside or outside the geometry.

  2. For each cell of the grid, check how many of its nodes are outside of the geometry boundary. Depending on a result of this check

    • skip a cell, if all its nodes are outside

    • skip a cell, if it is too small according to Size Threshold parameter

    • add a hexahedron in the mesh, if all nodes are inside

    • add a polyhedron or another cell type in the mesh, if some nodes are inside and some outside.

Viscous boundary layers are constructed as follows:

  • create an offset geometry with offset value equal to full layers thickness by using BRepOffset_MakeOffset OCCT class;

  • mesh the offset geometry with the Body Fitting algorithm;

  • create prisms of the layers by projecting boundary nodes of offset geometry onto the boundary of initial geometry.

Note

Viscous layers can’t be constructed on geometry with shared/internal faces.

To apply the algorithm when you define your mesh, select Body Fitting in the list of 3D algorithms and add Body Fitting Parameters hypothesis. The following dialog will appear:

Body Fitting Parameters hypothesis

_images/cartesian3D_hyp.png

Body Fitting Parameters hypothesis dialog

This dialog allows to define

  • Name of the algorithm.

  • Minimal size of a cell truncated by the geometry boundary. If the size of a truncated grid cell is Threshold times less than a initial cell size, then a mesh element is not created.

  • Implement Edges check-box activates incorporation of geometrical edges in the mesh.

    _images/cartesian_implement_edge.png

    Implement Edges switched off to the left and on to the right

  • Create Faces check-box activates creation on mesh faces.

  • Consider Shared and Internal Faces check-box activates treatment of faces shared by solids and internal. By default the algorithm considers only outer boundaries of the geometry.

  • Apply Threshold to Shared / Internal Faces check-box activates application of Threshold to cells cut by shared and internal faces, that can cause appearance of holes inside the mesh.

  • Definition mode allows choosing how Cartesian structured grid is defined. Location of nodes along each grid axis is defined individually:

    • You can specify the Coordinates of grid nodes. Insert button inserts a node at Step distance (negative or positive) from the selected node. Delete button removes the selected node. Double click on a coordinate in the list enables its edition. Note that node coordinates are measured along directions of axes that can differ from the directions of the Global Coordinate System.

    • You can define the Spacing of a grid as an algebraic formula f(t) where t is a position along a grid axis normalized at [0.0,1.0]. f(t) must be non-negative at 0. <= t <= 1. The whole extent of geometry can be divided into ranges with their own spacing formulas to apply; a t varies between 0.0 and 1.0 within each Range. Insert button divides a selected range into two. Delete button adds the selected sub-range to the previous one. Double click on a range in the list enables edition of its right boundary. Double click on a function in the list enables its edition.

Note

The actual grid spacing can slightly differ from the defined one. This is done for the best fitting of polyhedrons and for a better mesh quality on the interval boundaries. For example, if a constant Spacing is defined along an axis, the actual grid will fill the shape’s dimension L along this axis with round number of equal cells: Spacing_actual = L / round( L / Spacing_defined ).

  • Fixed Point group allows defining an exact location of a grid node in the direction defined by spacing. The following cases are possible:

    • If all three directions are defined by spacing, there will be a mesh node at the Fixed Point.

    • If two directions are defined by spacing, there will be at least a link between mesh nodes passing through the Fixed Point.

    • If only one direction is defined by spacing, there will be at least an element facet passing through the Fixed Point.

    • If no directions are defined by spacing, Fixed Point is disabled.

  • Directions of Axes group allows setting the directions of grid axes.

    • If Orthogonal Axes check-box is activated the axes remain orthogonal during their modification.

    • Selection buttons enable snapping corresponding axes to direction of a geometrical edge selected in the Object Browser. Edge direction is defined by coordinates of its end points.

    • Optimal Axes button runs an algorithm that tries to set the axes to maximize the number of generated hexahedra.

    • Reset button returns the axes in a default position parallel to the axes of the Global Coordinate System.

See Also a sample TUI Script of a Usage of Body Fitting algorithm.