A hybrid H1×H(curl)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ H ^1\times H (\mathrm {curl})$$\end{document} finite element formulation for a relaxed micromorphic continuum model of antiplane shear

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 H1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ H ^1$$\end{document}, such that standard nodal H1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ H ^1$$\end{document}-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 H1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ H ^1$$\end{document} and H(curl)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ H (\mathrm {curl})$$\end{document}, 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 the main 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.


Introduction
Materials with a pronounced microstructure such as metamaterials, see e.g. [2,3,7,19], porous media, composites etc., activate micro-motions which are not accounted for in classical continuum mechanics, where each material point is equipped with only three translational degrees of freedom. Therefore, several approaches to model such materials can be found in literature, such as multi-scale finite element methods [1,10,11] or generalized continuum theories. The latter can be classified into higher gradient theories [5,17,23,32] and so called micromorphic continuum theories [30,42]. These theories extend the kinematics of the material point. Depending on the extension one obtains for example micropolar [16,25,26], microstretch [38] or microstrain [13,15] theories. In its most general setting, as introduced by Eringen und Mindlin [12,22], a micromorphic continuum theory allows the material point to undergo an affine distortion independent of its macroscopic deformation arising from the displacement field. Consequently, in the micromorphic theory a material point is considered with 3 + 9 = 12 degrees of freedom, of which the microdistortion P encompasses 9. The various micromorphic theories differ in their proposition of the free energy functional. While classical theories incorporate the full gradient of the microdistortion ∇ P into the energy function [31], the relaxed micromorphic theory [20,31,33,34] considers only Curl P. The incorporation of the Curl of the microdistortion, formally known as the dislocation density, into the free energy functional relaxes the continuity assumptions on the microdistortion and enlarges the space of possible weak solutions, i.e. [H 1 ] 3 × H(Curl). Furthermore, the relaxed micromorphic theory aspires to capture the entire spectrum of mechanical behaviour between the macro and micro scale of the material. This is achieved via homogenization of the material parameters and the introduction of the characteristic length L c [19,29], which determines the influence of the dislocation density in the free energy functional. Specific analytical solutions to the full isotropic relaxed micromorphic model are presented in [35][36][37].
For non-trivial boundary value problems, solutions of continuum theories are approximated via the finite element method. While the standard Lagrange elements are well suited for solutions in H 1 , solutions in H(curl) may require a different class of elements, depending on the problem at hand. The lowest class of finite elements in H(curl), sometimes called edge elements, have been derived by Nédélec [27,28]. Extensions to higher order element formulations can be found in [8,9,41,44]. In this paper we consider finite element formulations employing either H 1 ×[H 1 ] 2 or H 1 × H(curl) and investigate their validity in correctly approximating results in the relaxed micromorphic continuum. Furthermore, we test both a primal and mixed formulation of the corresponding boundary problem for increasingly large values of the characteristic length L c . To that end, we consider a planar version of the relaxed micromorphic continuum, namely of antiplane shear [43]. More precisely, the matrix-Curl in 3D reduces to a scalar-curl of the microdistortion in 2D. However, the results of our investigation directly apply to the full three-dimensional version.
The paper is organized as follows: In the following section we introduce the planar relaxed micromorphic continuum. Section 3 is devoted to prove solvability of the primal and mixed problem and discussing properties in the limit case L c → ∞, in both the continuous and discrete settings, respectively. In Section 4 we present appropriate base functions for H(curl), the corresponding covariant Piola transformation for Nédélec finite elements and the resulting stiffness matrices. Finally, we present several numerical examples to confirm the theoretical results.
From now on, we consider the planar reduction of this continuum to antiplane shear, still capturing the main mathematical aspects of the three-dimensional version, namely the additional microdistortion and the curl where we employ the two-dimensional definitions of the curl and gradient operators In Eq. (2.3) we reduced the displacement to a scalar field u : ⊂ R 2 → R and the microdistortion P to a vector field ζ : ⊂ R 2 → R 2 . The displacement field u is now perpendicular to the plane of the domain. The elasticity tensors C e and C micro are replaced by the scalars μ e , μ micro > 0 and C c no longer appears.

