Chapter 1: Mathematical Foundations: 1.1-Tensors and Continuum Mechanics
Chapter 1: Mathematical Foundations: 1.1-Tensors and Continuum Mechanics
A tensor is a generalization of scalars, vectors and matrices and defined as n-dimensional array.
A Tensor may be classified in a particular coordinate system by a certain set of quantities, known as
its components. However, specified the components of a tensor in one coordinates system
determines the components in any other coordinates using the law of transformation.
Now why we use tensors, suppose we have a parallelepiped, where we have a force acting on the
center of gravity, this force have 3 components, x, y and z, then we cut the parallelepiped in the 3
dimensions (along xy, xz, yz) on each face of the cut we made, we have 3 faces, 2 which are shear
force (trying to tear the face) and 1 compression or tension.
On the left image you can see the cut, and on the right image you can see the forces acting on each
face, now we have a matrix of forces instead of vector, since these forces on each face have a
different effect on the object deformation. 2
The cartesian tensor uses the orthonormal system to be represented, otherwise it is general.
1
online reference of this information
2
a video to get the intuition of tensors
Therefore:
• tensors of zero order are specified by one component only and are called scalers (30 = 1)
• tensors of order one have 3 components and known as vectors (31 = 3)
𝑣1
⃗𝑉 = 𝑣1 𝑖 + 𝑣2 𝑗 + 𝑣3 𝑘⃗ , 𝑉
⃗ [𝑣2 ]
𝑣3
σ11 σ12 σ13
• 2
tensors of order 2 have 9 components (3 = 9) such as 3 × 3 matrix σ = [ 21 σ22 σ23 ]
σ
σ31 σ32 σ33
Tensors of order two (second order tensors), and higher order tensors are usually used in continuum
mechanics with the number of components is 3𝑁 in 2 (or 3) dimensional system (xoy system/ oxyz
system)
You can think of the rank as the number of basis vector needed to fully specify a component of the
tensor, for example a scaler which is represented by a tensor of rank 0, we need 0 basis vectors, also
a 3d vector is represented by a 3 basis vector, thus a tensor of rank 1.
The rank is defined by 3𝑛 , where n is the actual rank, and the number of basis vectors needed (if we
have 27 component, the rank of the tensor needed is thus 3 -since 27 = 33 , the tensor is not a
matrix.3
1.3- Indicial Notation: (The ideas below are not related to the definition of
tensor; the indicial notation is also called Einstein notation)
The components of a tensor of any order may be simplified and represented by the use of indicial
notation. meaning that the orthonormal base unit vectors {𝑖, 𝑗, 𝑘⃗ } are given in the form {𝑥
⃗⃗⃗⃗1 , ⃗⃗⃗⃗ 𝑥3 }
𝑥2 , ⃗⃗⃗⃗
In this notation, letter indices: may be used as subscripts (𝑥 ⃗⃗⃗𝑖 ). Typical examples
⃗⃗⃗𝑖 ), or superscripts (𝑥
3
a video to get the intuition of tensors
are: 𝑎𝑖 , 𝑇𝑖𝑗 , ∈𝑖𝑗𝑘 , 𝑏 𝑖 , 𝑅 𝑝𝑞 , with 𝑖, 𝑗, 𝑘, 𝑝, 𝑞 = 1, 2, 3 (1 for scalars, {1, 2} for 2 dimentional systems, {1,
2, 3} for 3 dimentional system)
The general idea of indicing is to shorthand long equations for example 𝑎11 + 𝑎22 + 𝑎33 may be
written 𝑎𝑖𝑖
i. Free Indices:
When an index occures unrepeated in a tensor, (appears only once), that index takes the
values of 1, 2, 3, … N (N range of index).
These indices are known as free indices
example:
𝑣1
⃗ |𝑣2 , 𝑣𝑖 represents one of these components
𝑉
𝑣3
• 𝑎𝑖 + 𝑏𝑖 represents one of the follwing 3 qunatities 𝑎1 + 𝑏1 𝑜𝑟 𝑎2 + 𝑏2 𝑜𝑟 𝑎3 + 𝑏3
• 𝑎𝑖𝑗 represents one of the follwing 9 terms
𝑎11 , 𝑎12, 𝑎13, 𝑎21, 𝑎22, 𝑎23, 𝑎31, 𝑎32, 𝑎33
The tensorial rank (order) of a give term is equal to the number of free indices appearing in
that term
A term is defined as a group of multiplied elements like a or ab, even if both can be assigned
different superscripts or subscripts. As long as the index is used once per term it is
considered a free index.4 The free index represents one value of a set of possible values (like
given in the example above).5
4
the first page of this document
5
a video on index notation
Example:
𝑎1
• 𝑎𝑖 + 𝑏𝑖𝑖 = (𝑎2 ) + (𝑏11 + 𝑏22 + 𝑏33 )
𝑎3 ↓
𝑠𝑐𝑎𝑙𝑎𝑟
𝑣𝑒𝑐𝑡𝑜𝑟
• 𝑎𝑖 . 𝑏𝑖𝑖
↓ ↓
𝑅 𝐷
It is possible to represent the dummy index by any letter and this will not change the value of the
tensor so we can write: 𝑎𝑖𝑖 = 𝑎𝑗𝑗 = 𝑎𝑘𝑘 where 𝑖 = 1 → 3
𝑗=1→3
𝑘=1→3
• 𝑥𝑖 = 𝐶𝑖 𝑗 . 𝑍𝑗
↓ ↓ ↓ ↓
𝐹 𝐹𝐷 𝐷
𝑥1 = 𝑐11 𝑧1 + 𝑐12 𝑧2 + 𝑐13 𝑧3 𝑥1 𝑐11 𝑐12 𝑐13 𝑧1
𝑥2 = 𝑐21 𝑧1 + 𝑐22 𝑧2 + 𝑐23 𝑧3 (𝑥2 ) = (𝑐21 𝑐22 𝑐23 ) (𝑧2 )
𝑥3 = 𝑐31 𝑧1 + 𝑐32 𝑧2 + 𝑐33 𝑧3 𝑥3 𝑐31 𝑐32 𝑐33 𝑧3
• 𝐴 𝑖 𝑗 = 𝐵i 𝑝 𝐶𝑖 𝑗 . 𝑍pq 𝑖, 𝑗 = 1,2
↓ ↓ ↓ ↓ ↓
𝐹𝐹 𝐷 𝐹𝐷 𝑝, 𝑞
𝐴11 = 𝐵11 𝐶11 𝐷11 + 𝐵11 𝐶12 𝐷12 + 𝐵12 𝐶11 𝐷21 + 𝐵12 𝐶12 𝐷22
𝐴12 = 𝐵11 𝐶21 𝐷11 + 𝐵11 𝐶22 𝐷12 + 𝐵12 𝐶21 𝐷21 + 𝐵12 𝐶22 𝐷22
𝐴21 =
𝐴21 =
(The last two values are left empty by the doctor)
𝛅𝒊𝒋 examples:
o 𝐶2 = 𝑎3 𝑏1 − 𝑎1 𝑏3
o 𝐶3 = 𝑎1 𝑏2 − 𝑎2 𝑏1
𝑖 𝑗 𝑘⃗
also 𝑎 ∧ 𝑏 = |𝑎1 𝑎2 𝑎3 | = (𝑎2 𝑏3 − 𝑎3 𝑏2 )𝑖 − (𝑎1 𝑏3 − 𝑎3 𝑏1 )𝑗 + (𝑎1 𝑏2 − 𝑎2 𝑏1 )𝑘⃗
⃗
𝑏1 𝑏2 𝑏3
• The box product 𝑎 ∧ 𝑏⃗. 𝑐 may be writer as:
𝑎 ∧ 𝑏⃗. 𝑐 = λ ⇒ ϵ𝑖𝑗𝑘 𝑎𝑖 𝑏𝑗 𝑐𝑘 = λ i, j, k= 1, 2, 3
𝑎1 𝑎2 𝑎3
⃗
𝑎 ∧ 𝑏. 𝑐 = |𝑏1 𝑏2 𝑏3 | = 𝑎1 (𝑏2 𝑐3 − 𝑐2 𝑏3 )
𝑐1 𝑐2 𝑐3
• The dual vector of an arbitrary second order tensor 𝑇𝑖𝑗 is defined by:
o 𝑣𝑖 = ϵ𝑖𝑗𝑘 𝑇𝑗𝑘
o 𝑣1 = 𝜖111 𝑇11 + 𝜖112 𝑇12 + 𝜖113 𝑇13 + 𝜖121 𝑇21 + 𝜖122 𝑇22 + 𝜖123 𝑇23 + 𝜖131 𝑇31 +
𝜖132 𝑇32 + 𝜖133 𝑇33 = 𝑇23 − 𝑇32
o 𝑣2 =
o 𝑣3 =
(the last two elements were left blank by the doctor)
Suppose for this explanation that 𝑓(𝑥1 , 𝑥2 , 𝑥3 ) = (𝑓1 (𝑥1 , 𝑥2 , 𝑥3 ), 𝑓2 (𝑥1 , 𝑥2 , 𝑥3 )), in this case “f” is a
vector-valued function 𝐑𝟑 → 𝐑𝟐
In the subscript of a function “f” (the place of 2 in 𝜕2 for example), we can introduce the comma “,”
where numbers before it represent the specific sub-function (for this example it is 𝑓1 & 𝑓2, and the
number after it represent the variable, we are deriving with respect to (in this example we have
𝑥1 , 𝑥2 , 𝑥3 ), if we put 2 numbers (as you will see in doctor’s examples) after the comma, then we have
double derivatives with respect to the variables indicated by the numbers (if you see the examples it
will be clear)
∂𝑓 ∂𝑓 ∂𝑓
→ = 𝑓,1 , = 𝑓,2 , … , = 𝑓,𝑖
∂𝑥1 ∂𝑥2 ∂𝑥𝑖
∂𝑓 ∂𝑓 ∂𝑓
→ 2 = 𝑓,11 , 2 = 𝑓,22 , … , 2 = 𝑓,𝑖𝑖
∂𝑥1 ∂𝑥1 ∂𝑥𝑖
∂𝑓
→ = 𝑓,𝑖𝑗
∂𝑥𝑖 𝑦𝑗
∂𝑓𝑖 ∂𝑓𝑖
→ = 𝑓𝑖,𝑗 , = 𝑓𝑖,𝑖
∂𝑥𝑗 ∂𝑥𝑖
since “𝑖” is found two times in the same term it is a dummy index, thus it represents a summation
rather than selection
1.7.1- Definitions: below you will find representation of some formulas using index notation
i. Divergence:
𝜕 𝜕 𝜕
⃗ ) = ∇. 𝑉(𝑥, 𝑦, 𝑧) =
𝑑𝑖𝑣(𝑉 (𝑉𝑥 ) + (𝑉𝑦 ) + (𝑉𝑧 )
𝜕𝑥 𝜕𝑦 𝜕𝑧
this the divergence formula you saw in Analysis II
∂𝑉
𝑑𝑖𝑣(𝑉⃗ ) = 𝑖 = 𝑉𝑖,𝑖
∂𝑥𝑖
∂𝑉1 ∂𝑉2 ∂𝑉3 ⃗
∂𝑉 ∂𝑉⃗ ⃗
∂𝑉
⃗)=
𝑑𝑖𝑣(𝑉 + + = ⃗⃗⃗
𝑉,𝑖 . ⃗𝑒⃗𝑖 = . 𝑒⃗⃗⃗1 + . 𝑒⃗⃗⃗2 + . 𝑒⃗⃗⃗
∂𝑥1 ∂𝑥2 ∂𝑥3 ∂𝑥1 ∂𝑥2 ∂𝑥3 3
where the dummy index 𝑖 = 1,2,3
ii. Gradient:
𝜕𝑓
ۊ 𝑥𝜕ۇ
𝜕𝑓
⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ (𝑓(𝑥, 𝑦, 𝑧)) = ۋ ۈ
𝑔𝑟𝑎𝑑 ۋ𝑦𝜕ۈ
ۋ ۈ
𝜕𝑓
ی 𝑧𝜕ۉ
this the gradient formula you saw in Analysis II
𝜕𝑓 𝜕𝑓 𝜕𝑓
⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝑔𝑟𝑎𝑑𝑓 = 𝑓,𝑖 ⃗𝑒⃗𝑖 = 𝑒⃗⃗⃗1 + 𝑒⃗⃗⃗2 + 𝑒⃗⃗⃗
𝜕𝑥1 𝜕𝑥2 𝜕𝑥3 3
iii. Rotation of ⃗𝑽:
∂
⃗)=∇
𝑅𝑜𝑡(𝑉 ⃗ ∧𝑉
⃗ ⃗ = ⃗𝑒⃗𝑖
where ∇ = ⃗𝑒⃗𝑖 ∂𝑖
∂𝑥𝑖
⃗
∂𝑉
= ⃗𝑒⃗𝑖 ∧ = ⃗𝑒⃗𝑖 ∧ ⃗⃗⃗⃗⃗
𝑉,𝑖
∂𝑥𝑖
⃗⃗⃗⃗⃗ ∧ = 𝑒⃗⃗⃗1
⃗⃗⃗⃗⃗ ∧ 𝑉,1 + 𝑒⃗⃗⃗2
⃗⃗⃗⃗⃗ ∧ 𝑉,2 + 𝑒⃗⃗⃗3
𝑉,3
𝑖 𝑗 ⃗
𝑘 ⃗𝑘 𝑖 𝑖 𝑗 ⃗𝑘
𝑗
⃗⃗⃗⃗⃗ ∧ 𝑖⃗𝑒⃗
| = 𝑖𝑉, 1 0 0 0 |+| 0
||+ 0 0 | 1
1
∂𝑉1 ∂𝑉2 ∂𝑉3 ∂𝑉3 ∂𝑉1 ∂𝑉2
∂𝑉1 ∂𝑉2 ∂𝑉3
∂𝑥1 ∂𝑥1 ∂𝑥1 ∂𝑥2 ∂𝑥2 ∂𝑥2
∂𝑥3 ∂𝑥3 ∂𝑥3
𝜕𝑉3 𝜕𝑉2
0 𝜕𝑉3 𝜕𝑉2 −
𝜕𝑉3 − ۇۊ 𝜕𝑥2 𝜕𝑥3
−ۇ ۊ 𝜕𝑥3ۇ ۊ 𝜕𝑥2ۇ ۊ ۋ 𝜕𝑉3 𝜕𝑉3
ۈ = ۋ 𝜕𝑉3ۈ +ۋ 0ۈ +ۋ 𝜕𝑥1ۈ =ۋ 𝑥𝜕 𝜕𝑥 −ۈ
𝜕𝑉2 𝜕𝑉2 𝜕𝑥1 1ۈ ۋ1
− 𝜕𝑉2 𝜕𝑉2
ی 0ۉ ی 𝜕𝑥1ۉ ی 𝜕𝑥1ۉ −
ی 𝜕𝑥1 𝜕𝑥1ۉ
Or:
𝜕𝑉3 𝜕𝑉2
−
𝑖 𝑗 ⃗𝑘 ۊ 𝜕𝑥2 𝜕𝑥3ۇ
∂ ۋ 𝜕𝑉3 𝜕𝑉3ۈ | ∂
𝑖⃗𝑒⃗ ∂ || = ⃗
𝑉∧
∂
ۈ=| −
𝑖𝑥∂ ۋ
∂𝑥1 ∂𝑥2 ∂𝑥3 ۋ 𝜕𝑥1 𝜕𝑥1ۈ
V_1 V_2 𝑉_3 𝜕𝑉2 𝜕𝑉2
−
ی 𝜕𝑥1 𝜕𝑥1ۉ
𝑉∂
𝑉 ∧ ∇⃗ = ) ⃗
𝑉(i.e., Rot ⃗⃗⃗⃗⃗ ∧ 𝑖⃗𝑒⃗ = ⃗
𝑗𝑉,𝑖 = ⃗𝑒⃗𝑖 ∂𝑖 ∧ 𝑣 = ϵ𝑖𝑗𝑘 ∂𝑗 𝑣𝑘 = ϵ𝑖𝑗𝑘 ∂𝑥𝑘 = ϵ𝑖𝑗𝑘 𝑉𝑘,
𝑖
1
𝑖⃗𝑒⃗ ∧ 𝑖⃗⃗𝑒) 𝑗= (𝑉𝑗,𝑖 − 𝑉𝑖,
2
In a given region of space, the vector function of position 𝐹 = 𝐹(𝑥) is defined at every
point on the smooth curve C. If the differential tangent vector to the curve at an arbitrary
⃗⃗⃗⃗ , the integral:
point P is 𝑑𝑥
𝑥𝑏
∫𝐹 . ⃗⃗⃗⃗
𝑑𝑥 = ∫ 𝐹 . ⃗⃗⃗⃗
𝑑𝑥
𝑐 𝑥𝑎
Where 𝑛⃗ is the outward unit normal on the position side of S, and dS is the differential
element of the surface S
In indicial form
i.e., the circulation of the vector 𝐹 along the curve C is equal to the rotational of 𝐹 over the
whole surface S limited by the curve C.
Example:
𝐹1 (𝑥1 , 𝑥) 𝑑𝑥1 0
The vector 𝐹 = |𝐹2 (𝑥1 , 𝑥2 ) ⃗⃗⃗⃗
𝑑𝑥 |𝑑𝑥2 𝑛⃗ |0
𝐹3 (𝑥1 , 𝑥2 ) ⏟ 0 1
The surface S is in the plane xoy
𝐹 . ⃗⃗⃗⃗
𝑑𝑥 = 𝐹1 𝑑𝑥1 + 𝐹2 𝑑𝑥2 + 𝐹3 𝑑𝑥3 = 𝐹1 𝑑𝑥1 + 𝐹2 𝑑𝑥2
𝜕𝑉3 𝜕𝑉2
−
𝑥𝜕ۇ2 𝜕𝑥3 ۊ
𝜕𝑉 𝜕𝑉
⃗⃗⃗⃗⃗⃗⃗ (𝐹 ) = ۈ3 − 3 ۋ
𝑅𝑜𝑡 𝑥𝜕ۈ ۋ
ۈ1 𝜕𝑥1 ۋ
𝜕𝑉2 𝜕𝑉2
−
𝑥𝜕ۉ1 𝜕𝑥1 ی
Applying Stoke’s theorem:
∂𝑓2 ∂𝑓1
∫ (𝐹1 𝑑𝑥1 + 𝐹2 𝑑𝑥2 ) = ∫ ( − ) 𝑑𝑥1 𝑑𝑥2
𝐶 𝑆 ∂𝑥1 ∂𝑥2