User:Luc.stpierre

The concept of objectivity in science means that qualitative and quantitative descriptions of physical phenomena remain unchanged independently of the point of view from which they are observed. All physical processes (e.g. material properties) must be invariant under changes of observers.

Euclidean tranformation
Physical processes can be described by an observer denoted by $$ O $$. In Euclidean three-dimensional space and time, an observer can measure relative positions of points in space and intervals of time.

Consider an event in Euclidean space characterized by the pairs $$ (x_0,t_0) $$ and $$ (x,t) $$ where $$x$$ is a position vector and $$t$$ is a scalar representing time. This pair is mapped to another one denoted by the $$*$$ superscript. This mapping is done with the orthogonal time-dependent second order tensor $$ Q(t) $$ in a way such that the distance between the pairs is kept the same. Therefore one can write:


 * $$\ x^*-x_0^*=Q(t)(x-x_0). $$

By introducing a vector $$ C(t) $$ and a real number $$ \alpha $$ denoting the time shift, the relationship between $$ x $$ and $$ x^* $$ can be expressed


 * $$\ x^*=c(t)+Q(t)x \quad \text{where} \quad c(t)=x_0^*-Q(t)x_0 \quad \text{and} \quad \alpha=t^*-t=t_0^*-t_0. $$

The one-to-one mapping connection of the pair $$ (x,t) $$ with its corresponding pair $$ (x^*, t^*) $$ is referred to as an Euclidean transformation.

Displacement
A physical quantity like displacement should be invariant relative to a change of observer. Consider one event recorded by two observers; for $$ O$$, point $$ x$$ moves to position $$y$$ whereas for $$ O^* $$, the same point $$ x^*$$ moves to $$y^*$$. For $$ O$$, the displacement is $$ u=y-x $$. On the other hand, for $$ O^*$$, one can write:


 * $$\ \begin{align} u^* &= y^*-x^* \\

&=C(t)+Q(t)y-C(t)-Q(t)x \\ &=Q(t)(y-x) \\ & =Q(t)u. \end{align} $$

Any spatial vector field $$ u $$ that transforms such that:


 * $$\ u^*=Q(t)u, $$

is said to be objective, since $$ |u^*|=|u| $$.

Velocity
Because $$ Q(t) $$ is a rotation matrix, $$ Q(t)^T Q(t) = I $$ where $$ I $$ is the identity matrix. Using this relation, the inverse of the Euclidean transformation can be written as:


 * $$\ x = Q(t)^T[x^*-c(t)]. $$

The velocity can be obtained by differentiating the above expression:


 * $$\ v(x,t) = \dot{x} = \dot{Q}(t)^T[x^*-c(t)] + Q(t)^T[v^*-\dot{c}(t)]. $$

By reorganizing the terms in the above equation, one can obtain:


 * $$\ \begin{align} v^*(x^*,t) & = Q(t)v + \dot{c}(t) - Q(t)\dot{Q}(t)^T[x^*-c(t)]\\

& = Q(t)v + \dot{c}(t) + \Omega(t)[x^*-c(t)], \end{align}$$ where


 * $$\ \Omega(t) = \dot{Q}(t)Q(t) = -\Omega(t)^T = -Q(t)\dot{Q}(t)^T, $$

is a skew tensor representing the spin of the reference frame of observer $$ O $$ relative to the reference frame of observer $$ O^* $$ (Holzapfel 2000). To simplify the mathematical notation, the arguments of functions will no longer be written.

From the above expression, one can conclude that velocity is not objective because of the presence of the extra terms $$ \dot{c} $$ and $$ \Omega[x^*-c] $$. Nevertheless, the velocity field can be made objective by constraining the change of observer to:


 * $$\ \dot{c} + \Omega(x^*-c) = 0, $$

A time-independant rigid transformation such as:


 * $$\ x^* = c_0 + Q_0x \quad \text{where} \quad \dot{Q}_0 = 0 \quad \text{and} \quad \dot{Q_0}=0,$$

respects this condition.

Acceleration
The material time derivative of the spatial velocity $$v$$ returns the spatial acceleration $$a$$. By differentiating the transformation law for the spatial velocity, one can obtain:


 * $$\ a^* = v^* = \dot{Q}v + Qa + \ddot{c} + \dot{\Omega}(x^*-c) + \Omega(v-\dot{c}), $$

which can be rewritten as follow:


 * $$\ a^* = Qa + \ddot{c} + (\dot{\Omega}-\Omega^2)(x^*-c) + 2\Omega(v-\dot{c}). $$

Just like the spatial velocity, the acceleration is not an objective quantity for a general change of observer (Holzapfel 2000). As for the spatial velocity, the acceleration can also be made objective by constraining the change of observer. One possibility would be to use the time-independent rigid transformation introduced above.