Remark 2.1
The simplification of the model to antiplane shear serves to facilitate the mathematical analysis of the model and allows for a thorough investigation of the numerical behaviour of finite element solutions in the relaxed micromorphic theory. Whether this reduced model can be applied to real-world metamaterials is unclear at this time. For applications of the full three-dimensional theory see [19,20].
In order to find functions minimizing the potential energy I we calculate the variations with respect to u and ζ Partial integration of Eq. (2.5a) and Eq. (2.5b) yields the strong form including boundary conditions (see Appendix A for more details) where τ and ν denote the outer tangent and normal vector on the boundary, see Fig. 1, and with u and ζ the displacement and microdistortion fields on u D and ζ D are prescribed. From a mathematical point of view, it is possible to prescribe the tangential components of the microdistortion ζ on the boundary ζ D . This is used to test our numerical formulation in Sect. 5. However, from the point of view of physics it is impossible to control the microdistortion of the continuum with no direct relation to the displacement u and as such, the consistent coupling condition ζ , τ = ∇ u, τ arises on the Dirichlet boundary, being common to both u and ζ , enforcing the condition ζ D ⊂ u D . Furthermore, Dirichlet boundary data for the microdistortion ζ are not required for the existence of a unique solution here, as coercivity in the appropriate spaces is still determined.

Continuous case
In this section we prove the existence and uniqueness of the weak form of the planar relaxed micromorphic continuum. Further, the corresponding mixed formulation is presented, whose coercivity constant is independent of L c . Finally, we study necessary and sufficient conditions such that ∇u = ζ is guaranteed in the limit L c → ∞. For simplicity, we assume homogeneous Dirichlet conditions on the entire boundary throughout this section, i.e., u = 0 and ζ , τ = 0 on u D = ζ D = ∂ , and mention that the proof can be readily adapted for inhomogeneous and mixed boundary conditions as long as the Dirichlet boundary for the displacements is non-trivial, | u D | > 0, [14]. We define the following Hilbert spaces and their respective norms which are based on the Lebesgue norm and space Further, we use the product space Proof Using Cauchy-Schwarz and triangle inequality yields the continuity of a(·, ·) for all {u, ζ }, {δu, δζ } ∈ X with the constant c 1 = max 2μ e , 2μ micro , μ macro L 2 c . By employing Young's 2 and Poincaré-Friedrich's 3 inequalities we show the bilinear form to be coercive when the constant ε is chosen as 1 > ε > μ e μ e + μ micro , which is possible for μ e , μ micro > 0. Consequently, the coercivity constant reads β = c 3 2 min 1, 1 This finishes the proof. 2 as a(·, ·) is then no longer coercive in this space because one cannot find a constant c > 0 such that ζ 2 2 , however, at the cost of sub-optimal convergence rates in the discretized setting. We present numerical examples, where the exact solution is in H(curl, ) but not in [H 1 ( )] 2 observing only slow convergence. If the exact solution is smooth, i.e. ζ is also in [H 1 ( )] 2 , optimal convergence is observed.

Remark 3.1 Note, that the proof fails when taking instead
An important aspect of the relaxed micromorphic continuum is its relation to the classical continuum theory (linear elasticity). This relation is governed by the material constants, where the characteristic length L c plays a significant role. We are therefore interested in robust computations with respect to L c .
The following result characterizes the conditions when a trivial solution with respect to L c is expected.  Considering the limit case L c = 0, the continuity of the bilinear form a(·, ·) follows automatically from Eq. (3.5). However, for coercivity to hold, the space for ζ must be changed to [L 2 ( )] 2 , i.e., the regularity of ζ is lost.
Proof The proof of existence and uniqueness follows exactly the same lines as the proof of Theorem 3.1. If ω = ∇r we can conclude as in the proof of Theorem 3.2 that ζ is a gradient field. Remark 3.2 Using Theorem 3.2 and assuming ω = 0, we can reformulate Eq. (2.6b) to retrieve ζ from the known field u ζ = ∇χ = μ e μ micro + μ e ∇u . (3.9) Furthermore, we can condensate Eq. (2.6a) into the Poisson equation where the homogenization of the material constants follows as in [29]. We notice, that Theorem 3.2 and Theorem 3.3 imply the field u is always independent of the microdistortion ζ in this setting. In the condensed state, the relation of the model with antiplane shear for membranes is apparent.

