Source author record

Patrizio Neff

Patrizio Neff appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

83works
13topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

83 published item(s)

preprint2023arXiv

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

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.

preprint2022arXiv

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

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.

preprint2022arXiv

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

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.

preprint2022arXiv

Cosserat micropolar elasticity: classical Eringen vs. dislocation form

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.

preprint2022arXiv

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

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.

preprint2022arXiv

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

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.

preprint2022arXiv

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

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.

preprint2022arXiv

Primal and mixed finite element formulations for the relaxed micromorphic model

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.

preprint2022arXiv

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

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.

preprint2021arXiv

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

One approach for the simulation of metamaterials is to extend an associated continuum theory concerning its kinematic equations, and the relaxed micromorphic continuum represents such a model. It incorporates the Curl of the nonsymmetric microdistortion in the free energy function. This suggests the existence of solutions not belonging to H1, 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.

preprint2021arXiv

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

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.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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.

preprint2020arXiv

Rank-one convexity vs. ellipticity for isotropic functions

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.

preprint2020arXiv

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

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.

preprint2020arXiv

The isotropic Cosserat shell model including terms up to $O(h^5)$. Part I: Derivation in matrix notation

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.

preprint2020arXiv

The isotropic Cosserat shell model including terms up to $O(h^5)$. Part II: Existence of minimizers

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$.

preprint2019arXiv

A note on non-homogeneous deformations with homogeneous Cauchy stress for a strictly rank-one convex energy in isotropic hyperelasticity

It has recently been shown that for a Cauchy stress response induced by a strictly rank-one convex hyperelastic energy potential, a homogeneous Cauchy stress tensor field cannot correspond to a non-homogeneous deformation if the deformation gradient has discrete values, i.e. if the deformation is piecewise affine linear and satisfies the Hadamard jump condition. In this note, we expand upon these results and show that they do not hold for arbitrary deformations by explicitly giving an example of a strictly rank-one convex energy and a non-homogeneous deformation such that the induced Cauchy stress tensor is constant. In the planar case, our example is related to another previous result concerning criteria for generalized convexity properties of conformally invariant energy functions, which we extend to the case of strict rank-one convexity.

preprint2019arXiv

Chirality in the plane

It is well-known that many three-dimensional chiral material models become non-chiral when reduced to two dimensions. Chiral properties of the two-dimensional model can then be restored by adding appropriate two-dimensional chiral terms. In this paper we show how to construct a three-dimensional chiral energy function which can achieve two-dimensional chirality induced already by a chiral three-dimensional model. The key ingredient to this approach is the consideration of a nonlinear chiral energy containing only rotational parts. After formulating an appropriate energy functional, we study the equations of motion and find explicit soliton solutions displaying two-dimensional chiral properties.

preprint2016arXiv

A new view on boundary conditions in the Grioli-Koiter-Mindlin-Toupin indeterminate couple stress model

In this paper we consider the Grioli-Koiter-Mindlin-Toupin linear isotropic indeterminate couple stress model. Our main aim is to show that, up to now, the boundary conditions have not been completely understood for this model. As it turns out, and to our own surprise, restricting the well known boundary conditions stemming from the strain gradient or second gradient models to the particular case of the indeterminate couple stress model, does not always reduce to the Grioli-Koiter-Mindlin-Toupin set of accepted boundary conditions. We present, therefore, a proof of the fact that when specific "mixed" kinematical and traction boundary conditions are assigned on the boundary, no "a priori" equivalence can be established between Mindlin's and our approach.

preprint2016arXiv

A panorama of dispersion curves for the weighted isotropic relaxed micromorphic model

We consider the weighted isotropic relaxed micromorphic model and provide an in depth investigation of the characteristic dispersion curves when the constitutive parameters of the model are varied. The weighted relaxed micromorphic model generalizes the classical relaxed micromorphic model previously introduced by the authors, since it features the Cartan-Lie decomposition of the tensors $P_{,t}$ and Curl $P$ in their dev, dev sym, skew and spheric part. It is shown that the split of the tensor $P_{,t}$ in the micro-inertia provide an independent control of the cut-offs of the optic benches. This is crucial for the future calibration of the relaxed micromorphic model on real band-gap metamaterials. Even if the physical interest of the introduction of the split of the tensor Curl $P$ is less evident than in the previous case, we discuss in detail which is its effect on the dispersion curves. Finally, we also provide a complete parametric study involving all the constitutive parameters of the introduced model, so giving rise to an exhaustive panorama of dispersion curves for the relaxed micromorphic model.

preprint2016arXiv

A review on wave propagation modeling in band-gap metamaterials via enriched continuum models

In the present contribution we show that the relaxed micromorphic model is the only non-local continuum model which is able to account for the description of band-gaps in metamaterials for which the kinetic energy accounts separately for micro and macro-motions without considering a micro-macro coupling. Moreover, we show that when adding a gradient inertia term which indeed allows for the description of the coupling of the vibrations of the microstructure to the macroscopic motion of the unit cell, other enriched continuum models of the micromorphic type may allow the description of the onset of band-gaps. Nevertheless, the relaxed micromorphic model proves to be yet the most effective enriched continuum model which is able to describe multiple band-gaps in non-local metamaterials.

preprint2016arXiv

Complete band gaps including non-local effects occur only in the relaxed micromorphic model

In this paper we substantiate the claim implicitly made in previous works that the relaxed micromorphic model is the only linear, isotropic, reversibly elastic, nonlocal generalized continuum model able to describe complete band-gaps on a phenomenological level. To this end, we recapitulate the response of the standard Mindlin-Eringen micromorphic model with the full micro-distortion gradient of P, the relaxed micromorphic model depending only on the Curl P of the micro-distortion P, and a variant of the standard micromorphic model in which the curvature depends only on the divergence Div P of the micro distortion. The Div-model has size-effects but the dispersion analysis for plane waves shows the incapability of that model to even produce a partial band gap. Combining the curvature to depend quadratically on Div P and Curl P shows that such a model is similar to the standard Mindlin-Eringen model which can eventually show only a partial band gap.

preprint2016arXiv

Existence result for a dislocation based model of single crystal gradient plasticity with isotropic or linear kinematic hardening

We consider a dislocation-based rate-independent model of single crystal gradient plasticity with isotropic or linear kinematic hardening. The model is weakly formulated through the so-called primal form of the flow rule as a variational inequality for which a result of existence and uniqueness is obtained using the functional analytical framework developed by Han-Reddy.

preprint2016arXiv

First evidence of non-locality in real band-gap metamaterials: determining parameters in the relaxed micromorphic model

In this paper we propose the first estimate of some elastic parameters of the relaxed micromorphic model on the basis of real experiments of transmission of longitudinal plane waves across an interface separating a classical Cauchy material (steel plate) and a phononic crystal (steel plate with fluid-filled holes). A procedure is set up in order to identify the parameters of our model by superimposing the experimentally-based profile of the reflection coefficient (plotted as function of the frequency of the traveling waves) with the analogous profile obtained via simulations based upon the relaxed micromorphic model. We end up with the determination of 5 out of 6 constitutive parameters which are featured by the relaxed micromorphic model in the isotropic case, plus the determination of the micro-inertia parameter. The sixth elastic parameter, namely the Cosserat couple modulus $μ_{c}$, still remains undetermined, since experimental data concerning the transmission properties of the considered interface for transverse incident waves are not yet available. A fundamental result of the present paper is the estimate of the non-locality intrinsically associated to the underlying microstructure of the metamaterial. As a matter of fact, we appraise that the characteristic length $L_{c}$ measuring the non-locality of the considered phononic crystal is of the order of $1/3$ of the diameter of the considered fluid-filled holes.