Objectivity for higher-order tensor fields
A tensor field of order $$ n $$ and denoted $$ u_1 \otimes \dots \otimes u_n $$ is objective if, during a general change of observer, the transformation is given by:


 * $$\ (u_1 \otimes \dots \otimes u_n)^* = Qu_1 \otimes \dots \otimes Qu_n. $$

Example for a second order tensor
Introducing a second order tensor $$ A = u_1 \otimes u_2 $$, one can find with the above definition of objectivity that:


 * $$\ A^* = (u_1 \otimes u_2)^* = Qu_1 \otimes Qu_2 = Q(u_1 \otimes u_2)Q^T = QAQ^T.$$

Example for a scalar field
The general condition of objectivity for a tensor of order $$ n $$ can be applied to a scalar field $$ \Phi $$ for which $$ n=0 $$. The transformation would give:


 * $$\ \Phi ^* = \Phi. $$

Physically, this means that a scalar field is independent of the observer. Temperature is an example of scalar field and it is easy to understand that the temperature at a given point in a room and at a given time would have the same value for any observer.

Deformation gradient
The deformation gradient at point $$ x $$ and at its associated point $$ x^* $$ is a second order tensor given by:


 * $$\ F = \frac{\partial x}{\partial X} \qquad \textrm{and} \qquad F^* = \frac{\partial x^*}{\partial X}, $$

where $$ X $$ represents the material coordinates. Using the chain rule, one can write:


 * $$\ F^* = \frac{\partial x^*}{\partial x} \frac{\partial x}{\partial X} = QF. $$

From the above equation, one can conclude that the deformation gradient $$ F $$ is objective even though it transforms like a vector and not like a second order tensor. This is because one index of the tensor describes the material coordinates $$ X $$ which are independent of the observer (Holzapfel 2000).

Cauchy stress tensor
The Cauchy traction vector $$ t $$ is related to the Cauchy stress tensor $$ \sigma $$ at a given point $$ x $$ by the outward normal to the surface $$ n $$ such that: $$ t = \sigma n $$. The Cauchy traction vector for another observer can be written simply $$ t^*=\sigma ^* n^* $$, where $$ t $$ and $$ n $$ are both objective vectors. Knowing that, one can write:


 * $$\ \begin{array}{rrcl} & t^* &=& \sigma ^* n^*\\

\Rightarrow & Qt & = & \sigma ^* Qn \\ \Rightarrow & Q\sigma n &=& \sigma ^* Qn \\ \Rightarrow & \sigma ^* &=& Q\sigma Q^T. \end{array}$$

This demonstrates that the Cauchy stress tensor is objective.

Piola-Kirchhoff stress tensors
The first Piola-Kirchhoff stress tensor $$ P $$ is defined as:


 * $$\ PF^T = J\sigma, $$

where $$ J=\det(F) $$. It is also interesting to know that since $$ Q $$ is a rotation matrix:


 * $$\ J^* = \det(F^*) = \det(QF) = \det(Q)\det(F) = \det(F) = J. $$

Using identities developed previously, one can write:


 * $$\ \begin{array}{rrcl} & P^*(F^*)^T & = & J^*\sigma ^* \\

\Rightarrow & P^*(QF)^T & = & J Q \sigma Q^T \\ \Rightarrow & P^*F^T Q^T & = & Q J \sigma Q^T \\ \Rightarrow & P^*F^T Q^T & = & Q P F^T Q^T \\ \Rightarrow & P^* & = & Q P. \end{array}$$

This proves that the first Piola-Kirchhoff stress tensor is objective. Similarly to the deformation gradient, this second order tensor transforms like a vector.

The second Piola-Kirchhoff stress tensor $$ S=F^{-1}P $$ is also objective and transforms like a scalar field. This can be easily demonstrated:


 * $$\ S^* = (F^*)^{-1}P^* = (QF)^{-1}QP = F^{-1}Q^{-1}QP = F^{-1}P = S. $$

The three stress tensors studied here, $$ \sigma $$, $$ P $$ and $$ S $$, were all found to be objective therefore they are all suitable to describe the material response and develop constitutive laws because they are independent of the observer.

Objectivity rates
It was shown above that even if a displacement field is objective, the velocity field is not. An objective vector $$ u^*=Qu $$ and an objective tensor $$ A^* = QAQ^T $$ usually do not conserve their objectivity through time differentiation as demonstrated below:


 * $$\ u^* = \dot{Q}u + Q\dot{u} \quad \text{and} \quad A^* = \dot{Q}AQ^T + Q\dot{A}Q^T + QA\dot{Q}^T.$$

Objectivity rates are modified material derivatives that allows to have an objective time differentiation. Before presenting some examples of ojectivity rates, certain other quantities need to be introduced. First, the spatial velocity gradient $$l$$ is defined as:


 * $$\ l = \dot{F}F^{-1} = d+w,$$