Remark 3.3
We note that the previous result does not hold in the full three-dimensional relaxed micromorphic continuum, i.e. the absence of external moments does not automatically imply P = ∇χ for χ ∈ [H 1 ( )] 3 .
Having considered the limit of the characteristic length L c → 0, we reformulate Problem 3.4 as an equivalent mixed formulation in order to examine its limit for L c → ∞. We start by introducing the new variable and constructing a new bilinear form by multiplying it with a test function The restriction to m ∈ L 2 0 ( ) follows from the Stoke's theorem We introduce the (bi-)linear forms and the resulting mixed formulation reads: where the Lagrange multiplier m has the physical meaning of a moment stress tensor. The limit case lim L c → ∞ of Eq. (3.15) is well-defined, resulting in the problem: (3.16b) Consequently, at the limit lim L c → ∞ we have curl 2D ζ = 0. We now show existence and uniqueness of both mixed problems and that in the limit case L c → ∞ the solution of Eq. (3.15) converges to the solution of Eq. (3.16) with quadratic convergence rate in L c .
where c 2 does not depend on L c .

Remark 3.5
In the full micromorphic continuum, where the gradient takes the place of the curl of the microdistortion existence and uniqueness follow similarly with the space However, the limit case L c → ∞ yields ∇ζ = 0 and consequently ζ = const., for which nontrivial boundary conditions cannot be considered, compare also Section 5.7 for a numerical example.
To conclude this section we investigate the necessary and sufficient conditions such that in the limit L c → ∞ the solution satisfies ∇u = ζ . This state represents a zoom into the microstructure in the three-dimensional theory with microscopic stiffness given by μ micro [29]. In Theorem 3.2 we found sufficient conditions to obtain a gradient field for the microdistortion, which, however, does not have to be ∇u.
The following theorem states that only for a zero right-hand side f , but arbitrary ω, the desired behaviour is achieved.
if and only if f = 0, where c does not depend on L c .
Remark 3. 6 We can weaken the assumptions of Theorem 3.5 to u D = ζ D = ∅. Further, also non-homogeneous Dirichlet data can be considered, provided the consistent coupling condition ζ , τ = ∇ u, τ on ζ D holds. From the proof of Theorem 3.5 we obtain from the existence of a potential such that ζ = ∇ . Thus, ζ is expected to be in H(curl, ) as in general ∇ / ∈ [H 1 ( )] 2 for ∈ H 1 ( ).

Discrete case
Motivated by the de' Rham complex (see Fig. 2) we formulate a finite element combining base functions from both H 1 ( ) and H(curl, ) (and L 2 ( ) for the mixed formulation) setting Throughout this work we will use meshes consisting of quadrilaterals. On each element we denote the set of quadrilateral polynomials by Q n,m = span{x k y j | 0 ≤ k ≤ n, 0 ≤ j ≤ m}, compare also Eq. (4.5), and further the set of Nédélec ansatz functions by We start with the Lax-Milgram setting by defining X h = V h × U h . We note that solvability of the discretized problem follows directly from the continuous one as X h ⊂ X. Using Cea's lemma for the quasi-best approximation we can generate convergence estimates a priori.
(3.28) Fig. 2 The de' Rham complex in two dimensions depicting Hilbert spaces and approximation spaces connected by differential and interpolation operators. The kernel of one differential operator is exactly the range of the previous differential operator on its space and the differential and projection operators commute 123 Proof By inserting the interpolation operators associated through the commuting diagram we find where | · | H k denotes the standard Sobolev semi-norm.
Note that the constant c in Eq. (3.28) depends on L c . One may prove robust estimates in this setting. We, however, test for robustness with respect to L c in the context of mixed methods and use the equivalence of both.
In general the solvability of the discretized mixed problem does not follow from the continuous one. However, thanks to the commuting property of the de' Rham complex, the discrete kernel coercivity and the LBB condition follow immediately. Thus, we obtain the quasi-best approximation error where c is independent of L c . Inequality Eq. (3.32) states that, as long as the discretization error is not reached, we have quadratic convergence to the limit case lim L c → ∞. Due to the equivalence of the primal formulation Problem 3.4 and the mixed Eq. (3.15) we can deduce that the solution of Problem 3.4 is also robust with respect to L c . As we will see in the numerical examples, the mixed formulation is better suited for extremely large values of L c due to rounding errors.