preprint2016arXiv

Geometry of logarithmic strain measures in solid mechanics

We consider the two logarithmic strain measures\[ω_{\rm iso}=\|\mathrm{dev}_n\log U\|=\|\mathrm{dev}_n\log \sqrt{F^TF}\|\quad\text{ and }\quad ω_{\rm vol}=|\mathrm{tr}(\log U)|=|\mathrm{tr}(\log\sqrt{F^TF})|\,,\]which are isotropic invariants of the Hencky strain tensor $\log U$, and show that they can be uniquely characterized by purely geometric methods based on the geodesic distance on the general linear group $\mathrm{GL}(n)$. Here, $F$ is the deformation gradient, $U=\sqrt{F^TF}$ is the right Biot-stretch tensor, $\log$ denotes the principal matrix logarithm, $\|.\|$ is the Frobenius matrix norm, $\mathrm{tr}$ is the trace operator and $\mathrm{dev}_n X$ is the $n$-dimensional deviator of $X\in\mathbb{R}^{n\times n}$. This characterization identifies the Hencky (or true) strain tensor as the natural nonlinear extension of the linear (infinitesimal) strain tensor $\varepsilon=\mathrm{sym}\nabla u$, which is the symmetric part of the displacement gradient $\nabla u$, and reveals a close geometric relation between the classical quadratic isotropic energy potential \[μ\,\|\mathrm{dev}_n\mathrm{sym}\nabla u\|^2+\fracκ{2}\,[\mathrm{tr}(\mathrm{sym}\nabla u)]^2=μ\,\|\mathrm{dev}_n\varepsilon\|^2+\fracκ{2}\,[\mathrm{tr}(\varepsilon)]^2\]in linear elasticity and the geometrically nonlinear quadratic isotropic Hencky energy\[μ\,\|\mathrm{dev}_n\log U\|^2+\fracκ{2}\,[\mathrm{tr}(\log U)]^2=μ\,ω_{\rm iso}^2+\frac\kappa2\,ω_{\rm vol}^2\,,\]where $μ$ is the shear modulus and $κ$ denotes the bulk modulus. Our deduction involves a new fundamental logarithmic minimization property of the orthogonal polar factor $R$, where $F=R\,U$ is the polar decomposition of $F$. We also contrast our approach with prior attempts to establish the logarithmic Hencky strain tensor directly as the preferred strain tensor in nonlinear isotropic elasticity.

preprint2016arXiv

Minimal geodesics on GL(n) for left-invariant, right-O(n)-invariant Riemannian metrics

We provide an easy approach to the geodesic distance on the general linear group GL(n) for left-invariant Riemannian metrics which are also right-O(n)-invariant. The parametrization of geodesic curves and the global existence of length minimizing geodesics are deduced using simple methods based on the calculus of variations and classical analysis only. The geodesic distance is discussed for some special cases and applications towards the theory of nonlinear elasticity are indicated.

preprint2016arXiv

Modeling real phononic crystals via the weighted relaxed micromorphic model with free and gradient micro-inertia

In this paper the relaxed micromorphic continuum model with weighted free and gradient micro-inertia is used to describe the dynamical behavior of a real two-dimensional phononic crystal for a wide range of wavelengths. In particular, a periodic structure with specific micro-structural topology and mechanical properties, capable of opening a phononic band-gap, is chosen with the criterion of showing a low degree of anisotropy (the band-gap is almost independent of the direction of propagation of the traveling wave). A Bloch wave analysis is performed to obtain the dispersion curves and the corresponding vibrational modes of the periodic structure. A linear-elastic, isotropic, relaxed micromorphic model including both a free micro-inertia (related to free vibrations of the microstructures) and a gradient micro-inertia (related to the motions of the microstructure which are coupled to the macro-deformation of the unit cell) is introduced and particularized to the case of plane wave propagation. The parameters of the relaxed model, which are independent of frequency, are then calibrated on the dispersion curves of the phononic crystal showing an excellent agreement in terms of both dispersion curves and vibrational modes. Almost all the homogenized elastic parameters of the relaxed micromorphic model result to be determined. This opens the way to the design of morphologically complex meta-structures which make use of the chosen phononic structure as the basic building block and which preserve its ability of "stopping" elastic wave propagation at the scale of the structure.

preprint2016arXiv

On the dislocation density tensor in the Cosserat theory of elastic shells

We consider the Cosserat continuum in its finite strain setting and discuss the dislocation density tensor as a possible alternative curvature strain measure in three-dimensional Cosserat models and in Cosserat shell models. We establish a close relationship (one-to-one correspondence) between the new shell dislocation density tensor and the bending-curvature tensor of 6-parameter shells.

preprint2016arXiv

Reflection and transmission of elastic waves at interfaces embedded in non-local band-gap metamaterials: a comprehensive study via the relaxed micromorphic model

In this paper we derive, by means of a suitable least action principle, the duality jump conditions to be imposed at surfaces of discontinuity of the material properties in non-dissipative, linear-elastic, isotropic, Mindlin's and relaxed micromorphic media, respectively. The introduced theoretical framework allows the transparent set-up of different types of micro-macro connections which are intrinsically compatible with the governing bulk equations. To illustrate the interest of the many introduced jump conditions, we focus on the case of an interface between a classical Cauchy continuum on one side and a relaxed micromorphic one on the other side. As expected, we find a complete reflection in the frequency intervals for which band-gaps are known to occur in the relaxed micromorphic continuum and precise microstructure-related reflective patterns are identified. We repeat a similar study for analogous connections between a classical Cauchy continuum and a Mindlin's micromorphic one and we show that the reflective properties of the considered interfaces are drastically modified due to the fact that band-gaps are not allowed in standard Mindlin's micromorphic media. The present work opens the way towards the possibility of conceiving complex metastructures in which band-gap metamaterials and classical materials are coupled together to produce structures with completely new and unorthodox properties with respect to wave propagation, transmission and reflection. Last, but not least, indirect measurements of the material coefficients of the relaxed micromorphic model based upon real experiments of reflection and transmission in band-gap metamaterials are uncovered by the present work which makes them finally realizable in the short term.

preprint2016arXiv

Rotational invariance conditions in elasticity, gradient elasticity and its connection to isotropy

For homogeneous higher gradient elasticity models we discuss frame-indifference and isotropy requirements. To this end, we introduce the notions of local versus global SO(3)-invariance and identify frame-indifference (traditionally) with global left SO(3)-invariance and isotropy with global right SO(3)-invariance. For specific restricted representations, the energy may also be local left SO(3)-invariant as well as local right SO(3)-invariant. Then we turn to linear models and consider a consequence of frame-indifference together with isotropy in nonlinear elasticity and apply this joint invariance condition to some specific linear models. The interesting point is the appearance of finite rotations in transformations of a geometrically linear model. It is shown that when starting with a linear model defined already in the infinitesimal symmetric strain $\varepsilon = {\rm sym} \, {\rm Grad}[u]$, the new invariance condition is equivalent to isotropy of the linear formulation. Therefore, it may be used also in higher gradient elasticity models for a simple check of isotropy and for extensions to anisotropy. In this respect we consider in more detail variational formulations of the linear indeterminate couple stress model, a new variant of it with symmetric force stresses and general linear gradient elasticity.

preprint2016arXiv

The exponentiated Hencky strain energy in modelling tire derived material for moderately large deformations