where $$d$$ is a symmetric tensor and $$w$$ is a skew tensor called the spin tensor. For a given $$l$$, $$d$$ and $$w$$ are uniquely defined. The Euclidean transformation for the spatial velocity gradient can be written as:


 * $$\ \begin{align} l^* & = \dot{F}^*(F^*)^{-1} \\

& = (\dot{Q}F+Q\dot{F})(QF)^{-1} \\ & = (\dot{Q}F+Q\dot{F})F^{-1}Q^T \\ & = \dot{Q}FF^{-1}Q^T + Q\dot{F}F^{-1}Q^{-1} \\ & = \dot{Q}Q^T + QlQ^{-1} \\ & = \Omega + QlQ^{-1}. \end{align}$$

Substituting $$l=d+w$$ in the above equation, one can obtain the two following relations:


 * $$\ \dot{Q} = w^*Q - Qw \quad \text{and} \quad \dot{Q}^T = -Q^Tw^* + wQ^T $$

Substituting the above result in the previously obtained equation for the rate of an objective vector one can write:


 * $$\ \begin{array}{rrcl} & u^* &=& \dot{Q}u + Q\dot{u} \\

\Rightarrow & u^* &=& (w^*Q - Qw)u + Q\dot{u} \\ \Rightarrow & u^* &=& w^*u^* - Qwu + Q\dot{u} \\ \Rightarrow & (u-wu)^* &=& Q(\dot{u}-wu) \\ \Rightarrow & \bar{u}^* &=& Q\bar{u}, \end{array}$$

where the co-rotational rate of the objective vector field $$u$$ is defined as:
 * $$\ \bar{u}=\dot{u}-wu, $$

and represents an objective quantity. Similarly, using the above equations, one can obtain the co-rotational rate of the objective second-order tensor field $$A$$:


 * $$\ \begin{array}{rrcl} & (\dot{A}-wA+Aw)^* &=& Q(\dot{A}-wA+Aw)Q^T \\

\Rightarrow & \bar{A}^* &=& Q\bar{A}^*AQ^T. \end{array}$$

This co-rotational rate second order tensor is defined as:


 * $$\ \bar{A} = \dot{A}-wA+Aw. $$

This objective rate is know as the Jaumann-Zaremba rate and it is often used in plasticity theory. Many different objective rates can be developed and the reader is referred to Holzapfel 2000 for more details.

Invariance of material response
The principal of material invariance basically means that the material properties are independent of the observer. In this section it will be shown how this principle adds constraints to constitutive laws.

Cauchy-elastic materials
A Cauchy-elastic material is one for which the stress field at a given time depends only on the current state of deformation (Holzapfel 2000). In other words, the material is independent of the deformation path and consequently, also independent of time.

Neglecting the effect of temperature and assuming the body to be homogeneous, a constitutive equation for the Cauchy stress tensor can be formulated based on the deformation gradient:


 * $$\ \sigma = G(F).$$

This constitutive equation for another arbitrary observer can be written $$ \sigma^* = G(F^*) $$. Knowing that the Cauchy stress tensor $$ \sigma $$ and the deformation gradient $$ F $$ are objective quantities, one can write:


 * $$\ \begin{array}{rrcl} & \sigma^* &=& G(F^*) \\

\Rightarrow & Q\sigma Q^T &=& G(QF) \\ \Rightarrow & QG(F) Q^T &=& G(QF). \end{array}$$

The above is a condition that the constitutive law $$ G $$ has to respect to make sure that the response of the material will be independent of the observer. Similar conditions can be derived for constitutive laws relating the deformation gradient to the first or second Piola-Kirchhoff stress tensors.

Isotropic Cauchy-elastic materials
Here, it will be assumed that the Cauchy stress tensor $$ \sigma $$ is a function of the left Cauchy-Green tensor $$ b=FF^T $$. The constitutive equation may be written:


 * $$\ \sigma = h(b). $$

In order to find the restriction on $$ h $$ that will ensure that the principle of material frame-indifference, one can write:


 * $$\ \begin{array}{rrcl} & \sigma^* &=& h(b^*) \\

\Rightarrow & Q \sigma Q^T &=& h(F^*(F^*)^T) \\ \Rightarrow & Q h(b) Q^T &=& h(QFF^TQ^T) \\ \Rightarrow & Q h(b) Q^T &=& h(QbQ^T). \end{array}$$

A constitutive equation that respects the above condition is said to be isotropic (Holzapfel 2000). Physically, this characteristic means that the material has no preferential direction. Wood and most fibre-reinforced composites are generally stronger in the direction of their fibres therefore they are not isotropic materials (they are qualified as anisotropic).