Abstract:The first part of the cumulative thesis contains the numerical analysis of different $hp$-finite element discretizations related to two different weak formulations of a model problem in elastoplasticity with linearly kinematic hardening. Thereby, the weak formulation either takes the form of a variational inequality of the second kind, including a non-differentiable plasticity functional, or represents a mixed formulation, in which the non-smooth plasticity functional is resolved by a Lagrange multiplier. As the non-differentiability of the plasticity functional causes many difficulties in the numerical analysis and the computation of a discrete solution it seems advantageous to consider discretizations of the mixed formulation. In a first work, an a priori error analysis of an higher-order finite element discretization of the mixed formulation (explicitly including the discretization of the Lagrange multiplier) is presented. The relations between the three different $hp$-discretizations are studied in a second work where also a reliable a posteriori error estimator that also satisfies some (local) efficiency estimates is derived. In a third work, an efficient semi-smooth Newton solver is proposed, which is obtained by reformulating a discretization of the mixed formulation as a system of decoupled nonlinear equations. The second part of the thesis introduces a new $hp$-adaptive algorithm for solving variational equations, in which the automatic mesh refinement does not rely on the use of an a posteriori error estimator or smoothness indicators but is based on comparing locally predicted error reductions.

What problem does this paper attempt to address?

The main problem that this paper attempts to solve is how to solve boundary value problems efficiently and accurately by the hp - adaptive finite element method (hp - FEM) in elastoplastic material modeling. Specifically, the paper focuses on the following aspects: 1. **Discretization of Weak Forms**: For the linear kinematic hardening model in elastoplastic problems, two different discretization methods of weak forms are studied: - One is the second form of variational inequality, which contains a non - differentiable plastic function. - The other is the mixed form, in which the non - differentiable plastic function is resolved by Lagrange multipliers. 2. **Numerical Analysis and Error Estimation**: - Numerical analysis is carried out for different hp - finite element discretizations of the above two weak forms. - A reliable a posteriori error estimator is derived, which is applicable to any discretization conforming to the displacement field and plastic strain and satisfies certain local efficiency estimates. 3. **Efficient Solution Algorithms**: - An efficient semi - smooth Newton solver is proposed. By reformulating the discretization of the mixed form and converting it into a set of decoupled nonlinear equations, superlinear convergence is achieved. 4. **hp - Adaptive Strategies**: - A new hp - adaptive algorithm based on local error reduction prediction is introduced. This algorithm does not rely on traditional a posteriori error estimates or smoothness indicators but performs automatic mesh refinement by comparing locally predicted error reductions. ### Specific Contributions - **Part I: hp - Finite Element Method for Elastoplasticity** - The application of different hp - finite element discretization methods in elastoplastic problems is studied, including the discretizations of variational inequalities and mixed forms. - A reliable a posteriori error estimator is derived and its local efficiency is proven. - An efficient semi - smooth Newton solver is proposed and its robustness and effectiveness in different situations are demonstrated. - **Part II: hp - Adaptivity Based on Local Error Reductions** - A new hp - adaptive strategy is introduced, which is based on local error reduction prediction rather than relying on traditional a posteriori error estimates or smoothness indicators. - Error reduction is predicted by solving small - scale linear problems with lower computational cost, thereby guiding adaptive mesh refinement. - This concept is proposed in the framework of abstract Hilbert space and applied to hp - finite element discretizations, demonstrating its effectiveness and applicability. ### Summary The core objective of this paper is to improve the accuracy and efficiency of solving boundary value problems in elastoplastic problems by introducing and analyzing hp - finite element discretization methods and their adaptive strategies. In particular, through improved numerical analysis, efficient solution algorithms, and innovative adaptive strategies, difficulties in traditional methods are solved, such as numerical challenges brought by non - differentiable plastic functions and the efficient approximation of solutions with low regularity.

$hp$-FEM for Elastoplasticity & $hp$-Adaptivity Based on Local Error Reductions

A Posteriori Error Estimates for $hp$-FE Discretizations in Elastoplasticity

An adaptive non-symmetric finite volume and boundary element coupling method for a fluid mechanics interface problem

An $hp$-adaptive strategy based on locally predicted error reductions

Residual-Based a posteriori error estimation for hp-adaptive finite element methods for the Stokes equations

Mixed Finite Elements of Higher-Order in Elastoplasticity

hp adaptive Discontinuous Galerkin strategies driven by a posteriori error estimation with application to aeronautical flow problems

Adaptive Refinement for $hp$-Version Trefftz Discontinuous Galerkin Methods for the Homogeneous Helmholtz Problem

An implementation of hp-FEM for the fractional Laplacian

Convergence and Optimality of hp-AFEM

An $hp$-adaptive discontinuous Galerkin discretization of a static anti-plane shear crack model

$hp$-error analysis of mixed-order hybrid high-order methods for elliptic problems on simplicial meshes

A sparse hierarchical $hp$-finite element method on disks and annuli

Unifying a posteriori error analysis of five piecewise quadratic discretisations for the biharmonic equation

Implementation of $hp$-adaptive discontinuous finite element methods in Dune-Fem

An arbitrarily high order unfitted finite element method for elliptic interface problems with automatic mesh generation, Part II. Piecewise-smooth interfaces

Unified $hp$-HDG Frameworks for Friedrichs' PDE systems

Mixed finite element methods for nonlinear equations : a priori and a posteriori error estimates

Convergence of Adaptive Mixed Finite Element Methods for the Hodge Laplacian Equation: without harmonic forms

Matrix-Free Higher-Order Finite Element Methods for Hyperelasticity

Adaptive hp-Finite Element Method for Electromagnetic Field Logging Problems