This work presents a hyper-viscoelastic model based on the Hencky-logarithmic strain tensor to model the response of a Tire Derived Material (TDM) undergoing moderately large deformations. TDM is a composite made by cold forging a mix of rubber fibers and grains, obtained by grinding scrap tires, and polyurethane binder. The mechanical properties are highly influenced by the presence of voids associated with the granular composition and low tensile strength due to the weak connection at the grain-matrix interface. For these reasons, TDM use is restricted to applications concerning a limited range of deformations. Experimental tests show that a central feature of the response is connected to highly nonlinear behavior of the material under volumetric deformation which conventional hyperelastic models fail in predicting. The strain energy function presented here is a variant of the exponentiated Hencky strain energy proposed by Neff et al., which for moderate strains is as good as the quadratic Hencky model and in the large strain region improves several important features from a mathematical point of view. The proposed form of the exponentiated Hencky energy possesses a set of parameters uniquely determined in the infinitesimal strain regime and an orthogonal set of parameters to determine the nonlinear response. The hyperelastic model is additionally incorporated in a finite deformation viscoelasticity framework that accounts for the two main dissipation mechanisms in TDMs, one at the microscale level and one at the macroscale level. The model is capable of predicting different deformation modes in a certain range of frequency and amplitude with a unique set of parameters with most of them having a clear physical meaning. Moreover, by comparing the predictions from the proposed constitutive model with experimental data we conclude that the new constitutive model gives accurate prediction.

preprint2016arXiv

Well-posedness for the microcurl model in both single and polycrystal gradient plasticity

We consider the recently introduced microcurl model which is a variant of strain gradient plasticity in which the curl of the plastic distortion is coupled to an additional micromorphic-type field. For both single crystal and polycrystal cases, we formulate the model and show its well-posedness in the rate-independent case provided some local hardening (isotropic or linear kinematic) is taken into account. To this end, we use the functional analytical framework developed by Han-Reddy. We also compare the model to the relaxed micromorphic model as well as to a dislocation-based gradient plasticity model.

preprint2015arXiv

A variant of the linear isotropic indeterminate couple stress model with symmetric local force-stress, symmetric nonlocal force-stress, symmetric couple-stresses and complete traction boundary conditions

In this paper we venture a new look at the linear isotropic indeterminate couple stress model in the general framework of second gradient elasticity and we propose a new alternative formulation which obeys Cauchy-Boltzmann's axiom of the symmetry of the force stress tensor. For this model we prove the existence of solutions for the equilibrium problem. Relations with other gradient elastic theories and the possibility to switch from a {4th order} (gradient elastic) problem to a 2nd order micromorphic model are also discussed with a view of obtaining symmetric force-stress tensors. It is shown that the indeterminate couple stress model can be written entirely with symmetric force-stress and symmetric couple-stress. The difference of the alternative models rests in specifying traction boundary conditions of either rotational type or strain type. If rotational type boundary conditions are used in the partial integration, the classical anti-symmetric nonlocal force stress tensor formulation is obtained. Otherwise, the difference in both formulations is only a divergence--free second order stress field such that the field equations are the same, but the traction boundary conditions are different. For these results we employ a novel integrability condition, connecting the infinitesimal continuum rotation and the infinitesimal continuum strain. Moreover, we provide the complete, consistent traction boundary conditions for both models.

preprint2015arXiv

An ellipticity domain for the distortional Hencky-logarithmic strain energy

We describe ellipticity domains for the isochoric elastic energy $ F\mapsto \|{\rm dev}_n\log U\|^2=\bigg\|\log \frac{\sqrt{F^TF}}{(\det F)^{1/n}}\bigg\|^2 =\frac{1}{4}\,\bigg\|\log \frac{C}{({\rm det} C)^{1/n}}\bigg\|^2 $ for $n=2,3$, where $C=F^TF$ for $F\in {\rm GL}^+(n)$. Here, ${\rm dev}_n\log {U} =\log {U}-\frac{1}{n}\, {\rm tr}(\log {U})\cdot 1\!\!1$ is the deviatoric part of the logarithmic strain tensor $\log U$. For $n=2$ we identify the maximal ellipticity domain, while for $n=3$ we show that the energy is Legendre-Hadamard elliptic in the set $\mathcal{E}_3\bigg(W_{_{\rm H}}^{\rm iso}, {\rm LH}, U, \frac{2}{3}\bigg)\,:=\,\bigg\{U\in{\rm PSym}(3) \;\Big|\, \|{\rm dev}_3\log U\|^2\leq \frac{2}{3}\bigg\}$, which is similar to the von-Mises-Huber-Hencky maximum distortion strain energy criterion. Our results complement the characterization of ellipticity domains for the quadratic Hencky energy $ W_{_{\rm H}}(F)=μ\,\|{\rm dev}_3\log U\|^2+ \fracκ{2}\,[{\rm tr} (\log U)]^2 $, $U=\sqrt{F^TF}$ with $μ>0$ and $κ>\frac{2}{3}\, μ$, previously obtained by Bruhns et al.

preprint2015arXiv

Correct traction boundary conditions in the indeterminate couple stress model

In this paper we consider the Grioli-Koiter-Mindlin-Toupin indeterminatecouple stress model. The main aim is to show that the traction boundary conditions were not yet completely deduced. As it turns out, and to our own surprise, restricting the boundary condition framework from the strain gradient models to the couple stress model does not reduce to Mindlin's set of accepted boundary conditions. We present therefore, for the first time the complete, consistent set of traction boundary conditions.

preprint2015arXiv

Existence results in dislocation based rate-independent isotropic gradient plasticity with kinematical hardening and plastic spin: The case with symmetric local backstress

In this paper we use convex analysis and variational inequality methods to establish an existence result for a model of infinitesimal rate-independent gradient plasticity with kinematic hardening and plastic spin, in which the local backstress tensor remains symmetric. The model features a defect energy contribution which is quadratic in the dislocation density tensor Curl p, giving rise to nonlocal non-symmetric kinematic hardening. Use is made of a recently established Korn's type inequality for incompatible tensor fields. The solution space for the non-symmetric plastic distortion is naturally H(Curl) together with suitable tangential boundary conditions on the plastic distortion. Connections to other models are established as well.

preprint2015arXiv

Integrability conditions between the first and second Cosserat deformation tensor in geometrically nonlinear micropolar models and existence of minimizers

In this note we extend integrability conditions for the symmetric stretch tensor $U$ in the polar decomposition of the deformation gradient $\nablaφ=F=R\,U$ to the non-symmetric case. In doing so we recover integrability conditions for the first Cosserat deformation tensor. Let $F=\bar R\,\bar U$ with $\bar R:Ω\subset\mathbb{R}^3\longrightarrow\mathrm{SO}(3)$ and $\bar U:Ω\subset\mathbb{R}^3\longrightarrow \mathrm{GL}(3)$. Then $\mathfrak{K}:={\bar R}^T\mathrm{Grad}\,{\bar R}=\mathrm{Anti}\Big( \frac{1}{\mathrm{det} \bar U}\Big[\bar U(\mathrm{Curl} \bar U)^T-\frac{1}{2} \mathrm{tr}(\bar U(\mathrm{Curl} \bar U)^T) 1\!\!1 \Big]\bar U\Big),$ giving a connection between the first Cosserat deformation tensor $\bar U$ and the second Cosserat tensor ${\mathfrak{K}}$. (Here, Anti denotes an isomorphism between $\mathbb{R}^{3\times 3}$ and $\mathfrak{So}(3):=\{\,\mathfrak{A}\in\mathbb{R}^{3\times 3\times 3}\,|\,\mathfrak{A}.u\in\mathfrak{so}(3)\;\forall u\in \mathbb{R}^3\}$.) The formula shows that it is not possible to prescribe $\bar U$ and $\mathfrak{K}$ independent from each other. We also propose a new energy formulation of geometrically nonlinear Cosserat models which completely separate the effects of nonsymmetric straining and curvature. For very weak constitutive assumptions (no direct boundary condition on rotations, zero Cosserat couple modulus, quadratic curvature energy) we show existence of minimizers in Sobolev-spaces.

