Discretization Analysis in FEM Models

MATEC Web of Conferences, Jan 2016

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.

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:

https://www.matec-conferences.org/articles/matecconf/pdf/2016/16/matecconf_spbwosce2016_01063.pdf

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)


This is a preview of a remote PDF: https://www.matec-conferences.org/articles/matecconf/pdf/2016/16/matecconf_spbwosce2016_01063.pdf

Giulia Martello. Discretization Analysis in FEM Models, MATEC Web of Conferences, 2016, 53, DOI: 10.1051/matecconf/20165301063