4.7 Article

A hybrid H1 x H(curl) finite element formulation for a relaxed micromorphic continuum model of antiplane shear

Journal

COMPUTATIONAL MECHANICS
Volume 68, Issue 1, Pages 1-24

Publisher

SPRINGER
DOI: 10.1007/s00466-021-02002-8

Keywords

Relaxed micromorphic continuum; Edge elements; Nedelec elements; Curl based energy; Mixed formulation; Combined Hilbert spaces; Metamaterials

Funding

  1. Austrian Science Fund (FWF) [W 1245]
  2. DFG [440935806, Neff 902/10-1]

Ask authors/readers for more resources

The text discusses a simulation approach for metamaterials by extending the associated continuum theory concerning its kinematic equations using a relaxed micromorphic continuum model. By incorporating the Curl of the nonsymmetric microdistortion in the free energy function, solutions not belonging to H-1 are suggested, highlighting potential limitations of standard nodal H-1-finite elements. The use of base functions from Hilbert spaces H-1 and H(curl) is proposed as a method to address convergence issues in this class of problems, with a focus on a reduced two-dimensional relaxed micromorphic continuum for simplicity.
One approach for the simulation of metamaterials is to extend an associated continuum theory concerning its kinematic equations, and the relaxed micromorphic continuum represents such a model. It incorporates the Curl of the nonsymmetric microdistortion in the free energy function. This suggests the existence of solutions not belonging to H-1, such that standard nodal H-1-finite elements yield unsatisfactory convergence rates and might be incapable of finding the exact solution. Our approach is to use base functions stemming from both Hilbert spaces H-1 and H(curl), demonstrating the central role of such combinations for this class of problems. For simplicity, a reduced two-dimensional relaxed micromorphic continuum describing antiplane shear is introduced, preserving themain computational traits of the three-dimensional version. This model is then used for the formulation and a multi step investigation of a viable finite element solution, encompassing examinations of existence and uniqueness of both standard and mixed formulations and their respective convergence rates.

Authors

I am an author on this paper
Click your name to claim this paper and add it to your profile.

Reviews

Primary Rating

4.7
Not enough ratings

Secondary Ratings

Novelty
-
Significance
-
Scientific rigor
-
Rate this paper

Recommended

Article Engineering, Multidisciplinary

Avoiding membrane locking with Regge interpolation

Michael Neunteufel, Joachim Schoeberl

Summary: In this paper, a method to overcome membrane locking of thin shells is proposed by inserting an interpolation operator into the membrane energy term to weaken implicitly given kernel constraints. The approach decreases the number of constraints on triangular meshes compared to reduced integration techniques, and can be easily incorporated into existing shell elements. The performance of the method is demonstrated through several benchmark examples.

COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING (2021)

Article Computer Science, Interdisciplinary Applications

Fluid-structure interaction with H(div)-conforming finite elements

Michael Neunteufel, Joachim Schoberl

Summary: This study presents a novel application of the (high-order) H(div)-conforming Hybrid Discontinuous Galerkin finite element method for monolithic fluid-structure interaction (FSI), which yields exact divergence free fluid velocity solutions by introducing the Piola transformation. With the use of hp-refinement strategies, singularities and boundary layers are overcome leading to optimal spatial convergence rates. Copyright (C) 2020 Elsevier Ltd.

COMPUTERS & STRUCTURES (2021)

Article Engineering, Multidisciplinary

A finite-strain model for incomplete damage in elastoplastic materials

David Melching, Michael Neunteufel, Joachim Schoeberl, Ulisse Stefanelli

Summary: This study introduces a three-dimensional model capable of describing coupled damage and plastic effects in solids at finite strains, with rate-independent, associative, and unidirectional damage evolution features. Numerical simulations are conducted using multiphysics finite element software for various 2D and 3D settings under different choices of boundary conditions and possibly in presence of pre-damaged regions.

COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING (2021)

Article Engineering, Multidisciplinary

Three-field mixed finite element methods for nonlinear elasticity

Michael Neunteufel, Astrid S. Pechstein, Joachim Schoeberl

Summary: This paper extends the TDNNS method to nonlinear elasticity by lifting the distributional derivatives of the displacement vector to a regular strain tensor using the Hu-Washizu principle. Three different methods are introduced using either the deformation gradient, the Cauchy-Green strain tensor, or both as independent variables. The good performance and accuracy of the methods are demonstrated through numerical examples where stress and strain variables are locally eliminated in linear sub-problems, resulting in an equation system solely in terms of displacement variables.

COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING (2021)

Article Engineering, Multidisciplinary

Primal and mixed finite element formulations for the relaxed micromorphic model

Adam Sky, Michael Neunteufel, Ingo Muench, Joachim Schoeberl, Patrizio Neff

Summary: The classical continuum theory fails to capture the mechanical behavior of materials with pronounced microstructure, but the relaxed micromorphic continuum provides an alternative method. By enriching the kinematics of the mathematical model, this theory introduces a microdistortion field with nine extra degrees of freedom for each material point. We describe the construction of appropriate finite elements using Nedelec and Raviart-Thomas subspaces and explore the numerical behavior of the relaxed micromorphic model.

COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING (2022)

Article Engineering, Multidisciplinary

On [H1]3x3, [H(curl)]3 and H(sym Curl) finite elements for matrix-valued Curl problems

Adam Sky, Ingo Muench, Patrizio Neff

Summary: In this study, we test the numerical behavior of matrix-valued fields approximated by finite element subspaces of [H-1](3x3), [H(curl)](3), and H(sym Curl) for a linear abstract variational problem related to the relaxed micromorphic model. We introduce the formulation of the corresponding finite elements, perform numerical benchmarks, and present our conclusions. The introduction of the newly appropriate space H(sym Curl) reduces the continuity assumptions of the classical micromorphic model and provides a larger solution space for the microdistortion.

JOURNAL OF ENGINEERING MATHEMATICS (2022)

Article Engineering, Multidisciplinary

On pressure robustness and independent determination of displacement and pressure in incompressible linear elasticity

Adam Zdunek, Michael Neunteufel, Waldemar Rachowicz

Summary: We investigate the possibility to determine the divergence-free displacement independently from the pressure reaction for a class of boundary-value problems in incompressible linear elasticity. For convex domains, there is only one variational boundary value problem that allows the independent determination. The weakly divergence-free displacement can be computed pressure robustly by decomposing the total body force using the Helmholtz decomposition.

COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING (2023)

Article Engineering, Multidisciplinary

A Reissner-Mindlin plate formulation using symmetric Hu-Zhang elements via polytopal transformations

Adam Sky, Michael Neunteufel, Jack S. Hale, Andreas Zilian

Summary: In this work, new finite element discretisations of the shear-deformable Reissner-Mindlin plate problem are developed based on the Hellinger-Reissner principle. Conforming Hu-Zhang elements are used to discretise the bending moments and the rotation field, resulting in highly accurate approximations and satisfaction of the Kirchhoff-Love constraint. The formulation is extended using Raviart-Thomas elements for the shear stress to preserve optimal convergence rates for small thicknesses.

COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING (2023)

Article Computer Science, Interdisciplinary Applications

Numerical shape optimization of the Canham-Helfrich-Evans bending energy

Michael Neunteufel, Joachim Schoeberl, Kevin Sturm

Summary: In this paper, a novel numerical scheme is proposed for the Canham-Helfrich-Evans bending energy. It employs a three-field lifting procedure to convert the distributional shape operator to an auxiliary mean curvature field. By introducing the energetic conjugate scalar stress field as a Lagrange multiplier, the resulting fourth order problem is transformed into a mixed saddle point problem involving second order differential operators. The authors also derive analytical first variations and discuss the computation of shape derivatives in the finite element software NGSolve. Numerical simulations are provided to demonstrate the effectiveness of the proposed scheme and method.

JOURNAL OF COMPUTATIONAL PHYSICS (2023)

Article Computer Science, Theory & Methods

A flexible sparse matrix data format and parallel algorithms for the assembly of finite element matrices on shared memory systems

Adam Sky, Cesar Polindara, Ingo Muench, Carolin Birk

Summary: This paper focuses on the assembly process of the global stiffness matrix in finite element methods, exploring different algorithms and their efficiency on shared memory systems using C++. The use of atomic synchronization primitives for data-race free algorithms and data structures is a key aspect of the investigation. Additionally, a new flexible storage format for sparse matrices is proposed and its performance is compared with the compressed row storage format using abstract benchmarks based on common characteristics of finite element problems.

PARALLEL COMPUTING (2023)

Article Mathematics, Applied

Higher order Bernstein-Bézier and Nédélec finite elements for the relaxed micromorphic model

Adam Sky, Ingo Muench, Gianluca Rizzi, Patrizio Neff

Summary: This work applies the polytopal template methodology for the construction of Nedelec elements in the relaxed micromorphic model. Dual numbers are used in conjunction with Bernstein-Bezier polynomials to compute hp-FEM solutions of the model. The Nedelec base functions inherit the optimal complexity of the underlying Bernstein-Bezier basis through the polytopal template methodology.

JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS (2024)

Article Computer Science, Interdisciplinary Applications

Fully and semi-automated shape differentiation in NGSolve

Peter Gangl, Kevin Sturm, Michael Neunteufel, Joachim Schoeberl

Summary: This paper presents a framework for automated shape differentiation in the finite element software NGSolve, combining the Lagrangian approach and automated differentiation capabilities. Users can customize the level of automation needed, with automatic generation of first- and second-order shape derivatives discussed for both constrained and unconstrained model problems. Numerical experiments validate the accuracy of computed derivatives, and first- and second-order shape optimization algorithms are demonstrated for various numerical optimization examples.

STRUCTURAL AND MULTIDISCIPLINARY OPTIMIZATION (2021)

No Data Available