Details
| Original language | English |
|---|---|
| Pages (from-to) | 31-40 |
| Number of pages | 10 |
| Journal | Finite Elements in Analysis and Design |
| Volume | 165 |
| Early online date | 20 Aug 2019 |
| Publication status | Published - 1 Nov 2019 |
Abstract
For the accurate imposition of surface loads using the Finite Element Method normally the load is discretized at the surface facets. Therefore appropriate surface shape functions are needed. If the surface contains kinks, which occurs often in contact cases, the imposition is more complicated. In order to simplify the imposition of surface loads, an alternative approach is presented. This formulation is purely based on volume contributions of the discretized elements of the body. The surface nodes are automatically identified and the linearization is straightforward. No specific surface information is necessary. The imposition is simply based on the application of the divergence theorem. The only prerequisite is the fulfillment of the integration constraint, a necessary requirement for Galerkin solution schemes. With this, boundary nodes are directly identified by possessing a non zero normal vector whereas for inner nodes this vector is identically zero. Moreover, the normal vectors at the boundary nodes contain all the information of the surface and correspond to resultant nodal normal vectors. This is especially advantageous in the case of surfaces that contain kinks. It also simplifies the search algorithm in computational contact mechanics to find the closest distance to other bodies. This approach also works for any kind of shape functions. The nodal force vectors are always determined accurately. The advantages and the simple handling of this approach are demonstrated by means of several examples including follower loads and contact cases. Additionally, the influence of the isoparametric concept on the integration constraint is investigated by evaluating the behavior of different shape functions on an irregular grid. Although this new approach is only demonstrated within the context of the Finite Element Method, due to its generic derivation it can be applied to any Galerkin solution scheme which fulfills the integration constraint.
Keywords
- Divergence theorem, Finite Element Method, Integration constraint, Surface loads
ASJC Scopus subject areas
- Mathematics(all)
- Analysis
- Engineering(all)
- General Engineering
- Computer Science(all)
- Computer Graphics and Computer-Aided Design
- Mathematics(all)
- Applied Mathematics
Cite this
- Standard
- Harvard
- Apa
- Vancouver
- BibTeX
- RIS
In: Finite Elements in Analysis and Design, Vol. 165, 01.11.2019, p. 31-40.
Research output: Contribution to journal › Article › Research › peer review
}
TY - JOUR
T1 - Direct nodal imposition of surface loads using the divergence theorem
AU - Weißenfels, Christian
PY - 2019/11/1
Y1 - 2019/11/1
N2 - For the accurate imposition of surface loads using the Finite Element Method normally the load is discretized at the surface facets. Therefore appropriate surface shape functions are needed. If the surface contains kinks, which occurs often in contact cases, the imposition is more complicated. In order to simplify the imposition of surface loads, an alternative approach is presented. This formulation is purely based on volume contributions of the discretized elements of the body. The surface nodes are automatically identified and the linearization is straightforward. No specific surface information is necessary. The imposition is simply based on the application of the divergence theorem. The only prerequisite is the fulfillment of the integration constraint, a necessary requirement for Galerkin solution schemes. With this, boundary nodes are directly identified by possessing a non zero normal vector whereas for inner nodes this vector is identically zero. Moreover, the normal vectors at the boundary nodes contain all the information of the surface and correspond to resultant nodal normal vectors. This is especially advantageous in the case of surfaces that contain kinks. It also simplifies the search algorithm in computational contact mechanics to find the closest distance to other bodies. This approach also works for any kind of shape functions. The nodal force vectors are always determined accurately. The advantages and the simple handling of this approach are demonstrated by means of several examples including follower loads and contact cases. Additionally, the influence of the isoparametric concept on the integration constraint is investigated by evaluating the behavior of different shape functions on an irregular grid. Although this new approach is only demonstrated within the context of the Finite Element Method, due to its generic derivation it can be applied to any Galerkin solution scheme which fulfills the integration constraint.
AB - For the accurate imposition of surface loads using the Finite Element Method normally the load is discretized at the surface facets. Therefore appropriate surface shape functions are needed. If the surface contains kinks, which occurs often in contact cases, the imposition is more complicated. In order to simplify the imposition of surface loads, an alternative approach is presented. This formulation is purely based on volume contributions of the discretized elements of the body. The surface nodes are automatically identified and the linearization is straightforward. No specific surface information is necessary. The imposition is simply based on the application of the divergence theorem. The only prerequisite is the fulfillment of the integration constraint, a necessary requirement for Galerkin solution schemes. With this, boundary nodes are directly identified by possessing a non zero normal vector whereas for inner nodes this vector is identically zero. Moreover, the normal vectors at the boundary nodes contain all the information of the surface and correspond to resultant nodal normal vectors. This is especially advantageous in the case of surfaces that contain kinks. It also simplifies the search algorithm in computational contact mechanics to find the closest distance to other bodies. This approach also works for any kind of shape functions. The nodal force vectors are always determined accurately. The advantages and the simple handling of this approach are demonstrated by means of several examples including follower loads and contact cases. Additionally, the influence of the isoparametric concept on the integration constraint is investigated by evaluating the behavior of different shape functions on an irregular grid. Although this new approach is only demonstrated within the context of the Finite Element Method, due to its generic derivation it can be applied to any Galerkin solution scheme which fulfills the integration constraint.
KW - Divergence theorem
KW - Finite Element Method
KW - Integration constraint
KW - Surface loads
UR - http://www.scopus.com/inward/record.url?scp=85070802454&partnerID=8YFLogxK
U2 - 10.1016/j.finel.2019.06.009
DO - 10.1016/j.finel.2019.06.009
M3 - Article
AN - SCOPUS:85070802454
VL - 165
SP - 31
EP - 40
JO - Finite Elements in Analysis and Design
JF - Finite Elements in Analysis and Design
SN - 0168-874X
ER -