-
The Biot stress -- right stretch relation for the compressible Neo-Hooke-Ciarlet-Geymonat model and Rivlin's cube problem
Authors:
Ionel-Dumitrel Ghiba,
Franz Gmeineder,
Sebastian Holthausen,
Robert J. Martin,
Patrizio Neff
Abstract:
The aim of the paper is to recall the importance of the study of invertibility and monotonicity of stress-strain relations for investigating the non-uniqueness and bifurcation of homogeneous solutions of the equilibrium problem of a hyperelastic cube subjected to equiaxial tensile forces. In other words, we reconsider a remarkable possibility in this nonlinear scenario: Does symmetric loading lead…
▽ More
The aim of the paper is to recall the importance of the study of invertibility and monotonicity of stress-strain relations for investigating the non-uniqueness and bifurcation of homogeneous solutions of the equilibrium problem of a hyperelastic cube subjected to equiaxial tensile forces. In other words, we reconsider a remarkable possibility in this nonlinear scenario: Does symmetric loading lead only to symmetric deformations or also to asymmetric deformations? If so, what can we say about monotonicity for these homogeneous solutions, a property which is less restrictive than the energetic stability criteria of homogeneous solutions for Rivlin's cube problem. For the Neo-Hooke type materials we establish what properties the volumetric function $h$ depending on ${\rm det}\, F$ must have to ensure the existence of a unique radial solution (i.e. the cube must continue to remain a cube) for any magnitude of radial stress acting on the cube. The function $h$ proposed by Ciarlet and Geymonat satisfies these conditions. However, discontinuous equilibrium trajectories may occur, characterized by abruptly appearing non-symmetric deformations with increasing load, and a cube can instantaneously become a parallelepiped. Up to the load value for which the bifurcation in the radial solution is realized local monotonicity holds true. However, after exceeding this value, monotonicity no longer occurs on homogeneous deformations which, in turn, preserve the cube shape.
△ Less
Submitted 7 August, 2024;
originally announced August 2024.
-
Intrinsic mixed-dimensional beam-shell-solid couplings in linear Cosserat continua via tangential differential calculus
Authors:
Adam Sky,
Jack S. Hale,
Andreas Zilian,
Stéphane P. A. Bordas,
Patrizio Neff
Abstract:
We present an approach to the coupling of mixed-dimensional continua by employing the mathematically enriched linear Cosserat micropolar model. The kinematical reduction of the model to lower dimensional domains leaves its fundamental degrees of freedom intact. Consequently, the degrees of freedom intrinsically agree even at the interface with a domain of a different dimensionality. Thus, this app…
▽ More
We present an approach to the coupling of mixed-dimensional continua by employing the mathematically enriched linear Cosserat micropolar model. The kinematical reduction of the model to lower dimensional domains leaves its fundamental degrees of freedom intact. Consequently, the degrees of freedom intrinsically agree even at the interface with a domain of a different dimensionality. Thus, this approach circumvents the need for intermediate finite elements or mortar methods. We introduce the derivations of all models of various dimensions using tangential differential calculus. The coupling itself is then achieved by defining a mixed-dimensional action functional with consistent Sobolev trace operators. Finally, we present numerical examples involving a three-dimensional silicone-rubber block reinforced with a curved graphite shell on its lower surface, a three-dimensional silver block reinforced with a graphite plate and beams, and lastly, intersecting silver shells reinforced with graphite beams.
△ Less
Submitted 17 July, 2024;
originally announced July 2024.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
The Legendre-Hadamard condition in Cosserat elasticity theory
Authors:
Milad Shirani,
David J. Steigmann,
Patrizio Neff
Abstract:
The Legendre-Hadamard necessary condition for energy minimizers is derived in the framework of Cosserat elasticity theory.
The Legendre-Hadamard necessary condition for energy minimizers is derived in the framework of Cosserat elasticity theory.
△ Less
Submitted 21 September, 2020; v1 submitted 2 June, 2020;
originally announced June 2020.
-
$L^p$-trace-free generalized Korn inequalities for incompatible tensor fields in three space dimensions
Authors:
Peter Lewintan,
Patrizio Neff
Abstract:
For $1<p<\infty$ we prove an $L^p$-version of the generalized trace-free Korn inequality for incompatible tensor fields $P$ in $ W^{1,\,p}_0(\operatorname{Curl}; Ω,\mathbb{R}^{3\times3})$. More precisely, let $Ω\subset\mathbb{R}^3$ be a bounded Lipschitz domain. Then there exists a constant $c>0$ such that \[ \|{ P }\|_{L^p(Ω,\mathbb{R}^{3\times3})}\leq c\,\left(\|{\operatorname{dev} \operatorname…
▽ More
For $1<p<\infty$ we prove an $L^p$-version of the generalized trace-free Korn inequality for incompatible tensor fields $P$ in $ W^{1,\,p}_0(\operatorname{Curl}; Ω,\mathbb{R}^{3\times3})$. More precisely, let $Ω\subset\mathbb{R}^3$ be a bounded Lipschitz domain. Then there exists a constant $c>0$ such that \[ \|{ P }\|_{L^p(Ω,\mathbb{R}^{3\times3})}\leq c\,\left(\|{\operatorname{dev} \operatorname{sym} P }\|_{L^p(Ω,\mathbb{R}^{3\times3})} + \|{ \operatorname{dev} \operatorname{Curl} P }\|_{L^p(Ω,\mathbb{R}^{3\times3})}\right) \] holds for all tensor fields $P\in W^{1,\,p}_0(\operatorname{Curl}; Ω,\mathbb{R}^{3\times3})$, i.e., for all $P\in W^{1,\,p}(\operatorname{Curl}; Ω,\mathbb{R}^{3\times3})$ with vanishing tangential trace $ P\times ν=0 $ on $ \partialΩ$ where $ν$ denotes the outward unit normal vector field to $\partialΩ$ and $\operatorname{dev} P := P -\frac13 \operatorname{tr}(P)\,\mathbb{1}_3$ denotes the deviatoric (trace-free) part of $P$. We also show the norm equivalence \[ \|{ P }\|_{L^p(Ω,\mathbb{R}^{3\times3})}+\|{\operatorname{Curl} P }\|_{L^p(Ω,\mathbb{R}^{3\times3})}\leq c\,\left(\|{\operatorname{dev} \operatorname{sym} P }\|_{L^p(Ω,\mathbb{R}^{3\times3})} + \|{ \operatorname{dev}\operatorname{Curl} P }\|_{L^p(Ω,\mathbb{R}^{3\times3})}\right) \] for tensor fields $P\in W^{1,\,p}_0(\operatorname{Curl}; Ω,\mathbb{R}^{3\times3})$. These estimates also hold true for tensor fields with vanishing tangential trace only on a relatively open (non-empty) subset $Γ\subseteq \partialΩ$ of the boundary.
△ Less
Submitted 13 October, 2021; v1 submitted 10 April, 2020;
originally announced April 2020.
-
The isotropic Cosserat shell model including terms up to $O(h^5)$. Part II: Existence of minimizers
Authors:
Ionel-Dumitrel Ghiba,
Mircea Bîrsan,
Peter Lewintan,
Patrizio Neff
Abstract:
We show the existence of global minimizers for a geometrically nonlinear isotropic elastic Cosserat 6-parameter shell model. The proof of the main theorem is based on the direct methods of the calculus of variations using essentially the convexity of the energy in the nonlinear strain and curvature measures. We first show the existence of the solution for the theory including $O(h^5)$ terms. The e…
▽ More
We show the existence of global minimizers for a geometrically nonlinear isotropic elastic Cosserat 6-parameter shell model. The proof of the main theorem is based on the direct methods of the calculus of variations using essentially the convexity of the energy in the nonlinear strain and curvature measures. We first show the existence of the solution for the theory including $O(h^5)$ terms. The energy allows us to show the coercivity for terms up to order $O(h^5)$ and the convexity of the energy. Secondly, we consider only that part of the energy including $O(h^3)$ terms. In this case the obtained minimization problem is not the same as those previously considered in the literature, since the influence of the curved initial shell configuration appears explicitly in the expression of the coefficients of the energies for the reduced two-dimensional variational problem and additional mixed bending-curvature and curvature terms are present. While in the theory including $O(h^5)$ the conditions on the thickness $h$ are those considered in the modelling process and they are independent of the constitutive parameter, in the $O(h^3)$-case the coercivity is proven on some more restrictive conditions under the thickness $h$.
△ Less
Submitted 15 September, 2020; v1 submitted 19 March, 2020;
originally announced March 2020.
-
The isotropic Cosserat shell model including terms up to $O(h^5)$. Part I: Derivation in matrix notation
Authors:
Ionel-Dumitrel Ghiba,
Mircea Bîrsan,
Peter Lewintan,
Patrizio Neff
Abstract:
We present a new geometrically nonlinear Cosserat shell model incorporating effects up to order $O(h^5)$ in the shell thickness $h$. The method that we follow is an educated 8-parameter ansatz for the three-dimensional elastic shell deformation with attendant analytical thickness integration, which leads us to obtain completely two-dimensional sets of equations in variational form. We give an expl…
▽ More
We present a new geometrically nonlinear Cosserat shell model incorporating effects up to order $O(h^5)$ in the shell thickness $h$. The method that we follow is an educated 8-parameter ansatz for the three-dimensional elastic shell deformation with attendant analytical thickness integration, which leads us to obtain completely two-dimensional sets of equations in variational form. We give an explicit form of the curvature energy using the orthogonal Cartan-decomposition of the wryness tensor. Moreover, we consider the matrix representation of all tensors in the derivation of the variational formulation, because this is convenient when the problem of existence is considered, and it is also preferential for numerical simulations. The step by step construction allows us to give a transparent approximation of the three-dimensional parental problem. The resulting 6-parameter isotropic shell model combines membrane, bending and curvature effects at the same time. The Cosserat shell model naturally includes a frame of orthogonal directors, the last of which does not necessarily coincide with the normal of the surface. This rotation-field is coupled to the shell-deformation and augments the well-known Reissner-Mindlin kinematics (one independent director) with so-called in-plane drill rotations, the inclusion of which is a decisive for subsequent numerical treatment and existence proofs. As a major novelty, we determine the constitutive coefficients of the Cosserat shell model in dependence on the geometry of the shell which are otherwise difficult to guess.
△ Less
Submitted 14 September, 2020; v1 submitted 1 March, 2020;
originally announced March 2020.
-
$L^p$-versions of generalized Korn inequalities for incompatible tensor fields in arbitrary dimensions with $p$-integrable exterior derivative
Authors:
Peter Lewintan,
Patrizio Neff
Abstract:
For $n\ge2$ and $1<p<\infty$ we prove an $L^p$-version of the generalized Korn-type inequality for incompatible, $p$-integrable tensor fields $P:Ω\to \mathbb{R}^{n\times n}$ having $p$-integrable generalized $\underline{\operatorname{Curl}}$ and generalized vanishing tangential trace $P\,τ_l=0$ on $\partial Ω$, denoting by $\{τ_l\}_{l=1,\ldots, n-1}$ a moving tangent frame on $\partialΩ$, more pre…
▽ More
For $n\ge2$ and $1<p<\infty$ we prove an $L^p$-version of the generalized Korn-type inequality for incompatible, $p$-integrable tensor fields $P:Ω\to \mathbb{R}^{n\times n}$ having $p$-integrable generalized $\underline{\operatorname{Curl}}$ and generalized vanishing tangential trace $P\,τ_l=0$ on $\partial Ω$, denoting by $\{τ_l\}_{l=1,\ldots, n-1}$ a moving tangent frame on $\partialΩ$, more precisely we have: $$\| P \|_{L^p(Ω,\mathbb{R}^{n\times n})}\leq c\,(\| \operatorname{sym} P\|_{L^p(Ω,\mathbb{R}^{n \times n})} + \|\underline{\operatorname{Curl}} P \|_{L^p(Ω,(\mathfrak{so}(n))^n)} ),$$ where the generalized $\underline{\operatorname{Curl}}$ is given by $ (\underline{\operatorname{Curl}})_{ijk} :=\partial_i P_{kj}-\partial_j P_{ki}$ and $c=c(n,p,Ω)>0$.
△ Less
Submitted 6 September, 2021; v1 submitted 18 December, 2019;
originally announced December 2019.
-
Nečas-Lions lemma revisited: An $L^p$-version of the generalized Korn inequality for incompatible tensor fields
Authors:
Peter Lewintan,
Patrizio Neff
Abstract:
For $1<p<\infty$ we prove an $L^p$-version of the generalized Korn inequality for incompatible tensor fields $P$ in $ W^{1,\,p}_0(\operatorname{Curl}; Ω,\mathbb{R}^{3\times3})$. More precisely, let $Ω\subset\mathbb{R}^3$ be a bounded Lipschitz domain. Then there exists a constant $c>0$ such that \begin{equation*} \| P\|_{L^p(Ω,\mathbb{R}^{3\times3})}\leq c\,\left( \|\operatorname{sym} P\|_{L^p(Ω,\…
▽ More
For $1<p<\infty$ we prove an $L^p$-version of the generalized Korn inequality for incompatible tensor fields $P$ in $ W^{1,\,p}_0(\operatorname{Curl}; Ω,\mathbb{R}^{3\times3})$. More precisely, let $Ω\subset\mathbb{R}^3$ be a bounded Lipschitz domain. Then there exists a constant $c>0$ such that \begin{equation*} \| P\|_{L^p(Ω,\mathbb{R}^{3\times3})}\leq c\,\left( \|\operatorname{sym} P\|_{L^p(Ω,\mathbb{R}^{3\times3})} + \| \operatorname{Curl}P \|_{L^p(Ω, \mathbb{R}^{3\times3})}\right)\end{equation*} holds for all tensor fields $P\in W^{1,\,p}_0(\operatorname{Curl}; Ω,\mathbb{R}^{3\times3})$, i.e., for all $P\in W^{1,\,p}(\operatorname{Curl}; Ω,\mathbb{R}^{3\times3})$ with vanishing tangential trace $ P\times ν=0 $ on $ \partialΩ$ where $ν$ denotes the outward unit normal vector field to $\partialΩ$. For compatible $P=D u$ this recovers an $L^p$-version of the classical Korn's first inequality $$ \|D u \|_{L^p(Ω,\mathbb{R}^{3\times 3})} \le c\, \|\operatorname{sym}D u\|_{L^p(Ω,\mathbb{R}^{3\times3})} \quad \text{with }D u \times ν= 0 \quad \text{on $\partial Ω$}, $$ and for skew-symmetric $P=A\in\mathfrak{so}(3)$ an $L^p$-version of the Poincaré inequality $$ \|A\|_{L^p(Ω,\mathfrak{so}(3))}\le c\, \|\operatorname{Curl} A\|_{L^p(Ω,\mathbb{R}^{3\times3})} \quad \text{with } A \times ν= 0 \ \Leftrightarrow \ A=0 \quad \text{on $\partial Ω$}. $$
△ Less
Submitted 15 May, 2021; v1 submitted 18 December, 2019;
originally announced December 2019.
-
Existence results for non-homogeneous boundary conditions in the relaxed micromorphic model
Authors:
Ionel-Dumitrel Ghiba,
Patrizio Neff,
Sebastian Owczarek
Abstract:
In this paper we use a property of the extension operator from the space of tangential traces of ${\rm H}({\rm curl};Ω)$ in the context of the linear relaxed micromorphic model, a theory which is recently used to describe the behaviour of some metamaterials showing unorthodox behaviors with respect to elastic wave propagation. We show that the new property is important for existence results of str…
▽ More
In this paper we use a property of the extension operator from the space of tangential traces of ${\rm H}({\rm curl};Ω)$ in the context of the linear relaxed micromorphic model, a theory which is recently used to describe the behaviour of some metamaterials showing unorthodox behaviors with respect to elastic wave propagation. We show that the new property is important for existence results of strong solution for non-homogeneous boundary condition in both the dynamic and the static case.
△ Less
Submitted 31 October, 2019;
originally announced October 2019.
-
Refined dimensional reduction for isotropic elastic Cosserat shells with initial curvature
Authors:
Mircea Birsan,
Ionel-Dumitrel Ghiba,
Robert J. Martin,
Patrizio Neff
Abstract:
Using a geometrically motivated 8-parameter ansatz through the thickness, we reduce a three-dimensional shell-like geometrically nonlinear Cosserat material to a fully two-dimensional shell model. Curvature effects are fully taken into account. For elastic isotropic Cosserat materials, the integration through the thickness can be performed analytically and a generalized plane stress condition allo…
▽ More
Using a geometrically motivated 8-parameter ansatz through the thickness, we reduce a three-dimensional shell-like geometrically nonlinear Cosserat material to a fully two-dimensional shell model. Curvature effects are fully taken into account. For elastic isotropic Cosserat materials, the integration through the thickness can be performed analytically and a generalized plane stress condition allows for a closed-form expression of the thickness stretch and the nonsymmetric shift of the midsurface in bending. We obtain an explicit form of the elastic strain energy density for Cosserat shells, including terms up to order $O(h^5)$ in the shell thickness $h$. This energy density is expressed as a quadratic function of the nonlinear elastic shell strain tensor and the bending-curvature tensor, with coefficients depending on the initial curvature of the shell.
△ Less
Submitted 27 September, 2019;
originally announced September 2019.
-
The axiomatic introduction of arbitrary strain tensors by Hans Richter -- a commented translation of "Strain tensor, strain deviator and stress tensor for finite deformations"
Authors:
Patrizio Neff,
Kai Graban,
Eva Schweickert,
Robert J. Martin
Abstract:
We provide a faithful translation of Hans Richter's important 1949 paper "Verzerrungstensor, Verzerrungsdeviator und Spannungstensor bei endlichen Formänderungen" from its original German version into English, complemented by an introduction summarizing Richter's achievements.
We provide a faithful translation of Hans Richter's important 1949 paper "Verzerrungstensor, Verzerrungsdeviator und Spannungstensor bei endlichen Formänderungen" from its original German version into English, complemented by an introduction summarizing Richter's achievements.
△ Less
Submitted 12 September, 2019;
originally announced September 2019.