# Constitutive models based on a quadri-dimensionnal formalism for numerical simulations of finite transformations

Emmanuelle Rouhaud, Benoît Panicaud, Arjen Roos, Najoua Mhenni,
Richard Kerner

► To cite this version:

Emmanuelle Rouhaud, Benoît Panicaud, Arjen Roos, Najoua Mhenni, Richard Kerner. Constitutive models based on a quadri-dimensionnal formalism for numerical simulations of finite transformations. 11e colloque national en calcul des structures, CSMA, May 2013, Giens, France. hal-01717080

HAL Id: hal-01717080

https://hal.science/hal-01717080

Submitted on 25 Feb 2018

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

Public Domain
# Constitutive models based on a quadri-dimensionnal formalism for numerical simulations of finite transformations

Emmanuelle ROUHAUD<fnref n="1" />, Benoît PANICAUD<fnref n="2" />, Arjen ROOS<fnref n="3" />, Najoua BEN MHENNI<fnref n="4" />, Richard KERNER<fnref n="5" />

**Résumé** — La description correcte des non-linéarités lors des transformations finies d'un milieu continu est une nécessité pour accrocher un comportement cinématique réaliste [1]. Nous proposons ainsi d'utiliser les outils mathématiques de la géométrie différentielle dans le cadre d'un formalisme quadridimensionnel et curviligne de l'espace-temps pour écrire des modèles de comportement. Cette approche garantie une description dite "covariante" des transformations subies par la matière, c'est-à-dire valable pour tous les systèmes de référence [2]. Cette approche a fait ses preuves en physique, en particulier en relativité générale [3].

Le but est de résoudre les difficultés qui se posent encore avec la notion d'objectivité matérielle, ambigüe en 3D [4]. Le principe de covariance généralisé est donc appliqué à des transformations qui restent celles de la mécanique classique des milieux continus (c'est-à-dire sans aucun effet relativiste et en particulier pour des vitesses qui restent toujours bien en deçà de la vitesse de la lumière). Ainsi, appliqué aux modèles de comportement purement macroscopique, il est montré comment les projections de la 4D sur l'espace 3D aboutissent à la construction de relations de comportement, par un formalisme variationnel. Ces relations sont ensuite déclinées pour des comportements particuliers (élasticité linéaire ou non-linéaire, pour différents choix de potentiels). En parallèle, une approche incrémentale est également proposée où le formalisme quadri-dimensionnel apporte des réponses uniques sur le choix de la dérivée "objective", comparé à la multitude des choix possibles en 3D. Les calculs éléments finis sur différentes structures montrent la comparaison entre ces différentes approchent pour différentes sollicitations. Une approche micromécanique pourra être finalement envisagée dans le cadre de l'élasticité pour en montrer à la fois les forces et les faiblesses.

**Mots clés** — transformations finies, géométrie différentielle, objectivité matérielle, élastoplasticité

# 1 Introduction

To say that a constitutive model has to verify the principle of material objectivity to ensure its frame-indifference is presented as common wisdom in the area of continuum mechanics. This problematic has been extensively studied, detailed and discussed. Among the points that have been subjected to debates are the principle itself [5], a (re)definition of the change of observers [5, 6] with inertial/non inertial consideration [7] or its link with the superposition of a rigid body motion [8, 9, 1]. An important issue concerns the definition of material frame-indifferent time derivatives to represent, objectively, the variations of a tensor with respect to time. These derivatives appear in incremental forms of constitutive models when the invariance under superposed rigid body motions is indeed observed. But the time derivative of an objective tensor is in general not objective. To work out this problem, 3D objective transport operators are defined; they are also referred to as "objective rates" or "invariant time fluxes". They are applied in particular to the Cauchy stress tensor. The result of each of these transports is proven to be invariant with respect to superposition of rigid body motions [10, 11, 1]. The difficulty resides in the fact that there are "infinitely many possible objective time fluxes that may be used" [1].

In the end, the validity of such an "objective" approach, and the restriction that it imposes on the