preprint2015arXiv

Loss of ellipticity for non-coaxial plastic deformations in additive logarithmic finite strain plasticity

In this paper we consider the additive logarithmic finite strain plasticity formulation from the view point of loss of ellipticity in elastic unloading. We prove that even if an elastic energy $F\mapsto W(F)=\hat{W}(\log U)$ defined in terms of logarithmic strain $\log U$, where $U=\sqrt{F^T\, F}$, is everywhere rank-one convex as a function of $F$, the new function $F\mapsto \widetilde{W}(F)=\hat{W}(\log U-\log U_p)$ need not remain rank-one convex at some given plastic stretch $U_p$ (viz. $E_p^{\log}:=\log U_p$). This is in complete contrast to multiplicative plasticity in which $F\mapsto W(F\, F_p^{-1})$ remains rank-one convex at every plastic distortion $F_p$ if $F\mapsto W(F)$ is rank-one convex. We show this disturbing feature with the help of a recently considered family of exponentiated Hencky energies.

preprint2015arXiv

Numerical Treatment of a Geometrically Nonlinear Planar Cosserat Shell Model

We present a new way to discretize a geometrically nonlinear elastic planar Cosserat shell. The kinematical model is similar to the general 6-parameter resultant shell model with drilling rotations. The discretization uses geodesic finite elements, which leads to an objective discrete model which naturally allows arbitrarily large rotations. Finite elements of any approximation order can be constructed. The resulting algebraic problem is a minimization problem posed on a nonlinear finite-dimensional Riemannian manifold. We solve this problem using a Riemannian trust-region method, which is a generalization of Newton's method that converges globally without intermediate loading steps. We present the continuous model and the discretization, discuss the properties of the discrete model, and show several numerical examples, including wrinkles of thin elastic sheets in shear.

preprint2015arXiv

On some fundamental misunderstandings in the indeterminate couple stress model. A comment on recent papers of A.R. Hadjesfandiari and G.F. Dargush

In a series of papers which are either published [A.R. Hadjesfandiari and G.F. Dargush, Couple stress theory for solids, Int. J. Solids Struct. 48, 2496-2510, 2011; A.R. Hadjesfandiari and G.F. Dargush, Fundamental solutions for isotropic size-dependent couple stress elasticity, Int. J. Solids Struct. 50, 1253-1265, 2013] or available as preprints Hadjesfandiari and Dargush have reconsidered the linear indeterminate couple stress model. They are postulating a certain physically plausible split in the virtual work principle. Based on this postulate they claim that the second-order couple stress tensor must always be skew-symmetric. Since they use an incomplete set of boundary conditions in their virtual work principle their statement contains unrecoverable errors. This is shown by specifying their development to the isotropic case. However, their choice of constitutive parameters is mathematically possible and still yields a well-posed boundary value problem.

preprint2015arXiv

On the convexity of nonlinear elastic energies in the right Cauchy-Green tensor

We present a sufficient condition under which a weak solution of the Euler-Lagrange equations in nonlinear elasticity is already a global minimizer of the corresponding elastic energy functional. This criterion is applicable to energies $W(F)=\widehat{W}(F^TF)=\widehat{W}(C)$ which are convex with respect to the right Cauchy-Green tensor $C=F^TF$, where $F$ denotes the gradient of deformation. Examples of such energies exhibiting a blow up for $\det F\to0$ are given.

preprint2015arXiv

On the generalized sum of squared logarithms inequality

Assume $n\geq 2$. Consider the elementary symmetric polynomials $e_k(y_1,y_2,\ldots, y_n)$ and denote by $E_0,E_1,\ldots,E_{n-1}$ the elementary symmetric polynomials in reverse order \begin{align*} E_k(y_1,y_2,\ldots,y_n):=e_{n-k}(y_1,y_2,\ldots,y_n)=\sum_{i_1<\ldots<i_{n-k}} y_{i_1}y_{i_2}\ldots y_{i_{n-k}}\, , \quad k\in \{0,1,\ldots,n{-}1 \}\, . \end{align*} Let moreover $S$ be a nonempty subset of $\{0,1,\ldots,n{-}1\}$. We investigate necessary and sufficient conditions on the function $f\colon\,I\to\mathbb{R}$, where $I\subset\mathbb{R}$ is an interval, such that the inequality \begin{align} \label{abstract_inequality} f(a_1)+f(a_2)+\ldots+f(a_n)\leq f(b_1)+f(b_2)+\ldots+f(b_n) \tag{*} \end{align} holds for all $a=(a_1,a_2,\ldots,a_n)\in I^n$ and $b=(b_1,b_2,\ldots,b_n)\in I^n$ satisfying $$E_k(a)< E_k(b) \ \hbox{for } k\in S\quad \hbox{and} \quad E_k(a)=E_k(b) \ \hbox{for } k\in \{0,1,\ldots,n{-}1 \}\setminus S\, .$$ As a corollary, we obtain \eqref{abstract_inequality} if $2\leq n\leq 4$, $f(x)=\log^2x$ and $S=\{1,\dotsc,n-1\}$, which is the sum of squared logarithms inequality previously known for $2\le n\le 3$.

preprint2015arXiv

On the sum of squared logarithms inequality and related inequalities

We consider the sum of squared logarithms inequality and investigate possible connections with the theory of majorization. We also discuss alternative sufficient conditions on two sets of vectors $a,b\in\mathbb{R}_+^n$ so that $\sum_{i=1}^n(\log a_i)^2\ \leq\ \sum_{i=1}^n(\log b_i)^2\,.\notag $ Generalizations of some inequalities from information theory are obtained, including a generalized information inequality and a generalized log sum inequality, which states for $a,b\in\mathbb{R}_+^n$ and $k_1,...,k_n\in [0,\infty)$: $ \sum_{i=1}^na_i\,\log\prod_{s=1}^m(\frac{a_i}{b_i} + k_s)\ \geq\ \log\prod_{s=1}^m(1+k_s)\,.\notag $

preprint2015arXiv

Soliton-like solutions based on geometrically nonlinear Cosserat micropolar elasticity

The Cosserat model generalises an elastic material taking into account the possible microstructure of the elements of the material continuum. In particular, within the Cosserat model the structured material point is rigid and can only experience microrotation, which is also known as micropolar elasticity. We present the geometrically nonlinear theory taking into account all possible interaction terms between the elastic and microelastic structure. This is achieved by considering the irreducible pieces of the deformation gradient and of the dislocation curvature tensor. In addition we also consider the so-called Cosserat coupling term. In this setting we seek soliton type solutions assuming small elastic displacements, however, we allow the material points to experience full rotations which are not assumed to be small. By choosing a particular ansatz we are able to reduce the system of equations to a Sine-Gordon type equation which is known to have soliton solutions.

preprint2015arXiv

The exponentiated Hencky-logarithmic strain energy. Improvement of planar polyconvexity

