11. Principal Axis and Invariants

Initialization

Needs["TensorCalculus4`Tensorial`"]

Needs["TContinuumMechanics2`TContinuumMechanics`"]

Needs["DrawGraphics`DrawingMaster`"]

numequ = 1 ;

oldflavors = IndexFlavors ;

ClearIndexFlavor/@oldflavors ;

DeclareIndexFlavor[{black, Black}, {red, Red}, {green, ForestGreen}, {star, SuperStar}, {blue, Blue}, {hat, OverHat}, {tilde, OverTilde}, {bar, OverBar}]

Base2d = {1, 2} ;

Base3d = {1, 2, 3} ;

DeclareBaseIndices[Base3d] ; {NDim, BaseIndices} ;

TensorLabelFormat[tb, OverBar[t]]

TensorLabelFormat[σb, OverBar[σ]]

DefineTensorShortcuts[{{d, dF, , , x}, 1}, {{g, t, tb, β, δ, σ, σb, ℰ, b}, 2}] ;

Principal Direction

We look to area elements such that the force dF_ i^i==d_ i^i σ_ (i  j)^(i  j) acting on it is normal to the element. This means that the covariant components of the force

res = dF_j^j == d_i^i σ_ (ij)^(ij) == σ d_j^j

dF_j^j == d_i^i σ_ (ij)^(ij) == σ d_j^j

res[[2]] == (res[[3]]/.j→i) g_ (ij)^(ij)

res1 = %//Simplify

d_i^i σ_ (ij)^(ij) == σ d_i^i g_ (ij)^(ij)

d_i^i (σ g_ (ij)^(ij) - σ_ (ij)^(ij)) == 0

A non-trivial solution exists if,