Appropriate base functions
In the following we demonstrate the construction of the hybrid element in the linear case. The finite elements for the mixed formulation are employed directly using the open source finite element library NETGEN/NGSolve 5 [39,40].
For the mapping of x and y, see Fig. 3, we make use of linear quadrilateral Lagrange nodal base functions Hy e , where n is the number of finite elements in the mesh. As shown in Fig. 3, the elements are mapped via However, for ζ we make use of linear Nédélec base functions of the first type for quadrilaterals [4,21,27,44]. These functions are built around approximations of the curl operator. The corresponding spaces are those of quadrilateral polynomials (4.5) The weak form of the curl in the 2D space is formulated via Greens' formula 6 Therefore, the curl in is fully determined by its interface and inner rotation field. Consequently, we can decompose the two terms, such that the elements' dofs determine the interpolated field completely. This can be confirmed by setting all dofs to zero, checking for a vanishing field. The corresponding dofs and degrees of the polynomial spaces have been defined by Nédélec [27]. The element's boundary has been decomposed as ∂ = 1 ∪ 2 ∪ 3 ∪ 4 . The dofs read 4k edge dofs: Applying the dofs along all edges with the variable basis we find our base functions (4.10) The factor 1 /2 is chosen instead of the resulting 1 /4 as to simplify prescription on the Dirichlet boundary. The functions are depicted in Fig. 4. For the mixed formulation involving m ∈ L 2 ( ) the corresponding finite element space is given by piece-wise constants, N 0 (ξ, η) = 1. To enforce zero mean value, i.e. m ∈ L 2 0 ( ), a Lagrange multiplier λ ∈ R has to be used, leading to one additional equation in the final system.

123
Using higher polynomial orders, we can achieve faster convergence rates and better approximations. The finite element software NGSolve offers the use of hierarchical high order base functions for H 1 , H(curl), and L 2 spaces [44]. We employ NGSolve in our investigation of the mixed formulation with higher order base functions.

Covariant Piola transformation
In the previous section we formulated our base functions for the curl in the parametric space. In order to preserve the properties of the base function ϑ j acting on the curve's tangents ς (see Fig. 3 where θ j is the base function in the physical space and ∂ e = 1 ∪ 2 ∪ 3 ∪ 4 , the so called covariant Piola transformation is required [24]. The transformation is achieved by considering the push forward of the boundaries' normal vectors where J is the Jacobi matrix of the element mappings. In two dimensions the normal vectors on the element boundary and ν are the 90 • rotation of the tangent vectors given by finally yielding the definition of a transformation preserving integration along the tangent The transformation in Eq. (4.15) alone cannot guarantee the aligned orientation of base functions on the edges of neighbouring elements [44]. In order to achieve conformity we introduce a topological correction function ψ j based on the global orientation of edges given by node collections as demonstrated in Fig. 5. The drawings in Fig. 5 show the different roles of the mapping functions: 1. The covariant Piola transformation scales the projection onto the edge tangent.
2. The topological correction function sets a consistent orientation.
Thus, the final form of our edge base functions reads (4.16) Using Eq. (4.16) for the approximation of the microdistortion ζ yields For vectors undergoing a covariant Piola transformation, the transformation of the curl operator simplifies to

Element stiffness matrices
where K e , K micro and K macro are the element stiffness matrices employing the base function matrices H and according to Eq. (4.17) and Eq. (4.2), respectively  In order to compare our formulation, we also derive a nodal H 1 -finite element In contrast to the hybrid element, the approach in Eq. (4.23) requires 8 dofs per element for the microdistortion. Using Eq. (4.23) we obtain the following stiffness matrices for the nodal element In conclusion, we compare the hybrid element having 8 degrees of freedom in total with the nodal element having 12 degrees of freedom. The difference in the overall degrees of freedom results from the vectorial approach to the microdistortion in the hybrid element.

Numerical examples
In following examples we construct analytical solutions by imposing predefined displacement and microdistortion fields and calculating the resulting right-hand side.