<footnote n="1">¹ Université de Technologie de Troyes, rouhaud@utt.fr</footnote>

<footnote n="2">² Université de Technologie de Troyes, benoit.panicaud@utt.fr</footnote>

<footnote n="3">³ ONERA, arjen.roos@onera.fr</footnote>

<footnote n="4">⁴ Université de Technologie de Troyes, najoua.ben_mhenni@utt.fr</footnote>

<footnote n="5">⁵ Université Pierre et Marie Curie, richard.kerner@upmc.fr</footnote>
constitutive models are often questioned and reconsidered; see for example [12, 13, 14, 6, 15, 7, 9, 16, 5, 17, 18, 19, 20, 21]. Thus, if everyone agrees on the necessity to define frame-indifferent entities and equations, the usage, applications and consequences of the notion of objectivity still constitute an open subject of debate.

Differential geometry, also known as Ricci-calculus, [22, 23], offers a convenient framework to model the finite transformations of a material continuum. It proposes a mathematical formalism for the use of tensors expressed in arbitrary coordinate systems within a differentiable manifold. The physical entities are then represented with tensor fields of the considered manifold. It is recognized as a formalism of choice to describe the straining motion of material continua within a classical three-dimensional (3D) context. Within such a formalism, tensors are independent of any arbitrary change of coordinate systems and it is possible to define derivative operators that are tensors of the considered space. Within its four-dimensional (4D) context, differential geometry has found a major and essential application in physics with the theory of General Relativity, which has shown its ability to deal properly with space-time transformations. The covariance principle guarantees the proper expression of all physical relations in any kind of systems of reference, as developed and reviewed by Landau and Lifshitz [3]. Indeed, as written by Eringen [24]: "Attempts to secure the invariance of the physical relations of motion from the observer have produced one of the great triumphs of twentieth-century physics. (...) Attempts to free the principles of classical mechanics from the motion of an observer were resolved by Einstein in his general theory of relativity...". The only appropriate way to define frame independence is thus to consider four-dimensional quantities, i.e. to construct physics and material mechanics within the scope of 4D physics. The present work thus proposes to use the covariance principle to describe a continuum, in order to guaranty the frame-indifference of tensors, equations and models. A major interest consists in using such a method to construct material constitutive relations. The objective of the present paper is to illustrate the method and use the 4D formalism to construct models for elastic transformations.

## 2 Material objectivity

In three dimensions, a rigid frame of reference is classically associated with an observer and offers the possibility to parametrize the position and instant of an occurring event. The frame of reference is constituted by a 3D Cartesian coordinate system to specify any positions, to which an origin for time and a Galilean chronology are associated. If the same event is now described by two observers with the respective coordinates and time $(x^i, t)$ and $(\tilde{x}^i, t)$, then an orthogonal matrix $Q$ and a vector $\delta$ can always be found such that:

$$ \tilde{x}^i = Q_j^i(t)x^j + \delta^i(t) \quad (1) $$

where $Q$ and $\delta$ correspond respectively to the instantaneous rotation and to the instantaneous translation of one frame with respect to the other. Then a second rank tensor $T$ is said to be *frame-indifferent or objective* if it verifies the transformation rule:

$$ \tilde{T}^{ij} = Q_k^i(t)Q_l^j(t)T^{kl} \quad (2) $$

for any changes of observers. Similar equations may be written for other component types and tensor ranks.

Truesdell and Noll [1] define the principle of objectivity as: "it is a fundamental principle of classical physics that material properties are indifferent, i.e., independent of the frame of reference or observer". Nemat Nasser [11] defines it as: "Constitutive relations must remain invariant under any rigid-body rotation of the reference coordinate system. This is called objectivity or the material frame indifference." As stated by Liu [25], "the principle of material frame indifference plays an important role in the development of continuum mechanics, by delivering restrictions on the formulation of the constitutive functions of material bodies. It is embedded in the idea that material properties should be independent of observations made by different observers. Since different observers are related by a time-dependent rigid transformation, known as a Euclidean transformation, material frame-indifference is sometimes interpreted as invariance under superposed rigid body motions". The two following physical notions could be distinguished indeed:
- On one hand, the invariance under superposed rigid body motions. For instance, when an unstrained body is animated with a rigid body motion, no stress is generated in the material. This property is related to the material and the observed mechanical quantities. It concerns most materials and mechanical phenomena.

