Lamination-based efficient treatment of weak discontinuities for non-conforming finite element meshes

Jedrzej Dobrzanski,Kajetan Wojtacki,Stanislaw Stupkiewicz

DOI: https://doi.org/10.1016/j.compstruc.2023.107209

2024-06-05

Abstract:When modelling discontinuities (interfaces) using the finite element method, the standard approach is to use a conforming finite-element mesh in which the mesh matches the interfaces. However, this approach can prove cumbersome if the geometry is complex, in particular in 3D. In this work, we develop an efficient technique for a non-conforming finite-element treatment of weak discontinuities by using laminated microstructures. The approach is inspired by the so-called composite voxel technique that has been developed for FFT-based spectral solvers in computational homogenization. The idea behind the method is rather simple. Each finite element that is cut by an interface is treated as a simple laminate with the volume fraction of the phases and the lamination orientation determined in terms of the actual geometrical arrangement of the interface within the element. The approach is illustrated by several computational examples relevant to the micromechanics of heterogeneous materials. Elastic and elastic-plastic materials at small and finite strain are considered in the examples. The performance of the proposed method is compared to two alternative, simple methods showing that the new approach is in most cases superior to them while maintaining the simplicity.

Computational Engineering, Finance, and Science

What problem does this paper attempt to address?

This paper focuses on the treatment of weak discontinuities in the finite element method (FEM), especially in the case of complex geometric shapes, particularly in 3D. The traditional approach is to use a matched finite element mesh that conforms to the interface, but this can be difficult and time-consuming when dealing with complex geometries. The paper proposes a new non-matching finite element technique called the Laminated Element Technique (LET) to efficiently handle weak discontinuities by using laminated microstructures.

The inspiration for this method comes from the composite voxel technique in spectralsolvers based on the Fast Fourier Transform (FFT). Each finite element severed by the interface is treated as a simple laminated structure, where the phase volume fraction and lamination direction are determined by the actual interface's geometric arrangement within the element. This approach does not require additional global degrees of freedom and is implemented only at the element level, improving efficiency in handling complex geometries while maintaining simplicity.

The paper demonstrates the method of LET through several computational examples related to the micromechanics of multiphase materials, and compares its performance with two simple non-matching methods (element-level allocation and Gaussian point-level allocation). Although LET may not be as accurate as more complex methods such as Extended Finite Element Method (XFEM) in terms of accuracy, its advantage lies in its simplicity and outperforms the reference methods in most cases.

In addition, the paper discusses other methods for handling discontinuities such as XFEM, Immersed Interface Finite Element Method (IIFEM), and Adaptive Finite Element Method (AFEM), and points out the advantages and disadvantages of these methods. The main contribution of LET is to provide a simplified approach for handling weak discontinuities in non-matching grids, especially for 3D models, without significantly increasing computational costs.