Benchmark for an imposed vanishing microdistortion
We impose the predefined fields In order to constrain the numerical solution to that of our proposed fields in Eq. (5.1), we set the following Dirichlet boundary conditions In the following example we set for simplicity and extract the resulting force and moment (the right-hand side) with irregular meshes under h-refinement, as shown in Fig. 7. Both element formulations converge towards the analytical solution, see Fig. 6. The microdistortion field ζ displayed in Fig. 8 approaches zero with each refinement, satisfying the imposed field. We notice faster convergence in the hybrid element.

Benchmark for a non-vanishing imposed microdistortion
In the following step in our investigation we test our finite element formulations for a non-vanishing microdistortion field ζ , specifically a rotation field, as to determine the convergence behaviour of the nodal element with respect to the curl stiffness. We set = [−4, 4] × [ −4, 4], μ e = μ macro = μ micro = L c = 1 and the fields x( with the corresponding Dirichlet boundary conditions The following force and moment are extracted, for details see Appendix B, Consequently, the curl term is neither explicitly nor implicitly omitted. We compare the displacement u and the error ζ −ζ L 2 for both element formulations on an irregular mesh undergoing refinement, see Figs. 9, 10 and 11.
As shown in Fig. 10, both elements converge towards the analytical solution. However, we notice differences in the convergence rates, namely the nodal element converges faster in ζ .

Solutions in H(curl)
where ζ follows from Eqs. (2.6a) and (2.6b). Note that the boundary data of ζ jumps and is therefore not in We test both elements on an irregular mesh undergoing refinement Fig. 12. We note the hybrid element finds the exact solution immediately with a coarse mesh, whereas the nodal element requires a much higher level of refinement in order to deliver a viable approximation. The nodal element localizes the error due to the discontinuity further with each refinement as seen in Fig. 13. The convergence graph in Fig. 14 depicts the slow sub-optimal convergence of the nodal element, compare Eq. (3.28). Note, the error in the hybrid element for the same meshes is always at a factor 10 −15 for both u and ζ . Due to the higher continuity conditions of the nodal element, it could never find the analytical solution, but would converge further towards it with each refinement.
We present a second example allowing us to compare the convergence rates for both formulations. Let    H(curl, ), not in [H 1 ( )] 2 . We use structured quadrilateral meshes (see Fig. 16) resolving the interface at x = 0.5, where the normal component of the exact solution of ζ jumps, with linear, quadratic and cubic polynomials for the nodal elements. We observe that higher polynomial degrees do not increase the convergence rate and only sub-optimal root-convergence is achieved (see Fig. 15). For linear and quadratic ansatz functions in the primal H(curl) method we observe optimal convergence rates.

Convergence for L c → 0
As mentioned in Section 3, the characteristic length L c represents an important term in the relaxed micromorphic theory. This scalar governs the relation of the relaxed micromorphic continuum to the standard Cauchy continuum. In the previous examples we have been able to generate stable results for the case L c = 1. In this example we consider the limit L c → 0, which can be interpreted as a highly homogenous material. In the L c = 0 setting, the relaxed We use u to recover the analytical solution for ζ ζ = μ e μ micro + μ e ∇ u = −2 cos(x) sin(x) −2 cos(y) sin(y) , (5.14) and the resulting right-hand side

Robustness in L c
The upper limit of the characteristic length L c is defined to be infinity.
from which we can extract the resulting force fields according to Eqs. (2.6a) and (2.6b). We test for convergence using linear elements.
As expected from the theory, we observe uniform convergence up to the point where rounding errors occur in the primal method for very large L c terms. The convergences of the mixed formulation remains stable for all values of L c as it is not affected by rounding errors, cf. Fig. 19. Using lowest order linear nodal elements for ζ leads to non-robust behaviour in L c in terms of immense locking. Considering quadratic Lagrange elements overcomes this locking phenomena, however, at the cost of more dofs. To test the convergence depending on L c , Eq. (3.18), for the case L c → ∞ we use quadratic elements -i.e., quadratic H 1 and Nédélec elements, and linear L 2 elements for m in the mixed formulation -in NGSolve and four different structured grids. The same domain as in the previous example is considered and for the limit solution Eq. (5.16) is used, with L c → ∞ in Eq. (5.16b). Again, the primal methods suffers for large values of L c from rounding errors, whereas for the mixed method we observe the expected quadratic convergence rate up to the discretization error, compare Eq.

Convergence for L c → ∞
We prove the theoretical result of Theorem 3.5, with the same domain, boundary conditions, and material constants as in the previous example, by setting the external force and moments , ω = ∇r + D curl ( ) , (5.17) and testing for convergence ∇u − ζ H(curl) = O(L −2 c ) for L c → ∞ using NGSolve with linear base functions.
The results are computed using the primal method. By staying within the rounding precision bounds retrieved from our investigation of the robustness in L c , we are able to find results converging quadratically to the previously derived expectations, see Fig. 21.

The consistent coupling condition
We conclude our investigation by considering the consistent coupling condition on both the full and relaxed micromorphic continuum models using NGSolve with the primal method. We set the domain = [−4, 4] × [−4, 4] with the material parameters μ e , μ micro , μ macro = 1, the boundary conditions (5.18) and the external forces and test for convergence in both micromorphic formulations with increasing characteristic lengths L c . As observed in Fig. 22, the relaxed micromorphic continuum converges towards a finite energy, whereas the non-trivial boundary conditions on the full micromorphic continuum lead to boundary-layers and consequently, everincreasing energy for L c → ∞. The result is consistent with the problematic mentioned in Remark 3.5.

Conclusions and outlook
The relaxed micromorphic continuum theory introduces the Curl operator in the formulation of the free energy functional. As a result, the solution of the weak form lies in the combined space H 1 × H(curl). The Lax-Milgram theorem confirms this result by assuring existence and uniqueness for the combined space. Our benchmarks with a completely nodal finite element show its capacity to approximate solutions in the combined space. However, the tests also show its inability to find the exact solution for discontinuous microdistortion fields and the corresponding sub-optimal convergence. A comparison between the linear nodal and hybrid element formulations also reveals the difference in the arising elemental stiffness matrices, namely K nodal ∈ R 12×12 and K hybrid ∈ R 8×8 , resulting in slower computation times for the nodal element. In contrast, the hybrid element yields stable approximations and convergence rates for all tested scenarios, being capable of finding the exact solution also for discontinuous microdistortion fields. The relaxed micromorphic theory aims to capture the mechanical behaviour of metamaterials, highly homogeneous materials and the entire spectrum in between. To that end, the characteristic length L c takes the role of a weighting parameter, determining the influence of the energy from the dislocation density (the energy depending on the curl operator). The range of the characteristic length L c is an open topic of research into metamaterials. However, from a theoretical point of view, it may vary between zero and infinity. Our tests reveal the arising instability of convergence where increasingly large L c parameters are concerned and emergence of locking effects if linear nodal elements are chosen to approximate the microdistortion. For the case of the hybrid element, lost precision can be recovered via the formulation of the corresponding mixed problem. Locking effects in the nodal version of the microdistortion can be alleviated via higher order polynomials at the cost of increased dofs. In addition, also in L c = 0 setting, where the external moment ω vanishes, we recognize the optimality of using H(curl)-elements for the computation of the microdistortion, seeing as it is in fact the natural space for the microdistortion in this setting. Lastly, we recognize the advantage of the relaxed micromorphic continuum with regard to its ability to generate finite energies as L c → ∞ for arbitrary boundary conditions. These findings build the basis for the extension of the formulation to the fully three-dimensional or a statically condensed two-dimensional version of the full relaxed micromorphic continuum.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/.

A Derivation of the strong form
In order to find the strong form of the Euler-Lagrange equations to Problem 3.4 we start with the most general setting The Dirichlet and Neumann boundary parts of u and ζ and assume that | u D | > 0 (for Lax-Milgram solvability). We assume smooth fields such that we can integrate by parts. Using the Green identity where ν is the normal vector on the boundary, and splitting the boundary terms of the first weak form Eq. (2.5a), we find 2μ e (∇u − ζ ), ∇δu − δu, f dX As the Dirichlet data is directly incorporated into the space we have δu = 0 on u D and thus, for given Dirichlet data u, we obtain the strong form

B Constructing analytical solutions
The predefined fields are given by u and ζ . We redefine the variables of the strong form u * = u − u and ζ * = ζ − ζ and insert them into the partial differential equation