- On the other hand, the independence with respect to the change of observers. For instance, when an observer moves around an unstrained material, no stress is generated in the material. This is a fundamental principle of classical physics.

Within a classical three-dimensional formalism, as stated by Liu, both properties result in the same mathematical condition described by Eq. 1 such that the term "objective" represents, ambiguously, both characteristics. It is further important to note that the 3D classical definition of observer-independence refers to invariance over Euclidian transformations and this is not verified by all 3D tensors.

### 3 The 4D formalism and the principle of covariance

As opposed to classical mechanics, the space is four-dimensional (4D); the fourth dimension references time as a coordinate. A curvilinear coordinate system $x^\mu$ references events such that:

$$x^\mu = (x^1, x^2, x^3, x^4) = (x^i, ct) \qquad (3)$$

Here, the indices $\mu$ varies from 1 to 4, while the indices $i$ varies from 1 to 3. Note that $x^i$ corresponds to the classical 3D coordinates. Greek indices refer to 4D entities while Latin indices refer only to their spatial parts. The fourth coordinate $x^4$ is the time $t$ multiplied by reference speed $c$ to be homogeneous with a length dimension as the other three spatial coordinates. Let $z^\mu$ represent the 4D Cartesian inertial coordinates. This coordinate system is said to be Galilean/Minkowskian. An interval $ds$ is then defined in the 4D space-time such that:

$$(ds)^2 = (dz^4)^2 - (dz^1)^2 - (dz^2)^2 - (dz^3)^2 = g_{\mu\nu}dx^\mu dx^\nu \qquad (4)$$

where $g_{\mu\nu}$ are the covariant components of the metric tensor of the coordinate system $x^\mu$.

In this space-time, physical entities are represented by 4D tensors. Consider now the same event parametrized with two different curvilinear coordinates $x^\mu$ and $\tilde{x}^\mu$. Tensors are indifferent to arbitrary changes of coordinate systems. Thus, the components of the second rank tensor $\alpha$ always transform through a change of coordinates from $x^\mu$ to $\tilde{x}^\mu$ as [24, 22, 26, 27]:

$$\tilde{\alpha}^{\mu\nu} = \frac{\partial \tilde{x}^{\mu}}{\partial x^{\lambda}} \frac{\partial \tilde{x}^{\nu}}{\partial x^{\kappa}} \alpha^{\lambda\kappa} \qquad (5)$$

It is possible to write similar equations for the components of tensors of any rank and for any components types: covariant, contravariant or mixed.

The objective of this work is to describe the finite transformations of a material body as it is traditionally studied in continuum mechanics. Consider an observer who measures the physical entities characterizing the state of the material body for each instant of time and for each point of space. In any case, the absolute velocity is always small with respect to the speed of light anywhere in the matter. Further, Newtonian mechanics is embedded in General Relativity, as a limiting theory concerning phenomena for which the absolute speed of each material point is negligible compared to the speed of light [28, 29, 30]. The physical hypotheses that are retained correspond to the classical hypotheses of Newtonian physics. It is thus supposed through out this work that the observed transformation is a classical continuum transformation encountered in mechanics: in particular, the speed of any point is negligible compared to the speed of light and the transformation does not change the definition of time. This is why we choose to always define time with respect to its absolute value in any frames of reference. It should be stressed however, that introducing the 4D-formalism means automatically passing to the realm of Relativity and introducing the Lorentz invariance. The velocity $c$ serves in our case, uniquely as useful dimensional constant unifying space and time units.

