Showing 1–50 of 166 results for author: Neff, P

Yet another best approximation isotropic elasticity tensor in plane strain

Authors: Jendrik Voss, Panos Gourgiotis, Peter Lewintan, Adam Sky, Patrizio Neff

Abstract: For plane strain linear elasticity, given any anisotropic elasticity tensor $\mathbb{C}_{\rm aniso}$, we determine a best approximating isotropic counterpart $\mathbb{C}_{\rm iso}$. This is not done by using a distance measure on the space of positive definite elasticity tensors (Euclidean or logarithmic distance) but by considering two simple isotropic analytic solutions (center of dilatation and… ▽ More For plane strain linear elasticity, given any anisotropic elasticity tensor $\mathbb{C}_{\rm aniso}$, we determine a best approximating isotropic counterpart $\mathbb{C}_{\rm iso}$. This is not done by using a distance measure on the space of positive definite elasticity tensors (Euclidean or logarithmic distance) but by considering two simple isotropic analytic solutions (center of dilatation and concentrated couple) and best fitting these radial solutions to the numerical anisotropic solution based on $\mathbb{C}_{\rm aniso}$. The numerical solution is done via a finite element calculation, and the fitting via a subsequent quadratic error minimization. Thus, we obtain the two Lamé-moduli $μ$, $λ$ (or $μ$ and the bulk-modulus $κ$) of $\mathbb{C}_{\rm aniso}$. We observe that our so-determined isotropic tensor $\mathbb{C}_{\rm iso}$ coincides with neither the best logarithmic fit of Norris nor the best Euclidean fit. Our result calls into question the very notion of a best-fit isotropic elasticity tensor to a given anisotropic material. △ Less

Submitted 28 June, 2024; originally announced June 2024.

MSC Class: 74A10; 74B05; 74M25

Novel $H^\mathrm{dev}(\mathrm{Curl})$-conforming elements on regular triangulations and Clough--Tocher splits for the planar relaxed micromorphic model

Authors: Adam Sky, Michael Neunteufel, Peter Lewintan, Panos Gourgiotis, Andreas Zilian, Patrizio Neff

Abstract: In this work we present a consistent reduction of the relaxed micromorphic model to its corresponding two-dimensional planar model, such that its capacity to capture discontinuous dilatation fields is preserved. As a direct consequence of our approach, new conforming finite elements for $H^\mathrm{dev}(\mathrm{Curl},A)$ become necessary. We present two novel $H^\mathrm{dev}(\mathrm{Curl},A)$-confo… ▽ More In this work we present a consistent reduction of the relaxed micromorphic model to its corresponding two-dimensional planar model, such that its capacity to capture discontinuous dilatation fields is preserved. As a direct consequence of our approach, new conforming finite elements for $H^\mathrm{dev}(\mathrm{Curl},A)$ become necessary. We present two novel $H^\mathrm{dev}(\mathrm{Curl},A)$-conforming finite element spaces, of which one is a macro element based on Clough--Tocher splits, as well as primal and mixed variational formulations of the planar relaxed micromorphic model. Finally, we demonstrate the effectiveness of our approach with two numerical examples. △ Less

Submitted 23 May, 2024; originally announced May 2024.

Global regularity for a physically nonlinear version of the relaxed micromorphic model on Lipschitz domains

Authors: Dorothee Knees, Sebastian Owczarek, Patrizio Neff

Abstract: In this paper, we investigate the global higher regularity properties of weak solutions for a linear elliptic system coupled with a nonlinear Maxwell-type system defined on Lipschitz domains. The regularity result is established using a modified finite difference approach. These adjusted finite differences involve inner variations in conjunction with a Piola-type transformation to preserve the cur… ▽ More In this paper, we investigate the global higher regularity properties of weak solutions for a linear elliptic system coupled with a nonlinear Maxwell-type system defined on Lipschitz domains. The regularity result is established using a modified finite difference approach. These adjusted finite differences involve inner variations in conjunction with a Piola-type transformation to preserve the curl-structure within the matrix Maxwell system. The proposed method is further applied to the linear relaxed micromorphic model. As a result, for a physically nonlinear version of the relaxed micromorphic model, we demonstrate that for arbitrary $ε> 0$, the displacement vector $u$ belongs to $H^{\frac{3}{2}-ε}(Ω)$, and the microdistortion tensor $P$ belongs to $H^{\frac{1}{2}-ε}(Ω)$ while $\Curl P$ belongs to $H^{\frac{1}{2}-ε}(Ω)$. △ Less

Submitted 26 March, 2024; originally announced March 2024.

A computational approach to identify the material parameters of the relaxed micromorphic model

Authors: Mohammad Sarhil, Lisa Scheunemann, Peter Lewintan, Jörg Schröder, Patrizio Neff

Abstract: We determine the material parameters in the relaxed micromorphic generalized continuum model for a given periodic microstructure in this work. This is achieved through a least squares fitting of the total energy of the relaxed micromorphic homogeneous continuum to the total energy of the fully-resolved heterogeneous microstructure, governed by classical linear elasticity. The relaxed micromorphic… ▽ More We determine the material parameters in the relaxed micromorphic generalized continuum model for a given periodic microstructure in this work. This is achieved through a least squares fitting of the total energy of the relaxed micromorphic homogeneous continuum to the total energy of the fully-resolved heterogeneous microstructure, governed by classical linear elasticity. The relaxed micromorphic model is a generalized continuum that utilizes the $\Curl$ of a micro-distortion field instead of its full gradient as in the classical micromorphic theory, leading to several advantages and differences. The most crucial advantage is that it operates between two well-defined scales. These scales are determined by linear elasticity with microscopic and macroscopic elasticity tensors, which respectively bound the stiffness of the relaxed micromorphic continuum from above and below. While the macroscopic elasticity tensor is established a priori through standard periodic first-order homogenization, the microscopic elasticity tensor remains to be determined. Additionally, the characteristic length parameter, associated with curvature measurement, controls the transition between the micro- and macro-scales. Both the microscopic elasticity tensor and the characteristic length parameter are here determined using a computational approach based on the least squares fitting of energies. This process involves the consideration of an adequate number of quadratic deformation modes and different specimen sizes. We conduct a comparative analysis between the least square fitting results of the relaxed micromorphic model, the fitting of a skew-symmetric micro-distortion field (Cosserat-micropolar model), and the fitting of the classical micromorphic model with two different formulations for the curvature... △ Less

Submitted 28 January, 2024; originally announced January 2024.

Journal ref: Computer Methods in Applied Mechanics and Engineering 2024

On the representation of fourth and higher order anisotropic elasticity tensors in generalized continuum models

Authors: Marco Valerio d'Agostino, Robert J. Martin, Peter Lewintan, Davide Bernardini, Alexandre Danescu, Patrizio Neff

Abstract: The classification of all fourth order anisotropic tensor classes for classical linear elasticity is well known. We review the mathematical procedure behind this classification and extend it to fourth order elasticity tensors acting on non-symmetric matrices. These tensors naturally appear in generalised continuum models. Based on tensor-symmetrization we provide the most general form of these ten… ▽ More The classification of all fourth order anisotropic tensor classes for classical linear elasticity is well known. We review the mathematical procedure behind this classification and extend it to fourth order elasticity tensors acting on non-symmetric matrices. These tensors naturally appear in generalised continuum models. Based on tensor-symmetrization we provide the most general form of these tensors for the case of orthotropic, transversely isotropic, cubic and isotropic materials. We give a self-contained presentation and provide a detailed calculation for simple examples. △ Less

Submitted 8 January, 2024; originally announced January 2024.

An essay on deformation measures in isotropic thin shell theories. Bending versus curvature

Authors: Ionel-Dumitrel Ghiba, Peter Lewintan, Adam Sky, Patrizio Neff

Abstract: It has become commonplace for the stored energy function of any realistic shell model to align ``within first order" with the classical Koiter membrane-bending (flexural) shell model. In this paper, we assess whether certain extended Cosserat shell models are consistent with the classical linear Koiter model. In doing this, we observe that there are numerous reasons why a modified version of the c… ▽ More It has become commonplace for the stored energy function of any realistic shell model to align ``within first order" with the classical Koiter membrane-bending (flexural) shell model. In this paper, we assess whether certain extended Cosserat shell models are consistent with the classical linear Koiter model. In doing this, we observe that there are numerous reasons why a modified version of the classical Koiter model should be considered, a consensus reached not only by Koiter himself but also by Sanders and Budiansky, who independently developed the same theory during the same period. To provide a comprehensive overview of the strain measures employed in our Cosserat shell models, this paper presents them in a unified manner and compares them with the strain measures previously utilized in the literature. We show that all our new strain tensors either generalize (in the case of nonlinear constrained or unconstrained models) or coincide (in the case of the linear constrained model) with the strain tensors recognized as the ``best" or those possessing a well-defined geometric interpretation connected to bending or curvature. △ Less

Submitted 19 December, 2023; v1 submitted 17 December, 2023; originally announced December 2023.

Green's functions for the isotropic planar relaxed micromorphic model -- concentrated force and concentrated couple

Authors: Panos Gourgiotis, Gianluca Rizzi, Peter Lewintan, Davide Bernardini, Adam Sky, Angela Madeo, Patrizio Neff

Abstract: We derive the Green's functions (concentrated force and couple in an infinite space) for the isotropic planar relaxed micromorphic model. Since the relaxed micromorphic model particularises into the microstretch, Cosserat (micropolar), couple-stress, and linear elasticity model for certain choices of material parameters, we recover the fundamental solutions in all these cases. We derive the Green's functions (concentrated force and couple in an infinite space) for the isotropic planar relaxed micromorphic model. Since the relaxed micromorphic model particularises into the microstretch, Cosserat (micropolar), couple-stress, and linear elasticity model for certain choices of material parameters, we recover the fundamental solutions in all these cases. △ Less

Submitted 22 February, 2024; v1 submitted 18 October, 2023; originally announced November 2023.

Rayleigh waves in isotropic elastic materials with micro-voids

Authors: Emilian Bulgariu, Ionel-Dumitrel Ghiba, Hassam Khan, Patrizio Neff

Abstract: In this paper, we show that a general method introduced by Fu and Mielke allows to give a complete answer on the existence and uniqueness of a subsonic solution describing the propagation of surface waves in an isotropic half space modelled with the linear theory of isotropic elastic materials with micro-voids. Our result is valid for the entire class of materials admitting real wave propagation w… ▽ More In this paper, we show that a general method introduced by Fu and Mielke allows to give a complete answer on the existence and uniqueness of a subsonic solution describing the propagation of surface waves in an isotropic half space modelled with the linear theory of isotropic elastic materials with micro-voids. Our result is valid for the entire class of materials admitting real wave propagation which include auxetic materials (negative Poisson's ration) and composite materials with negative-stiffness inclusions (negative Young's modulus). Moreover, the used method allows to formulate a simple and complete numerical strategy for the computation of the solution. △ Less

Submitted 27 September, 2023; originally announced September 2023.

Comments: arXiv admin note: text overlap with arXiv:2104.13143

From frequency-dependent models to frequency-independent enriched continua for mechanical metamaterials

Authors: Gianluca Rizzi, Marco Valerio d'Agostino, Jendrik Voss, Davide Bernardini, Patrizio Neff, Angela Madeo

Abstract: Mechanical metamaterials have recently gathered increasing attention for their uncommon mechanical responses enabling unprecedented applications for elastic wave control. To model the mechanical response of large metamaterials' samples made up of base unit cells, so-called homogenization or upscaling techniques come into play trying to establish an equivalent continuum model describing these macro… ▽ More Mechanical metamaterials have recently gathered increasing attention for their uncommon mechanical responses enabling unprecedented applications for elastic wave control. To model the mechanical response of large metamaterials' samples made up of base unit cells, so-called homogenization or upscaling techniques come into play trying to establish an equivalent continuum model describing these macroscopic metamaterials' characteristics. A common approach is to assume a priori that the target continuum model is a classical linear Cauchy continuum featuring the macroscopic displacement as the only kinematical field. This implies that the parameters of such continuum models (density and/or elasticity tensors) must be considered to be frequency-dependent to capture the complex metamaterials' response in the frequency domain. These frequency-dependent models can be useful to describe some of the aforementioned macroscopic metamaterials' properties, yet, they suffer some drawbacks such as featuring negative masses and/or elastic coefficients in some frequency ranges. More than being counter-intuitive, this implies that the considered Cauchy continuum is not positive-definite for all the considered frequencies. In this paper, we present a procedure, based on the definition of extra kinematical variables (with respect to displacement alone) and the use of the inverse Fourier transform in time, to convert a frequency-dependent model into an enriched continuum model of the micromorphic type. All the parameters of the associated enriched model are constant (i.e., frequency-independent) and the model itself remains positive-definite for all the considered frequency ranges. The response of the frequency-dependent model and the associated micromorphic model coincide in the frequency domain, in particular when looking at the dispersion curves. △ Less

Submitted 15 February, 2024; v1 submitted 17 September, 2023; originally announced September 2023.

Comments: 37 pages, 10 figures

Explicit formula for the Gamma-convergence homogenized quadratic curvature energy in isotropic Cosserat shell models

Authors: Maryam Mohammadi Saem, Emilian Bulgariu, Ionel-Dumitrel Ghiba, Patrizio Neff

Abstract: We show how to explicitly compute the homogenized curvature energy appearing in the isotropic $Γ$-limit for flat and for curved initial configuration Cosserat shell models, when a parental three-dimensional minimization problem on $Ω\subset \mathbb{R}^3$ for a Cosserat energy based on the second order dislocation density tensor… ▽ More We show how to explicitly compute the homogenized curvature energy appearing in the isotropic $Γ$-limit for flat and for curved initial configuration Cosserat shell models, when a parental three-dimensional minimization problem on $Ω\subset \mathbb{R}^3$ for a Cosserat energy based on the second order dislocation density tensor $α:=\overline{R} ^T {\rm Curl}\,\overline{R} \in \mathbb{R}^{3\times 3}$, $\overline{R}\in {\rm SO}(3)$ is used. △ Less

Submitted 12 September, 2023; originally announced September 2023.

Novel $H(\mathrm{sym} \mathrm{Curl})$-conforming finite elements for the relaxed micromorphic sequence

Authors: Adam Sky, Michael Neunteufel, Peter Lewintan, Andreas Zilian, Patrizio Neff

Abstract: In this work we construct novel $H(\mathrm{sym} \mathrm{Curl})$-conforming finite elements for the recently introduced relaxed micromorphic sequence, which can be considered as the completion of the $\mathrm{div} \mathrm{Div}$-sequence with respect to the $H(\mathrm{sym} \mathrm{Curl})$-space. The elements respect $H(\mathrm{Curl})$-regularity and their lowest order versions converge optimally for… ▽ More In this work we construct novel $H(\mathrm{sym} \mathrm{Curl})$-conforming finite elements for the recently introduced relaxed micromorphic sequence, which can be considered as the completion of the $\mathrm{div} \mathrm{Div}$-sequence with respect to the $H(\mathrm{sym} \mathrm{Curl})$-space. The elements respect $H(\mathrm{Curl})$-regularity and their lowest order versions converge optimally for $[H(\mathrm{sym} \mathrm{Curl}) \setminus H(\mathrm{Curl})]$-fields. This work introduces a detailed construction, proofs of linear independence and conformity of the basis, and numerical examples. Further, we demonstrate an application to the computation of metamaterials with the relaxed micromorphic model. △ Less

Submitted 15 August, 2023; originally announced August 2023.

A global higher regularity result for the static relaxed micromorphic model on smooth domains

Authors: Dorothee Knees, Sebastian Owczarek, Patrizio Neff

Abstract: We derive a global higher regularity result for weak solutions of the linear relaxed micromorphic model on smooth domains. The governing equations consist of a linear elliptic system of partial differential equations that is coupled with a system of Maxwell-type. The result is obtained by combining a Helmholtz decomposition argument with regularity results for linear elliptic systems and the class… ▽ More We derive a global higher regularity result for weak solutions of the linear relaxed micromorphic model on smooth domains. The governing equations consist of a linear elliptic system of partial differential equations that is coupled with a system of Maxwell-type. The result is obtained by combining a Helmholtz decomposition argument with regularity results for linear elliptic systems and the classical embedding of $H(\mathrm{div};Ω)\cap H_0(\mathrm{curl};Ω)$ into $H^1(Ω)$. △ Less

Submitted 5 July, 2023; originally announced July 2023.

Comments: 14 pages

MSC Class: 35Q74; 35B65; 49N60; 74A35; 74G40

A geometrically nonlinear Cosserat shell model for orientable and non-orientable surfaces: Discretization with geometric finite elements

Authors: Lisa Julia Nebel, Oliver Sander, Mircea Bîrsan, Patrizio Neff

Abstract: We investigate discretizations of a geometrically nonlinear elastic Cosserat shell with nonplanar reference configuration originally introduced by Bîrsan, Ghiba, Martin, and Neff in 2019. The shell model includes curvature terms up to order 5 in the shell thickness, which are crucial to reliably simulate high-curvature deformations such as near-folds or creases. The original model is generalized t… ▽ More We investigate discretizations of a geometrically nonlinear elastic Cosserat shell with nonplanar reference configuration originally introduced by Bîrsan, Ghiba, Martin, and Neff in 2019. The shell model includes curvature terms up to order 5 in the shell thickness, which are crucial to reliably simulate high-curvature deformations such as near-folds or creases. The original model is generalized to shells that are not homeomorphic to a subset of $\mathbb{R}^2$. For this, we replace the originally planar parameter domain by an abstract two-dimensional manifold, and verify that the hyperelastic shell energy and three-dimensional reconstruction are invariant under changes of the local coordinate systems. This general approach allows to determine the elastic response for even non-orientable surfaces like the Möbius strip and the Klein bottle. We discretize the model with a geometric finite element method and, using that geometric finite elements are $H^1$-conforming, prove that the discrete shell model has a solution. Numerical tests then show the general performance and versatility of the model and discretization method. △ Less

Submitted 1 August, 2023; v1 submitted 22 April, 2023; originally announced April 2023.

The Gamma limit of the simple shear problem in nonlinear Cosserat elasticity

Authors: Thomas Blesgen, Patrizio Neff

Abstract: The zero and first order Gamma-limit of vanishing internal length scale are studied for the mechanical energy of a shear problem in geometrically nonlinear Cosserat elasticity. The convergence of the minimizers is shown and the limit functionals are characterized. The zero and first order Gamma-limit of vanishing internal length scale are studied for the mechanical energy of a shear problem in geometrically nonlinear Cosserat elasticity. The convergence of the minimizers is shown and the limit functionals are characterized. △ Less

Submitted 18 February, 2023; originally announced February 2023.

Comments: 27 pages, 6 figures

MSC Class: 35A15; 74B20

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

Authors: Adam Sky, Ingo Muench, Gianluca Rizzi, Patrizio Neff

Abstract: The relaxed micromorphic model is a generalized continuum model that is well-posed in the space $X = [H^1]^3 \times [H(\textrm{curl})]^3$. Consequently, finite element formulations of the model rely on $H^1$-conforming subspaces and Nédélec elements for discrete solutions of the corresponding variational problem. This work applies the recently introduced polytopal template methodology for the cons… ▽ More The relaxed micromorphic model is a generalized continuum model that is well-posed in the space $X = [H^1]^3 \times [H(\textrm{curl})]^3$. Consequently, finite element formulations of the model rely on $H^1$-conforming subspaces and Nédélec elements for discrete solutions of the corresponding variational problem. This work applies the recently introduced polytopal template methodology for the construction of Nédélec elements. This is done in conjunction with Bernstein-Bézier polynomials and dual numbers in order to compute hp-FEM solutions of the model. Bernstein-Bézier polynomials allow for optimal complexity in the assembly procedure due to their natural factorization into univariate Bernstein base functions. In this work, this characteristic is further augmented by the use of dual numbers in order to compute their values and their derivatives simultaneously. The application of the polytopal template methodology for the construction of the Nédélec base functions allows them to directly inherit the optimal complexity of the underlying Bernstein-Bézier basis. We introduce the Bernstein-Bézier basis along with its factorization to univariate Bernstein base functions, the principle of automatic differentiation via dual numbers and a detailed construction of Nédélec elements based on Bernstein-Bézier polynomials with the polytopal template methodology. This is complemented with a corresponding technique to embed Dirichlet boundary conditions, with emphasis on the consistent coupling condition. The performance of the elements is shown in examples of the relaxed micromorphic model. △ Less

Submitted 4 January, 2023; originally announced January 2023.

Korn-Maxwell-Sobolev inequalities for general incompatibilities

Authors: Franz Gmeineder, Peter Lewintan, Patrizio Neff

Abstract: We establish a family of coercive Korn-type inequalities for generalised incompatible fields in the superlinear growth regime under sharp criteria. This extends and unifies several previously known inequalities that are pivotal to the existence theory for a multitude of models in continuum mechanics in an optimal way. Different from our preceding work (ArXiv 2206.10373), where we focussed on the c… ▽ More We establish a family of coercive Korn-type inequalities for generalised incompatible fields in the superlinear growth regime under sharp criteria. This extends and unifies several previously known inequalities that are pivotal to the existence theory for a multitude of models in continuum mechanics in an optimal way. Different from our preceding work (ArXiv 2206.10373), where we focussed on the case $p=1$ and incompatibilities governed by the matrix curl, the case $p>1$ considered in the present paper gives us access to substantially stronger results from harmonic analysis but conversely deals with more general incompatibilities. Especially, we obtain sharp generalisations of recently proved inequalities by the last two authors and Müller (Calc. Var. PDE 60 (2021), 150) in the realm of incompatible Korn-type inequalities with conformally invariant dislocation energy. However, being applicable to higher order scenarios as well, our approach equally gives the first and sharp inequalities involving Kröner's incompability tensor $\mathrm{inc}$. △ Less

Submitted 7 December, 2023; v1 submitted 26 December, 2022; originally announced December 2022.

Comments: 33 pages, 3 figures

MSC Class: 35A23; 26D10; 35Q74/35Q75; 46E35

Linear constrained Cosserat-shell models including terms up to ${O}(h^5)$. Conditional and unconditional existence and uniqueness

Authors: Ionel-Dumitrel Ghiba, Patrizio Neff

Abstract: In this paper we linearise the recently introduced geometrically nonlinear constrained Cosserat-shell model. In the framework of the linear constrained Cosserat-shell model, we provide a comparison of our linear models with the classical linear Koiter shell model and the "best" first order shell model. For all proposed linear models we show existence and uniqueness based on a Korn's inequality for… ▽ More In this paper we linearise the recently introduced geometrically nonlinear constrained Cosserat-shell model. In the framework of the linear constrained Cosserat-shell model, we provide a comparison of our linear models with the classical linear Koiter shell model and the "best" first order shell model. For all proposed linear models we show existence and uniqueness based on a Korn's inequality for surfaces. △ Less

Submitted 22 November, 2022; originally announced November 2022.

Comments: arXiv admin note: text overlap with arXiv:2208.04574, arXiv:2010.14308

Regularity for a geometrically nonlinear flat Cosserat micropolar membrane shell with curvature

Authors: Andreas Gastel, Patrizio Neff

Abstract: We consider the rigorously derived thin shell membrane $Γ$-limit of a three-dimensional isotropic geometrically nonlinear Cosserat micropolar model and deduce full interior regularity of both the midsurface deformation $m:ω\subset{\mathbb R}^2\to{\mathbb R}^3$ and the orthogonal microrotation tensor field $R:ω\subset{\mathbb R}^2\to SO(3)$. The only further structural assumption is that the curvat… ▽ More We consider the rigorously derived thin shell membrane $Γ$-limit of a three-dimensional isotropic geometrically nonlinear Cosserat micropolar model and deduce full interior regularity of both the midsurface deformation $m:ω\subset{\mathbb R}^2\to{\mathbb R}^3$ and the orthogonal microrotation tensor field $R:ω\subset{\mathbb R}^2\to SO(3)$. The only further structural assumption is that the curvature energy depends solely on the uni-constant isotropic Dirichlet type energy term $|DR|^2$. We use Rivière's regularity techniques of harmonic map type systems for our system which couples harmonic maps to $SO(3)$ with a linear equation for $m$. The additional coupling term in the harmonic map equation is of critical integrability and can only be handled because of its special structure. △ Less

Submitted 19 November, 2022; originally announced November 2022.

MSC Class: 58E20 (Primary) 74G40; 74B20; 35J35 (Secondary)

Size-effects of metamaterial beams subjected to pure bending: on boundary conditions and parameter identification in the relaxed micromorphic model

Authors: Mohammad Sarhil, Lisa Scheunemann, Jörg Schröder, Patrizio Neff

Abstract: In this paper we model the size-effects of metamaterial beams under bending with the aid of the relaxed micromorphic continuum. We analyze first the size-dependent bending stiffness of heterogeneous fully discretized metamaterial beams subjected to pure bending loads. Two equivalent loading schemes are introduced which lead to a constant moment along the beam length with no shear force. The relaxe… ▽ More In this paper we model the size-effects of metamaterial beams under bending with the aid of the relaxed micromorphic continuum. We analyze first the size-dependent bending stiffness of heterogeneous fully discretized metamaterial beams subjected to pure bending loads. Two equivalent loading schemes are introduced which lead to a constant moment along the beam length with no shear force. The relaxed micromorphic model is employed then to retrieve the size-effects. We present a procedure for the determination of the material parameters of the relaxed micromorphic model based on the fact that the model operates between two well-defined scales. These scales are given by linear elasticity with micro and macro elasticity tensors which bound the relaxed micromorphic continuum from above and below, respectively. The micro elasticity tensor is specified as the maximum possible stiffness that is exhibited by the assumed metamaterial while the macro elasticity tensor is given by standard periodic first-order homogenization. For the identification of the micro elasticity tensor, two different approaches are shown which rely on affine and non-affine Dirichlet boundary conditions of candidate unit cell variants with the possible stiffest response. The consistent coupling condition is shown to allow the model to act on the whole intended range between macro and micro elasticity tensors for both loading cases. We fit the relaxed micromorphic model against the fully resolved metamaterial solution by controlling the curvature magnitude after linking it with the specimen's size. The obtained parameters of the relaxed micromorphic model are tested for two additional loading scenarios. △ Less

Submitted 28 March, 2023; v1 submitted 31 October, 2022; originally announced October 2022.

Journal ref: Computational Mechanics 2023

A local regularity result for the relaxed micromorphic model based on inner variations

Authors: Dorothee Knees, Sebastian Owczarek, Patrizio Neff

Abstract: In this paper we study local higher regularity properties of a linear elliptic system that is coupled with a system of Maxwell-type. The regularity result is proved by means of a modified finite difference argument. These modified finite differences are based on inner variations combined with a Piola-type transformation in order to preserve the $\curl$-structure in the Maxwell system. The result i… ▽ More In this paper we study local higher regularity properties of a linear elliptic system that is coupled with a system of Maxwell-type. The regularity result is proved by means of a modified finite difference argument. These modified finite differences are based on inner variations combined with a Piola-type transformation in order to preserve the $\curl$-structure in the Maxwell system. The result is applied to the relaxed micromorphic model. △ Less

Submitted 9 August, 2022; originally announced August 2022.

Comments: 12 pages

MSC Class: 35Q74; 35B65; 49N60; 74A35; 74G40

A linear isotropic Cosserat shell model including terms up to $O(h^5)$. Existence and uniqueness

Authors: Ionel-Dumitrel Ghiba, Mircea Birsan, Patrizio Neff

Abstract: In this paper we derive the linear elastic Cosserat shell model incorporating effects up to order $O(h^5)$ in the shell thickness $h$ as a particular case of the recently introduced geometrically nonlinear elastic Cosserat shell model. The existence and uniqueness of the solution is proven in suitable admissible sets. To this end, inequalities of Korn-type for shells are established which allow to… ▽ More In this paper we derive the linear elastic Cosserat shell model incorporating effects up to order $O(h^5)$ in the shell thickness $h$ as a particular case of the recently introduced geometrically nonlinear elastic Cosserat shell model. The existence and uniqueness of the solution is proven in suitable admissible sets. To this end, inequalities of Korn-type for shells are established which allow to show coercivity in the Lax-Milgram theorem. We are also showing an existence and uniqueness result for a truncated $O(h^3)$ model. Main issue is the suitable treatment of the curved reference configuration of the shell. Some connections to the classical Koiter membrane-bending model are highlighted. △ Less

Submitted 9 August, 2022; originally announced August 2022.

Comments: arXiv admin note: text overlap with arXiv:2010.14308, arXiv:2003.08594

A geometrically nonlinear Cosserat (micropolar) curvy shell model via Gamma convergence

Authors: Maryam Mohammadi Saem, Ionel-Dumitrel Ghiba, Patrizio Neff

Abstract: Using $Γ$-convergence arguments, we construct a nonlinear membrane-like Cosserat shell model on a curvy reference configuration starting from a geometrically nonlinear, physically linear three-dimensional isotropic Cosserat model. Even if the theory is of order $O(h)$ in the shell thickness $h$, by comparison to the membrane shell models proposed in classical nonlinear elasticity, beside the chang… ▽ More Using $Γ$-convergence arguments, we construct a nonlinear membrane-like Cosserat shell model on a curvy reference configuration starting from a geometrically nonlinear, physically linear three-dimensional isotropic Cosserat model. Even if the theory is of order $O(h)$ in the shell thickness $h$, by comparison to the membrane shell models proposed in classical nonlinear elasticity, beside the change of metric, the membrane-like Cosserat shell model is still capable to capture the transverse shear deformation and the {Cosserat}-curvature due to remaining Cosserat effects. We formulate the limit problem by scaling both unknowns, the deformation and the microrotation tensor, and by expressing the parental three-dimensional Cosserat energy with respect to a fictitious flat configuration. The model obtained via $Γ$-convergence is similar to the membrane {(no $O(h^3)$ flexural terms, but still depending on the Cosserat-curvature)} Cosserat shell model derived via a derivation approach but these two models do not coincide. Comparisons to other shell models are also included. △ Less

Submitted 18 July, 2022; originally announced July 2022.

Optimal incompatible Korn-Maxwell-Sobolev inequalities in all dimensions

Authors: Franz Gmeineder, Peter Lewintan, Patrizio Neff

Abstract: We characterise all linear maps $\mathcal{A}\colon\mathbb{R}^{n\times n}\to\mathbb{R}^{n\times n}$ such that, for $1\leq p<n$, \begin{align*} \|P\|_{L^{p^{*}}(\mathbb{R}^{n})}\leq c\,\Big(\|\mathcal{A}[P]\|_{L^{p^{*}}(\mathbb{R}^{n})}+\|\mathrm{Curl} P\|_{L^{p}(\mathbb{R}^{n})} \Big) \end{align*} holds for all compactly supported $P\in C_{c}^{\infty}(\mathbb{R}^{n};\mathbb{R}^{n\times n})$, where… ▽ More We characterise all linear maps $\mathcal{A}\colon\mathbb{R}^{n\times n}\to\mathbb{R}^{n\times n}$ such that, for $1\leq p<n$, \begin{align*} \|P\|_{L^{p^{*}}(\mathbb{R}^{n})}\leq c\,\Big(\|\mathcal{A}[P]\|_{L^{p^{*}}(\mathbb{R}^{n})}+\|\mathrm{Curl} P\|_{L^{p}(\mathbb{R}^{n})} \Big) \end{align*} holds for all compactly supported $P\in C_{c}^{\infty}(\mathbb{R}^{n};\mathbb{R}^{n\times n})$, where $\mathrm{Curl} P$ displays the matrix curl. Being applicable to incompatible, that is, non-gradient matrix fields as well, such inequalities generalise the usual Korn-type inequalities used e.g. in linear elasticity. Different from previous contributions, the results gathered in this paper are applicable to all dimensions and optimal. This particularly necessitates the distinction of different constellations between the ellipticities of $\mathcal{A}$, the integrability $p$ and the underlying space dimensions $n$, especially requiring a finer analysis in the two-dimensional situation. △ Less

Submitted 29 June, 2023; v1 submitted 17 June, 2022; originally announced June 2022.

MSC Class: 35A23; 26D10; 35Q74/35Q75; 46E35

Journal ref: CalcVar PDE 62 (2023)

Simple shear in nonlinear Cosserat micropolar elasticity: Existence of minimizers, numerical simulations and occurrence of microstructure

Authors: Thomas Blesgen, Patrizio Neff

Abstract: Deformation microstructure is studied for a 1D-shear problem in geometrically nonlinear Cosserat elasticity. Microstructure solutions are described analytically and numerically for zero characteristic length scale. Deformation microstructure is studied for a 1D-shear problem in geometrically nonlinear Cosserat elasticity. Microstructure solutions are described analytically and numerically for zero characteristic length scale. △ Less

Submitted 13 September, 2022; v1 submitted 16 June, 2022; originally announced June 2022.

Cosserat micropolar elasticity: classical Eringen vs. dislocation form

Authors: Ionel-Dumitrel Ghiba, Gianluca Rizzi, Angela Madeo, Patrizio Neff

Abstract: In this paper we do a comparative presentation of the linear isotropic Cosserat elastic model from two perspectives: the classical Mindlin-Eringen-Nowacki description in terms of a microrotation vector and a new formulation in terms of a skew-symmetric matrix and a curvature energy in dislocation form. We provide the reader with an alternative representation of the energy for the isotropic Cossera… ▽ More In this paper we do a comparative presentation of the linear isotropic Cosserat elastic model from two perspectives: the classical Mindlin-Eringen-Nowacki description in terms of a microrotation vector and a new formulation in terms of a skew-symmetric matrix and a curvature energy in dislocation form. We provide the reader with an alternative representation of the energy for the isotropic Cosserat model to ease the comparison with the relaxed micromorphic model and the geometrically nonlinear Cosserat elastic model. △ Less

Submitted 6 June, 2022; originally announced June 2022.

Unfolding engineering metamaterials design: relaxed micromorphic modeling of large-scale acoustic meta-structures

Authors: F. Demore, G. Rizzi, M. Collet, P. Neff, A. Madeo

Abstract: In this paper, we present a unit cell showing a band-gap in the lower acoustic domain. The corresponding metamaterial is made up of a periodic arrangement of this unit cell. We rigorously show that the relaxed micromorphic model can be used for metamaterials' design at large scales as soon as suficiently large specimens are considered. We manufacture the metamaterial via metal etching procedures a… ▽ More In this paper, we present a unit cell showing a band-gap in the lower acoustic domain. The corresponding metamaterial is made up of a periodic arrangement of this unit cell. We rigorously show that the relaxed micromorphic model can be used for metamaterials' design at large scales as soon as suficiently large specimens are considered. We manufacture the metamaterial via metal etching procedures applied to a titanium plate so as to show that its production for realistic applications is viable. Experimental tests are also carried out confirming that the metamaterials' response is in good agreement with the theoretical design. In order to show that our micromorphic model opens unprecedented possibilities in metastructural design, we conceive a finite-size structure that is able to focus elastic energy in a confined region, thus enabling its possible subsequent re-use. Indeed, thanks to the introduction of a well-posed set of micromorphic boundary conditions, we can combine different metamaterials and classical Cauchy materials in such a way that the elastic energy produced by a source of vibrations is focused in specific collection points. The design of this structure would have not been otherwise possible (via e.g., direct simulations), due to the large dimensions of the metastructure, couting hundreds of unit cells. △ Less

Submitted 14 April, 2022; originally announced April 2022.

Modeling a labyrinthine acoustic metamaterial through an inertia-augmented relaxed micromorphic approach

Authors: Jendrik Voss, Gianluca Rizzi, Patrizio Neff, Angela Madeo

Abstract: We present an inertia-augmented relaxed micromorphic model that enriches the relaxed micromorphic model previously introduced by the authors via a term $\text{Curl}\dot{P}$ in the kinetic energy density. This enriched model allows us to obtain a good overall fitting of the dispersion curves while introducing the new possibility of describing modes with negative group velocity that are known to tri… ▽ More We present an inertia-augmented relaxed micromorphic model that enriches the relaxed micromorphic model previously introduced by the authors via a term $\text{Curl}\dot{P}$ in the kinetic energy density. This enriched model allows us to obtain a good overall fitting of the dispersion curves while introducing the new possibility of describing modes with negative group velocity that are known to trigger negative refraction effects. The inertia-augmented model also allows for more freedom on the values of the asymptotes corresponding to the cut-offs. In the previous version of the relaxed micromorphic model, the asymptote of one curve (pressure or shear) is always bounded by the cut-off of the following curve of the same type. This constraint does not hold anymore in the enhanced version of the model. While the obtained curves' fitting is of good quality overall, a perfect quantitative agreement must still be reached for very small wavelengths that are close to the size of the unit cell. △ Less

Submitted 5 April, 2022; originally announced April 2022.

MSC Class: 74A10; 74B05; 74J05; 74M25

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

Authors: Adam Sky, Ingo Muench, Patrizio Neff

Abstract: In this work we test the numerical behaviour of matrix-valued fields approximated by finite element subspaces of $[\mathit{H}^1]^{3\times 3}$, $[\mathit{H}(\mathrm{curl})]^3$ and $\mathit{H}(\mathrm{sym}\mathrm{Curl})$ for a linear abstract variational problem connected to the relaxed micromorphic model. The formulation of the corresponding finite elements is introduced, followed by numerical benc… ▽ More In this work we test the numerical behaviour of matrix-valued fields approximated by finite element subspaces of $[\mathit{H}^1]^{3\times 3}$, $[\mathit{H}(\mathrm{curl})]^3$ and $\mathit{H}(\mathrm{sym}\mathrm{Curl})$ for a linear abstract variational problem connected to the relaxed micromorphic model. The formulation of the corresponding finite elements is introduced, followed by numerical benchmarks and our conclusions. The relaxed micromorphic continuum model reduces the continuity assumptions of the classical micromorphic model by replacing the full gradient of the microdistortion in the free energy functional with the Curl. This results in a larger solution space for the microdistortion, namely $[\mathit{H}(\mathrm{curl})]^3$ in place of the classical $[\mathit{H}^1]^{3\times 3}$. The continuity conditions on the microdistortion can be further weakened by taking only the symmetric part of the Curl. As shown in recent works, the new appropriate space for the microdistortion is then $\mathit{H}(\mathrm{sym}\mathrm{Curl})$. The newly introduced space gives rise to a new differential complex for the relaxed micromorphic continuum theory. △ Less

Submitted 17 February, 2022; originally announced February 2022.

Primal and mixed finite element formulations for the relaxed micromorphic model

Authors: Adam Sky, Michael Neunteufel, Ingo Muench, Joachim Schöberl, Patrizio Neff

Abstract: The classical Cauchy continuum theory is suitable to model highly homogeneous materials. However, many materials, such as porous media or metamaterials, exhibit a pronounced microstructure. As a result, the classical continuum theory cannot capture their mechanical behaviour without fully resolving the underlying microstructure. In terms of finite element computations, this can be done by modellin… ▽ More The classical Cauchy continuum theory is suitable to model highly homogeneous materials. However, many materials, such as porous media or metamaterials, exhibit a pronounced microstructure. As a result, the classical continuum theory cannot capture their mechanical behaviour without fully resolving the underlying microstructure. In terms of finite element computations, this can be done by modelling the entire body, including every interior cell. The relaxed micromorphic continuum offers an alternative method by instead enriching the kinematics of the mathematical model. The theory introduces a microdistortion field, encompassing nine extra degrees of freedom for each material point. The corresponding elastic energy functional contains the gradient of the displacement field, the microdistortion field and its Curl (the micro-dislocation). Therefore, the natural spaces of the fields are $[\mathit{H}^1]^3$ for the displacement and $[\mathit{H}(\mathrm{curl})]^3$ for the microdistortion, leading to unusual finite element formulations. In this work we describe the construction of appropriate finite elements using Nédélec and Raviart-Thomas subspaces, encompassing solutions to the orientation problem and the discrete consistent coupling condition. Further, we explore the numerical behaviour of the relaxed micromorphic model for both a primal and a mixed formulation. The focus of our benchmarks lies in the influence of the characteristic length $L_\mathrm{c}$ and the correlation to the classical Cauchy continuum theory. △ Less

Submitted 17 February, 2022; originally announced February 2022.

Numerical approaches for investigating quasiconvexity in the context of Morrey's conjecture

Authors: Jendrik Voss, Robert J. Martin, Oliver Sander, Siddhant Kumar, Dennis M. Kochmann, Patrizio Neff

Abstract: Deciding whether a given function is quasiconvex is generally a difficult task. Here, we discuss a number of numerical approaches that can be used in the search for a counterexample to the quasiconvexity of a given function $W$. We will demonstrate these methods using the planar isotropic rank-one convex function \[ W_{\rm magic}^+(F)=\frac{λ_{\rm max}}{λ_{\rm min}}-\log\frac{λ_{\rm max}}{λ_{\rm… ▽ More Deciding whether a given function is quasiconvex is generally a difficult task. Here, we discuss a number of numerical approaches that can be used in the search for a counterexample to the quasiconvexity of a given function $W$. We will demonstrate these methods using the planar isotropic rank-one convex function \[ W_{\rm magic}^+(F)=\frac{λ_{\rm max}}{λ_{\rm min}}-\log\frac{λ_{\rm max}}{λ_{\rm min}}+\log\det F=\frac{λ_{\rm max}}{λ_{\rm min}}+2\logλ_{\rm min}\,, \] where $λ_{\rm max}\geqλ_{\rm min}$ are the singular values of $F$, as our main example. In a previous contribution, we have shown that quasiconvexity of this function would imply quasiconvexity for all rank-one convex isotropic planar energies $W:\operatorname{GL}^+(2)\rightarrow\mathbb{R}$ with an additive volumetric-isochoric split of the form \[ W(F)=W_{\rm iso}(F)+W_{\rm vol}(\det F)=\widetilde W_{\rm iso}\bigg(\frac{F}{\sqrt{\det F}}\bigg)+W_{\rm vol}(\det F) \] with a concave volumetric part. This example is therefore of particular interest with regard to Morrey's open question whether or not rank-one convexity implies quasiconvexity in the planar case. △ Less

Submitted 17 January, 2022; originally announced January 2022.

MSC Class: 74B20; 74A10; 26B25

The consistent coupling boundary condition for the classical micromorphic model: existence, uniqueness and interpretation of parameters

Authors: Marco Valerio d'Agostino, Gianluca Rizzi, Hassam Khan, Peter Lewintan, Angela Madeo, Patrizio Neff

Abstract: We consider the classical Mindlin-Eringen linear micromorphic model with a new strictly weaker set of displacement boundary conditions. The new consistent coupling condition aims at minimizing spurious influences from arbitrary boundary prescription for the additional microdistortion field P. In effect, P is now only required to match the tangential derivative of the classical displacement u which… ▽ More We consider the classical Mindlin-Eringen linear micromorphic model with a new strictly weaker set of displacement boundary conditions. The new consistent coupling condition aims at minimizing spurious influences from arbitrary boundary prescription for the additional microdistortion field P. In effect, P is now only required to match the tangential derivative of the classical displacement u which is known at the Dirichlet-part of the boundary. We derive the full boundary condition, in adding the missing Neumann condition on the Dirichlet-part. We show existence and uniqueness of the static problem for this weaker boundary condition. These results are based on new coercive inequalities for incompatible tensor fields with prescribed tangential part. Finally, we show that compared to classical Dirichlet conditions on u and P, the new boundary condition modifies the interpretation of the constitutive parameters. △ Less

Submitted 22 December, 2021; originally announced December 2021.

Lagrange and $H(\operatorname{curl},{\cal B})$ based Finite Element formulations for the relaxed micromorphic model

Authors: Jörg Schröder, Mohammad Sarhil, Lisa Scheunemann, Patrizio Neff

Abstract: Modeling the unusual mechanical properties of metamaterials is a challenging topic for the mechanics community and enriched continuum theories are promising computational tools for such materials. The so-called relaxed micromorphic model has shown many advantages in this field. In this contribution, we present the significant aspects related to the relaxed micromorphic model realization with the f… ▽ More Modeling the unusual mechanical properties of metamaterials is a challenging topic for the mechanics community and enriched continuum theories are promising computational tools for such materials. The so-called relaxed micromorphic model has shown many advantages in this field. In this contribution, we present the significant aspects related to the relaxed micromorphic model realization with the finite element method. The variational problem is derived and different FEM-formulations for the two-dimensional case are presented. These are a nodal standard formulation $H^1({\cal B}) \times H^1({\cal B})$ and a nodal-edge formulation $H^1({\cal B}) \times H(\operatorname{curl}, {\cal B})$, where the latter employs the Nédélec space. However, the implementation of higher-order Nédélec elements is not trivial and requires some technicalities which are demonstrated. We discuss the convergence behavior of Lagrange-type and tangential-conforming finite element discretizations. Moreover, we analyze the characteristic length effect on the different components of the model and reveal how the size-effect property is captured via this characteristic length. △ Less

Submitted 1 December, 2021; originally announced December 2021.

Journal ref: Computational Mechanics 2022

Metamaterial shields for inner protection and outer tuning through a relaxed micromorphic approach

Authors: Gianluca Rizzi, Patrizio Neff, Angela Madeo

Abstract: In this paper, a coherent boundary value problem to model metamaterials' behavior based on the relaxed micromorphic model is established. This boundary value problem includes well-posed boundary conditions, thus disclosing the possibility of exploring the scattering patterns of finite-size metamaterials' specimens. Thanks to the simplified model's structure (few frequency- and angle-independent pa… ▽ More In this paper, a coherent boundary value problem to model metamaterials' behavior based on the relaxed micromorphic model is established. This boundary value problem includes well-posed boundary conditions, thus disclosing the possibility of exploring the scattering patterns of finite-size metamaterials' specimens. Thanks to the simplified model's structure (few frequency- and angle-independent parameters), we are able to unveil the scattering metamaterial's response for a wide range of frequencies and angles of propagation of the incident wave. These results are an important stepping stone towards the conception of more complex large-scale meta-structures that can control elastic waves and recover energy. △ Less

Submitted 23 November, 2021; originally announced November 2021.

Comments: 12 pages, 9 figures

Analytical solution of the uniaxial extension problem for the relaxed micromorphic continuum and other generalized continua (including full derivations)

Authors: Gianluca Rizzi, Hassam Khan, Ionel-Dumitrel Ghiba, Angela Madeo, Patrizio Neff

Abstract: We derive analytical solutions for the uniaxial extension problem for the relaxed micromorphic continuum and other generalized continua. These solutions may help in the identification of material parameters of generalized continua which are able to disclose size-effects. We derive analytical solutions for the uniaxial extension problem for the relaxed micromorphic continuum and other generalized continua. These solutions may help in the identification of material parameters of generalized continua which are able to disclose size-effects. △ Less

Submitted 19 July, 2021; originally announced July 2021.

Comments: 16 pages, 5 figures. arXiv admin note: substantial text overlap with arXiv:2104.11322

Polyconvex anisotropic hyperelasticity with neural networks

Authors: Dominik K. Klein, Mauricio Fernández, Robert J. Martin, Patrizio Neff, Oliver Weeger

Abstract: In the present work, two machine learning based constitutive models for finite deformations are proposed. Using input convex neural networks, the models are hyperelastic, anisotropic and fulfill the polyconvexity condition, which implies ellipticity and thus ensures material stability. The first constitutive model is based on a set of polyconvex, anisotropic and objective invariants. The second ap… ▽ More In the present work, two machine learning based constitutive models for finite deformations are proposed. Using input convex neural networks, the models are hyperelastic, anisotropic and fulfill the polyconvexity condition, which implies ellipticity and thus ensures material stability. The first constitutive model is based on a set of polyconvex, anisotropic and objective invariants. The second approach is formulated in terms of the deformation gradient, its cofactor and determinant, uses group symmetrization to fulfill the material symmetry condition, and data augmentation to fulfill objectivity approximately. The extension of the dataset for the data augmentation approach is based on mechanical considerations and does not require additional experimental or simulation data. The models are calibrated with highly challenging simulation data of cubic lattice metamaterials, including finite deformations and lattice instabilities. A moderate amount of calibration data is used, based on deformations which are commonly applied in experimental investigations. While the invariant-based model shows drawbacks for several deformation modes, the model based on the deformation gradient alone is able to reproduce and predict the effective material behavior very well and exhibits excellent generalization capabilities. In addition, the models are calibrated with transversely isotropic data, generated with an analytical polyconvex potential. For this case, both models show excellent results, demonstrating the straightforward applicability of the polyconvex neural network constitutive models to other symmetry groups. △ Less

Submitted 25 November, 2021; v1 submitted 20 June, 2021; originally announced June 2021.

Morrey's conjecture for the planar volumetric-isochoric split. Part I: least convex energy functions

Authors: Jendrik Voss, Robert J. Martin, Ionel-Dumitrel Ghiba, Patrizio Neff

Abstract: We consider Morrey's open question whether rank-one convexity already implies quasiconvexity in the planar case. For some specific families of energies, there are precise conditions known under which rank-one convexity even implies polyconvexity. We will extend some of these findings to the more general family of energies $W:\operatorname{GL}^+(n)\rightarrow\mathbb{R}$ with an additive volumetric-… ▽ More We consider Morrey's open question whether rank-one convexity already implies quasiconvexity in the planar case. For some specific families of energies, there are precise conditions known under which rank-one convexity even implies polyconvexity. We will extend some of these findings to the more general family of energies $W:\operatorname{GL}^+(n)\rightarrow\mathbb{R}$ with an additive volumetric-isochoric split, i.e. \[ W(F)=W_{\rm iso}(F)+W_{\rm vol}(\det F)=\widetilde W_{\rm iso}\bigg(\frac{F}{\sqrt{\det F}}\bigg)+W_{\rm vol}(\det F)\,, \] which is the natural finite extension of isotropic linear elasticity. Our approach is based on a condition for rank-one convexity which was recently derived from the classical two-dimensional criterion by Knowles and Sternberg and consists of a family of one-dimensional coupled differential inequalities. We identify a number of \enquote{least} rank-one convex energies and, in particular, show that for planar volumetric-isochorically split energies with a concave volumetric part, the question of whether rank-one convexity implies quasiconvexity can be reduced to the open question of whether the rank-one convex energy function \[ W_{\rm magic}^+(F)=\frac{λ_{\rm max}}{λ_{\rm min}}-\log\frac{λ_{\rm max}}{λ_{\rm min}}+\log\det F=\frac{λ_{\rm max}}{λ_{\rm min}}-2\logλ_{\rm min} \] is quasiconvex. In addition, we demonstrate that under affine boundary conditions, $W_{\rm magic}^+(F)$ allows for non-trivial inhomogeneous deformations with the same energy level as the homogeneous solution, and show a surprising connection to the work of Burkholder and Iwaniec in the field of complex analysis. △ Less

Submitted 30 June, 2021; v1 submitted 22 June, 2021; originally announced June 2021.

MSC Class: 74B20; 74A10; 26B25

Boundary and interface conditions in the relaxed micromorphic model: exploring finite-size metastructures for elastic wave control

Authors: Gianluca Rizzi, Marco Valerio d'Agostino, Patrizio Neff, Angela Madeo

Abstract: In this paper, we establish well-posed boundary and interface conditions for the relaxed micromorphic model that are able to unveil the scattering response of fully finite-size metamaterials' samples. The resulting relaxed micromorphic boundary value problem is implemented in finite element simulations describing the scattering of a square metamaterial's sample whose side counts 9 unit cells. The… ▽ More In this paper, we establish well-posed boundary and interface conditions for the relaxed micromorphic model that are able to unveil the scattering response of fully finite-size metamaterials' samples. The resulting relaxed micromorphic boundary value problem is implemented in finite element simulations describing the scattering of a square metamaterial's sample whose side counts 9 unit cells. The results are validated against a direct finite element simulation encoding all the details of the underlying metamaterial's microstructure. The relaxed micromorphic model can recover the scattering metamaterial's behavior for a wide range of frequencies and for all possible angles of incidence, thus showing that it is suitable to describe dynamic anisotropy. Finally, thanks to the model's computational performances, we can design a metastructure combining metamaterials and classical materials in such a way that it acts as a protection device while providing energy focusing in specific collection points. These results open important perspectives for the short-term design of sustainable structures that can control elastic waves and recover energy. △ Less

Submitted 3 May, 2021; originally announced May 2021.

Comments: 16 pages, 11 figures

Existence and uniqueness of Rayleigh waves in isotropic elastic Cosserat materials and algorithmic aspects

Authors: Hassam Khan, Ionel-Dumitrel Ghiba, Angela Madeo, Patrizio Neff

Abstract: We discuss the propagation of surface waves in an isotropic half space modelled with the linear Cosserat theory of isotropic elastic materials. To this aim we use a method based on the algebraic analysis of the surface impedance matrix and on the algebraic Riccati equation, and which is independent of the common Stroh formalism. Due to this method, a new algorithm which determines the amplitudes a… ▽ More We discuss the propagation of surface waves in an isotropic half space modelled with the linear Cosserat theory of isotropic elastic materials. To this aim we use a method based on the algebraic analysis of the surface impedance matrix and on the algebraic Riccati equation, and which is independent of the common Stroh formalism. Due to this method, a new algorithm which determines the amplitudes and the wave speed in the theory of isotropic elastic Cosserat materials is described. Moreover, the method allows to prove the existence and uniqueness of a subsonic solution of the secular equation, a problem which remains unsolved in almost all generalised linear theories of elastic materials. Since the results are suitable to be used for numerical implementations, we propose two numerical algorithms which are viable for any elastic material. Explicit numerical calculations are made for alumunium-epoxy in the context of the Cosserat model. Since the novel form of the secular equation for isotropic elastic material has not been explicitly derived elsewhere, we establish it in this paper, too. △ Less

Submitted 1 February, 2022; v1 submitted 27 April, 2021; originally announced April 2021.

MSC Class: 74J15; 74M25; 74H05; 74J05

Analytical solution of the cylindrical torsion problem for the relaxed micromorphic continuum and other generalized continua (including full derivations)

Authors: Gianluca Rizzi, Geralf Hütter, Hassam Khan, Ionel Dumitrel Ghiba, Angela Madeo, Patrizio Neff

Abstract: We solve the St.Venant torsion problem for an infinite cylindrical rod whose behaviour is described by a family of isotropic generalized continua, including the relaxed micromorphic and classical micromorphic model. The results can be used to determine the material parameters of these models. Special attention is given to the possible nonphysical stiffness singularity for a vanishing rod diameter,… ▽ More We solve the St.Venant torsion problem for an infinite cylindrical rod whose behaviour is described by a family of isotropic generalized continua, including the relaxed micromorphic and classical micromorphic model. The results can be used to determine the material parameters of these models. Special attention is given to the possible nonphysical stiffness singularity for a vanishing rod diameter, since slender specimens are in general described as stiffer. △ Less

Submitted 7 April, 2021; originally announced April 2021.

Comments: 41 pages, 29 figures

On in-plane drill rotations for Cosserat surfaces

Authors: Maryam Mohammadi Saem, Peter Lewintan, Patrizio Neff

Abstract: We show under some natural smoothness assumptions that pure in-plane drill rotations as deformation mappings of a $C^2$-smooth regular shell surface to another one parametrized over the same domain are impossible provided that the rotations are fixed at a portion of the boundary. Put otherwise, if the tangent vectors of the new surface are obtained locally by only rotating the given tangent vector… ▽ More We show under some natural smoothness assumptions that pure in-plane drill rotations as deformation mappings of a $C^2$-smooth regular shell surface to another one parametrized over the same domain are impossible provided that the rotations are fixed at a portion of the boundary. Put otherwise, if the tangent vectors of the new surface are obtained locally by only rotating the given tangent vectors, and if these rotations have a rotation axis which coincides everywhere with the normal of the initial surface, then the two surfaces are equal provided they coincide at a portion of the boundary. In the language of differential geometry of surfaces we show that any isometry which leaves normals invariant and which coincides with the given surface at a portion of the boundary, is the identity mapping. △ Less

Submitted 16 August, 2021; v1 submitted 20 February, 2021; originally announced February 2021.

MSC Class: 74B20; 74A20; 74A10; 74K25; 53A05; 53A10

Journal ref: Proceedings of the Royal Society A, Mathematical, Physical and Engineering Sciences (2021)

Towards the conception of complex engineering meta-structures: relaxed-micromorphic modelling of low-frequency mechanical diodes/high-frequency screens

Authors: Gianluca Rizzi, Domenico Tallarico, Patrizio Neff, Angela Madeo

Abstract: In this paper we show that an enriched continuum model of the micromorphic type (Relaxed Micromorphic Model) can be used to model metamaterials' response in view of their use for meta-structural design. We focus on the fact that the reduced model's structure, coupled with the introduction of well-posed interface conditions, allows us to easily test different combinations of metamaterials' and clas… ▽ More In this paper we show that an enriched continuum model of the micromorphic type (Relaxed Micromorphic Model) can be used to model metamaterials' response in view of their use for meta-structural design. We focus on the fact that the reduced model's structure, coupled with the introduction of well-posed interface conditions, allows us to easily test different combinations of metamaterials' and classical-materials bricks, so that we can eventually end-up with the conception of a meta-structure acting as a mechanical diode for low/medium frequencies and as a total screen for higher frequencies. Thanks to the reduced model's structure, we are also able to optimize this meta-structure so that the diode-behaviour is enhanced for both "pressure" and "shear" incident waves and for all possible angles of incidence. △ Less

Submitted 26 May, 2022; v1 submitted 21 December, 2020; originally announced December 2020.

Comments: 21 pages, 20 figures (47 pictures). arXiv admin note: text overlap with arXiv:2007.14942

Analytical solutions of the cylindrical bending problem for the relaxed micromorphic continuum and other generalized continua (including full derivations)

Authors: Gianluca Rizzi, Geralf Hütter, Angela Madeo, Patrizio Neff

Abstract: We consider the cylindrical bending problem for an infinite plate as modelled with a family of generalized continuum models, including the micromorphic approach. The models allow to describe length scale effects in the sense that thinner specimens are comparatively stiffer. We provide the analytical solution for each case and exhibit the predicted bending stiffness. The relaxed micromorphic contin… ▽ More We consider the cylindrical bending problem for an infinite plate as modelled with a family of generalized continuum models, including the micromorphic approach. The models allow to describe length scale effects in the sense that thinner specimens are comparatively stiffer. We provide the analytical solution for each case and exhibit the predicted bending stiffness. The relaxed micromorphic continuum shows bounded bending stiffness for arbitrary thin specimens, while classical micromorphic continuum or gradient elasticity as well as Cosserat models [35] exhibit unphysical unbounded bending stiffness for arbitrary thin specimens. This finding highlights the advantage of using the relaxed micromorphic model, which has a definite limit stiffness for small samples and which aids in identifying the relevant material parameters. △ Less

Submitted 6 December, 2020; originally announced December 2020.

Comments: 50 pages, 28 figures (38 pictures)

Korn inequalities for incompatible tensor fields in three space dimensions with conformally invariant dislocation energy

Authors: Peter Lewintan, Stefan Müller, Patrizio Neff

Abstract: Let $Ω\subset \mathbb{R}^3$ be an open and bounded set with Lipschitz boundary and outward unit normal $ν$. For $1<p<\infty$ we establish an improved version of the generalized $L^p$-Korn inequality for incompatible tensor fields $P$ in the new Banach space $… ▽ More Let $Ω\subset \mathbb{R}^3$ be an open and bounded set with Lipschitz boundary and outward unit normal $ν$. For $1<p<\infty$ we establish an improved version of the generalized $L^p$-Korn inequality for incompatible tensor fields $P$ in the new Banach space $$ W^{1,\,p,\,r}_0(\operatorname{dev}\operatorname{sym}\operatorname{Curl}; Ω,\mathbb R^{3\times3}) = \{ P \in L^p(Ω,\mathbb R^{3\times3}) \mid \operatorname{dev} \operatorname{sym} \operatorname{Curl} P \in L^r(Ω,\mathbb R^{3\times3}),\ \operatorname{dev} \operatorname{sym} (P \times ν) = 0 \text{ on $\partial Ω$}\} $$ where $$ r \in [1, \infty), \qquad \frac1r \le \frac1p + \frac13, \qquad r >1 \quad \text{if $p = \frac32$.}$$ Specifically, there exists a constant $c=c(p,Ω,r)>0$ such that the inequality \[ \|P \|_{L^p}\leq c\,\left(\|\operatorname{sym} P \|_{L^p} + \|\operatorname{dev}\operatorname{sym} \operatorname{Curl} P \|_{L^{r}}\right) \] holds for all tensor fields $P\in W^{1,\,p, \, r}_0(\operatorname{dev}\operatorname{sym}\operatorname{Curl})$. Here, $\operatorname{dev} X := X -\frac13 \operatorname{tr}(X)\,\mathbb{1}$ denotes the deviatoric (trace-free) part of a $3 \times 3$ matrix $X$ and the boundary condition is understood in a suitable weak sense. △ Less

Submitted 5 July, 2021; v1 submitted 20 November, 2020; originally announced November 2020.

Comments: Ref.numbers are now visible

MSC Class: Primary: 35A23; Secondary: 35B45; 35Q74; 46E35

Journal ref: Calculus of Variations and Partial Differential Equations (2021)

A constrained Cosserat shell model up to order $O(h^5)$: Modelling, existence of minimizers, relations to classical shell models and scaling invariance of the bending tensor

Authors: Ionel-Dumitrel Ghiba, Mircea Bîrsan, Peter Lewintan, Patrizio Neff

Abstract: We consider a recently introduced geometrically nonlinear elastic Cosserat shell model incorporating effects up to order $O(h^5)$ in the shell thickness $h$. We develop the corresponding geometrically nonlinear constrained Cosserat shell model, we show the existence of minimizers for the $O(h^5)$ and $O(h^3)$ case and we draw some connections to existing models and classical shell strain measures.… ▽ More We consider a recently introduced geometrically nonlinear elastic Cosserat shell model incorporating effects up to order $O(h^5)$ in the shell thickness $h$. We develop the corresponding geometrically nonlinear constrained Cosserat shell model, we show the existence of minimizers for the $O(h^5)$ and $O(h^3)$ case and we draw some connections to existing models and classical shell strain measures. Notably, the role of the appearing new bending tensor is highlighted and investigated with respect to an invariance condition of Acharya [Int. J. Solids and Struct., 2000] which will be further strengthened. △ Less

Submitted 28 October, 2020; v1 submitted 27 October, 2020; originally announced October 2020.

Journal ref: Journal of Elasticity (2021)

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

Authors: Adam Sky, Michael Neunteufel, Ingo Münch, Joachim Schoeberl, Patrizio Neff

Abstract: 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, such that standard nodal H1-finite elements yield unsat… ▽ More 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, such that standard nodal H1-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 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 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. △ Less

Submitted 2 March, 2021; v1 submitted 28 September, 2020; originally announced September 2020.

A rank-one convex, non-polyconvex isotropic function on $\operatorname{GL}^+(2)$ with compact connected sublevel sets

Authors: Jendrik Voss, Ionel-Dumitrel Ghiba, Robert J. Martin, Patrizio Neff

Abstract: According to a 2002 theorem by Cardaliaguet and Tahraoui, an isotropic, compact and connected subset of the group $\operatorname{GL}^+(2)$ of invertible $2\times2-\,$matrices is rank-one convex if and only if it is polyconvex. In a 2005 Journal of Convex Analysis article by Alexander~Mielke, it has been conjectured that the equivalence of rank-one convexity and polyconvexity holds for isotropic fu… ▽ More According to a 2002 theorem by Cardaliaguet and Tahraoui, an isotropic, compact and connected subset of the group $\operatorname{GL}^+(2)$ of invertible $2\times2-\,$matrices is rank-one convex if and only if it is polyconvex. In a 2005 Journal of Convex Analysis article by Alexander~Mielke, it has been conjectured that the equivalence of rank-one convexity and polyconvexity holds for isotropic functions on $\operatorname{GL}^+(2)$ as well, provided their sublevel sets satisfy the corresponding requirements. We negatively answer this conjecture by giving an explicit example of a function $W:\operatorname{GL}^+\to\mathbb{R}$ which is not polyconvex, but rank-one convex as well as isotropic with compact and connected sublevel sets. △ Less

Submitted 21 September, 2020; originally announced September 2020.

MSC Class: 26B25; 26A51; 74B20

Rank-one convexity vs. ellipticity for isotropic functions

Authors: Robert J. Martin, Jendrik Voss, Ionel-Dumitrel Ghiba, Patrizio Neff

Abstract: It is well known that a twice-differentiable real-valued function $W:\operatorname{GL}^+(n)\rightarrow\mathbb{R}$ on the group $\operatorname{GL}^+(n)$ of invertible $n\times n-$matrices with positive determinant is rank-one convex if and only if it is Legendre-Hadamard elliptic. Many energy functions arising from interesting applications in isotropic nonlinear elasticity, however, are not necessa… ▽ More It is well known that a twice-differentiable real-valued function $W:\operatorname{GL}^+(n)\rightarrow\mathbb{R}$ on the group $\operatorname{GL}^+(n)$ of invertible $n\times n-$matrices with positive determinant is rank-one convex if and only if it is Legendre-Hadamard elliptic. Many energy functions arising from interesting applications in isotropic nonlinear elasticity, however, are not necessarily twice differentiable everywhere on $\operatorname{GL}^+(n)$, especially at points with non-simple singular values. Here, we show that if an isotropic function $W$ on $\operatorname{GL}^+(n)$ is twice differentiable at each $F\in\operatorname{GL}^+(n)$ with simple singular values and Legendre-Hadamard elliptic at each such $F$, then $W$ is already rank-one convex under strongly reduced regularity assumptions. In particular, this generalization makes (local) ellipticity criteria accessible as criteria for (global) rank-one convexity to a wider class of elastic energy potentials expressed in terms of ordered singular values. Our results are also directly applicable to so-called conformally invariant energy functions. We also discuss a classical ellipticity criterion for the planar case by Knowles and Sternberg which has often been used in the literature as a criterion for global rank-one convexity and show that for this purpose, it is still applicable under weakened regularity assumptions. △ Less

Submitted 26 August, 2020; originally announced August 2020.

MSC Class: 74B20; 26B25

Sharp rank-one convexity conditions in planar isotropic elasticity for the additive volumetric-isochoric split

Authors: Jendrik Voss, Ionel-Dumitrel Ghiba, Robert J. Martin, Patrizio Neff

Abstract: We consider the volumetric-isochoric split in planar isotropic hyperelasticity and give a precise analysis of rank-one convexity criteria for this case, showing that the Legendre-Hadamard ellipticity condition separates and simplifies in a suitable sense. Starting from the classical two-dimensional criterion by Knowles and Sternberg, we can reduce the conditions for rank-one convexity to a family… ▽ More We consider the volumetric-isochoric split in planar isotropic hyperelasticity and give a precise analysis of rank-one convexity criteria for this case, showing that the Legendre-Hadamard ellipticity condition separates and simplifies in a suitable sense. Starting from the classical two-dimensional criterion by Knowles and Sternberg, we can reduce the conditions for rank-one convexity to a family of one-dimensional coupled differential inequalities. In particular, this allows us to derive a simple rank-one convexity classification for generalized Hadamard energies of the type $W(F)=\fracμ{2}\frac{\lVert F\rVert^2}{\det F}+f(\det F)$; such an energy is rank-one convex if and only if the function $f$ is convex. △ Less

Submitted 11 August, 2020; v1 submitted 10 August, 2020; originally announced August 2020.

MSC Class: 26B25; 26A51; 74A10; 74B20

Exploring metamaterials' structures through the relaxed micromorphic model: switching an acoustic screen into an acoustic absorber

Authors: Gianluca Rizzi, Manuel Collet, Félix Demore, Bernhard Eidel, Patrizio Neff, Angela Madeo

Abstract: While the design of always new metamaterials with exotic static and dynamic properties is attracting deep attention in the last decades, little effort is made to explore their interactions with other materials. This prevents the conception of (meta-)structures that can enhance metamaterials' unorthodox behaviours and that can be employed in real engineering applications. In this paper, we give a f… ▽ More While the design of always new metamaterials with exotic static and dynamic properties is attracting deep attention in the last decades, little effort is made to explore their interactions with other materials. This prevents the conception of (meta-)structures that can enhance metamaterials' unorthodox behaviours and that can be employed in real engineering applications. In this paper, we give a first answer to this challenging problem by showing that the relaxed micromorphic model with zero static characteristic length can be usefully applied to describe the refractive properties of simple meta-structures for extended frequency ranges and for any direction of propagation of the incident wave. Thanks to the simplified model's structure, we are able to efficiently explore different configurations and to show that a given meta-structure can drastically change its overall refractive behaviour when varying the elastic properties of specific meta-structural elements. In some cases, changing the stiffness of a homogeneous material which is in contact with a metamaterial's slab, reverses the structure's refractive behaviour by switching it from an acoustic screen (total reflection) into an acoustic absorber (total transmission). The present paper clearly indicates that, while the study and enhancement of the intrinsic metamaterials' properties is certainly of great importance, it is even more challenging to enable the conception of meta-structures that can eventually boost the use of metamaterials in real-case applications. △ Less

Submitted 29 July, 2020; originally announced July 2020.

Comments: 21 pages, 14 figures

A note on local higher regularity in the dynamic linear relaxed micromorphic model

Authors: Sebastian Owczarek, Ionel-Dumitrel Ghiba, Patrizio Neff

Abstract: We consider the regularity question of solutions for the dynamic initial-boundary value problem for the linear relaxed micromorphic model. This generalized continuum model couples a wave-type equation for the displacement with a generalized Maxwell-type wave equation for the micro-distortion. Naturally solutions are found in ${\rm H}^1$ for the displacement $u$ and ${\rm H}({\rm Curl})$ for the mi… ▽ More We consider the regularity question of solutions for the dynamic initial-boundary value problem for the linear relaxed micromorphic model. This generalized continuum model couples a wave-type equation for the displacement with a generalized Maxwell-type wave equation for the micro-distortion. Naturally solutions are found in ${\rm H}^1$ for the displacement $u$ and ${\rm H}({\rm Curl})$ for the microdistortion $P$. Using energy estimates for difference quotients, we improve this regularity. We show ${\rm H}^1_{\rm loc}$-regularity for the displacement field, ${\rm H}^1_{\rm loc}$-regularity for the micro-distortion tensor $P$ and that ${\rm Curl}\,P$ is ${\rm H}^1$-regular if the data is sufficiently smooth. △ Less

Submitted 14 May, 2021; v1 submitted 9 June, 2020; originally announced June 2020.