March 18th, 2019
approx reading time
In general, we can decompose errors in FEA—finite element analysis—in three main groups:
Furthermore, it must be kept in mind that the finite element method, in addition to the other numerical methods (FVM, FDM, BEM) consists of approximations.
The finite element description is a boundary value problem (BVP), which means there is a differential equation with a number of constraints called boundary conditions.
The solution to a BVP must resolve the differential equation and satisfy the boundary condition (which I will elaborate in more detail in another article).1
Errors of this type can include:
• the wrong geometric description: for instance, we use axial symmetry or rotational symmetry, but often forget that we are dealing with an antisymmetric load
• the wrong definition of the material: for example, the limit of Poisson’s ratio in isotropic materials
• the wrong definition of the load: trying to simplify complex load states or a number of loads with one load (depends on the case). But often “the difference between the effects of two different but statically equivalent loads becomes very small at sufficiently large distances from load” as Saint-Venant’s Principle expresses
• the wrong boundary conditions: if you forget to fix model (rigid body motion) or set boundary conditions, you might inadvertently falsify the results
• the wrong type of analysis
Concerning the analysis types, there are several options:
The purpose of a helmet is to protect the person who wears it from a head injury during impact. In this project, the impact of a human skull with and without a helmet was simulated with a nonlinear dynamic analysis. Download this case study for free.
A finite element model will sometimes contain a so-called singularity. This means there are points in your model where values tend toward an infinite value.4
If you are a newbie to finite element analysis, you might assume that a singularity is a term derived from a science-fiction movie like Star Trek.
The image below shows you what a singularity is:5
A singularity, as shown in the picture above, is basically the same in your FEA model. For example, in both cases, the infinite density and gravity are equivalent to an infinite stress at a sharp corner.
Singularities can be confusing because they cause an accuracy problem inside the model, which implies a problem with visualization because singularities extend the range of the stresses. This means that smaller stresses can appear to be negligible.
But what causes singularities? If you some online research, you will find a multitude of causes and explanations. However, summa summarum, we can confidently identify boundary conditions as the major cause.
Another predestinated problem implies that singularity problems are cracks because they can be seen as a 180° re-entrant corner. Luckily, the crack tip does not have to be studied every time, rather the focus can be placed on the stress intensity factor, using the J-Integral, or even looking at the energy dissipated during fracture (strain energy release rate).
Additionally, it is important to pay attention if you apply a force to a single node. This will give you infinite stresses! In the real world, forces are never applied on a single node due to the Saint-Venant’s principle, which tells us that “the difference between the effects of two different but statically equivalent loads becomes very small at sufficiently large distances from the load.” I also encountered this problem when I modeled a Hertzian contact in ANSYS and applied my force to a single node on the top of my cylinder. It took me two days to realize my mistake.
Second-order elements can assume either a concave or convex shape, meaning that they can display better results and deformation by easily mapping to curvilinear geometry. This is a big advantage when handling nonlinear material, i.e., elastoplastic or hyperelastic material. But everything comes at a price. The computational effort and time consumption are higher than using first-order tetrahedral elements.
h-method: reducing the size of the mesh
p-method: increase of the polynomial order in the element, which is good for regions with a low-stress gradient
r-method: relocates the position of a node
or combinations of the methods mentioned: hr, which is good for regions with large stress gradients, and hp or hpr
A nice mnemonic for these methods:
h-method = MesH
p-method = Polynomial order
r-method = Relocate node
General warnings to avoid errors in FEA:
Useful information: Exact integrations provide a structure that is too stiff!2
“Science walks forward on two feet, namely theory and experiment… but continuous progress is only made by the use of both.” – Robert A. Millikan
I hope this article gave you some useful tips on how to avoid errors in FEA and better understand singularities. If you’d like to learn more about the SimScale cloud-based simulation platform and its features, download this overview.
Alternatively, watch the recording of our “More Durable Consumer Products with FEA” webinar. All you need to do is fill out this short form and it will play automatically.
Subscribe to the SimScale Blog
Your hub for everything you need to know about
simulation and the world of CAE