(*11.1*)ResultFrame[(res1101 = S == Det[res1[[1, 2]]] == det ( res1[[1, 2]]//ArrayExpansion[red @ i, red @ j]//Kronecker[g]//MatrixForm) == 0)]

This equation has three eigenvalues (one is necessarily real, but the other two may be either real or complex conjugate) σ_ (m)={σ_ (1),σ_ (2),σ_ (3)}associated to three area elements (d_m) _ i^i=={(d_m) _ 1^1,(d_m) _ 2^2,(d_m) _ 3^3}

Unsymmetric tensor

If we consider a general second-order tensor t,

res = t_ (ij)^(ij) v_i^i == T v_j^j

res[[1]] == (res[[2]]/.j→i) g_ (ij)^(ij) ; (*11.2*)

ResultFrame[res1102 = %//Simplify]

t_ (ij)^(ij) v_i^i == T v_j^j

       (T g_ (ij)^(ij) - t_ (ij)^(ij)) v_i^i == 0      (11.2)

res2 = Collect[-Det[ res1102[[1, 1]]//ArrayExpansion[red @ i, red @ j]//Kronecker[g]], T] == 0

or

which is called the characteristic equation of the tensor t. If we change the frame _ i^i, the three eigenvalues of (12.2) remain unchanged. This suppose that the coefficients are themselves invariants. Hence we have here three invariants built with t :

(*11.04*)ResultFrame[(res1104 = {t1[i_] = tI == t_ (ii)^(ii), t2[i_, j_] = tII == (-t_ (ij)^(ij) t_ (ji)^(ji) + t_ (ii)^(ii) t_ (jj)^(jj)), t3[i_, j_] = tIII == Det[t_ (ij)^(ij)]})//TableForm]

The invariants  tI,tII, and   tIII are respectively of order 1, 2 and  3 in t.  All combinaitions of  tI,tII, and   tIII are also invariants. An important invariant of degree 2 is 2 tI^2-3 tII, which can be expressed simply as a function of the zero trace tensor Overscript[t, _].

(*11.05*)ResultFrame[res1105 = Overscript[t, _] _ (ij)^(ij) == t_ (ij)^(ij) - 1/3 g_ (ij)^(ij) t_ (mm)^(mm)]

      Overscript[t, _] _ (ij)^(ij) == t_ (ij)^(ij) - 1/3 g_ (ij)^(ij) t_ (mm)^(mm)       (11.5)

(*11.06*)ResultFrame[res1106 = 2t1[i][[1]] t1[j][[1]] - 3t2[i, j][[1]] == 3Overscript[t, _] _ (ij)^(ij) Overscript[t, _] _ (ji)^(ji)//ToFlavor[red]]

      2 tI^2 - 3 tII == 3 Overscript[t, _] _ (ij)^(ij) Overscript[t, _] _ (ji)^(ji)       (11.6)

Proof :

tb§_ (i_j_)^(i_j_) := Module[{m}, t_ (ij)^(ij) - 1/3 g_ (ij)^(ij) t_ (mm)^(mm)//ToFlavor[Head[i]]]

2t1[i][[1]] t1[j][[1]] - 3t2[i, j][[1]] == (2t1[i][[2]] t1[j][[2]] - 3t2[i, j][[2]]//Simplify) == 3Overscript[t, _] _ (ij)^(ij) Overscript[t, _] _ (ji)^(ji)//ToFlavor[red]

%/.tb→tb§ ;

(EinsteinSum[]/@%[[{2, 3}]])//Kronecker[g]//FullSimplify

2 tI^2 - 3 tII == 3 t_ (ij)^(ij) t_ (ji)^(ji) - t_ (ii)^(ii) t_ (jj)^(jj) == 3 Overscript[t, _] _ (ij)^(ij) Overscript[t, _] _ (ji)^(ji)

True

For a symmetrical tensor t the eigenvectors (v_m) _ i^i and (v_n) _ i^i associated to two different eigenvalues T_mand T_nare orthogonal. We write

Tensor[t, {red[i], Void}, {Void, red[j]}]*Tensor[Subscript[v,"(m)"], {Void},
   {red[i]}] == Subscript[T,"(m)"]*Tensor[Subscript[v,"(m)"], {Void}, {red[j]}]
Tensor[t, {red[j], Void}, {Void, red[l]}]*Tensor[Subscript[v,"(n)"], {Void},
   {red[j]}] == Subscript[T,"(n)"]*Tensor[Subscript[v,"(n)"], {Void}, {red[l]}]

t_ (ij)^(ij) v_(m)_i^i == T_(m) v_(m)_j^j

t_ (jl)^(jl) v_(n)_j^j == T_(n) v_(n)_l^l

the second equation is multiplied by g_ (i  l)^(i  l)

g_ (il)^(il) t_ (jl)^(jl) v_"(n)"_j^j == g_ (il)^(il) T_"(n)" v_"(n)"_l^l//MetricSimplify[g]

%//UpDownSwap[red @ j]

t_ (ji)^(ji) v_(n)_j^j == T_(n) v_(n)_i^i

t_ (ji)^(ji) v_(n)_j^j == T_(n) v_(n)_i^i

and we multiply the first by v_(n)_j^j, and the new (second) equation by v_(m)_i^i

t_ (ij)^(ij) v_"(m)"_i^iv_"(n)"_j^j == T_"(m)" v_"(m)"_j^jv_"(n)"_j^j

t_ (ji)^(ji) v_"(n)"_j^j v_"(m)"_i^i == T_"(n)" v_"(n)"_i^i v_"(m)"_i^i

t_ (ij)^(ij) v_(m)_i^i v_(n)_j^j == T_(m) v_(m)_j^j v_(n)_j^j

t_ (ji)^(ji) v_(m)_i^i v_(n)_j^j == T_(n) v_(m)_i^i v_(n)_i^i

The left hand side are equal if and only if t_ (i  j)^(i  j)==t_ (j  i)^(j  i), that is to say if t_ (i  j)^(i  j)is symmetrical.

(-t_ (ji)^(ji) + t_ (ij)^(ij)) v_(m)_i^i v_(n)_j^j == (T_(m)- T_(n)) v_(m)_i^i v_(n)_i^i

In this case the right hand side are equal and

(*11.07*)ResultFrame[res1107 = v_"(m)"_i^i v_"(n)"_i^i == 0]

      v_(m)_i^i v_(n)_i^i == 0      (11.7)

if  m ≠ n. If the equation (11.2) has a double root, it is always possible to choose two orthogonal eigenvectors in the degenerate subspace. It can also be easily shown that all the eigenvalues of a symmetrical tensor must be real.

New base vectors

The eigenvectors orthogonal _(m),

_"(m)" == Tensor[v_"(m)", {red @ i}, {Void}] d[red @ i]

_(m)== _i^i v_(m)_i^i

can be chosen as a new base (blue for instance),

d[blue @ m] == βud[red @ i, blue @ m] d[red @ i]

_m^m == _i^i β_ (im)^(im)

and we notice that their contravariant components(v_m) _ i^iare identical to the basis transformation β_ (i  m)^(i  m).
From (12.5) we find,

"for m ≠ n :"

res1107/.v_"(m)"_i^i→g_ (ij)^(ij) v_"(m)"_j^j

%/.v_"(m)"_i_^i_→β_ (im)^(im)/.v_"(n)"_i_^i_→β_ (in)^(in)

%//KroneckerAbsorb[β]

for m ≠ n :

g_ (ij)^(ij) v_(m)_j^j v_(n)_i^i == 0

g_ (ij)^(ij) β_ (in)^(in) β_ (jm)^(jm) == 0

g_ (nm)^(nm) == 0

In addition, the blue basis can be taken orthonormal, if we choose the eigenvectors _ m^m as unit vectors :

gdd[blue @ m, blue @ n] == δdd[blue @ m, blue @ n]

g_ (mn)^(mn) == δ_ (mn)^(mn)

The reference frame is cartesian. The tensor t is diagonal in this basis.

tud[red @ i, red @ j] Tensor[v_"(m)", {Void}, {red @ i}] == T_mTensor[v_"(m)", {Void}, {red @ j}]

(%//UpDownSwap[red @ i])/.Tensor[v_"(m)", {Void}, {red @ j}] →Tensor[v_"(m)", {red @ i}, {Void}] gdd[red @ i, red @ j]

%/.Tensor[v_"(m)", {red @ i_}, {Void}] →βud[red @ i, blue @ m]

Times[#, βud[red @ j, blue @ n]] &/@%

%//KroneckerAbsorb[β]

t_ (ij)^(ij) v_(m)_i^i == T_m v_(m)_j^j

t_ (ij)^(ij) v_(m)_i^i == T_m g_ (ij)^(ij) v_(m)_i^i

t_ (ij)^(ij) β_ (im)^(im) == T_m g_ (ij)^(ij) β_ (im)^(im)

t_ (ij)^(ij) β_ (im)^(im) β_ (jn)^(jn) == T_m g_ (ij)^(ij) β_ (im)^(im) β_ (jn)^(jn)

t_ (mn)^(mn) == T_m g_ (mn)^(mn)

On the diagonal of the associated matrix, we have the eigenvalues T_m.

Tensors of Stress and Strain

We come back to the area element with normal d_m,

d_m == Tensor[d_m, {Void}, {red[j]}] u[red @ j]

d_m == _j^j d_m_j^j

pointing in the principal direction. The force is given by,

d == Tensor[d_m, {Void}, {red[i]}] u[red @ j] σud[red[i], red[j]] == σ_m  d_m

d == _j^j σ_ (ij)^(ij) d_m_i^i == d_m σ_m

σ_m is the principal stress. It is a normal stress, root of

res1101[[2]] == Det[σuu[red @ i, red @ j] - σ guu[red @ i, red @ j]] == 0

Det[σ g_ (ij)^(ij) - σ_ (ij)^(ij)] == Det[-σ g_ (ij)^(ij) + σ_ (ij)^(ij)] == 0

To a general stress field σ_ (i  j)^(i  j)(x_ k^k) corresponds its stress trajectories (tangent to the σ_m) to which we can associate a curvilinear coordinate system _ m^m.   

Two invariants σI and σII built as shown above, have already been encountered in chapter 4, equations (4.14),

s == 1/3σud[red @ i, red @ i] == σI/3

s == σ_ (ii)^(ii)/3 == σI/3

and (4.21), using (12.4) for σ,

rul = Solve[res1105//ToFlavor[red], t_ (ij)^(ij)][[1]]/.{t→σ, tb→σb}

{σ_ (ij)^(ij) →1/3 (g_ (ij)^(ij) σ_ (mm)^(mm) + 3 Overscript[σ, _] _ (ij)^(ij))}

(rul/.{i→j, j→i})

{σ_ (ji)^(ji) →1/3 (g_ (ji)^(ji) σ_ (mm)^(mm) + 3 Overscript[σ, _] _ (ji)^(ji))}

The strain tensor ℰ_ (i  j)^(i  j):
It is also a symmetric tensor which has three principal axis. From (4.11 or 12) applied to the reference frame _ m^m, we see that in an isotropic elastic medium, the principal axis of stress and strain coincide.

res412 = σ_ (ip)^(ip) == 2 μ ℰ_ (ip)^(ip) + λ g_ (ip)^(ip) ℰ_ (mm)^(mm)

σ_ (ip)^(ip) == 2 μ ℰ_ (ip)^(ip) + λ g_ (ip)^(ip) ℰ_ (mm)^(mm)

Curvature

We have considered in Chapter 8, the curvature tensor b_ (α  β)^(α  β), and its two invariants the mean curvature b_ (α  α)^(α  α), and the Gaussian curvature {b_ (α  β)^(α  β)}(corresponding to the two invariants bI and bIII in the above notation). The invariant bII appears to be a combition of the two invariants bI and bIII .
b_ (α  β)^(α  β)is symmetrical, and has two orthogonal principal directions _ λ^λfor which b_ (1  2)^(1  2)==b_ (2  1)^(2  1)==0. b_ (1  1)^(1  1)and b_ (2  2)^(2  2)are the principal curvatures. The lines of the curvature field are called the lines of principal curvature. They form an orthogonal set. This orthonormal set is not necessarily the best to use due to the boundary conditions which in general do not coincide with these lines.

ClearTensorShortcuts[{{d, dF, , , x}, 1}, {{g, t, tb, β, δ, σ, σb, ℰ, b}, 2}] ;


Created by Mathematica  (November 27, 2007) Valid XHTML 1.1!