Discretization Analysis in FEM Models
MATEC Web of Conferences
Discretization Analysis in FEM Models
Giulia Martello 0
0 Università degli Studi di Brescia, Mechanical and Industrial Engineering Department , Brescia , Italy
Discretization allows to switch from an infinite number of elements in a finite number and thanks to that the infinite degrees of freedom of the system become a finite number. The discretization of the domain consists of a subdivision into subdomains, which form a lattice or mesh. This work deals on different types of mesh and propose a choice for a specific problem.
1 Introduction
The discretisation [1-4] of the solution domain produces a numerical description of the computational
domain [5-9]. The space is divided into a finite number of discrete regions or cells [10, 11], which
form a mesh [12, 13].
This allows to switch from an infinite number of elements in a finite number and thanks to that the
infinite degrees of freedom of the system become a finite number [14].
The discretization with Finite Element Method (FEM) [15-19] is based on the principle which is
called with the latin expression “divide et impera” (thus “divide and rule”). This means that a complex
problem is decomposed into many simpler problems, whose solutions lead to the solution of the
general problem [20-23]. In this case the model of difficult solution is decomposed into many
elements that can be solved more easily.
Alongside spatial discretization, also temporal discretization is, for transient simulations, fundamental
as well, thus the time interval is also split into a finite number of time-steps [24-33]. This leads to less
computational efforts and an analysis of smaller deformation of the general problem.
The mesh generation is a fundamental step in the simulation because it is an integral part of the
numerical solution and must satisfy certain criteria to ensure a valid, and hence accurate, solution
[3439]. The mesh size [40-42] (defined for each direction and identified by the characteristic length h)
has to be optimized [43-48], indeed it should be fine enough to make the problem easily resolvable
and at the same time it has not to be too fine because the higher number of cells, the higher
computational costs. Therefore the mesh density [49] depends also on the machine which is used for
computations. After an accurate observation of the model shape and of gradients involved, a proper
mesh can be obtained by refining it in the regions with physical needs, for example close to
geometrical details that imply high gradients of the variable evaluated.
In the beam considered the variable which has a fundamental role is the stress component S12,
parallel to the y axis.
2 Types of meshing
Typically, cells used in FEM are simple geometric figures (in the one-dimensional problem segments,
triangles or quadrilaterals in 2D problems, tetrahedrons or hexagonal in 3D problems). The topology
is defined by special points called nodes.
The main types of cells that one can choose in a three-dimensional model follow.
• HEX: this type of elements are polyhedrons with 6 faces and quadrilateral faces, therefore these
elements have 8 knots. An example of a mesh constructed by using these elements is shown in Figure
2.
• HEX-dominated: HEX elements dominate, nevertheless it is also allowed to use some triangular
base prisms (wedge) in the transition regions. An example in Figure 2.4.
• TET: composition of elements only tetrahedral, thus with 4 knots, as shown in figure 4
• WEDGE: composition of wedge-shaped elements with 6 knots. More precisely, it refers to prisms
with a triangular base, as shown in figure 5.
3 Meshing effects and selection
Given these four possible forms of item, one must choose the most appropriate geometry treated. It
can be noted that HEX-dominated and HEX types differ only if the geometry is not completely
regular; for example a triangular base prisms may appear if one uses this type of elements for a square
plate with smoothed corner or with a hole in the middle. Since the geometry considered is extremely
regular, this option is not taken into account. Also the tetragonal form of elements can be discarded,
since these elements can be used to simplify a complex geometry but if the starting shape is already
regular the use of these elements can lead to a more complicated problem. The more plausible
possibilities are therefore HEX and WEDGE.
Concerning elements with triangular base, the concept explained regarding case tetrahedral is still
true, even if to a lesser degree. Indeed, also WEDGE elements introduce some irregularity in a
regular model.
HEX mesh, or in cubes or generic parallelepipeds, is the simplest and most intuitive mesh which
perfectly reflects the geometry of the model. Nevertheless before approving it permanently also
boundary conditions and loads must be taken into account, since they could make irregular stresses
even if within a regular geometry.
As seen, the beam is analyzed by using a fixed joint at one end (...truncated)