Physical observations lead to the conclusion that the presence, nature and number of the observers do not change the physical phenomena undergone by the matter. Observers should thus agree on the evaluation of these physical phenomena. This corresponds to the principle of general covariance first
formulated by A. Einstein [28]. It states that a physical equation must be expressed using a mathematical form that is strictly the same for each observer : the model should be independent of any arbitrary change of observers. This is possible with a 4D covariant description of physics. Within this formalism, all tensors are by construction independent of arbitrary changes of (possibly accelerated and deformed) observers. Because the fourth coordinate represents time, the observers are completely defined once the 4D coordinate system is chosen and 4D coordinate transformations describe changes of observers. We thus propose to describe the finite transformations of materials with a 4D covariant formalism.

## 4 Covariance and invariance under superposition of rigid body motion

We now propose to write the 3D change of observers defined by Eq. 1 using the 4D formalism to see the implications of the covariance principle. Such a change of observer is described with :

$$ \left\{ \begin{array}{l} \tilde{x}^i = Q^i_j x^j + \delta^i \\ \tilde{x}^4 = x^4 = ct \end{array} \right. \qquad (6) $$

This corresponds to a change of coordinates (from $x^\mu$ to $\tilde{x}^\mu$) of the 4D space-time. The 4D change of coordinates of Eq. 11 is then equivalent to the 3D change of frame represented by Eq. 1. Because the mapping of the events with the 4D coordinate system includes time, there is thus no mathematical difference between a change of 4D coordinate systems and a change of frame of reference (as opposed to the 3D approach). In other words, any change of 3D observers can be represented by a 4D change of coordinates. The Jacobian matrix of this change of coordinates is then :

$$ \frac{\partial \tilde{x}^\mu}{\partial x^\nu} = \begin{pmatrix} Q^i_j & \frac{1}{c} (\dot{\delta}^i + \dot{Q}^i_j x^j) \\ 0, 0, 0 & 1 \end{pmatrix} \qquad (7) $$

where $\dot{\delta}^i \equiv \frac{d\delta^i}{dt}$ and $\dot{Q}^i_j \equiv \frac{dQ^i_j}{dt}$. The factor $1/c$ comes from the differentiation with respect to $x^4 = ct$. It should be noted that $\tilde{x}^\mu$ is a curvilinear coordinate system. Then, with $T^{\mu\nu}$ and $\tilde{T}^{\mu\nu}$ the contravariant components of the second rank 4D tensor $T$ in the respective coordinate systems, it can be written (applying Eq. ?? with Eq. 7) :

$$ \left\{ \begin{aligned} \tilde{T}^{ij} &= Q^i_k Q^j_l T^{kl} + \frac{1}{c} (\dot{\delta}^i + \dot{Q}^i_m x^m) Q^j_l T^{4l} + \frac{1}{c} (\dot{\delta}^j + \dot{Q}^j_m x^m) Q^i_k T^{k4} \\ &\quad + \frac{1}{c^2} (\dot{\delta}^i + \dot{Q}^i_m x^m) (\dot{\delta}^j + \dot{Q}^j_m x^m) T^{44} \\ \tilde{T}^{4i} &= Q^i_j T^{4j} + \frac{1}{c} (\dot{\delta}^i + \dot{Q}^i_j x^j) T^{44} \\ \tilde{T}^{44} &= T^{44} \end{aligned} \right. \qquad (8) $$

In any cases, 4D tensors verify equations such as Eq. 8 because they are frame-indifferent. This equation is equivalent to Eq. 2 only if :

$$ \frac{1}{c} (\dot{\delta}^i + \dot{Q}^i_m x^m) Q^j_l T^{4l} + \frac{1}{c} (\dot{\delta}^j + \dot{Q}^j_m x^m) Q^i_k T^{k4} + \frac{1}{c^2} (\dot{\delta}^i + \dot{Q}^i_m x^m) (\dot{\delta}^j + \dot{Q}^j_m x^m) T^{44} \to 0 \qquad (9) $$

This is possible, letting $c$ tend to infinity in our non-relativistic context and provided $T^{4i}, T^{j4}$ and $T^{44}$ are small compared to $c$. It is thus demonstrated that if an entity is a tensor of the 4D space-time and if its fourth components have an adapted limit when $c$ tends to infinity such that Eq. 9 is verified, then this quantity is invariant with respect to the superposition of rigid body motions. But the terms of Eq 9 do not always vanish, mainly because the fourth components of a tensor could be proportional to $c$ (see the example of the four-velocity in the following Section), such that a 4D tensor is frame-indifferent in any case but may or may not be invariant under rigid body motion superposition. We thus have shown that the 4D formulation enables to make a clear distinction between

- the property of frame indifference : by construction verified by a all 4D tensors with Eq. 5,

- the property of invariance under rigid body motion superposition : a tensor has to verify Eq. 9.

Because these properties are now clearly distinguished, it is necessary to clarify what is meant by "objective" : to prevent possible ambiguities, in this work we use the expressions "frame-indifference" and
"invariance under rigid body motion superposition". Using the 4D formalism, a law, and thus a constitutive model, is frame-indifferent as soon as 4D tensors are used. It is a principle a physics that has to be verified. It is nevertheless possible to construct a model that depends on the superposition of rigid body motion.

# 5 Material transformations within a 4D context

Consider the finite transformation of a material body. A particle of matter is defined as an elementary 3D volume of matter. Similarly to the 3D case, we define configurations of the matter corresponding in the space-time formulation to space-like hyper-surfaces (3D volumes) of the space-time manifold. The motion of the particles corresponds to a set of worldlines (trajectories) that spans a connected open domain of the 4D manifold. The coordinates of the events undergone by a given particle are described with the specification of bijective and differentiable functions noted $\xi^\mu$ within the inertial system $z^\mu$ and $\xi^\mu$ within a coordinate system $x^\mu$. The four-velocity of a particle of matter, within the hypotheses of this work may be defined as :

$$u^\mu = \frac{d\xi^\mu}{dt} = \left( \frac{d\xi^i}{dt}, c \right) \qquad (10)$$

The definition of the four-velocity can be understood as relative to the chosen coordinate system.

General definitions for 4D deformation and strain tensors have been proposed by Lamoureux-Brousse [31] for General Relativity applications. They are applied here to the 4D non-relativistic limit. To take a full advantage of the 4D formalism, this author has proposed to compare two different motions $\xi^i$ and $\xi^{i'}$ of the same material body. These two motions are described within two different coordinate systems respectively $x^\mu$ and $x'^\mu$, associated to the two respective covariant components of the metric tensor $g_{\mu\nu}$ and $g'_{\mu\nu}$. To quantify the difference between one motion and the other, a 4D generalization of the deformation gradient is proposed such that :

$$F_{\nu}^{\mu}(\xi^{\prime\mu}, \xi^{\nu}) \equiv \frac{\partial \xi^{\prime\mu}}{\partial \xi^{\nu}} \qquad (11)$$

The 4D generalization of the Cauchy's deformation tensor is defined as

$$b_{\mu\nu} \equiv F_{\mu}^{\alpha} F_{\nu}^{\beta} g'_{\alpha\beta} \qquad (12)$$

The second rank energy-momentum four-tensor of the matter corresponds to the 4D generalization of the 3D Cauchy stress tensor; it is possible to show that its spatial components reduce to the 3D Cauchy stress components in the non-relativistic limit [3, 29, 32]. This tensor is a second rank tensor density and can be shown to have a weight equal to one [33, 26]. The use of a tensor density to describe the stress tensor is not common in continuum mechanics, although it is mentioned in [24, 26]. It is necessary to take the weight of the tensor into account because it appears in the definitions of the derivative operators below.

# 6 Derivatives in a 4D context

The 4D formalism offers the possibility to define derivatives and time derivative operators that are 4D tensors, in other words that are frame-indifferent. The covariant derivative operator is first defined such that :

$$d\alpha^{\mu\nu} = \nabla_{\lambda}\alpha^{\mu\nu} dx^{\lambda} = \left( \frac{\partial\alpha^{\mu\nu}}{\partial x^{\lambda}} + \Gamma^{\mu}_{\kappa\lambda} \alpha^{\kappa\nu} + \Gamma^{\nu}_{\kappa\lambda} \alpha^{\mu\kappa} \right) dx^{\lambda} \qquad (13)$$

for a second rank tensor and where the $\Gamma^\nu_{\kappa\lambda}$ are the coefficients of the affine connection. A covariant rate corresponding to an increment of the tensor as seen by a point of the manifold and following the direction of the motion may be defined. It is a tensor of the 4D space-time [3]. Within the hypotheses of this work, it becomes the operator noted $\frac{D}{Dt}$ in this work, such that :

$$\lim_{t' \to t} \frac{\alpha(\xi^i, t') - \alpha(\xi^i, t)}{t' - t} = \frac{d\xi^\lambda}{dt} \nabla_\lambda \alpha = \frac{D\alpha}{Dt} \qquad (14)$$
The Lie derivative, 4D operator corresponding to the variation of the tensor as seen by the particle within its motion is then defined as :

$$ \lim_{t' \to t} \frac{\alpha \left( \xi^{\mu} + \frac{d\xi^{\mu}}{dt}(t' - t) \right) - \alpha(\xi^{\mu})}{t' - t} = \mathcal{L}_u(\alpha) \quad (15) $$

For a second rank tensor of weight W, one has

$$ \mathcal{L}_u(\alpha)^{\mu\nu} = u^{\lambda} \frac{\partial \alpha^{\mu\nu}}{\partial x^{\lambda}} - \alpha^{\lambda\nu} \frac{\partial u^{\mu}}{\partial x^{\lambda}} - \alpha^{\mu\lambda} \frac{\partial u^{\nu}}{\partial x^{\lambda}} + W \alpha^{\mu\nu} \frac{\partial u^{\lambda}}{\partial x^{\lambda}} \quad (16) $$

Both, the covariant rate and the Lie derivative are frame-indifferent. The Lie derivative can be shown to be invariant under rigid body motion superposition. The Lie derivative of the Cauchy stress tensor within an inertial 3D Cartesian frame corresponds to Truesdell's transport. Truesdell's rate is thus the only objective transport that corresponds to a time derivative of the stress tensor.

Using the derivatives above it is possible to evaluate the four acceleration as the variation of the four-velocity as seen by a point of the manifold following the direction of the motion. The covariant derivative can be applied to the velocity to define the 4D rate of deformation d as :

$$ d_{\nu}^{\mu} = \frac{1}{2} (\nabla_{\nu} u^{\mu} + \nabla_{\mu} u^{\nu}) \quad (17) $$

It is also possible to define a spin tensor. It can be shown that d, the spin tensor and the four acceleration are 4D tensors and are thus *frame-indifferent* [34, 35, 30, 3].

# 7 Construction of hyperelastic and hypoelastic models

Using the tensor definitions given in the previous sections it is possible to construct isotropic hyper-elastic models such that [36] :

$$ \sigma^{\mu\nu} = \sqrt{|b_{\beta}^{\alpha}|} (a_1 b^{\kappa\lambda} (g_{\kappa\lambda} - b_{\kappa\lambda}) + a_0 \Delta T) b^{\mu\nu} + \sqrt{|b_{\beta}^{\alpha}|} a_2 b^{\mu}_{\kappa} (b^{\kappa\nu} - b^{\kappa}_{\lambda} b^{\lambda\nu}) \quad (18) $$

where the $a_i$ are scalar constants. It is possible to further detail the different terms :

- $a_0 \Delta T b^{\mu\nu}$ is the thermal expansion due to a non-isothermal transformation;

- $a_1 b^{\kappa\lambda} (g_{\kappa\lambda} - b_{\kappa\lambda}) b^{\mu\nu}$ corresponds to the "linear", homogeneous and isotropic elastic behaviour;

- $a_2 b^{\mu}_{\kappa} (b^{\kappa\nu} - b^{\kappa}_{\lambda} b^{\lambda\nu})$ are non-linear quadratic and cubic terms due to elastic behaviour.

Hypo-elastic models may also be constructed, like for example :

$$ \mathcal{L}_u(\sigma)^{\mu\nu} = 2\mu d^{\mu\nu} + \lambda d^{\kappa\lambda} g_{\kappa\lambda} g^{\mu\nu} \quad (19) $$

where the constants are here chosen to be the classical Lamé elastic constants.

We consider now a particular kinematics to illustrate how finite transformations can be handled with the chosen formalism and we compare this 4D approach with the classical 3D relations. When necessary, formal computations have been performed using the software Mathematica©. The following simple isothermal and pure shear transformation is considered

$$ \begin{pmatrix} \zeta^1 = Z^1 + \gamma(t)Z^2 \\ \zeta^2 = Z^2 \\ \zeta^3 = Z^3 \\ \zeta^4 = ct \end{pmatrix} \quad (20) $$

where $\gamma$ depends on the fourth coordinate. To respect the hypotheses of this work, the geometry and shear speed of the considered motion are such that the velocities for all the events are negligible compared to the speed of light. The results are given in Figure 1 for the hyperelastic model and Figure 2 for the hypoelastic model. In the second case, the 4D model is compared with solutions that would be obtained using classical 3D objective stress transports instead of the Lie derivative.
<figure><img src="image_1.png" /><figcaption>Fig. 1 – Hyper-elastic model. Shear stress $\sigma_{12}$ versus slip $\gamma$, for different sets of parameters. Lamé constant $a_2$ is 75GPa. The ratio $r$ describe a kinematics non-linearity ($r \equiv 1 + \frac{a_1}{a_2}$) and takes the different values listed in the caption.</figcaption></figure>

<figure><img src="image_2.png" /><figcaption>Fig. 2 – Hypo-elastic model. Shear stress $\sigma_{12}$ versus time, for different stress transports.</figcaption></figure>

# 8 Conclusion

The application of the covariance principle within the 4D differential geometry formalism ensures the frame-indifference of any physical relations. The present approach enables to obtain constitutive equations that have not necessarily a direct matching with a 3D classical relation. We have proposed a four-dimensional formulation for hyper and hypoelastic models corresponding. The general 4D differential geometry formalism thus appears to be a promising tool to compute other kind of material constitutive relations.

# Références

[1] C. Truesdell and W. Noll. *The Non-Linear Field Theories of Mechanics*. Springer, third edition, 2003.

[2] L. J. Kempers. The principle of material frame-indifference and the covariance principle. *Il Nuovo Cimento B*, 103 :227, 1989.

[3] L.D. Landau and E.M. Lifshitz. *The Classical Theory of Fields*. Elsevier, fourth edition, 1975.

[4] M. Frewer. More clarity on the concept of material frame-indifference in classical continuum mechanics. *Acta Mechanica*, 202 :213–246, 2009.

[5] P. Rougée. Kinematics of finite deformations. *Arch. Mech.*, 44 :527–556, 1992.

[6] J. Garrigues. *Fondements de la Mécanique des Milieux Continus*. Hermès Science Publications, 2007.

[7] A. I. Murdoch. On material frame-indifference, intrinsic spin, and certain constitutive relations motivated by the kinetic theory of gases. *Archive for Rational Mechanics and Analysis*, 83 :185–194, 1983.

[8] A.I. Murdoch. Some primitive concepts in continuum mechanics regarded in terms of objective space-time molecular averaging : The key role played by inertial observers. *Journal of Elasticity*, 84(1) :69–97, 2006.

[9] W. Muschik and L. Restuccia. Systematic remarks on objectivity and frame-indifference, liquid crystal theory as an example. *Arch. Appl. Mech.*, 78 :837–854, 2008.
[10] J. E. Marsden and J. R. Hugues. *Mathematical foundation of elasticity*. Prentice Hall, 1983.

[11] S. Nemat-Nasser. *Plasticity, A Treatise on Finite Deformation of Heterogeneous Inelastic Materials*. Cambridge University Press, 2004.

[12] J. Besson, G. Cailletaud, J.-L. Chaboche, and S. Forest. *Non-Linear Mechanics of Materials*. Springer, 2009.

[13] M.K. Duszek and P. Perzyna. The localization of plastic deformation in thermoplastic solids. *Int. J. Solids Struc.*, 27 :1419–1443, 1991.

[14] Z. Fiala. Time derivative obtained by applying the riemannian manifold of riemannian metrics to kinematics of continua. *C. R. Mecanique*, 332 :97–102, 2004.

[15] A. Meyers, P. Schieße, and O. T Bruhns. Some comments on objectivity rates of symmetric eulerian tensors with application to eulerian strain rates. *Acta Mechanica*, 139 :91–103, 2000.

[16] W. Muschik. Is the heat flux density really non-objective? a glance back, 40 years later. *Continuum Mech. Thermodyn.*, 2012.

[17] B. Schieck and H. Stumpf. Deformation analysis for finite elastic-plastic strains in a lagrangean-type description. *Int. J. Solids. Struc.*, 30 :2639–2660, 1993.

[18] J. C. Simo and M. Ortiz. A unified approach to finite deformation elastoplastic analysis based on the use of hyperelastic constitutive equations. *Computer Methods in Applied Mechanics and Engineering*, 49 :221–245, 1985.

[19] H. Stumpf and U. Hoppe. The application of tensor algebra on manifolds to nonlinear continuum mechanics – invited survey article. *Math. Mech.*, 77 :327–339, 1997.

[20] B. Svendsen and A. Bertram. On frame-indifference and form-invariance in constitutive theory. *Acta Mechanica*, 132 :195–207, 1999.

[21] C. K. Valanis. Elasticity of space-time : basis of newton’s 2nd law of motion. *Journal of Engineering Mechanics*, 129 :1039–1047, 2003.

[22] T. Levi-Civita, E. Persico, and Long M. *The absolute differential calculus*. Dover Phoenix Editions, 2005.

[23] J.A. Schouten. *Ricci-calculus : An Introduction to Tensor Analysis and Its Geometrical Applications*. Springer-Verlag, 1954.

[24] A. C. Eringen. *Nonlinear Theory of Continuous Media*. McGraw-Hill, 1962.

[25] I.-S. Liu. On euclidean objectivity and the principle of material frame-indifference. *Continuum Mech. Thermodyn.*, 16 :177–183, 2004.

[26] J. G. Oldroyd. On the formulation of rheological equations of state. *Proc. R. Soc. Lond. A*, 200(1063) :523–541, 1950.

[27] wikipedia. Tensor density http://en.wikipedia.org/wiki/tensor_density, 13 April 2012.

[28] P. Havas. Four-dimensional formulations of newtonian mechanics and their relation to the special and the general theory of relativity. *Reviews of Modern Physics*, 1964.

[29] S. Weinberg. *Gravitation and Cosmology : Principles and Applications of the General Theory of Relativity*. John Wiley & Sons Inc., 1972.

[30] M. Boratav and R. Kerner. *Relativité*. Ellipses, 1991.

[31] L. Lamoureux-Brousse. Infinitesimal deformations of finite conjugacies in non-linear classical or general relativistic theory of elasticity. *Physica D.*, 35 :203–219, 1989.

[32] D. N Williams. The elastic energy momentum tensor in special relativity. *Annals of Physics*, 196 :345–380, 1989.

[33] C.W. Misner, K.S. Thorne, and J.A. Wheeler. *Gravitation*. W. H. Freeman and Co Ltd., 1973.

[34] A. Bressan. *Relativistic Theories of Materials*. Springer-Verlag, Berlin, 1978.

[35] R.A. Grot and A.C. Eringen. Relativistic continuum mechanics part i mechanics and thermodynamics. *Int. J. Engng Sci*, 4 :611–638, 1966.

[36] B. Panicaud and E. Rouhaud. A frame-indifferent model for a thermo-elastic material - beyond the three-dimensional eulerian and lagrangean descriptions. *Journal of Continuum Mechanics and Thermodynamics*, 2013.