In this paper we improve the result about the polyconvexity of the energies from the family of isotropic volumetric-isochoric decoupled strain exponentiated Hencky energies defined in the first part of this series, i.e. $$ W_{_{\rm eH}}(F)= \left\{\begin{array}{lll} \fracμ{k}\,e^{k\,\|{\rm dev}_n\log U\|^2}+\fracκ{2\,\widehat{k}}\,e^{\widehat{k}\,[(\log {\rm det} U)]^2}&\text{if}& {\rm det}\, F>0,\\ +\infty &\text{if} &{\rm det} F\leq 0\,, \end{array}\right. $$ where $F=\nabla φ$ is the gradient of deformation, $U=\sqrt{F^T F}$ is the right stretch tensor and ${\rm dev}_n\log {U}$ is the deviatoric part of the strain tensor $\log U$. The main result in this paper is that in plane elastostatics, i.e. for $n=2$, the energies of this family are polyconvex for $k\geq \frac{1}{4}$, $\widehat{k}\geq \frac{1}{8}$, extending a previous result which proves polyconvexity for $k\geq \frac{1}{3}$, $\widehat{k}\geq \frac{1}{8}$. This leads immediately to an extension of the existence result.

preprint2015arXiv

The geometrically nonlinear Cosserat micropolar shear-stretch energy. Part I: A general parameter reduction formula and energy-minimizing microrotations in 2D

In any geometrically nonlinear quadratic Cosserat-micropolar extended continuum model formulated in the deformation gradient field $F := \nablaφ: Ω\to \mathrm{GL}^+(n)$ and the microrotation field $R: Ω\to \mathrm{SO}(n)$, the shear-stretch energy is necessarily of the form \begin{equation*} W_{μ,μ_c}(R\,;F) := μ\,\left\lVert{\mathrm{sym}(R^T F - \boldsymbol{1})}\right\rVert^2 + μ_c\,\left\lVert{\mathrm{skew}(R^T F - \boldsymbol{1})}\right\rVert^2\;, \end{equation*} where $μ> 0$ is the Lamé shear modulus and $μ_c \geq 0$ is the Cosserat couple modulus. In the present contribution, we work towards explicit characterizations of the set of optimal Cosserat microrotations $\mathrm{argmin}_{R\,\in\,\mathrm{SO}(n)}{W_{μ,μ_c}(R\,;F)}$ as a function of $F \in \mathrm{GL}^+(n)$ and weights $μ> 0$ and $μ_c \geq 0$. For $n \geq 2$, we prove a parameter reduction lemma which reduces the optimality problem to two limit cases: $(μ, μ_c) = (1,1)$ and $(μ,μ_c) = (1,0)$. In contrast to Grioli's theorem, we derive non-classical minimizers for the parameter range $μ> μ_c \geq 0$ in dimension $n\!=\!2$. Currently, optimality results for $n \geq 3$ are out of reach for us, but we contribute explicit representations for $n\!=\!2$ which we name $\mathrm{rpolar}^{\pm}_{μ,μ_c}(F) \in \mathrm{SO}(2)$ and which arise for $n\!=\!3$ by fixing the rotation axis a priori. Further, we compute the associated reduced energy levels and study the non-classical optimal Cosserat rotations $\mathrm{rpolar}^\pm_{μ,μ_c}(F_γ)$ for simple planar shear.

preprint2015arXiv

The geometrically nonlinear Cosserat micropolar shear-stretch energy. Part II: Non-classical energy-minimizing microrotations in 3D and their computational validation

In any geometrically nonlinear, isotropic and quadratic Cosserat micropolar extended continuum model formulated in the deformation gradient field $F = \nablaφ: Ω\to GL^+(n)$ and the microrotation field $R: Ω\to SO(n)$, the shear-stretch energy is necessarily of the form $$W_{μ,μ_c}(R;F) = μ\, \| sym(R^T F - 1) \|^2 + μ_c\, \| skew(R^T F - 1) \|^2 .$$ We aim at the derivation of closed form expressions for the minimizers of $W(R;F)$ in $SO(3)$, i.e., for the set of optimal Cosserat microrotations in dimension $n = 3$, as a function of $F \in GL^+(n)$. In a previous contribution (Part I), we have first shown that, for all $n \geq 2$, the full range of weights $μ> 0$ and $μ_c \geq 0$ can be reduced to either a classical or a non-classical limit case. We have then derived the associated closed form expressions for the optimal planar rotations in $SO(2)$ and proved their global optimality. In the present contribution (Part II), we characterize the non-classical optimal rotations in dimension n = 3. After a lift of the minimization problem to the unit quaternions, the Euler-Lagrange equations can be symbolically solved by the computer algebra system Mathematica. Among the symbolic expressions for the critical points, we single out two candidates $rpolar^{\pm}_{μ,μ_c}(F) \in SO(3)$ which we analyze and for which we can computationally validate their global optimality by Monte Carlo statistical sampling of $SO(3)$. Geometrically, our proposed optimal Cosserat rotations $rpolar^{\pm}_{μ,μ_c}(F)$ act in the "plane of maximal strain" and our previously obtained explicit formulae for planar optimal Cosserat rotations in $SO(2)$ reveal themselves as a simple special case. Further, we derive the associated reduced energy levels of the Cosserat shear--stretch energy and criteria for the existence of non-classical optimal rotations.

preprint2015arXiv

The sum of squared logarithms inequality in arbitrary dimensions

We prove the \emph{sum of squared logarithms inequality} (SSLI) which states that for nonnegative vectors $x, y \in \mathbb{R}^n$ whose elementary symmetric polynomials satisfy $e_k(x)\le e_k(y)$ (for $1\le k < n$) and $e_n(x)=e_n(y)$, the inequality $\sum_i (\log x_i)^2 \le \sum_i (\log y_i)^2$ holds. Our proof of this inequality follows by a suitable extension to the complex plane. In particular, we show that the function $f\colon M\subseteq \mathbb{C}^n\to \mathbb{R}$ with $f(z)=\sum_i(\log z_i)^2$ has nonnegative partial derivatives with respect to the elementary symmetric polynomials of $z$. This property leads to our proof. We conclude by providing applications and wider connections of the SSLI.

preprint2015arXiv

Wave propagation in pantographic 2D lattices with internal discontinuities

In the present paper we consider a 2D pantographic structure composed by two orthogonal families of Euler beams. Pantographic rectangular 'long' waveguides are considered in which imposed boundary displacements can induce the onset of traveling (possibly non-linear) waves. We performed numerical simulations concerning a set of dynamically interesting cases. The system undergoes large rotations which may involve geometrical non-linearities, possibly opening the path to appealing phenomena such as propagation of solitary waves. Boundary conditions dramatically influence the transmission of the considered waves at discontinuity surfaces. The theoretical study of this kind of objects looks critical, as the concept of pantographic 2D sheets seems to have promising possible applications in a number of fields, e.g. acoustic filters, vascular prostheses and aeronautic/aerospace panels.

preprint2014arXiv

A first regularity result for the Armstrong-Frederick cyclic hardening plasticity model with Cosserat effects

The purpose of this article is to prove the Hölder continuity up to the boundary of the displacement vector and the microrotation matrix for the quasistatic, rate-independent Armstrong-Frederick cyclic hardening plasticity model with Cosserat effects. This model is of non-monotone and non-associated type. In the case of two space dimensions we use the hole-filling technique of Widman and the Morrey's Dirichlet growth theorem.

preprint2014arXiv

A logarithmic minimization property of the unitary polar factor in the spectral norm and the Frobenius matrix norm

The unitary polar factor $Q=U$ in the polar decomposition of the matrix $Z=UH$ is the minimizer for both $\| \mathrm{Log}(Q^* Z)\|^2$ and its Hermitian part $\| \mathrm{sym Log}(Q^* Z)\|^2$ over both $\mathbb{R}$ and $\mathbb{C}$, for any given invertible matrix $Z$ in $\mathbb{C}^{n\times n}$ and any matrix logarithm $\mathrm{Log}$, not necessarily the principal logarithm $\mathrm{log}$. We prove this for the spectral matrix norm in any dimension and for the Frobenius matrix norm in two and three dimensions. The result shows that the unitary polar factor is the nearest orthogonal matrix to $Z$ not only in the normwise sense, but also in a geodesic distance. The derivation is based on Bhatia's generalization of Bernstein's trace inequality for the matrix exponential and a new sum of squared logarithms inequality.

preprint2014arXiv

Existence theorem for geometrically nonlinear Cosserat micropolar model under uniform convexity requirements

We reconsider the geometrically nonlinear Cosserat model for a uniformly convex elastic energy and write the equilibrium problem as a minimization problem. Applying the direct methods of the calculus of variations we show the existence of minimizers. We present a clear proof based on the coercivity of the elastically stored energy density and on the weak lower semi-continuity of the total energy functional. Use is made of the dislocation density tensor $\bar{\boldsymbol{K}}=\bar{\boldsymbol{R}}^T\,\mathrm{Curl}\,\bar{\boldsymbol{R}}$ as a suitable Cosserat curvature measure.

preprint2014arXiv

Poincare meets Korn via Maxwell: Extending Korn's First Inequality to Incompatible Tensor Fields

For a bounded three-dimensional domain with Lipschitz boundary we extend Korn's first inequality to incompatible tensor fields. For compatible tensor fields our estimate reduces to a non-standard variant of the well known Korn's first inequality. On the other hand, for skew-symmetric tensor fields our new estimate turns to Poincare's inequality. Therefore, our result may be viewed as a natural common generalization of Korn's first and Poincare's inequality. Decisive tools for this unexpected estimate are the classical Korn's first inequality, Helmholtz decompositions for mixed boundary conditions and the Maxwell estimate.

preprint2014arXiv

Rediscovering G.F. Becker's early axiomatic deduction of a multiaxial nonlinear stress-strain relation based on logarithmic strain

We discuss a completely forgotten work of the geologist G.F. Becker on the ideal isotropic nonlinear stress-strain function. In doing this we provide the original paper from 1893 newly typeset in LaTeX and with corrections of typographical errors as well as an updated notation. Due to the fact that the mathematical modelling of elastic deformations has evolved greatly since the original publication we give a modern reinterpretation of Becker's work, combining his approach with the current framework of the theory of nonlinear elasticity. Interestingly, Becker introduces a multiaxial constitutive law incorporating the logarithmic strain tensor, more than 35 years before the quadratic Hencky strain energy was introduced by Heinrich Hencky in 1929. Becker's deduction is purely axiomatic in nature. He considers the finite strain response to applied shear stresses and spherical stresses, formulated in terms of the principal strains and stresses, and postulates a principle of superposition for principal forces which leads, in a straightforward way, to a unique invertible constitutive relation. In the special case where Poisson's number is zero, the formulation is hyperelastic and the corresponding strain energy has the form of the maximum entropy function.

preprint2014arXiv

The axiomatic deduction of the quadratic Hencky strain energy by Heinrich Hencky

The introduction of the quadratic Hencky strain energy based on the logarithmic strain tensor log V is a milestone in the development of nonlinear elasticity theory in the first half of the 20th century. Since the original manuscripts are written in German, they are not easily accessible today. However, we believe that the deductive approach taken by Hencky deserves to be rediscovered today. In this work we have gathered parts of the original contributions "Über die Form des Elastizitätsgesetzes bei ideal elastischen Stoffen", "Welche Umstände bedingen die Verfestigung bei der bildsamen Verformung von festen isotropen Körpern?" and "Das Superpositionsgesetz eines endlich deformierten relaxationsfähigen elastischen Kontinuums und seine Bedeutung für eine exakte Ableitung der Gleichungen für die zähe Flüssigkeit in der Eulerschen Form" which center around this deductive approach. We tried to provide, for the first time, a faithful translation into English. All footnotes are our addition.

preprint2014arXiv

The exponentiated Hencky-logarithmic strain energy. Part I: Constitutive issues and rank-one convexity

We investigate a family of isotropic volumetric-isochoric decoupled strain energies $$ F\mapsto W_{_{\rm eH}}(F):=\widehat{W}_{_{\rm eH}}(U):=\left\{\begin{array}{lll} \fracμ{k}\,e^{k\,\|{\rm dev}_n\log {U}\|^2}+\fracκ{2\, {\widehat{k}}}\,e^{\widehat{k}\,[{ \rm tr}(\log U)]^2}&\text{if}& { \rm det} F>0,\\ +\infty &\text{if} &{ \rm det} F\leq 0, \end{array}\right.\quad $$ based on the Hencky-logarithmic (true, natural) strain tensor $\log U$, where $μ>0$ is the infinitesimal shear modulus, $κ=\frac{2μ+3λ}{3}>0$ is the infinitesimal bulk modulus with $λ$ the first Lamé constant, $k,\widehat{k}$ are dimensionless parameters, $F=\nabla φ$ is the gradient of deformation, $U=\sqrt{F^T F}$ is the right stretch tensor and ${\rm dev}_n\log {U} =\log {U}-\frac{1}{n} {\rm tr}(\log {U})\cdot 1\!\!1$ is the deviatoric part of the strain tensor $\log U$. For small elastic strains, $W_{_{\rm eH}}$ approximates the classical quadratic Hencky strain energy $$ F\mapsto W_{_{\rm H}}(F):=\widehat{W}_{_{\rm H}}(U):=μ\,\|{\rm dev}_n\log U\|^2+\fracκ{2}\,[{\rm tr}(\log U)]^2, $$ which is not everywhere rank-one convex. In plane elastostatics, i.e. $n=2$, we prove the everywhere rank-one convexity of the proposed family $W_{_{\rm eH}}$, for $k\geq \frac{1}{4}$ and $\widehat{k}\geq \frac{1}{8}$. Moreover, we show that the corresponding Cauchy (true)-stress-true-strain relation is invertible for $n=2,3$ and we show the monotonicity of the Cauchy (true) stress tensor as a function of the true strain tensor in a domain of bounded distortions. We also prove that the rank-one convexity of the energies belonging to the family $W_{_{\rm eH}}$ is not preserved in dimension $n=3$.

preprint2014arXiv

The exponentiated Hencky-logarithmic strain energy. Part II: Coercivity, planar polyconvexity and existence of minimizers

We consider a family of isotropic volumetric-isochoric decoupled strain energies $$ F\mapsto W_{\rm eH}(F):=\widehat{W}_{\rm eH}(U):=\left\{\begin{array}{lll} \fracμ{k}\,e^{k\,\|{\rm dev}_n\log {U}\|^2}+\fracκ{2\hat{k}}\,e^{\hat{k}\,[{\rm tr}(\log U)]^2}&\text{if}& {\rm det}\, F>0,\\ +\infty &\text{if} &{\rm det} F\leq 0, \end{array}\right.\quad $$ based on the Hencky-logarithmic (true, natural) strain tensor $\log U$, where $μ>0$ is the infinitesimal shear modulus, $κ=\frac{2μ+3λ}{3}>0$ is the infinitesimal bulk modulus with $λ$ the first Lamé constant, $k,\hat{k}$ are dimensionless parameters, $F=\nabla φ$ is the gradient of deformation, $U=\sqrt{F^T F}$ is the right stretch tensor and ${\rm dev}_n\log {U} =\log {U}-\frac{1}{n} {\rm tr}(\log {U})\cdot 1\!\!1$ is the deviatoric part (the projection onto the traceless tensors) of the strain tensor $\log U$. For small elastic strains the energies reduce to first order to the classical quadratic Hencky energy $$ F\mapsto W{_{\rm H}}(F):=\widehat{W}_{_{\rm H}}(U):=μ\,\|{\rm dev}_n\log U\|^2+\fracκ{2}\,[{\rm tr}(\log U)]^2, $$ which is known to be not rank-one convex. The main result in this paper is that in plane elastostatics the energies of the family $W_{_{\rm eH}}$ are polyconvex for $k\geq \frac{1}{3}$, $\widehat{k}\geq \frac{1}{8}$, extending a previous finding on its rank-one convexity. Our method uses a judicious application of Steigmann's polyconvexity criteria based on the representation of the energy in terms of the principal invariants of the stretch tensor $U$. These energies also satisfy suitable growth and coercivity conditions. We formulate the equilibrium equations and we prove the existence of minimizers by the direct methods of the calculus of variations.

preprint2014arXiv

The exponentiated Hencky-logarithmic strain energy. Part III: Coupling with idealized isotropic finite strain plasticity

We investigate an immediate application in finite strain multiplicative plasticity of the family of isotropic volumetric-isochoric decoupled strain energies \begin{align*} F\mapsto W_{_{\rm eH}}(F):=\hat{W}_{_{\rm eH}}(U):=\{\begin{array}{lll} \fracμ{k}\,e^{k\,\|{\rm dev}_n\log {U}\|^2}+\fracκ{\text{}{2\, {\hat{k}}}}\,e^{\hat{k}\,[{\rm tr}(\log U)]^2}&\text{if}& {\rm det}\, F>0,\\ +\infty &\text{if} &{\rm det} F\leq 0, \end{array}.\quad \end{align*} based on the Hencky-logarithmic (true, natural) strain tensor $\log U$. Here, $μ>0$ is the infinitesimal shear modulus, $κ=\frac{2μ+3λ}{3}>0$ is the infinitesimal bulk modulus with $λ$ the first Lamé constant, $k,\hat{k}$ are dimensionless fitting parameters, $F=\nabla φ$ is the gradient of deformation, $U=\sqrt{F^T F}$ is the right stretch tensor and ${\rm dev}_n\log {U} =\log {U}-\frac{1}{n}\, {\rm tr}(\log {U})\cdot 1\!\!1$ is the deviatoric part of the strain tensor $\log U$. Based on the multiplicative decomposition $F=F_e\, F_p$, we couple these energies with some isotropic elasto-plastic flow rules $F_p\,\frac{\rm d}{{\rm d} t}[F_p^{-1}]\in-\partial χ({\rm dev}_3 Σ_{e})$ defined in the plastic distortion $F_p$, where $\partial χ$ is the subdifferential of the indicator function $χ$ of the convex elastic domain $\mathcal{E}_{\rm e}(W_{\rm iso},{Σ_{e}},\frac{1}{3}{\boldsymbolσ}_{\!\mathbf{y}}^2)$ in the mixed-variant $Σ_{e}$-stress space and $Σ_{e}=F_e^T D_{F_e} W_{\rm iso}(F_e)$. While $W_{_{\rm eH}}$ may loose ellipticity, we show that loss of ellipticity is effectively prevented by the coupling with plasticity, since the ellipticity domain of $W_{_{\rm eH}}$ on the one hand, and the elastic domain in $Σ_{e}$-stress space on the other hand, are closely related.

preprint2014arXiv

The relaxed linear micromorphic continuum: well-posedness of the static problem and relations to the gauge theory of dislocations

In this paper we consider the equilibrium problem in the relaxed linear model of micromorphic elastic materials. The basic kinematical fields of this extended continuum model are the displacement $u\in \mathbb{R}^3$ and the non-symmetric micro-distortion density tensor $P\in \mathbb{R}^{3\times 3}$. In this relaxed theory a symmetric force-stress tensor arises despite the presence of microstructure and the curvature contribution depends solely on the micro-dislocation tensor ${\rm Curl}\, P$. However, the relaxed model is able to fully describe rotations of the microstructure and to predict non-polar size-effects. In contrast to classical linear micromorphic models, we allow the usual elasticity tensors to become positive-semidefinite. We prove that, nevertheless, the equilibrium problem has a unique weak solution in a suitable Hilbert space. The mathematical framework also settles the question of which boundary conditions to take for the micro-distortion. Similarities and differences between linear micromorphic elasticity and dislocation gauge theory are discussed and pointed out.

preprint2014arXiv

Well-posedness for dislocation based gradient visco-plasticity with isotropic hardening

In this work we establish the well-posedness for infinitesimal dislocation based gradient viscoplasticity with isotropic hardening for general gradient monotone plastic flows. We assume an additive split of the displacement gradient into non-symmetric elastic distortion and non-symmetric plastic distortion. The thermodynamic potential is augmented with a term taking the dislocation density tensor $\operatorname{Curl} p$ into account. The constitutive equations in the models we study are assumed to be of self-controlling type. Based on the generalized version of Korn's inequality for incompatible tensor fields (the non-symmetric plastic distortion) due to Neff/Pauly/Witsch the existence of solutions of quasi-static initial-boundary value problems under consideration is shown using a time-discretization technique and a monotone operator method.

preprint2013arXiv

A Riemannian approach to strain measures in nonlinear elasticity

The isotropic Hencky strain energy appears naturally as a distance measure of the deformation gradient to the set SO(n) of rigid rotations in the canonical left-invariant Riemannian metric on the general linear group GL(n). Objectivity requires the Riemannian metric to be left-GL(n)-invariant, isotropy requires the Riemannian metric to be right-O(n)-invariant. The latter two conditions are satisfied for a three-parameter family of Riemannian metrics on the tangent space of GL(n). Surprisingly, the final result is basically independent of the chosen parameters. In deriving the result, geodesics on GL(n) have to be parametrized and a novel minimization problem, involving the matrix logarithm for non-symmetric arguments, has to be solved.

preprint2013arXiv

A unifying perspective: the relaxed linear micromorphic continuum

We formulate a relaxed linear elastic micromorphic continuum model with symmetric Cauchy force-stresses and curvature contribution depending only on the micro-dislocation tensor. Our relaxed model is still able to fully describe rotation of the microstructure and to predict non-polar size-effects. It is intended for the homogenized description of highly heterogeneous, but non polar materials with microstructure liable to slip and fracture. In contrast to classical linear micromorphic models our free energy is not uniformly pointwise positive definite in the control of the independent constitutive variables. The new relaxed micromorphic model supports well-posedness results for the dynamic and static case. There, decisive use is made of new coercive inequalities recently proved by Neff, Pauly and Witsch and by Bauer, Neff, Pauly and Starke. The new relaxed micromorphic formulation can be related to dislocation dynamics, gradient plasticity and seismic processes of earthquakes. It unifies and simplifies the understanding of the linear micromorphic models.

preprint2013arXiv

Counterexamples in the theory of coerciveness for linear elliptic systems related to generalizations of Korn's second inequality

We show that the following generalized version of Korn's second inequality with nonconstant measurable matrix valued coefficients P ||DuP+(DuP)^T||_q+||u||_q >= c ||Du||_q for u in W_0^{1,q}(Ω;R^3), 1<q<{\infty} is in general false, even if P is in SO(3), while the Legendre-Hadamard condition and ellipticity on C^n for the quadratic form |Du P+(DuP)^T|^2 is satisfied. Thus Garding's inequality may be violated for formally positive quadratic forms.

preprint2013arXiv

Dev-Div- and DevSym-DevCurl-inequalities for incompatible square tensor fields with mixed boundary conditions

For an n-dimensional bounded domain we derive some inequalities bounding the norm of a square tensor field. Concerning the Div-Dev-inequality the bound is given by the trace-free part and the divergence and the tensor. In the case of the DevSym-Curl-inequality the bound is given by the trace-free and symmetric part and the curl of the tensor. For n=3 the bound is given by the trace-free symmetric part of the tensor and the trace-free part of the curl of the tensor. Some prototype applications are presented in which the new inequalities may be used to derive the coercivity of the models.

preprint2013arXiv

Existence of minimizers in the geometrically non-linear 6-parameter resultant shell theory with drilling rotations

The paper is concerned with the geometrically non-linear theory of 6-parametric elastic shells with drilling degrees of freedom. This theory establishes a general model for shells, which is characterized by two independent kinematic fields: the translation vector and the rotation tensor. Thus, the kinematical structure of 6-parameter shells is identical to that of Cosserat shells. We show the existence of global minimizers for the geometrically non-linear 2D equations of elastic shells. The proof of the existence theorem is based on the direct methods of the calculus of variations using essentially the convexity of the energy in the strain and curvature measures. Since our result is valid for general anisotropic shells, we analyze separately the particular cases of isotropic shells, orthotropic shells, and composite shells.

preprint2013arXiv

Existence theorems in the geometrically non-linear 6-parametric theory of elastic plates

In this paper we show the existence of global minimizers for the geometrically exact, non-linear equations of elastic plates, in the framework of the general 6-parametric shell theory. A characteristic feature of this model for shells is the appearance of two independent kinematic fields: the translation vector field and the rotation tensor field (representing in total 6 independent scalar kinematic variables). For isotropic plates, we prove the existence theorem by applying the direct methods of the calculus of variations. Then, we generalize our existence result to the case of anisotropic plates. We also present a detailed comparison with a previously established Cosserat plate model.

preprint2013arXiv

On an Extension of Korn's First Inequality to Incompatible Tensor Fields on Domains of Arbitrary Dimensions

For a bounded N-dimensional domain with Lipschitz boundary we extend Korn's first inequality to incompatible tensor fields. For compatible tensor fields our estimate reduces to a non-standard variant of the well known Korn's first inequality. On the other hand, for skew-symmetric tensor fields our new estimate turns to Poincare's inequality. Therefore, our result may be viewed as a natural common generalization of Korn's first and Poincare's inequality. Decisive tools for this unexpected estimate are the classical Korn's first inequality, Helmholtz decompositions for mixed boundary conditions and the Maxwell estimate.

preprint2013arXiv

On Grioli's minimum property and its relation to Cauchy's polar decomposition

The unitary polar factor of a matrix F is the unitary matrix Q realizing the minimum of the norm of F-Q over all unitary matrices Q. Tracing back the development on the optimality of the polar factor to its presumable roots, in this paper we present a commented translation of a note by G. Grioli from 1940 (Boll.Un.Math.Ital.2,452-455(1940)) showing the minimization property for the Frobenius norm in dimension 3; in his words: "We show that for the homogeneous displacement, tangent to any finite displacement, there exists a minimum property completely analogous to a well known property valid for infinitesimal displacements." Keywords: polar decomposition, optimality of the polar factor, Euclidean distance, geodesic distance, Euclidean movement

preprint2013arXiv

Sum of squared logarithms - An inequality relating positive definite matrices and their matrix logarithm

Let y1, y2, y3, a1, a2, a3 > 0 be such that y1 y2 y3 = a1 a2 a3 and y1 + y2 + y3 >= a1 + a2 + a3, y1 y2 + y2 y3 + y1 y3 >= a1 a2 + a2 a3 + a1 a3. Then the following inequality holds (log y1)^2 + (log y2)^2 + (log y3)^2 >= (log a1)^2 + (log a2)^2 + (log a3)^2. This can also be stated in terms of real positive definite 3x3-matrices P1, P2: If their determinants are equal det P1 = det P2, then tr P1 >= tr P2 and tr Cof P1 >= tr Cof P2 implies norm(log P1) >= norm(log P2), where log is the principal matrix logarithm and norm(P) denotes the Frobenius matrix norm. Applications in matrix analysis and nonlinear elasticity are indicated.

preprint2013arXiv

The Armstrong-Frederick cyclic hardening plasticity model with Cosserat effects

We propose an extension of the cyclic hardening plasticity model formulated by Armstrong and Frederick which includes micropolar effects. Our micropolar extension establishes coercivity of the model which is otherwise not present. We study then existence of solutions to the quasistatic, rate-independent Armstrong-Frederick model with Cosserat effects which is, however, still of non-monotone, non-associated type. In order to do this, we need to relax the pointwise definition of the flow rule into a suitable weak energy-type inequality. It is shown that the limit in the Yosida approximation process satisfies this new solution concept. The limit functions have a better regularity than previously known in the literature, where the original Armstrong-Frederick model has been studied.

preprint2013arXiv

The minimization of matrix logarithms - on a fundamental property of the unitary polar factor

We show that the unitary factor U in the polar decomposition of a nonsingular matrix Z = U H is the minimizer for both ||Log(Q^* Z)|| and ||sym (Log(Q^*Z))|| over unitary Q, for any given invertible complex n-times-n matrix Z, for any unitarily invariant norm and any n. We prove that U is the unique matrix with this property. As important tools we use a generalized Bernstein trace inequality and the theory of majorization.

preprint2013arXiv

The relaxed linear micromorphic continuum: existence, uniqueness and continuous dependence in dynamics

We study well-posedness for the relaxed linear elastic micromorphic continuum model with symmetric Cauchy force-stresses and curvature contribution depending only on the micro-dislocation tensor. In contrast to classical micromorphic models our free energy is not uniformly pointwise positive definite in the control of the independent constitutive variables. Another interesting feature concerns the prescription of boundary values for the micro-distortion field: only tangential traces may be determined which are weaker than the usual strong anchoring boundary condition. There, decisive use is made of new coercive inequalities recently proved by Neff, Pauly and Witsch and by Bauer, Neff, Pauly and Starke. The new relaxed micromorphic formulation can be related to dislocation dynamics, gradient plasticity and seismic processes of earthquakes.

preprint2013arXiv

Wave propagation in relaxed micromorphic continua: modelling metamaterials with frequency band-gaps

In this paper the relaxed micromorphic model proposed in [Patrizio Neff, Ionel-Dumitrel Ghiba, Angela Madeo, Luca Placidi, Giuseppe Rosi. A unifying perspective: the relaxed linear micromorphic continuum, submitted, 2013, arXiv:1308.3219; and Ionel-Dumitrel Ghiba, Patrizio Neff, Angela Madeo, Luca Placidi, Giuseppe Rosi. The relaxed linear micromorphic continuum: existence, uniqueness and continuous dependence in dynamics, submitted, 2013, arXiv:1308.3762] has been used to study wave propagation in unbounded continua with microstructure. By studying dispersion relations for the considered relaxed medium, we are able to disclose precise frequency ranges (band-gaps) for which propagation of waves cannot occur. These dispersion relations are strongly nonlinear so giving rise to a macroscopic dispersive behavior of the considered medium. We prove that the presence of band-gaps is related to a unique elastic coefficient, the so-called Cosserat couple modulus $μ_{c}$, which is also responsible for the loss of symmetry of the Cauchy force stress tensor. This parameter can be seen as the trigger of a bifurcation phenomenon since the fact of slightly changing its value around a given threshold drastically changes the observed response of the material with respect to wave propagation. We finally show that band-gaps cannot be accounted for by classical micromorphic models as well as by Cosserat and second gradient ones. The potential fields of application of the proposed relaxed model are manifold, above all for what concerns the conception of new engineering materials to be used for vibration control and stealth technology.