# A CAS Approach to Handle the Anisotropic Hooke’s Law for Cancellous Bone and Wood

Authors: Sandeep Kumar
Journal: Chinese Journal of Engineering
 (2014)

---

## Abstract

The present research entirely relies on the Computer Algebric Systems (CAS) to develop techniques for the data analysis of the sets of elastic constant data measurements. In particular, this study deals with the development of some appropriate programming codes that favor the data analysis of known values of elastic constants for cancellous bone, hardwoods, and softwood species. More precisely, a “Mathematica” code, which has an ability to unfold a fourth-order elasticity tensor is discussed. Also, an effort towards the fabrication of an appropriate “MAPLE” code has been exposed, that can calculate not only the eigenvalues and eigenvectors for cancellous bone, hardwoods, and softwood species, but also computes the nominal average of eigenvectors, average eigenvectors, average eigenvalues, and the average elasticity matrices for these materials. Further, using such a MAPLE code, the histograms corresponding to average elasticity matrices of 15 hardwood species have been plotted and the graphs for I, II, III, IV, V, and VI eigenvalues of each hardwood species against their apparent densities are also drawn.


---

## Body

## 1. Introduction

The study of material symmetry of 3-dimensional space is of great interest due to having crucial theoretical as well as practical significance. This is because a symmetrical space includes crystals and all homogeneous fields without exceptions: electric, magnetic, gravitational, and so forth.The variation of material properties with respect to direction at a stagnant point in a material is called material symmetry; for instance, if the material properties are same in all directions at some fixed point, they are called isotropic, whereas if the material properties show variation at the same point, they are called anisotropic [1]. Of course, the familiarity with material symmetries is the best way to categorize the materials. However, according to [2], many materials are anisotropic and inhomogeneous due to the varying composition of their constituents. In such materials, it becomes ticklish to identify the symmetries or more particularly, the elastic symmetries. The variable composition method to identify material’s elastic symmetry becomes complicated and hence to overcome this difficulty, an approach is developed by [3], called “averaging anisotropic elastic constant data.” In this approach, the identification of elastic symmetries and method of variable composition are analyzed separately. With the aid of this fabulous approach, a sheer volume of research towards material symmetry has been put forward by various researchers, for example, estimation of the effective transversely isotropic elastic constants of a material from known values of the material’s orthotropic elastic constants [4], anisotropic Hooke’s law for cancellous bone and wood [2], the validation of generalized Hooke’s law for coronary arteries [5], constitutive relationships of fabric, density, and elastic properties in cancellous bone architecture [6], analysis of elastic anisotropy of wood material for engineering applications [7], a multidimensional anisotropic strength criterion based on Kelvin’s modes [8], spectral decomposition of a 4th-order covariance tensor: application to diffusion tensor MRI [9], a normal distribution of tensor valued random variables: applications to diffusion tensor MRI [10], and so forth.Concerning computer assisted analysis of bone material, a few of the articles, like direct mechanical assessment of elastic symmetries and properties of trabecular bone architecture [11], relationships between morphology and bone elastic properties can be accurately quantified using high-resolution computer reconstructions [12], and intrinsic mechanical properties of trabecular calcaneus determined by finite element models using 3D synchrotron microtomography [13], have explored the subject material up to some great extent.Now, as far as the present proposed research concerned, we turn our concentration towards cancellous bone and wood (soft and hard). According to [4], the type of material symmetry exposed by bone tissue is often ambiguous. So, in order to attain accuracy, we have to measure the elastic constants of bone specimen, which may be represented as orthotropic. However, the orthotropy of the bone specimen may or may not be of higher degree in all the directions with respect to a fixed point in the specimen. Thereby, for all practical purposes the bone specimen can be represented as transversely isotropic or even isotropic.Now, the question arises is what is the need of studying the mechanical properties of cancellous bone and wood? The rigorous answer is that the knowledge about mechanical properties of cancellous bone is essential for the determination of bone fracture risk in osteoporosis and other pathological conditions involving impaired bone strength [6]. Also, the mechanical properties of cancellous bone are determined by the properties of its bone tissue and its architecture [4, 14–21]. Further, the wood is a cellulosic, semicrystalline, cellular material and the mechanical properties (e.g., elasticity, strength and rheology, etc.) are generally higher along the bole of a tree than across the bole [7].Mechanically, clear wood obeys the law of elastic orthotropic material like the bone does. In general the woody substance also exposes the high toughness and stiffness properties and these properties vary according to the type of wood and the direction in which the woody substance is examined, as the woody substance shows a higher degree of anisotropy.In addition, to deal with anisotropic Hooke’s law, one should be familiar with the algebra of the 4th-order tensor. We would like to emphasize that the 4th-order tensor algebra is not only involved in the study of material symmetry, but its crucial appearance has been set up in the field of diffusion tensor MRI [9, 10], visualization and processing of tensor fields [22], biomechanics [13], tissue mechanics [1], geophysics [23], and so forth.Though an adequate amount of research work has been dedicated towards the algebra of 4th-order tensors, then also, in the following section, we shall present a brief digest regarding this issue, as this issue hits the present study up to some great level.


## 2. The Algebra of Fourth-Order Tensor

The splendid word “Tensor” has been derived from the Latin word “Tensus”, which is the past participle of “tenděre” and stands for “Stretch.” This word was used in anatomy in the early 1704 to represent muscle and stretches. However, in Mathematics, it existed in 1846 when William Rowan Hamilton has explored his quaternion Algebra. Even though the Hamilton sense regarding tensor did not survive. The recent meaning of tensor is due to Voigt, who used the term “Tensortriple” in crystal elasticity around 1899.Likewise the second-order tensor, a fourth-order tensor𝔸 (say), is a linear map from Lin
V to Lin
V. That is, it assigns to each second-order tensor A the second-order tensor 𝔸
U; that is,
(1)
𝔸
:
Lin
(
V
)
⟶
Lin
(
V
)
,
U
⟼
𝔸
U
,
∀
U
∈
Lin
(
V
)
,
where V is a vector space and Lin
(
V
) is a linear space of dimension n
2 and is equivalent to
(2)
Lin
(
V
)
≑
A
:
V
⟶
V
,
u
⟼
A
u
,
∀
u
∈
V
.
Also, any 4th-order tensor 𝔸 has the representation
(3)
𝔸
=
A
i
j
k
l
e
i
⊗
e
j
⊗
e
k
⊗
e
l
,
where
(4)
A
i
j
k
l
=
〈
e
i
⊗
e
j
,
𝔸
(
e
k
⊗
e
k
)
〉
,
and the set of all four vectors, that is, the set of all linear maps from Lin
(
V
) to Lin
(
V
) is a linear space of dimension n
4 and is denoted by 𝕃
in
(
V
).The outer (or open) product of the two 4th-order tensors is given by the composition(5)
(
𝔸
𝔹
)
U
=
𝔸
(
𝔹
U
)
,
∀
U
∈
Lin
(
V
)
.
In abstract index notations, we have
(6)
(
𝔸
𝔹
)
i
j
k
l
=
A
i
j
p
q
B
p
p
q
k
l
.
Moreover, the transpose of 𝔸
∈
𝕃
in
(
V
) is a 4th-order tensor 𝔸
T and is defined by
(7)
〈
U
,
𝔸
T
V
〉
=
〈
V
,
𝔸
U
〉
,
∀
U
,
V
∈
Lin
(
V
)
.
On 𝕃
in
(
V
), we can define the inner product as
(8)
〈
(
A
)
,
𝔹
〉
≑
tr
⁡
(
𝔸
T
𝔹
)
=
A
i
j
k
l
B
i
j
k
l
.
The induced norm is
(9)
∥
𝔸
∥
=
[
tr
⁡
(
𝔸
T
𝔸
)
]
1
/
2
,
and the metric is
(10)
d
E
(
𝔸
,
𝔹
)
=
∥
𝔸
-
𝔹
∥
.
Hence, with respect to the inner product given by (8), the set {
e
i
⊗
e
j
⊗
e
k
⊗
e
l
}
{
1
≤
i
,
j
,
k
,
l
≤
n
} forms an orthonormal basis of 𝕃
in
(
V
).


### 2.1. Symmetries of a 4th-Order Tensor

In case of a second-order tensor, the only known symmetry represents an invariance under the mutual rotation of two indices, while for the 4th-order tensor, there are several notions of symmetry that represent invariance under exchanging pair of indices. Thus, for the 4th-order tensor generally we have three types of symmetries, namely, major, minor, and total symmetries, and these notions of symmetries widely play an important role in the theory of elasticity [24].


#### 2.1.1. Major Symmetry of a Fourth-Rank Tensor

For𝔸
∈
𝕃
in
(
V
), we say that 𝔸 is symmetric if 𝔸
=
𝔸
T, or in component form A
i
j
k
l
=
A
k
l
i
j and we say that 𝔸 is skew-symmetric if 𝔸
=
-
𝔸
T or A
i
j
k
l
=
-
A
k
l
i
j.Moreover, we have the decomposition(11)
𝔸
=
1
2
(
𝔸
+
𝔸
T
)
+
1
2
(
𝔸
-
𝔸
T
)
.
Then the set of all symmetric fourth-order tensors is
(12)
𝕊
ym
(
V
)
=
{
𝔸
∈
𝕃
in
(
V
)
∣
𝔸
=
𝔸
T
}
is a subspace of 𝕃
in
(
V
) of dimension n
2
(
n
2
+
1
)
/
2.In the theory of elasticity, this symmetry is often evoked as “major symmetry.” Since the major symmetry represents invariance under exchanging the pair(
i
,
j
) and (
k
,
l
), then if 𝔸
∈
𝕊
ym
(
V
), then ∃
{
λ
i
}
i
=
1,2
,
…
,
n
2
∈
ℝ and {
U
i
}
i
=
1,2
,
…
,
n
2 is an orthonormal basis of Lin
(
V
), such that
(13)
𝔸
=
∑
i
=
1
n
2
‍
λ
i
U
i
⊗
U
i
.
The real number λ
i is called the eigenvalues of 𝔸 associated with the eigentensor U
i.Such a representation is called the spectral decomposition of𝔸. The trace and the determinant of the fourth-order tensor are defined, respectively, as
(14)
tr
⁡
(
𝔸
)
=
∑
i
=
1
n
2
‍
λ
i
,
det
⁡
(
𝔸
)
=
∏
i
=
1
n
2
‍
λ
i
.


#### 2.1.2. Minor Symmetry

The second type of symmetry is called “minor symmetry” and is defined by(15)
〈
U
,
𝔸
V
〉
=
〈
U
T
,
𝔸
V
T
〉
,
∀
U
,
V
∈
Lin
(
V
)
.
In component form, A
i
j
k
l
=
A
j
i
k
l
=
A
i
j
l
k, 1
≤
i
,
j
,
k
,
l
≤
n.Now, the invariance under the exchange of first pair of the indices is called the “first minor symmetry” and the invariance under the exchange of second pair of indices is called “second minor symmetry.”The set of all 4th-order tensors that satisfy the minor symmetry is denoted by(16)
𝕊
ym
(
V
)
:
=
{
𝔸
∈
𝕃
in
(
V
)
∣
𝔸
satisfies
minor
symmetry
}
.


#### 2.1.3. Total Symmetry

A fourth-order tensor𝔸 is said to have total symmetry if in addition to satisfying the major and minor symmetries, it also satisfies
(17)
𝔸
u
⊗
U
⊗
v
=
𝔸
u
⊗
U
T
⊗
v
,
∀
U
∈
Lin
(
V
)
,
u
,
v
∈
V
,
or in abstract index notations, A
i
j
k
l
=
A
σ
(
i
)
σ
(
j
)
σ
(
j
)
σ
(
l
) for some permutation σ of [
1
:
n
].The set of the fourth-rank tensors bearing total symmetry is denoted by(18)
𝕊
y
m
total
(
V
)
=
{
𝔸
∈
𝕃
in
(
V
)
∣
𝔸
satisfies
total
symmetry
}
.
Thus any fourth-order tensor can be decomposed into its totally symmetric part 𝔸
s and its totally antisymmetric part 𝔸
a, that is, 𝔸
=
𝔸
s
+
𝔸
a.The components of totally symmetric and antisymmetric parts of a fourth-order tensor are described as [25]
(19)
A
i
j
k
l
s
=
1
3
(
A
i
j
k
l
+
A
i
k
j
l
+
A
i
l
k
j
)
,
A
i
j
k
l
a
=
1
3
(
A
i
j
k
l
-
A
i
k
j
l
-
A
i
l
j
k
)
.
Obviously, a fourth-order tensor is totally symmetric if 𝔸
=
𝔸
s or equivalently 𝔸
a
=
𝕆, where 𝕆 stands for a fourth-order null tensor. Likewise 𝔸 is antisymmetric if 𝔸
=
𝔸
a or equivalently 𝔸
s
=
𝕆.Now, it is straightforward to show that if𝔸 is totally symmetric and 𝔹 is a symmetric tensor, then tr
⁡
(
𝔸
𝔹
)
=
0.Reference [25] has also shown that there is an isomorphism (i.e., a 1-1 linear map) between totally symmetric fourth-rank tensor and a homogeneous polynomial of degree four. Also, if a totally symmetric fourth-order tensor is traceless, that is, tr
⁡
(
𝔸
)
=
0, then it is isomorphic to a harmonic polynomial of degree four. For this reason, a totally symmetric and traceless tensor is called harmonic tensor.To meet the objectives of proposed research, let us turn our attention to the flattening (sometimes called unfolding) of the fourth-order 3-dimensional tensor with minor symmetries, as it is customary in the 3-dimensional elasticity theory.In the following section, we discuss the notion of 9-dimensional representation of a 3-dimensional fourth-order tensor and then particularize this representation to a fourth-order elasticity tensor having minor symmetries due to the well-known anisotropic Hooke’s law.


## 2.1. Symmetries of a 4th-Order Tensor

In case of a second-order tensor, the only known symmetry represents an invariance under the mutual rotation of two indices, while for the 4th-order tensor, there are several notions of symmetry that represent invariance under exchanging pair of indices. Thus, for the 4th-order tensor generally we have three types of symmetries, namely, major, minor, and total symmetries, and these notions of symmetries widely play an important role in the theory of elasticity [24].


### 2.1.1. Major Symmetry of a Fourth-Rank Tensor

For𝔸
∈
𝕃
in
(
V
), we say that 𝔸 is symmetric if 𝔸
=
𝔸
T, or in component form A
i
j
k
l
=
A
k
l
i
j and we say that 𝔸 is skew-symmetric if 𝔸
=
-
𝔸
T or A
i
j
k
l
=
-
A
k
l
i
j.Moreover, we have the decomposition(11)
𝔸
=
1
2
(
𝔸
+
𝔸
T
)
+
1
2
(
𝔸
-
𝔸
T
)
.
Then the set of all symmetric fourth-order tensors is
(12)
𝕊
ym
(
V
)
=
{
𝔸
∈
𝕃
in
(
V
)
∣
𝔸
=
𝔸
T
}
is a subspace of 𝕃
in
(
V
) of dimension n
2
(
n
2
+
1
)
/
2.In the theory of elasticity, this symmetry is often evoked as “major symmetry.” Since the major symmetry represents invariance under exchanging the pair(
i
,
j
) and (
k
,
l
), then if 𝔸
∈
𝕊
ym
(
V
), then ∃
{
λ
i
}
i
=
1,2
,
…
,
n
2
∈
ℝ and {
U
i
}
i
=
1,2
,
…
,
n
2 is an orthonormal basis of Lin
(
V
), such that
(13)
𝔸
=
∑
i
=
1
n
2
‍
λ
i
U
i
⊗
U
i
.
The real number λ
i is called the eigenvalues of 𝔸 associated with the eigentensor U
i.Such a representation is called the spectral decomposition of𝔸. The trace and the determinant of the fourth-order tensor are defined, respectively, as
(14)
tr
⁡
(
𝔸
)
=
∑
i
=
1
n
2
‍
λ
i
,
det
⁡
(
𝔸
)
=
∏
i
=
1
n
2
‍
λ
i
.


### 2.1.2. Minor Symmetry

The second type of symmetry is called “minor symmetry” and is defined by(15)
〈
U
,
𝔸
V
〉
=
〈
U
T
,
𝔸
V
T
〉
,
∀
U
,
V
∈
Lin
(
V
)
.
In component form, A
i
j
k
l
=
A
j
i
k
l
=
A
i
j
l
k, 1
≤
i
,
j
,
k
,
l
≤
n.Now, the invariance under the exchange of first pair of the indices is called the “first minor symmetry” and the invariance under the exchange of second pair of indices is called “second minor symmetry.”The set of all 4th-order tensors that satisfy the minor symmetry is denoted by(16)
𝕊
ym
(
V
)
:
=
{
𝔸
∈
𝕃
in
(
V
)
∣
𝔸
satisfies
minor
symmetry
}
.


### 2.1.3. Total Symmetry

A fourth-order tensor𝔸 is said to have total symmetry if in addition to satisfying the major and minor symmetries, it also satisfies
(17)
𝔸
u
⊗
U
⊗
v
=
𝔸
u
⊗
U
T
⊗
v
,
∀
U
∈
Lin
(
V
)
,
u
,
v
∈
V
,
or in abstract index notations, A
i
j
k
l
=
A
σ
(
i
)
σ
(
j
)
σ
(
j
)
σ
(
l
) for some permutation σ of [
1
:
n
].The set of the fourth-rank tensors bearing total symmetry is denoted by(18)
𝕊
y
m
total
(
V
)
=
{
𝔸
∈
𝕃
in
(
V
)
∣
𝔸
satisfies
total
symmetry
}
.
Thus any fourth-order tensor can be decomposed into its totally symmetric part 𝔸
s and its totally antisymmetric part 𝔸
a, that is, 𝔸
=
𝔸
s
+
𝔸
a.The components of totally symmetric and antisymmetric parts of a fourth-order tensor are described as [25]
(19)
A
i
j
k
l
s
=
1
3
(
A
i
j
k
l
+
A
i
k
j
l
+
A
i
l
k
j
)
,
A
i
j
k
l
a
=
1
3
(
A
i
j
k
l
-
A
i
k
j
l
-
A
i
l
j
k
)
.
Obviously, a fourth-order tensor is totally symmetric if 𝔸
=
𝔸
s or equivalently 𝔸
a
=
𝕆, where 𝕆 stands for a fourth-order null tensor. Likewise 𝔸 is antisymmetric if 𝔸
=
𝔸
a or equivalently 𝔸
s
=
𝕆.Now, it is straightforward to show that if𝔸 is totally symmetric and 𝔹 is a symmetric tensor, then tr
⁡
(
𝔸
𝔹
)
=
0.Reference [25] has also shown that there is an isomorphism (i.e., a 1-1 linear map) between totally symmetric fourth-rank tensor and a homogeneous polynomial of degree four. Also, if a totally symmetric fourth-order tensor is traceless, that is, tr
⁡
(
𝔸
)
=
0, then it is isomorphic to a harmonic polynomial of degree four. For this reason, a totally symmetric and traceless tensor is called harmonic tensor.To meet the objectives of proposed research, let us turn our attention to the flattening (sometimes called unfolding) of the fourth-order 3-dimensional tensor with minor symmetries, as it is customary in the 3-dimensional elasticity theory.In the following section, we discuss the notion of 9-dimensional representation of a 3-dimensional fourth-order tensor and then particularize this representation to a fourth-order elasticity tensor having minor symmetries due to the well-known anisotropic Hooke’s law.


## 2.1.1. Major Symmetry of a Fourth-Rank Tensor

For𝔸
∈
𝕃
in
(
V
), we say that 𝔸 is symmetric if 𝔸
=
𝔸
T, or in component form A
i
j
k
l
=
A
k
l
i
j and we say that 𝔸 is skew-symmetric if 𝔸
=
-
𝔸
T or A
i
j
k
l
=
-
A
k
l
i
j.Moreover, we have the decomposition(11)
𝔸
=
1
2
(
𝔸
+
𝔸
T
)
+
1
2
(
𝔸
-
𝔸
T
)
.
Then the set of all symmetric fourth-order tensors is
(12)
𝕊
ym
(
V
)
=
{
𝔸
∈
𝕃
in
(
V
)
∣
𝔸
=
𝔸
T
}
is a subspace of 𝕃
in
(
V
) of dimension n
2
(
n
2
+
1
)
/
2.In the theory of elasticity, this symmetry is often evoked as “major symmetry.” Since the major symmetry represents invariance under exchanging the pair(
i
,
j
) and (
k
,
l
), then if 𝔸
∈
𝕊
ym
(
V
), then ∃
{
λ
i
}
i
=
1,2
,
…
,
n
2
∈
ℝ and {
U
i
}
i
=
1,2
,
…
,
n
2 is an orthonormal basis of Lin
(
V
), such that
(13)
𝔸
=
∑
i
=
1
n
2
‍
λ
i
U
i
⊗
U
i
.
The real number λ
i is called the eigenvalues of 𝔸 associated with the eigentensor U
i.Such a representation is called the spectral decomposition of𝔸. The trace and the determinant of the fourth-order tensor are defined, respectively, as
(14)
tr
⁡
(
𝔸
)
=
∑
i
=
1
n
2
‍
λ
i
,
det
⁡
(
𝔸
)
=
∏
i
=
1
n
2
‍
λ
i
.


## 2.1.2. Minor Symmetry

The second type of symmetry is called “minor symmetry” and is defined by(15)
〈
U
,
𝔸
V
〉
=
〈
U
T
,
𝔸
V
T
〉
,
∀
U
,
V
∈
Lin
(
V
)
.
In component form, A
i
j
k
l
=
A
j
i
k
l
=
A
i
j
l
k, 1
≤
i
,
j
,
k
,
l
≤
n.Now, the invariance under the exchange of first pair of the indices is called the “first minor symmetry” and the invariance under the exchange of second pair of indices is called “second minor symmetry.”The set of all 4th-order tensors that satisfy the minor symmetry is denoted by(16)
𝕊
ym
(
V
)
:
=
{
𝔸
∈
𝕃
in
(
V
)
∣
𝔸
satisfies
minor
symmetry
}
.


## 2.1.3. Total Symmetry

A fourth-order tensor𝔸 is said to have total symmetry if in addition to satisfying the major and minor symmetries, it also satisfies
(17)
𝔸
u
⊗
U
⊗
v
=
𝔸
u
⊗
U
T
⊗
v
,
∀
U
∈
Lin
(
V
)
,
u
,
v
∈
V
,
or in abstract index notations, A
i
j
k
l
=
A
σ
(
i
)
σ
(
j
)
σ
(
j
)
σ
(
l
) for some permutation σ of [
1
:
n
].The set of the fourth-rank tensors bearing total symmetry is denoted by(18)
𝕊
y
m
total
(
V
)
=
{
𝔸
∈
𝕃
in
(
V
)
∣
𝔸
satisfies
total
symmetry
}
.
Thus any fourth-order tensor can be decomposed into its totally symmetric part 𝔸
s and its totally antisymmetric part 𝔸
a, that is, 𝔸
=
𝔸
s
+
𝔸
a.The components of totally symmetric and antisymmetric parts of a fourth-order tensor are described as [25]
(19)
A
i
j
k
l
s
=
1
3
(
A
i
j
k
l
+
A
i
k
j
l
+
A
i
l
k
j
)
,
A
i
j
k
l
a
=
1
3
(
A
i
j
k
l
-
A
i
k
j
l
-
A
i
l
j
k
)
.
Obviously, a fourth-order tensor is totally symmetric if 𝔸
=
𝔸
s or equivalently 𝔸
a
=
𝕆, where 𝕆 stands for a fourth-order null tensor. Likewise 𝔸 is antisymmetric if 𝔸
=
𝔸
a or equivalently 𝔸
s
=
𝕆.Now, it is straightforward to show that if𝔸 is totally symmetric and 𝔹 is a symmetric tensor, then tr
⁡
(
𝔸
𝔹
)
=
0.Reference [25] has also shown that there is an isomorphism (i.e., a 1-1 linear map) between totally symmetric fourth-rank tensor and a homogeneous polynomial of degree four. Also, if a totally symmetric fourth-order tensor is traceless, that is, tr
⁡
(
𝔸
)
=
0, then it is isomorphic to a harmonic polynomial of degree four. For this reason, a totally symmetric and traceless tensor is called harmonic tensor.To meet the objectives of proposed research, let us turn our attention to the flattening (sometimes called unfolding) of the fourth-order 3-dimensional tensor with minor symmetries, as it is customary in the 3-dimensional elasticity theory.In the following section, we discuss the notion of 9-dimensional representation of a 3-dimensional fourth-order tensor and then particularize this representation to a fourth-order elasticity tensor having minor symmetries due to the well-known anisotropic Hooke’s law.


## 3. A 9-Dimensional Exposition of a 3-Dimensional 4th-Order Tensor and the Anisotropic Hooke’s Law

In accordance with the evolution of tensor theory, the algebra of the 3-dimensional fourth- order tensor is still not fully developed. For instance, the techniques for calculating the latent roots an latent tensors of a tensor likeC
i
j
k
l, or the computation of its inverse (which is called a compliance S
i
j
k
l), are still not widely available. There are also some certain psychological constraints; that is, we are trained to tackle matrices but not trained to deal with the multiway arrays (sometimes called matrix of matrices). Though, recently Tamra [26–29] has developed a tensor toolbox that runs under MATLAB [30], in my view this tool still needs some enhancements to handle symbolic tensor data. Moreover, Constantinescu and Korsunsky [31] have developed some crucial packages, namely, “Tensor2Analysis.m,” “IntegrateStrain.m,” “ParametricMesh.m,” and “VectorAnalysis.m,” which are compatible with Mathematica and Mathematica programming; using these packages enables one to compute symbolic computations regarding elasticity tensors. Anyways, we shall come back to these CAS assisted issues in the subsequent sections, where we shall demonstrate the proposed objectives using CAS.In order to study material symmetries and anisotropic Hooke’s law, it is customary to transform the 3-dimensional 4th-order tensor as a 9 × 9 matrix, that is, a 9-dimensional representation of a 3-dimensional fourth-order tensor. For this purpose, there is a basic notion of representing a fourth-order tensor as a second-order tensor [22]. According to [22], if 1
≤
i, j
≤
n, then we set
(20)
e
~
ψ
(
i
,
j
)
=
e
i
⊗
e
j
.
Thus, {
e
~
α
⊗
e
~
β
}
1
≤
α
,
β
≤
n
2 is an orthonormal basis of Lin
(
V
~
)
≡
𝕃
in
(
V
).Therefore, any fourth-rank tensor𝔸
∈
𝕃
in
(
V
) can be assumed as an n-dimensional fourth rank tensor 𝔸
:
=
A
i
j
k
l
e
i
⊗
e
j
⊗
e
k
⊗
e
l, or as an n
2-dimensional second-order tensor 𝔸
↔
A
~
=
A
~
α
β
e
α
⊗
e
β, where A
~
ψ
(
i
,
j
)
ψ
k
,
l
=
A
i
j
k
l, 1
≤
i
,
j
,
k
,
l
≤
n.Again, we have from (8) that
(21)
tr
⁡
(
𝔸
T
𝔹
)
=
〈
𝔸
,
𝔹
〉
=
tr
⁡
(
A
~
T
B
~
T
)
,
∥
𝔸
∥
=
∥
A
~
∥
;
hence the two-way map 𝔸
↔
A
~ is an isometry.Let us apply this concept to anisotropic Hooke’s law by assuming that the anisotropic Hooke’s law given below is generalized one and also valid in the case where stress and strain tensors are not necessarily symmetric.The anisotropic Hooke’s law in abstract index notations is often depicted as(22)
σ
i
j
=
C
i
j
k
l
ϵ
k
l
,
or
in
index
free
notations
σ
=
C
ϵ
,
where C
i
j
k
l are the components of an elastic tensor. Writing the stress, strain, and elastic tensors in usual tensor bases, we have
(23)
σ
=
σ
i
j
e
i
⊗
e
j
,
ϵ
=
ϵ
k
l
e
k
⊗
e
l
,
C
=
C
i
j
k
l
e
i
⊗
e
k
⊗
e
k
⊗
e
l
.
Now, when working with orthonormal basis, one needs to introduce a new basis which should be composed of the three diagonal elements
(24)
E
1
≡
e
1
⊗
e
1
,
E
2
≡
e
2
⊗
e
2
,
E
3
≡
e
3
⊗
e
3
,
the three symmetric elements
(25)
E
4
≡
1
2
(
e
1
⊗
e
2
+
e
2
⊗
e
1
)
,
E
5
≡
1
2
(
e
2
⊗
e
3
+
e
3
⊗
e
2
)
,
E
6
≡
1
2
(
e
3
⊗
e
1
+
e
1
⊗
e
3
)
,
and the three asymmetric elements
(26)
E
7
≡
1
2
(
e
1
⊗
e
2
-
e
2
⊗
e
1
)
,
E
8
≡
1
2
(
e
2
⊗
e
3
-
e
3
⊗
e
2
)
,
E
9
≡
1
2
(
e
3
⊗
e
1
-
e
1
⊗
e
3
)
.
In this new system of bases, the components of stress, strain, and elastic tensors, respectively, are defined as
(27)
σ
=
S
A
E
A
,
ϵ
=
E
A
E
A
,
C
=
C
A
B
e
A
⊗
e
B
,
where all the implicit sums concerning the indices {
A
,
B
,
…
.
} range from 1 to 9 and S is the compliance tensor.Thus, for Hooke’s law and for eigenstiffness-eigenstrain equations, one can have the following equivalences:(28)
σ
i
j
=
C
i
j
k
l
ϵ
k
l
⟺
S
A
=
C
A
B
E
B
,
C
i
j
k
l
ϵ
k
l
=
λ
ϵ
i
j
⟺
C
A
B
E
B
=
λ
E
A
.
Using elementary algebra, we can have the components of the stress and strain tensors in two bases:
(29)
(
S
1
S
2
S
3
S
4
S
5
S
6
S
7
S
8
S
9
)
=
(
σ
11
σ
22
σ
33
σ
(
12
)
σ
(
23
)
σ
(
31
)
σ
[
12
]
σ
[
23
]
σ
[
31
]
)
;
(
E
1
E
2
E
3
E
4
E
5
E
6
E
7
E
8
E
9
)
=
(
ϵ
11
ϵ
22
ϵ
33
ϵ
(
12
)
ϵ
(
23
)
ϵ
(
31
)
ϵ
[
12
]
ϵ
[
23
]
ϵ
[
31
]
)
,
where we have used the following notations:
(30)
θ
(
i
j
)
≡
1
2
(
θ
i
j
+
θ
j
i
)
;
θ
[
i
j
]
≡
1
2
(
θ
i
j
-
θ
j
i
)
.Now in the new basisE
A, the new components of stiffness tensor C are(31)
(
C
11
C
12
C
13
C
14
C
15
C
16
C
17
C
18
C
19
C
21
C
22
C
23
C
24
C
25
C
26
C
27
C
28
C
29
C
31
C
32
C
33
C
34
C
35
C
36
C
37
C
38
C
39
C
41
C
42
C
43
C
44
C
45
C
46
C
47
C
48
C
49
C
51
C
52
C
53
C
54
C
55
C
56
C
57
C
58
C
59
C
61
C
62
C
63
C
64
C
65
C
66
C
67
C
68
C
69
C
71
C
72
C
73
C
74
C
75
C
76
C
77
C
78
C
79
C
81
C
82
C
83
C
84
C
85
C
86
C
87
C
88
C
89
C
91
C
92
C
93
C
94
C
95
C
96
C
97
C
98
C
99
)
=
C
1111
C
1122
C
1133
C
2211
C
2222
C
3333
C
3311
C
3322
C
3333
C
11
(
12
)
C
11
(
23
)
C
11
(
31
)
C
22
(
12
)
C
22
(
23
)
C
22
(
31
)
C
33
(
12
)
C
33
(
23
)
C
33
(
31
)
C
11
[
12
]
C
11
[
23
]
C
11
[
31
]
C
22
[
12
]
C
22
[
23
]
C
22
[
31
]
C
33
[
12
]
C
33
[
23
]
C
33
[
31
]
C
(
12
)
11
C
(
12
)
22
C
(
12
)
33
C
(
23
)
11
C
(
23
)
22
C
(
23
)
33
C
(
31
)
11
C
(
31
)
22
C
(
31
)
33
C
(
1212
)
C
(
1223
)
C
(
1231
)
C
(
2312
)
C
(
2323
)
C
(
2331
)
C
(
3112
)
C
(
3123
)
C
(
3131
)
C
(
12
)
[
12
]
C
(
12
)
[
23
]
C
(
12
)
[
31
]
C
(
23
)
[
12
]
C
(
23
)
[
23
]
C
(
23
)
[
31
]
C
(
31
)
[
12
]
C
(
31
)
[
23
]
C
(
31
)
[
31
]
C
[
12
]
11
C
[
12
]
22
C
[
12
]
33
C
[
23
]
11
C
[
23
]
22
C
[
123
]
33
C
[
21
]
11
C
[
31
]
22
C
[
31
]
33
C
[
12
]
(
12
)
C
[
12
]
(
23
)
C
[
12
]
(
21
)
C
[
23
]
(
12
)
C
[
23
]
(
23
)
C
[
23
]
(
21
)
C
[
31
]
(
12
)
C
[
31
]
(
23
)
C
[
31
]
(
21
)
C
[
1212
]
C
[
1223
]
C
[
1231
]
C
[
2312
]
C
[
2323
]
C
[
2331
]
C
[
3112
]
C
[
3123
]
C
[
3131
]
,where
(32)
C
i
j
(
k
l
)
≡
1
2
(
C
i
j
k
l
+
C
i
j
l
k
)
;
C
i
j
[
k
l
]
≡
1
2
(
C
i
j
k
l
-
C
i
j
l
k
)
,
C
(
i
j
)
k
l
≡
1
2
(
C
i
j
k
l
+
C
j
i
k
l
)
;
C
[
i
j
]
k
l
≡
1
2
(
C
i
j
k
l
-
C
j
i
k
l
)
,
C
(
i
j
k
l
)
≡
1
2
(
C
i
j
k
l
+
C
i
j
l
k
+
C
j
i
k
l
+
C
j
i
l
k
)
,
C
(
i
j
)
[
k
l
]
≡
1
2
(
C
i
j
k
l
-
C
i
j
l
k
+
C
j
i
k
l
-
C
j
i
l
k
)
,
C
[
i
j
]
(
k
l
)
≡
1
2
(
C
i
j
k
l
+
C
i
j
l
k
-
C
j
i
k
l
-
C
j
i
l
k
)
,
C
[
i
j
k
l
]
≡
1
2
(
C
i
j
k
l
-
C
j
i
k
l
-
C
i
j
l
k
+
C
j
i
l
k
)
.
In case if we impose symmetries of stress and strain tensors, let us see what happens with generalized Hooke’s law (22).In generalized Hooke’s law, when the stress-strain symmetries do not affect the stiffness tensorC, the number of components of stiffness tensor in 3-dimensional space is equal to 3
4
=
81. Now, if we impose the symmetries of stress-strain tensors, that is, σ
i
j
=
σ
j
i and ϵ
k
l
=
ϵ
l
k, the C
i
j
k
l will be like C
i
j
k
l
=
C
j
i
k
l
=
C
i
j
l
k.Moreover, imposing symmetrical connection, that is,C
i
j
k
l
=
C
k
l
i
j, we would have only 21 significant components out of 81. Thus, if we flatten the stiffness tensor under the notion of Hooke’s law, we will definitely have a 6 × 6 matrix having only 21 independent elastic coefficients instead of 9 × 9 matrix having 81 elastic coefficients.Here, we depict a figure (see Figure1) which delineates the component reduction process for the stiffness tensor.Figure 1
Hooke’s law (reduction process of elastic coefficients). Here, still there are 36 components seen in the stiffness table, but due to the components, for example,C
2313
=
C
1323, and so forth, the counting of 36 components will be reduced up to 21. Also, in the above figure, the components having gray background expose symmetry.Now, with 21 significant independent components, the stiffness tensorC
i
j
k
l can be mapped on a symmetric 6 × 6 matrix.As the elasticity of a material is described by a fourth-order tensor with 21 independent components as shown in (Figure1) and the mathematical description of elasticity tensor appears in Hooke’s law, now how to map this 3-dimensional fourth order tensor on a 6 × 6 matrix?We are fortunate to have a long series of research papers concerning this issue. For instance, [2, 3, 32–39] and many more.The very first approach to map 21 significant components of the elasticity tensor on a symmetric 6 × 6 matrix was introduced by Voigt [32] and after that various successors of Voigt have used his fabulous notions for flattening of fourth-order tensors. But Lord Kelvin [40] found the Voigt notations are inadequate from the perspective of tensorial nature and then introduced his own advanced notations, now known as Kelvin’s mapping. More recently, [39] have introduced some sophisticated methodology to unfold a fourth-rank tensor called “MCN” (Mehrabadi and Cowin’s notations).Let us briefly go through these three notions of tensor flattening one by one.


### 3.1. The Voigt Six-Dimensional Notations for Unfolding an Elasticity Tensor

It is well known that the Voigt mapping preserves the elastic energy density of the material and elastic stiffness and is given by(33)
2
·
E
energy
=
σ
i
j
ϵ
i
j
=
σ
p
ϵ
p
.
The Voigt mapping receives this relation by the mapping rules
(34)
p
=
i
δ
i
j
+
(
1
-
δ
i
j
)
(
9
-
i
-
j
)
,
q
=
k
δ
k
l
+
(
1
-
δ
k
l
)
(
9
-
k
-
l
)
.
Using this rule, we have σ
i
j
=
σ
p, C
i
j
k
l
=
C
p
q, and ϵ
q
=
(
2
-
δ
k
l
)
ϵ
k
l.This Voigt mapping can be visualized as shown in Table1.Table 1
Voigt’s mappings for stress (left) and strain tensors (right), respectively.Thus, in accordance with Voigt’s mappings, Hooke’s law (22) can be represented in matrix form as follows:
(35)
(
σ
1
σ
2
σ
3
σ
4
σ
5
σ
6
)
=
(
C
11
C
12
C
13
C
14
C
15
C
16
C
21
C
22
C
23
C
24
C
25
C
26
C
31
C
32
C
33
C
34
C
35
C
36
C
41
C
42
C
43
C
44
C
45
C
46
C
51
C
52
C
53
C
54
C
55
C
56
C
61
C
62
C
63
C
64
C
65
C
66
)
(
ϵ
1
ϵ
2
ϵ
3
ϵ
4
ϵ
5
ϵ
6
)
,
where the simple index conversion rule of Voigt (see Table 2) is applied.Table 2
Index conversion rule of Voigt.
i
,
j
11
22
33
23 or 32
13 or 31
12 or 21
p
1
2
3
4
5
6But in the Voigt notations, many disadvantages were noticed. For instance,(1)
theσ
i
j and ϵ
k
l are treated differently;
(2)
the norms ofσ
i
j, ϵ
k
l, and C
i
j
k
l are not preserved;
(3)
the entries in all the three Voigt arrays (see (35)) are not the tensor or the vector components and thus one may not be able to have advantages of tensor algebra like tensor law of transformation, rotation of coordinate system, and so forth.


### 3.2. Kelvin’s Mapping Rules

Likewise Voigt’s mapping rules, Kelvin’s mapping rules also preserve the elastic energy density of material under the following methodology:(36)
p
=
i
δ
i
j
+
(
1
-
δ
i
j
)
(
9
-
i
-
j
)
,
q
=
k
δ
k
l
+
(
1
-
δ
k
l
)
(
9
-
k
-
l
)
.
In this way,
(37)
σ
p
=
(
δ
i
j
+
2
(
1
-
δ
i
j
)
)
ϵ
i
j
,
ϵ
q
=
(
δ
k
l
+
2
(
1
-
δ
k
l
)
)
ϵ
k
l
,
C
p
q
=
(
δ
i
j
+
2
(
1
-
δ
i
j
)
)
(
δ
k
l
+
2
(
1
-
δ
k
l
)
)
C
i
j
k
l
.
Naturally, Kelvin’s mapping preserves the norms of the three tensors and the stress and strain are treated identically. Also, the mappings of stress, strain, and elasticity tensors have all the properties of tensor of first- and second-rank tensors, respectively, in 6-dimensional space.But the only disadvantage of this mapping is that the values of stiffness components are changed. However, there is a simple tool for conversion between Voigt’s and Kelvin’s notations. For this purpose, one just needed a single array (see Table3). Thus,
(38)
σ
kelvin
=
σ
viogt
ξ
,
σ
voigt
=
σ
kelvin
ξ
,
ϵ
kelvin
=
ϵ
voigt
ξ
,
ϵ
voigt
=
ϵ
kelvin
ξ
.Table 3
Voigt to Kelvin and vice versa.
ξ
=
1
1
1
2
2
2


### 3.3. Mehrabadi-Cowin Notations (MCN)

To preserve the tensor properties of Hooke’s law during the unfolding of fourth-order elasticity tensor, [39] have introduced a new notion which is nothing but a conversion of Voigt’s notations into Kelvin ones. This conversion is done using the following relation:
(39)
C
kelvin
=
(
1
1
1
2
2
2
1
1
1
2
2
2
1
1
1
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
)
C
voigt
.
Exploiting this new notion, (35) becomes
(40)
(
σ
11
σ
22
σ
33
2
σ
23
2
σ
13
2
σ
12
)
=
(
C
1111
C
1122
C
1133
2
C
1123
2
C
1113
2
C
1112
C
2211
C
2222
C
2333
2
C
2223
2
C
2213
2
C
2212
C
3311
C
3322
C
3333
2
C
3323
2
C
3313
2
C
3312
2
C
2311
2
C
2322
2
C
2333
2
C
2323
2
C
2313
2
C
2312
2
C
1311
2
C
1322
2
C
1333
2
C
1323
2
C
1313
2
C
1312
2
C
1211
2
C
1222
2
C
1233
2
C
1223
2
C
1213
2
C
1212
)
*
(
ϵ
11
ϵ
22
ϵ
33
2
ϵ
23
2
ϵ
13
2
ϵ
12
)
.
Further, to involve the features of tensor in the matrix representation (40), [39] have evoked that (40) can be rewritten as
(41)
σ
^
=
C
^
ϵ
^
,
where the new six-dimensional stress and strain vectors denoted by σ
^ and ϵ
^, respectively, are multiplied by the factors 2. Also, C
^ is a new six-by-six matrix [39]. The matrix form of (40) is given as
(42)
(
σ
11
σ
22
σ
33
2
σ
23
2
σ
13
2
σ
12
)
=
(
C
11
C
12
C
13
2
C
14
2
C
15
2
C
16
C
12
C
22
C
23
2
C
24
2
C
25
2
C
26
C
13
C
23
C
33
2
C
34
2
C
35
2
C
36
2
C
14
2
C
24
2
C
34
2
C
44
2
C
45
2
C
46
2
C
15
2
C
25
2
C
35
2
C
45
2
C
55
2
C
56
2
C
16
2
C
26
2
C
36
2
C
46
2
C
56
2
C
66
)
*
(
ϵ
11
ϵ
22
ϵ
33
2
ϵ
23
2
ϵ
13
2
ϵ
12
)
.
The representation (42) is further produced in MCN pattern like below:
(43)
(
σ
^
1
σ
^
2
σ
^
3
σ
^
4
σ
^
5
σ
^
6
)
=
(
C
^
11
C
^
12
C
^
13
C
^
14
C
^
15
C
^
16
C
^
12
C
^
22
C
^
23
C
^
24
C
^
25
C
^
26
C
^
13
C
^
23
C
^
33
C
^
34
C
^
35
C
^
36
C
^
14
C
^
24
C
^
34
C
^
44
C
^
45
C
^
46
C
^
15
C
^
25
C
^
35
C
^
45
C
^
55
C
^
56
C
^
16
C
^
26
C
^
36
C
^
46
C
^
56
C
^
66
)
*
(
ϵ
^
1
ϵ
^
2
ϵ
^
3
ϵ
^
4
ϵ
^
5
ϵ
^
6
)
.
To deal with all the above representations of a fourth-order elasticity tensor as a six-by-six matrix, a simple conversion table has been proposed by [39], which is depicted as in Table 4.Table 4
The stress, strain, elasticity, and compliance tensors components in Voigt’s, Kelvin’s, and MCN patterns.
Voigt's notations
Kelvin's notations
MCN
Stress
σ
11
σ
1
σ
^
1
σ
22
σ
2
σ
^
2
σ
33
σ
3
σ
^
3
σ
23
σ
4
σ
^
4
σ
13
σ
5
σ
^
5
σ
12
σ
5
σ
^
6
Strain
ϵ
11
ϵ
1
ϵ
^
1
ϵ
22
ϵ
2
ϵ
^
2
ϵ
33
ϵ
3
ϵ
^
3
ϵ
23
ϵ
4
ϵ
^
4
ϵ
13
ϵ
15
ϵ
^
5
ϵ
12
ϵ
6
ϵ
^
6
Elasticity
C
1111
C
11
C
^
11
C
2222
C
22
C
^
22
C
3333
C
33
C
^
33
C
1122
C
12
C
^
12
C
1133
C
13
C
^
13
C
2233
C
23
C
^
23
C
2323
C
44
1
2
C
^
44
C
1313
C
55
1
2
C
^
55
C
1212
C
66
1
2
C
^
66
C
1323
C
54
1
2
C
^
54
C
1312
C
56
1
2
C
^
56
C
1223
C
64
1
4
C
^
64
C
2311
C
41
1
2
C
^
41
C
1311
C
51
1
2
C
^
51
C
1211
C
61
1
2
C
^
61
C
2322
C
42
1
2
C
^
42
C
1322
C
52
1
2
C
^
52
C
1222
C
62
1
2
C
^
62
C
2333
C
43
1
2
C
^
43
C
1333
C
53
1
2
C
^
53
C
1233
C
63
1
2
C
^
63
Compliance
K
1111
S
11
S
^
11
K
2222
S
22
S
^
22
K
3333
S
33
S
^
33
K
1122
S
12
S
^
12
K
1133
S
13
S
^
13
K
2233
S
23
S
^
23
K
2323
1
4
S
44
1
2
S
^
44
K
1313
1
4
S
55
1
2
S
^
55
K
1212
1
4
S
66
1
2
S
^
66
K
1323
1
4
S
54
1
2
S
^
54
K
1312
1
4
S
56
1
2
S
^
56
K
1223
1
4
S
64
1
2
S
^
64
K
2311
1
2
S
41
1
2
S
^
41
K
1311
1
2
S
51
1
2
S
^
51
K
1211
1
2
S
61
1
2
S
^
61
K
2322
1
2
S
42
1
2
S
^
42
K
1322
1
2
S
52
1
2
S
^
52
K
1222
1
2
S
62
1
2
S
^
62
K
2333
1
2
S
43
1
2
S
^
43
K
1333
1
2
S
53
1
2
S
^
53
K
1233
1
2
S
63
1
2
S
^
63Even though the foregoing detailed description is interesting and important from the viewpoint of those who are beginners in the field of mathematical elasticity, the modern scenario of the literature of mathematical elasticity now involves the assistance of various computer algebraic systems (CAS), like MATLAB [30], Maple [41], Mathematica [42], wxMaxima [43], and some peculiar packages like Tensor Toolbox [26–29], GRTensor [44], Ricci [45], and so forth, to handle particular problems of the scientific literature.However, in the following section, we shall delineate a CAS approach for Section3.


## 3.1. The Voigt Six-Dimensional Notations for Unfolding an Elasticity Tensor

It is well known that the Voigt mapping preserves the elastic energy density of the material and elastic stiffness and is given by(33)
2
·
E
energy
=
σ
i
j
ϵ
i
j
=
σ
p
ϵ
p
.
The Voigt mapping receives this relation by the mapping rules
(34)
p
=
i
δ
i
j
+
(
1
-
δ
i
j
)
(
9
-
i
-
j
)
,
q
=
k
δ
k
l
+
(
1
-
δ
k
l
)
(
9
-
k
-
l
)
.
Using this rule, we have σ
i
j
=
σ
p, C
i
j
k
l
=
C
p
q, and ϵ
q
=
(
2
-
δ
k
l
)
ϵ
k
l.This Voigt mapping can be visualized as shown in Table1.Table 1
Voigt’s mappings for stress (left) and strain tensors (right), respectively.Thus, in accordance with Voigt’s mappings, Hooke’s law (22) can be represented in matrix form as follows:
(35)
(
σ
1
σ
2
σ
3
σ
4
σ
5
σ
6
)
=
(
C
11
C
12
C
13
C
14
C
15
C
16
C
21
C
22
C
23
C
24
C
25
C
26
C
31
C
32
C
33
C
34
C
35
C
36
C
41
C
42
C
43
C
44
C
45
C
46
C
51
C
52
C
53
C
54
C
55
C
56
C
61
C
62
C
63
C
64
C
65
C
66
)
(
ϵ
1
ϵ
2
ϵ
3
ϵ
4
ϵ
5
ϵ
6
)
,
where the simple index conversion rule of Voigt (see Table 2) is applied.Table 2
Index conversion rule of Voigt.
i
,
j
11
22
33
23 or 32
13 or 31
12 or 21
p
1
2
3
4
5
6But in the Voigt notations, many disadvantages were noticed. For instance,(1)
theσ
i
j and ϵ
k
l are treated differently;
(2)
the norms ofσ
i
j, ϵ
k
l, and C
i
j
k
l are not preserved;
(3)
the entries in all the three Voigt arrays (see (35)) are not the tensor or the vector components and thus one may not be able to have advantages of tensor algebra like tensor law of transformation, rotation of coordinate system, and so forth.


## 3.2. Kelvin’s Mapping Rules

Likewise Voigt’s mapping rules, Kelvin’s mapping rules also preserve the elastic energy density of material under the following methodology:(36)
p
=
i
δ
i
j
+
(
1
-
δ
i
j
)
(
9
-
i
-
j
)
,
q
=
k
δ
k
l
+
(
1
-
δ
k
l
)
(
9
-
k
-
l
)
.
In this way,
(37)
σ
p
=
(
δ
i
j
+
2
(
1
-
δ
i
j
)
)
ϵ
i
j
,
ϵ
q
=
(
δ
k
l
+
2
(
1
-
δ
k
l
)
)
ϵ
k
l
,
C
p
q
=
(
δ
i
j
+
2
(
1
-
δ
i
j
)
)
(
δ
k
l
+
2
(
1
-
δ
k
l
)
)
C
i
j
k
l
.
Naturally, Kelvin’s mapping preserves the norms of the three tensors and the stress and strain are treated identically. Also, the mappings of stress, strain, and elasticity tensors have all the properties of tensor of first- and second-rank tensors, respectively, in 6-dimensional space.But the only disadvantage of this mapping is that the values of stiffness components are changed. However, there is a simple tool for conversion between Voigt’s and Kelvin’s notations. For this purpose, one just needed a single array (see Table3). Thus,
(38)
σ
kelvin
=
σ
viogt
ξ
,
σ
voigt
=
σ
kelvin
ξ
,
ϵ
kelvin
=
ϵ
voigt
ξ
,
ϵ
voigt
=
ϵ
kelvin
ξ
.Table 3
Voigt to Kelvin and vice versa.
ξ
=
1
1
1
2
2
2


## 3.3. Mehrabadi-Cowin Notations (MCN)

To preserve the tensor properties of Hooke’s law during the unfolding of fourth-order elasticity tensor, [39] have introduced a new notion which is nothing but a conversion of Voigt’s notations into Kelvin ones. This conversion is done using the following relation:
(39)
C
kelvin
=
(
1
1
1
2
2
2
1
1
1
2
2
2
1
1
1
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
)
C
voigt
.
Exploiting this new notion, (35) becomes
(40)
(
σ
11
σ
22
σ
33
2
σ
23
2
σ
13
2
σ
12
)
=
(
C
1111
C
1122
C
1133
2
C
1123
2
C
1113
2
C
1112
C
2211
C
2222
C
2333
2
C
2223
2
C
2213
2
C
2212
C
3311
C
3322
C
3333
2
C
3323
2
C
3313
2
C
3312
2
C
2311
2
C
2322
2
C
2333
2
C
2323
2
C
2313
2
C
2312
2
C
1311
2
C
1322
2
C
1333
2
C
1323
2
C
1313
2
C
1312
2
C
1211
2
C
1222
2
C
1233
2
C
1223
2
C
1213
2
C
1212
)
*
(
ϵ
11
ϵ
22
ϵ
33
2
ϵ
23
2
ϵ
13
2
ϵ
12
)
.
Further, to involve the features of tensor in the matrix representation (40), [39] have evoked that (40) can be rewritten as
(41)
σ
^
=
C
^
ϵ
^
,
where the new six-dimensional stress and strain vectors denoted by σ
^ and ϵ
^, respectively, are multiplied by the factors 2. Also, C
^ is a new six-by-six matrix [39]. The matrix form of (40) is given as
(42)
(
σ
11
σ
22
σ
33
2
σ
23
2
σ
13
2
σ
12
)
=
(
C
11
C
12
C
13
2
C
14
2
C
15
2
C
16
C
12
C
22
C
23
2
C
24
2
C
25
2
C
26
C
13
C
23
C
33
2
C
34
2
C
35
2
C
36
2
C
14
2
C
24
2
C
34
2
C
44
2
C
45
2
C
46
2
C
15
2
C
25
2
C
35
2
C
45
2
C
55
2
C
56
2
C
16
2
C
26
2
C
36
2
C
46
2
C
56
2
C
66
)
*
(
ϵ
11
ϵ
22
ϵ
33
2
ϵ
23
2
ϵ
13
2
ϵ
12
)
.
The representation (42) is further produced in MCN pattern like below:
(43)
(
σ
^
1
σ
^
2
σ
^
3
σ
^
4
σ
^
5
σ
^
6
)
=
(
C
^
11
C
^
12
C
^
13
C
^
14
C
^
15
C
^
16
C
^
12
C
^
22
C
^
23
C
^
24
C
^
25
C
^
26
C
^
13
C
^
23
C
^
33
C
^
34
C
^
35
C
^
36
C
^
14
C
^
24
C
^
34
C
^
44
C
^
45
C
^
46
C
^
15
C
^
25
C
^
35
C
^
45
C
^
55
C
^
56
C
^
16
C
^
26
C
^
36
C
^
46
C
^
56
C
^
66
)
*
(
ϵ
^
1
ϵ
^
2
ϵ
^
3
ϵ
^
4
ϵ
^
5
ϵ
^
6
)
.
To deal with all the above representations of a fourth-order elasticity tensor as a six-by-six matrix, a simple conversion table has been proposed by [39], which is depicted as in Table 4.Table 4
The stress, strain, elasticity, and compliance tensors components in Voigt’s, Kelvin’s, and MCN patterns.
Voigt's notations
Kelvin's notations
MCN
Stress
σ
11
σ
1
σ
^
1
σ
22
σ
2
σ
^
2
σ
33
σ
3
σ
^
3
σ
23
σ
4
σ
^
4
σ
13
σ
5
σ
^
5
σ
12
σ
5
σ
^
6
Strain
ϵ
11
ϵ
1
ϵ
^
1
ϵ
22
ϵ
2
ϵ
^
2
ϵ
33
ϵ
3
ϵ
^
3
ϵ
23
ϵ
4
ϵ
^
4
ϵ
13
ϵ
15
ϵ
^
5
ϵ
12
ϵ
6
ϵ
^
6
Elasticity
C
1111
C
11
C
^
11
C
2222
C
22
C
^
22
C
3333
C
33
C
^
33
C
1122
C
12
C
^
12
C
1133
C
13
C
^
13
C
2233
C
23
C
^
23
C
2323
C
44
1
2
C
^
44
C
1313
C
55
1
2
C
^
55
C
1212
C
66
1
2
C
^
66
C
1323
C
54
1
2
C
^
54
C
1312
C
56
1
2
C
^
56
C
1223
C
64
1
4
C
^
64
C
2311
C
41
1
2
C
^
41
C
1311
C
51
1
2
C
^
51
C
1211
C
61
1
2
C
^
61
C
2322
C
42
1
2
C
^
42
C
1322
C
52
1
2
C
^
52
C
1222
C
62
1
2
C
^
62
C
2333
C
43
1
2
C
^
43
C
1333
C
53
1
2
C
^
53
C
1233
C
63
1
2
C
^
63
Compliance
K
1111
S
11
S
^
11
K
2222
S
22
S
^
22
K
3333
S
33
S
^
33
K
1122
S
12
S
^
12
K
1133
S
13
S
^
13
K
2233
S
23
S
^
23
K
2323
1
4
S
44
1
2
S
^
44
K
1313
1
4
S
55
1
2
S
^
55
K
1212
1
4
S
66
1
2
S
^
66
K
1323
1
4
S
54
1
2
S
^
54
K
1312
1
4
S
56
1
2
S
^
56
K
1223
1
4
S
64
1
2
S
^
64
K
2311
1
2
S
41
1
2
S
^
41
K
1311
1
2
S
51
1
2
S
^
51
K
1211
1
2
S
61
1
2
S
^
61
K
2322
1
2
S
42
1
2
S
^
42
K
1322
1
2
S
52
1
2
S
^
52
K
1222
1
2
S
62
1
2
S
^
62
K
2333
1
2
S
43
1
2
S
^
43
K
1333
1
2
S
53
1
2
S
^
53
K
1233
1
2
S
63
1
2
S
^
63Even though the foregoing detailed description is interesting and important from the viewpoint of those who are beginners in the field of mathematical elasticity, the modern scenario of the literature of mathematical elasticity now involves the assistance of various computer algebraic systems (CAS), like MATLAB [30], Maple [41], Mathematica [42], wxMaxima [43], and some peculiar packages like Tensor Toolbox [26–29], GRTensor [44], Ricci [45], and so forth, to handle particular problems of the scientific literature.However, in the following section, we shall delineate a CAS approach for Section3.


## 4. Unfolding the Anisotropic Hooke’s Law with the Aid of “Mathematica”

“Mathematica” is a general computing environment intimated with organizing algorithmic, visualization, and user-friendly interface capabilities. Moreover, many mathematical algorithms encapsulated by “Mathematica” make the computation easy and fast [42, 46].A fabulous book entitled “Elasticity with Mathematica” has been produced by [31]. This book is equipped with some great Mathematica packages, namely, VectorAnalysis.m, Displacement.m, IntegrateStrain.m, ParametricMesh.m, and Tensor2Analysis.m.Overall the effort of [31] is greatly appreciated for producing a masterpiece with Mathematica packages, notebooks, and worked examples. Moreover, all the aforementioned packages, notebooks, and worked examples are freely available to readers and can be downloaded from the publisher’s website http://www.cambridge.org/9780521842013.We consider the following code that easily meets the purpose discussed in Section3.Here, kindly note that only the input Mathematica codes are being exposed. For considering the entire processing, an Electronic (see Appendix  A in Supplementary Material available online athttp://dx.doi.org/10.1155/2014/487314) is being proposed to readers.First of all install the “Tensor2Analysis” package in Mathematica using the following commandIn[1]:= << Tensor2Analysis`We have loaded the above package like Loading the package “Tensor2Analysis.m” by setting the following path:SetDirectory[C:/Program Files/Wolfram/Research/Mathematica/8.0/AddOns/Packages].Next is concerning the creation of tensor. The code in Algorithm1 generates the stress and strain tensors using the index rules specified by [31].Algorithm 1
I
n[
2
]
≔  SymIndex[2,  dim_]  ≔
Apply[  (Sort[List[##]]  &),  Array[  List,  Array[dim  &,  2]],  {2}  ]
SymIndex[
2
,  3
]
SymIndex[
2
,  6
]
In[
3
]
≔
SymIndex[
4
]  =
Apply[  Flatten[  Sort[  Map[  Sort,  Partition[List[##],  2
]]]]  &,
Array[  List,  Array[
3  &,  4
]],  {4}]Use the MakeTensor command to create tensorial expressions having components with respect to symmetric conditions. This command combines all the previous commands to do so (see Algorithm2).Algorithm 2
I
n[
4
]
≔
MakeName[myexp__]  ≔
ToExpression  //@  (  StringJoin  @@  Map[  ToString[  #]  &,  List[  myexp]])
MakeTensor[  mystring_,  2,  dim_  ]  ≔
Apply[(MakeName[mystring,  ##]  &),  SymIndex[2,  dim],  {2}]
MakeTensor[  mystring_,  4  ]  ≔
Apply[(MakeName[mystring,  ##]  &),  SymIndex[
4
],  {4}]
In[
5
]
≔
(  stress  =  MakeTensor["sigma",  2,  3])  //  MatrixForm
In[
6
]
≔
(  strain  =  MakeTensor["epsilon",  2,  3])  //  MatrixForm
In[
7
]
≔
(  C4  =  MakeTensor["C",  4]  )  //  MatrixFormNow the next code deals with the index rules that is able to handle Hooke’s tensor conversion from the fourth-order to the second-order tensor or second-order Voigt form using indexrule (see Algorithm3).Algorithm 3
I
n[
8
]  ≔
indexrule2to1  =  {
{
1,  1}  ->  1,  {2,  2}  ->  2,  {3,  3}  ->  3,  {2,  3}  ->
4  ,  {3,  1}  ->  5,  {1,  2}  ->  6,  {3,  2}  ->  4  ,  {1,  3}  ->
5,  {2,  1}  ->  6}
indexrule1to2  =  Map[Rule[#[[
2
]],  #[[
1
]]  ]  &  ,  indexrule2to1]
In[
9
]  ≔
Index6[
1
]  =  Range[  6
]  /.  indexrule1to2Afterwards, we have a crucial stage that concerns with the transformation of Voigt’s notations into a fourth-rank tensor and vice versa. This stage includes the codes like HookeVto4, Hooke4toV. Also the codes HookeVto2 and Hook2toV transform the Voigt notations into the second-order tensor notations and vice versa. Finally the code in Algorithm4 provides us a splendid form of fourth-order elasticity tensor as a six-by-six matrix in MCN notations.Algorithm 4
I
n[
10
]
≔
HookeVto4[  myC_]  ≔
Array[  myC[[  {#1,  #2}  /.  indexrule2to1,  {#3,  #4}  /.
indexrule2to1]]  &,  Array[
3  &,  4
]]
In[
11
]
≔
(C2  =  MakeTensor[
"
C",  2,  6])  //  MatrixForm
In[
12
]
≔
C2
In[
13
]
≔
(C2to4  =  HookeVto4[C2]  )  //  MatrixForm
In[
14
]
≔
C2to4
In[
15
]
≔
Hooke4toV[  myC_]  ≔
Apply[  Part[  myC,  #
#]  &,
Array[  Join[  #1  /.  indexrule1to2  ,  #2  /.  indexrule1to2]  &,
Array[6  &,  2]]  ,  {2}]
In[
16
]
≔
C2back  =  Hooke4toV[C2to4]  )  //  MatrixForm
In[
17
]
≔
(C4to2  =  Hooke4toV[C4]  )  //  MatrixForm
In[
18
]
≔
(C4back  =  HookeVto4[C4to2])  //  MatrixForm
In[
19
]
≔
Hooke2toV[myc2_]  ≔  Table
[
Which
[
i  <=  3  &&  j  <=  3  ,  myc2[[i,  j]],
4  <=  i  &&  j  <=  3  ,  myc2[[i,  j]]  /  2
∧
(1/2),
i  <=  3  &&  4  <=  j  ,  myc2[[i,  j]]/  2
∧
(1/2),
4  <=  i  &&  4  <=  j  ,  myc2[[i,  j]]/  2
]
,  {i,  6},  {j,  6}]
HookeVto2[mycV_]  ≔  Table[
Which[
i  <=  3  &&  j  <=  3  ,  mycV[[i,  j]],
4  <=  i  &&  j  <=  3  ,  mycV[[i,  j]]  *  2
∧
(1/2),
i  <=  3  &&  4  <=  j  ,  mycV[[i,  j]]  *  2
∧
(1/2),
4  <=  i  &&  4  <=  j  ,  mycV[[i,  j]]  *  2
],  {i,  6},  {j,  6}]
Hooke4to2[myC_]  ≔  HookeVto2[Hooke4toV[  myC  ]]
Hooke2to4[myC_]  ≔  HookeVto4[Hooke2toV[  myC  ]]
In[
20
]
≔
(CV  =  Hooke2toV[  C2  ]  )  //  MatrixForm
In[
21
]
≔
(C2back  =  HookeVto2[  CV  ]  )  //  MatrixForm
In[
22
]
≔
(CV  =  Hooke4to2[C4])  //  MatrixForm
In[
23
]
≔
C4back  =  Hooke2to4[CV])  //  MatrixFormWe have presented here a minimal code of mathematica developed by [31]. However, a numerical demonstration of this code has also been presented by [31] and the reader(s), interested in understanding this code in full, are requested to refer to the website already mentioned above.Let us now proceed to Section5 which in our view is the soul of the present work. In this section, we shall deal with the eigenvalues, eigenvectors, the nominal average of eigenvectors and the average eigenvectors, and so forth, for different species of softwoods, hardwoods and some specimen of cancellous bone.


## 5. Analysis of Known Values of Elastic Constants of Woods and Cancellous Bone Using MAPLE

Of course, “MAPLE” is a sophisticated CAS and produced under the results of over 30 years of cutting-edge research and development, which assists us in analyzing, exploring, visualizing, and manipulating almost every mathematical problem. Having more than 500 functions, this CAS offers broadness, depth, and performance to handle every phenomenon of mathematics. In a nutshell, it is a CAS that offers high performance mathematics capabilities with integrated numeric and symbolic computations.However, in the present study, we attempt to explore how this CAS handles complicated analysis regarding elastic constant data.The elastic constant data of our interest, for hardwoods and softwoods, as calculated by Hearmon [47, 48] are tabulated in Tables 5 and 6.Table 5
The elastic constants data for hardwoods.
S. no.
Species
ρ
C
^
11
C
^
22
C
^
33
C
^
12
C
^
13
C
^
23
C
^
44
C
^
55
C
^
66
1
Quipo
0.1
0.045
0.251
1.075
0.027
0.033
0.025
0.226
0.118
0.078
2
Quipo
0.2
0.159
0.427
3.446
0.069
0.131
0.178
0.430
0.280
0.144
3
White
0.38
0.547
1.192
10.041
0.399
0.360
0.555
1.442
1.344
0.022
4
Khaya
0.44
0.631
1.381
10.725
0.389
0.520
0.662
1.800
1.196
0.420
5
Mahogany
0.50
0.952
1.575
11.996
0.571
0.682
0.790
1.960
1.498
0.638
6
Mahogany
0.53
0.765
1.538
13.010
0.655
0.631
0.841
1.218
0.938
0.300
7
S. Germ
0.54
0.772
1.772
12.240
0.558
0.530
0.871
2.318
1.582
0.540
8
Maple
0.58
1.451
2.565
11.492
1.197
1.267
1.818
2.460
2.194
0.584
9
Walnut
0.59
0.927
1.760
12.432
0.707
0.936
1.312
1.922
1.400
0.460
10
Birch
0.62
0.898
1.623
1.7173
0.671
0.714
1.075
2.346
1.816
0.372
11
Y. Birch
0.64
1.084
1.697
15.288
0.777
0.883
1.191
2.120
1.942
0.480
12
Oak
0.67
1.350
2.983
16.958
1.007
1.005
1.463
2.380
1.532
0.784
13
Ash
0.68
1.135
2.142
16.958
0.827
0.917
1.427
2.684
1.784
0.540
14
Ash
0.80
1.439
2.439
17.000
1.037
1.485
1.968
1.720
1.218
0.500
15
Beech
0.74
1.659
3.301
15.437
1.279
1.433
2.142
3.216
2.112
0.912
Source: this data is consulted from [47, 48]. The number of repetitions of the particular species in this table shows the number of measurements done for that particular wood species. The units of ρ (bulk or apparent density) are g/cm3 and the units of elastic constants are GPa = 1
0
10 dynes/cm2.Table 6
The elastic constants data for softwoods.
S. no.
Species
ρ
C
^
11
C
^
22
C
^
33
C
^
12
C
^
13
C
^
23
C
^
44
C
^
55
C
^
66
1
Blasa
0.2
0.127
0.360
6.380
0.086
0.091
0.154
0.624
0.406
0.066
2
Spruce
0.39
0.572
1.030
11.950
0.262
0.365
0.506
1.498
1.442
0.078
3
Spruce
0.43
0.594
1.106
14.055
0.346
0.476
0.686
1.442
1.0
0.064
4
Spruce
0.44
0.443
0.775
16.286
0.192
0.321
0.442
1.234
1.52
0.072
5
Spruce
0.50
0.755
0.963
17.221
0.333
0.549
0.548
1.25
1.706
0.07
6
Douglas Fir
0.45
0.929
1.173
16.095
0.409
0.539
0.539
1.767
1.766
0.176
7
Douglas Fir
0.59
1.226
1.775
17.004
0.753
0.747
0.941
2.348
1.816
0.160
8
Pine
0.54
0.721
1.405
16.929
0.454
0.535
0.857
3.484
1.344
0.132
Source: this data is consulted from [47, 48]. The number of repetitions of the particular species in this table shows the number of measurements done for that particular wood species. The units of ρ (bulk or apparent density) are g/cm3 and the units of elastic constants are GPa = 1
0
10 dynes/cm2.For cancellous bone, we have the elastic constants data available for three specimens [2, 11]. These data are tabulated in Table 7.Table 7
The elastic constants data for the three specimens of cancellous bone.
Elastic constants
Specimen 1
Specimen 2
Specimen 3
C
^
11
9.86
791.2
698.3
C
^
12
2.147
360.9
281.6
C
^
13
2.773
318.6
284.5
C
^
14
−0.162
1.1
8.0
C
^
15
−2.93
5.4
−0.7
C
^
16
−0.821
−0.6
11.4
C
^
22
9.35
852.9
896.7
C
^
23
2.346
328.1
322.7
C
^
24
−1.204
0.9
−162.3
C
^
25
−0.936
−2.5
−3.2
C
^
26
0.445
−4.2
16.0
C
^
33
5.952
1473.0
916.0
C
^
34
−1.309
−2.4
−162.2
C
^
35
−1.791
16.9
−7.2
C
^
36
0.54
0.6
2.4
C
^
44
5.61
358.7
348.9
C
^
45
0.798
0.5
10.5
C
^
46
−2.159
5.1
−9.8
C
^
55
6.032
344.8
242.6
C
^
56
−0.335
−0.7
−64.4
C
^
66
7.425
230.4
237.8
Source: the data for cancellous bone has been taken from [2, 11].Precisely speaking, to meet the requirements of proposed research, a MAPLE code consisting of almost 81 steps has been developed. This MAPLE code (in its minimal form) is appended to Appendix A, while its entire working mechanism is included in an electronic appendix (Appendix B) and can be accessed fromhttp://dx.doi.org/10.1155/2014/487314. Particularly, as far as our intention concerns, analysis of elastic constants data regarding hardwoods, softwoods, and cancellous bone using MAPLE is consisting of the following steps.(i)
Computating the eigenvalues and eigenvectors for all the 15 hardwood species, 8 softwood species, and 3 specimens of cancellous bone using MAPLE.
(ii)
Computing the nominal averages of eigenvectors and the average eigenvectors for all the 15 hardwood species.
(iii)
Computing the average eigenvalues for hardwood species.
(iv)
Computing the average elasticity matrices for all the 15 hardwood species.
(v)
Plotting the histograms for elasticity matrices of all the 15 hardwood species.
(vi)
Plotting the graphs for I, II, III, IV, V, and VI eigenvalues of 15 hardwood species against their apparent densities (see Figures14–16; also see Figure 17 for a combined view).


### 5.1. Computing the Eigenvalues and Eigenvectors for Hardwoods Species, Softwood Species, and Cancellous Bone Using MAPLE

Here we present the outcomes of MAPLE code (which is mentioned in Appendix A) regarding the computation of eigenvalues and eigenvectors of the stiffness matricesC
^ concerning 15 hardwood species (see Table 5).It is well known that the eigenvalues and eigenvectors of a stiffness matrixC
^ or its compliance matrix S
^ are determined from the equations [3]
(44)
(
C
^
-
Λ
I
)
N
^
=
0
or
(
S
^
-
1
Λ
I
)
N
^
=
0
,
where I is the 6 × 6 identity matrix and N
^ represents the normalized eigenvectors of C
^ or S
^. Naturally, C
^ or S
^ being positive definite, therefore there would be six positive eigenvalues and these values are known as Kelvin’s moduli. These Kelvin’s moduli are denoted by Λ
i
,
i
=
1,2
,
…
,
6 and if possible are ordered in such a way that Λ
1
≥
Λ
2
≥
⋯
≥
Λ
6
>
0. Also, the eigenvalues of S
^ can simply de computed by taking the inverse of the eigenvalues of C
^.A MAPLE code mentioned in Step 16 of Appendix A enables us to simultaneously compute the eigenvalues and eigenvectors for all the 15 hardwood species. Also, the results of this MAPLE code are summarized in Tables8 and 9. Further, by simply substituting the elasticity constants from Table 6 in Step 2 to Step 11 of Appendix A, and performing Step 16 again, one can have the eigenvalues and eigenvectors for the 8 softwood species. The computed results for softwood species are summarized in Table 10. Similarly, substitution of elastic constants data for cancellous bone from Table 7 into Step 2 to Step 11 and repeating Step 16 of the MAPLE code yields the desired results tabulated in Table 11.Table 8
The eigenvalues and eigenvectors for 15 hardwoods species calculated using MAPLE.Λ
[
i
] stand for the eigenvalues of ith softwood species, and V
[
i
] stand for the corresponding eigenvectors.
S. no.
Hardwood species
EigenvaluesΛ
[
i
]
EigenvectorsV
[
i
]
1
Quipo
[
1.076866026
0.04067345638
0.2534605191
0.2260000000
0.1180000000
0.07800000000
]
[
0.03276733390
0.9918780499
-
0.1228992897
0.0
0.0
0.0
0.03131140851
-
0.1239237163
-
0.9917976143
0.0
0.0
0.0
0.9989724208
-
0.02865041271
0.03511774899
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
]
2
Quipo
[
3.461963876
0.1399776173
0.4300584948
0.4300000000
0.2800000000
0.1440000000
]
[
0.04079957804
0.9756137601
-
0.2156691559
0.0
0.0
0.0
0.05942480245
-
0.2178361212
-
0.9741745817
0.0
0.0
0.0
0.9973986606
-
0.02692981555
0.06686327745
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
]
3
White
[
10.09116301
0.3556701722
1.333166815
1.442000000
1.344000000
0.02200000000
]
[
0.04028666644
0.9048796003
-
0.4237568827
0.0
0.0
0.0
0.06399313350
-
0.4255671399
-
0.9026613361
0.0
0.0
0.0
0.9971368323
-
0.009247686096
0.07505076253
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
]
4
Khaya
[
10.80102614
0.4606834511
1.475290392
1.800000000
1.196000000
0.4200000000
]
[
0.05368516527
0.9266208315
-
0.3721447810
0.0
0.0
0.0
0.07220768570
-
0.3753089731
-
0.9240829102
0.0
0.0
0.0
0.9959437502
-
0.02273783078
0.08705767064
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
]
5
Mahogany
[
12.10253321
0.6098853915
1.810581412
1.960000000
1.498000000
0.6380000000
]
[
-
0.06484967920
-
0.8666704601
-
0.4946481887
0.0
0.0
0.0
-
0.07817061387
0.4985803653
-
0.8633116316
0.0
0.0
0.0
-
0.9948285648
0.01731853005
0.1000809391
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
]
6
Mahogany
[
13.10854783
0.3895295651
1.814922632
1.218000000
0.9380000000
0.3000000000
]
[
-
0.05490172356
-
0.8720473797
-
0.4863323643
0.0
0.0
0.0
-
0.07547523076
0.4892979419
-
0.8688446422
0.0
0.0
0.0
-
0.9956351187
0.01099502014
0.09268127742
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
]
7
S. Germ
[
12.34053410
0.5212570563
1.922208822
2.318000000
1.582000000
0.5400000000
]
[
0.04967528691
0.9161715971
-
0.3976958243
0.0
0.0
0.0
0.08463937788
-
0.4006166011
-
0.9123280736
0.0
0.0
0.0
0.9951726186
-
0.01165943224
0.09744493938
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
]
8
MAPLE
[
12.05449768
0.6868279555
2.766674361
2.460000000
2.194000000
0.5840000000
]
[
-
0.1387533797
-
0.8476774112
-
0.5120454135
0.0
0.0
0.0
-
0.2031928413
0.5304148448
-
0.8230265872
0.0
0.0
0.0
-
0.9692575344
0.01015375725
0.2458388325
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
]Table 9
In continuation to Table8.
S. no.
Hardwood species
EigenvaluesΛ
[
i
]
EigenvectorsV
[
i
]
9
Walnut
[
12.67869546
0.5204134969
1.919891015
1.922000000
1.400000000
0.4600000000
]
[
0.08621257414
0.8755641188
-
0.4753471023
0.0
0.0
0.0
0.1243595697
-
0.4828494241
-
0.8668282006
0.0
0.0
0.0
0.9884847438
-
0.01561753067
0.1505124700
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
]
10
Birch
[
3.168908680
0.4747157903
0.5946755301
2.346000000
1.816000000
0.3720000000
]
[
0.3961520086
0.9160614623
0.06240981414
0.0
0.0
0.0
0.6338768023
-
0.2236797319
-
0.7403833993
0.0
0.0
0.0
0.6642768888
-
0.3328645006
0.6692812840
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
]
11
Y. Birch
[
15.45414040
0.5549651499
2.059894445
2.120000000
1.942000000
0.4800000000
]
[
-
0.06591789198
-
0.8289381135
-
0.5554425577
0.0
0.0
0.0
-
0.08975775227
0.5593224991
-
0.8240763851
0.0
0.0
0.0
-
0.9937798435
0.004466103673
0.1112729817
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
]
12
Oak
[
17.18666431
0.8648974218
3.239438239
2.380000000
1.532000000
0.7840000000
]
[
0.06975010637
0.9079000793
-
0.4133429180
0.0
0.0
0.0
0.1071019008
-
0.4187725641
-
0.9017531389
0.0
0.0
0.0
0.9917984199
-
0.01862756541
0.1264472547
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
]
13
Ash
[
17.15567967
0.6696042476
2.409716064
2.684000000
1.784000000
0.5400000000
]
[
0.06190313535
0.8746549514
-
0.4807772027
0.0
0.0
0.0
0.09781749768
-
0.4846986500
-
0.8691944279
0.0
0.0
0.0
0.9932772717
-
0.006777439445
0.1155609239
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
]
14
Ash
[
17.42363268
0.7844815431
2.669885744
1.720000000
1.218000000
0.5000000000
]
[
-
0.1004077048
-
0.8554189669
-
0.5081108976
0.0
0.0
0.0
-
0.1363859604
0.5177044741
-
0.8446188165
0.0
0.0
0.0
-
0.9855542417
0.01550704374
0.1686486565
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
]
15
Beech
[
15.99002755
0.9558575345
3.451114905
3.216000000
2.112000000
0.9120000000
]
[
0.1135176976
0.8848520771
-
0.4518302069
0.0
0.0
0.0
0.1764907831
-
0.4654963442
-
0.8672739791
0.0
0.0
0.0
0.9777344911
-
0.01870707734
0.2090103088
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
]Table 10
The eigenvalues and eigenvectors for 8 softwoods species calculated using MAPLE.Λ
[
i
] stand for the eigenvalues of ith softwood species, and  V
[
i
] stand for the corresponding eigenvectors.
S. no.
Softwood species
EigenvaluesΛ
[
i
]
EigenvectorsV
[
i
]
1
Blasa
[
6.385324208
0.09845837744
0.3832174064
0.6240000000
0.4060000000
0.06600000000
]
[
0.01488818113
0.9510211778
-
0.3087669978
0.0
0.0
0.0
0.02575997614
-
0.3090635474
-
0.9506924583
0.0
0.0
0.0
0.9995572851
-
0.006200252135
0.02909967665
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
]
2
Spruce
[
11.98583568
0.4516442734
1.114520065
1.498000000
1.442000000
0.07800000000
]
[
-
0.03300264531
-
0.9144925828
-
0.4032544375
0.0
0.0
0.0
-
0.04689865645
0.4044467539
-
0.9133582749
0.0
0.0
0.0
-
0.9983543167
0.01123114838
0.05623628701
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
]
3
Spruce
[
14.10927768
0.4185382987
1.227183961
1.442000000
1.0
0.06400000000
]
[
-
0.03651783317
-
0.8968065074
-
0.4409132956
0.0
0.0
0.0
-
0.05361649666
0.4423303564
-
0.8952480817
0.0
0.0
0.0
-
0.9978936411
0.009052295016
0.06423657043
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
]
4
Spruce
[
16.30529953
0.3544288592
0.8442715844
1.234000000
1.520000000
0.07200000000
]
[
-
0.02057139660
-
0.9124371483
-
0.4086994866
0.0
0.0
0.0
-
0.02869707087
0.4091564350
-
0.9120128765
0.0
0.0
0.0
-
0.9993764538
0.007032901380
0.03460119282
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
]
5
Spruce
[
17.25845208
0.5092652753
1.171282625
1.250000000
1.706000000
0.07000000000
]
[
-
0.03391880763
-
0.8104111857
-
0.5848788055
0.0
0.0
0.0
-
0.03428302033
0.5858146064
-
0.8097196604
0.0
0.0
0.0
-
0.9988364173
0.007413313465
0.04765347716
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
]
6
Douglas Fir
[
16.13459769
0.6233733771
1.439028943
1.767000000
1.766000000
0.1760000000
]
[
-
0.03639420861
-
0.8057068065
-
0.5911954018
0.0
0.0
0.0
-
0.03697194691
0.5922678885
-
0.8048924305
0.0
0.0
0.0
-
0.9986533619
0.007435777607
0.05134366998
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
]
7
Douglas Fir
[
17.10150156
0.6986859431
2.204812526
2.348000000
1.816000000
0.1600000000
]
[
-
0.04991838192
-
0.8212998538
-
0.5683086323
0.0
0.0
0.0
-
0.06364825251
0.5704773676
-
0.8188433771
0.0
0.0
0.0
-
0.9967231590
0.004703484281
0.08075160717
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
]
8
Pine
[
16.99539133
0.4940019684
1.565606760
3.484000000
1.344000000
0.1320000000
]
[
-
0.03436105770
-
0.8973128643
-
0.4400556096
0.0
0.0
0.0
-
0.05585205149
0.4413516031
-
0.8955943895
0.0
0.0
0.0
-
0.9978476167
0.006195562437
0.06528207193
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
]Table 11
The eigenvalues and eigenvectors for 3 specimens of cancellous bone calculated using MAPLE.Λ
[
i
] stand for the eigenvalues of ith specimen of cancellous bone, and V
[
i
] stand for the corresponding eigenvectors.
Cancellous bone
EigenvaluesΛ
[
i
]
EigenvectorsV
[
i
]
Specimen 1
[
13.34869880
-
0.2071167473
8.644210768
6.837676373
3.768806080
4.123724855
]
[
0.3316115964
0.8831238953
0.2356441950
0.2242668491
0.05354786914
0.03787816350
0.6475092156
-
0.1099136965
0.1556396475
-
0.6940063479
-
0.1278953723
0.2154647417
0.4800906963
-
0.2832250652
0.07648827487
0.2494278793
0.6410471445
-
0.4585742261
-
0.2737532692
-
0.06354940263
0.4518342532
-
0.06214920576
0.5836952694
0.6101669758
-
0.3706582070
0.3314035780
-
0.1276606259
-
0.6337030234
0.3639432801
-
0.4499474345
0.1671829224
0.1180163925
-
0.8330343502
0.02008124471
0.3109325968
0.4087706408
]
Specimen 2
[
1754.634171
777.3848354
155.0688791
230.1790452
344.5174379
358.9156342
]
[
-
0.3057290719
-
0.3005337771
-
0.9030341401
0.01584500766
-
0.02175914011
-
0.003741933916
-
0.4311831740
-
0.8021570593
0.4130146809
0.0001362645283
-
0.003854432451
0.005396441669
-
0.8488209390
0.5154657137
0.1152135966
-
0.004741957786
0.02229739182
-
0.002068264586
0.0009402467441
-
0.005358584867
0.003987518000
-
0.03989381331
0.03796849613
-
0.9984595017
-
0.01058157817
0.02100470582
0.02092543658
0.006230946613
-
0.9987520802
-
0.03826770492
0.0009823402498
0.006977574167
0.01484146906
0.9990475909
0.008196718814
-
0.03958286493
]
Specimen 3
[
1418.542670
583.8388363
99.59058231
343.3818493
302.4968125
175.7492470
]
[
0.3244960223
-
0.004831643783
-
0.8451713783
0.4062092722
0.1164524794
-
0.04239310198
0.6461152835
-
0.7125334236
0.2668138654
0.04131783689
-
0.02431040234
0.03665187401
0.6621315621
0.7009969173
0.2633438279
0.02836289512
-
0.001711614840
0.005269111592
-
0.1962401199
0.009159859564
0.3740839141
0.8850746598
0.1904577552
0.04284654369
-
0.008597323304
-
0.002526763452
0.01897677365
0.1738799329
-
0.6977939665
-
0.6945566457
0.01533190599
-
0.02803230310
0.06961922985
-
0.1374446216
0.6801869735
-
0.7159517722
]


### 5.2. Computing the Average Elasticity Tensors for Woods and Cancellous Bone Using MAPLE

In order to compute average elasticity tensors for the hardwoods, softwoods, and cancellous bone, let us go through a brief mathematical description regarding this issue [3], as this notion helps us developing an appropriate MAPLE code to have desired results.Suppose, we haveM measurements of the elastic constants data for some particular species or material, C
^
I
,
C
^
I
I
,
…
,
C
^
M and we want to construct a tensor C
AVG representing an average elasticity tensor for that particular material. Then to do so, we need the following key algebra [3].(i)
The eigenvaluesΛ
k
Y, k
=
1,2
,
…
,
6, Y
=
1,2
,
…
,
M and the eigenvectors N
^
k
Y, k
=
1,2
,
…
,
6, Y
=
1,2
,
…
,
M for each of the Yth measurement of the particular species or material.
(ii)
The nominal average (NA),N
^
k
NA of the eigenvectors N
^
k
Y associated with the particular species. The nominal average is given by
(45)
N
^
k
NA
≡
1
M
∑
Y
=
1
M
‍
N
^
k
Y
.
(iii)
The averageN
^
k
AVG, which is obtained by the following formula:
(46)
N
^
k
AVG
≡
J
^
NA
N
^
k
NA
,
whereJ
^
NA stands for the inverse square root of the positive definite tensor |
∑
k
=
1
6
‍
N
^
k
NA
⊗
N
^
k
NA
|; that is,
(47)
J
^
NA
≡
(
∑
k
=
1
6
‍
N
^
k
NA
⊗
N
^
k
NA
)
-
1
/
2
.
(iv)
Now, we proceed to average the eigenvalue. For this, we need to transform each set of eigenvalues corresponding to each eigenvector to the average eigenvectorN
^
k
AVG; that is,
(48)
Λ
q
AVG
≡
1
M
∑
Y
=
1
M
‍
∑
k
=
1
6
‍
Λ
k
Y
(
N
^
k
Y
·
N
^
q
AVG
)
2
.
It is obvious thatΛ
k
AVG
>
0 for all k
=
1,2
,
…
,
6 and Λ
k
AVG
≠
Λ
k
NA, where Λ
k
NA is the nominal average of eigenvalues of material and is defined by
(49)
Λ
k
NA
≡
1
M
∑
Y
=
1
M
‍
Λ
k
Y
.
(v)
Eventually, we can now calculate the average elasticity matrixC
^
AVG in such a way that
(50)
C
^
AVG
=
∑
k
=
1
6
‍
Λ
k
AVG
N
^
k
AVG
⊗
N
^
k
AVG
.Taking care of the above outlined steps, we have developed some MAPLE codes (See Step 17 to Step 81 of Appendix A) that enable us to calculate the nominal averages, average eigenvectors, average eigenvalues, and the average elasticity matrices for each of the 15 hardwood species tabulated in Table 5. In addition to this, Appendix A also retains few steps, for example, Steps 21, 26, 31, 36, 41, 46, 51, 56, 61, 66, 73, and 80, which are particularly developed to plot histograms of average elasticity matrices of each of the hardwood species as well as graphs between I, II, III, IV, V, and VI eigenvalues of all hardwood species and the apparent densities (see Table 5).We summarize the above claimed consequences of each of the hardwood species one by one ss follows.


#### 5.2.1. MAPLE Evaluations for Quipo

(
1
)
The nominal averages of eigenvectors for the two measurements (S. nos. 1 and 2 of Table5) of Quipo are
(51)
[
0.0367834559700000
0.0453681054800000
0.998185540700000
0.0
0.0
0.0
]
,
[
0.983745905000000
-
0.170879918750000
-
0.0277901141300000
0.0
0.0
0.0
]
,
[
-
0.169284222800000
-
0.982986098000000
0.0509905132200000
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
2
)
The average eigenvectors for Quipo are(52)
[
0.03679134179
0.04537783172
0.9983995367
0.0
0.0
0.0
]
,
[
0.9859858256
-
0.1712690003
-
0.02785339027
0.0
0.0
0.0
]
,
[
-
0.1697052873
-
0.9854311019
0.05111734310
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
3
)
The averages eigenvalue for the two measurements (S. nos. 1 and 2 of Table5) of Quipo is 3.339500000.
(
4
)
The average elasticity matrix for Quipo is(53)
[
3.34725276837330
0.000112116752981933
0.00198542359441849
0.0
0.0
0.0
0.000112116752981933
3.34773746006714
-
0.000991802322788501
0.0
0.0
0.0
0.00198542359441849
-
0.000991802322788501
3.34013593769726
0.0
0.0
0.0
0.0
0.0
0.0
3.33950000000000
0.0
0.0
0.0
0.0
0.0
0.0
3.33950000000000
0.0
0.0
0.0
0.0
0.0
0.0
3.33950000000000
]
.(
5
)
The histogram of the average elasticity matrix for Quipo is shown in Figure2.Figure 2
Histogram showing average elasticity matrix of Quipo.


#### 5.2.2. MAPLE Evaluations for White

(
1
)
The nominal averages of eigenvectors for the single measurements (S. no. 3 of Table5) of White are(54)
[
0.04028666644
0.06399313350
0.9971368323
0.0
0.0
0.0
]
,
[
0.9048796003
-
0.4255671399
-
0.009247686096
0.0
0.0
0.0
]
,
[
-
0.4237568827
-
0.9026613361
0.07505076253
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
2
)
The average eigenvectors for White are(55)
[
0.04028666644
0.06399313350
0.9971368323
0.0
0.0
0.0
]
,
[
0.9048795994
-
0.4255671395
-
0.009247686087
0.0
0.0
0.0
]
,
[
-
0.4237568827
-
0.9026613361
0.07505076253
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
3
)
The average eigenvalue for the single measurement (S. no. 3 of Table5) of White is 14.58799998.
(
4
)
The average elasticity matrix for White is given as(56)
[
14.58799998
0.00000001785571215
-
0.00000001009489597
0.0
0.0
0.0
0.00000001785571215
14.58799997
0.000000002407020197
0.0
0.0
0.0
-
0.00000001009489597
0.000000002407020197
14.58799997
0.0
0.0
0.0
0.0
0.0
0.0
14.58799998
0.0
0.0
0.0
0.0
0.0
0.0
14.58799998
0.0
0.0
0.0
0.0
0.0
0.0
14.58799998
]
.(
5
)
The histogram of the average elasticity matrix for White is shown in Figure3.Figure 3
Histogram showing average elasticity matrix of white.


#### 5.2.3. MAPLE Evaluations for Khaya

(
1
)
The nominal averages of eigenvectors for the single measurements (S. no. 4 of Table5) of Khaya are(57)
[
0.05368516527
0.07220768570
0.9959437502
0.0
0.0
0.0
]
,
[
0.9266208315
-
0.3753089731
-
0.02273783078
0.0
0.0
0.0
]
,
[
-
0.3721447810
-
0.9240829102
0.08705767064
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
.
(
2
)
The average eigenvectors for Khaya are(58)
[
0.05368516527
0.07220768570
0.9959437502
0.0
0.0
0.0
]
,
[
0.9266208315
-
0.3753089731
-
0.02273783078
0.0
0.0
0.0
]
,
[
-
0.3721447806
-
0.9240829093
0.08705767055
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
3
)
The average eigenvalue for the single measurement (S. no. 4 of Table5) of Khaya is 16.15299998.
(
4
)
The average elasticity matrix for Khaya is given as(59)
[
16.15299998
0.000000005508173090
-
0.000000008722620038
0.0
0.0
0.0
0.000000005508173090
16.15299996
-
0.000000004345156817
0.0
0.0
0.0
-
0.000000008722620038
-
0.000000004345156817
16.15300000
0.0
0.0
0.0
0.0
0.0
0.0
16.15299998
0.0
0.0
0.0
0.0
0.0
0.0
16.15299998
0.0
0.0
0.0
0.0
0.0
0.0
16.15299998
]
.(
5
)
The histogram of the average elasticity matrix for Khaya is shown in Figure4.Figure 4
Histogram showing average elasticity matrix of Khaya.


#### 5.2.4. MAPLE Evaluations for Mahogany

(
1
)
The nominal averages of eigenvectors for the two measurements (S. nos. 5 and 6 of Table5) of Mahogany are(60)
[
-
0.0598757013800000
-
0.0768229223150000
-
0.995231841750000
0.0
0.0
0.0
]
,
[
-
0.869358919900000
0.493939153600000
0.0141567750950000
0.0
0.0
0.0
]
,
[
-
0.490490276500000
-
0.866078136900000
0.0963811082600000
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.(
2
)
The average eigenvectors for Mahogany are(61)
[
-
0.05987730132
-
0.07682497510
-
0.9952584354
0.0
0.0
0.0
]
,
[
-
0.8693926232
0.4939583026
0.01415732393
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
3
)
The average eigenvalue for the two measurements (S. nos. 5 and 6 of Table5) of Mahogany is 18.19400002.
(
4
)
The average elasticity matrix for Mahogany is given as(62)
[
18.19451173
-
0.0001438038661
0.0001358556898
0.0
0.0
0.0
-
0.0001438038661
18.19484018
-
0.0004764690574
0.0
0.0
0.0
0.0001358556898
-
0.0004764690574
18.19454255
0.0
0.0
0.0
0.0
0.0
0.0
18.19400002
18.19400002
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
18.19400002
]
.(
5
)
The histogram of the average elasticity matrix for Mahogany is shown in Figure5.Figure 5
Histogram showing average elasticity matrix of Mahogany.


#### 5.2.5. MAPLE Evaluations for S. Germ

(
1
)
The nominal averages of eigenvectors for the single measurement (S. no. 7 of Table5) of S. Gem are
(63)
[
0.04967528691
0.08463937788
0.9951726186
0.0
0.0
0.0
]
,
[
0.9161715971
-
0.4006166011
-
0.01165943224
0.0
0.0
0.0
]
,
[
-
0.3976958243
-
0.9123280736
0.09744493938
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.(
2
)
The average eigenvectors for S. Germ are(64)
[
0.04967528696
0.08463937796
0.9951726196
0.0
0.0
0.0
]
,
[
0.9161715980
-
0.4006166015
-
0.01165943225
0.0
0.0
0.0
]
,
[
-
0.3976958247
-
0.9123280745
0.09744493948
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
3
)
The average eigenvalue for the single measurement (S. no. 7 of Table5) of S. Germ is 19.22399998.
(
4
)
The average elasticity matrix for S. Germ is given as(65)
[
19.22399998
-
0.00000001176508811
-
0.00000001537920006
0.0
0.0
0.0
-
0.00000001176508811
19.22400000
-
0.000000007631928036
0.0
0.0
0.0
-
0.00000001537920006
-
0.000000007631928036
19.22400000
0.0
0.0
0.0
0.0
0.0
0.0
19.22399998
19.22399998
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
19.22399998
]
.(
5
)
The histogram of the average elasticity matrix for S. Germ is shown in Figure6.Figure 6
Histogram showing average elasticity matrix of S. Germ.


#### 5.2.6. MAPLE Evaluations for maple

(
1
)
The nominal averages of eigenvectors for the single measurement (S. no. 8 of Table5) of maple are
(66)
[
0.0
0.0
0.0
0.0
0.0
1.0
]
,
[
-
0.1387533797
-
0.2031928413
-
0.9692575344
0.0
0.0
0.0
]
,
[
-
0.8476774112
0.5304148448
0.01015375725
0.0
0.0
0.0
]
,
[
-
0.5120454135
-
0.8230265872
0.2458388325
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
.(
2
)
The average eigenvectors for maple are(67)
[
0.0
0.0
0.0
0.0
0.0
1.0
]
,
[
-
0.1387533798
-
0.2031928415
-
0.9692575354
0.0
0.0
0.0
]
,
[
-
0.8476774112
0.5304148448
0.01015375725
0.0
0.0
0.0
]
,
[
-
0.5120454135
-
0.8230265872
0.2458388325
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
.
(
3
)
The average eigenvalue for the single measurement (S. no. 8 of Table5) of maple is 20.74600000.
(
4
)
The average elasticity matrix for maple is given as(68)
[
20.74599999
0.000000004356660361
0.00000003402343994
0.0
0.0
0.0
0.000000004356660361
20.74600004
0.00000002232269567
0.0
0.0
0.0
0.00000003402343994
0.00000002232269567
20.74600000
0.0
0.0
0.0
0.0
0.0
0.0
20.74600000
20.74600000
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
20.74600000
]
.(
5
)
The histogram of the average elasticity matrix for MAPLE is shown in Figure7.Figure 7
Histogram showing average elasticity matrix of maple.


#### 5.2.7. MAPLE Evaluations for Walnut

(
1
)
The nominal averages of eigenvectors for the single measurement (S. no. 9 of Table5) of Walnut are
(69)
[
0.08621257414
0.1243595697
0.9884847438
0.0
0.0
0.0
]
,
[
0.8755641188
-
0.4828494241
-
0.01561753067
0.0
0.0
0.0
]
,
[
-
0.4753471023
-
0.8668282006
0.1505124700
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.(
2
)
The average eigenvectors for Walnut are(70)
[
0.08621257423
0.1243595698
0.9884847448
0.0
0.0
0.0
]
,
[
0.8755641188
-
0.4828494241
-
0.01561753067
0.0
0.0
0.0
]
,
[
-
0.4753471018
-
0.8668281997
0.1505124698
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
3
)
The average eigenvalue for the single measurement (S. no. 9 of Table5) of walnut is 18.90099997.
(
4
)
The average elasticity matrix for Walnut is given as(71)
[
18.90099999
0.00000001209663993
-
0.00000002532733996
0.0
0.0
0.0
0.00000001209663993
18.90099991
0.000000001701090138
0.0
0.0
0.0
-
0.00000002532733996
0.000000001701090138
18.90100001
0.0
0.0
0.0
0.0
0.0
0.0
18.90099997
18.90099997
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
18.90099997
]
.(
5
)
The histogram of the average elasticity matrix for Walnut is shown in Figure8.Figure 8
Histogram showing average elasticity matrix of Walnut.


#### 5.2.8. MAPLE Evaluations for Birch

(
1
)
The nominal averages of eigenvectors for the single measurement (S. no. 10 of Table5) of Birch are
(72)
[
0.3961520086
0.6338768023
0.6642768888
0.0
0.0
0.0
]
,
[
0.9160614623
-
0.2236797319
-
0.3328645006
0.0
0.0
0.0
]
,
[
0.06240981414
-
0.7403833993
0.6692812840
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
2
)
The average eigenvectors for Birch are(73)
[
0.3961520086
0.6338768023
0.6642768888
0.0
0.0
0.0
]
,
[
0.9160614614
-
0.2236797317
-
0.3328645003
0.0
0.0
0.0
]
,
[
0.06240981426
-
0.7403834008
0.6692812853
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
3
)
The average eigenvalue for the single measurement (S. no. 10 of Table5) of Birch is 8.772299998.
(
4
)
The average elasticity matrix for Birch is given as
(74)
[
8.772299997
-
0.00000003535236899
0.00000003421196978
0.0
0.0
0.0
-
0.00000003535236899
8.772300024
-
0.00000001649192439
0.0
0.0
0.0
0.00000003421196978
-
0.00000001649192439
8.772299994
0.0
0.0
0.0
0.0
0.0
0.0
8.772299998
8.772299998
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
8.772299998
]
.(
5
)
The histogram of the average elasticity matrix for Birch is shown in Figure9.Figure 9
Histogram showing average elasticity matrix of Birch.


#### 5.2.9. MAPLE Evaluations for Y. Birch

(
1
)
The nominal averages of eigenvectors for the single measurement (S. no. 11 of Table5) of Y. Birch are
(75)
[
-
0.06591789198
-
0.08975775227
-
0.9937798435
0.0
0.0
0.0
]
,
[
-
0.8289381135
0.5593224991
0.004466103673
0.0
0.0
0.0
]
,
[
-
0.5554425577
-
0.8240763851
0.1112729817
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
2
)
The average eigenvectors for Y. Birch are(76)
[
-
0.1387533798
-
0.2031928415
-
0.9692575354
0.0
0.0
0.0
]
,
[
-
0.8476774112
0.5304148448
0.01015375725
0.0
0.0
0.0
]
,
[
-
0.5120454135
-
0.8230265872
0.2458388325
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
3
)
The average eigenvalue for the single measurement (S. no. 11 of Table5) of Y. Birch is 22.28057495.
(
4
)
The average elasticity matrix for Y. Birch is given as(77)
[
22.28057494
0.000000004678921127
0.00000003654014285
0.0
0.0
0.0
0.000000004678921127
22.28057499
0.00000002397389829
0.0
0.0
0.0
0.00000003654014285
0.00000002397389829
22.28057495
0.0
0.0
0.0
0.0
0.0
0.0
22.28057495
22.28057495
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
22.28057495
]
.(
5
)
The histogram of the average elasticity matrix for Y. Birch is shown in Figure10.Figure 10
Histogram showing average elasticity matrix of Y. Birch.


#### 5.2.10. MAPLE Evaluations for Oak

(
1
)
The nominal averages of eigenvectors for the single measurement (S. no. 12 of Table5) of Oak are(78)
[
0.06975010637
0.1071019008
0.9917984199
0.0
0.0
0.0
]
,
[
0.9079000793
-
0.4187725641
-
0.01862756541
0.0
0.0
0.0
]
,
[
-
0.4133429180
-
0.9017531389
0.1264472547
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.(
2
)
The average eigenvectors for Oak are(79)
[
0.06975010637
0.1071019008
0.9917984199
0.0
0.0
0.0
]
,
[
0.9079000793
-
0.4187725641
-
0.01862756541
0.0
0.0
0.0
]
,
[
-
0.4133429184
-
0.9017531398
0.1264472548
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
3
)
The average eigenvalue for the single measurement (S. no. 12 of Table5) of Oak is 25.98699998.
(
4
)
The average elasticity matrix for Oak is given as(80)
[
25.98699997
-
0.00000001889254939
-
0.000000003638179937
0.0
0.0
0.0
-
0.00000001889254939
25.98700006
0.000000008575709980
0.0
0.0
0.0
-
0.000000003638179937
0.000000008575709980
25.98699998
0.0
0.0
0.0
0.0
0.0
0.0
25.98699998
25.98699998
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
25.98699998
]
.(
5
)
The histogram of the average elasticity matrix for Oak is shown in Figure11.Figure 11
Histogram showing average elasticity matrix of Oak.


#### 5.2.11. MAPLE Evaluations for Ash

(
1
)
The nominal averages of eigenvectors for the two measurements (S. nos. 13 and 14 of Table5) of Ash are
(81)
[
-
0.0192522847250000
-
0.0192842313600000
0.00386151499999998
0.0
0.0
0.0
]
,
[
0.00961799224999998
0.0165029120500000
0.00436480214750000
0.0
0.0
0.0
]
,
[
-
0.494444050150000
-
0.856906622200000
0.142104790200000
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
2
)
The average eigenvectors for Ash are(82)
[
-
25.41745843
-
25.45963537
5.098090872
0.0
0.0
0.0
]
,
[
25.05315863
42.98714977
11.36953303
0.0
0.0
0.0
]
,
[
-
0.4949599718
-
0.8578007511
0.1422530677
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
3
)
The average eigenvalue for the two measurements (S. nos. 13 and 14 of Table5) of Ash is 24.77950014.
(
4
)
The average elasticity matrix for Ash is given as(83)
[
31567.9169478799
42732.4383657779
3845.57503470365
0.0
0.0
0.0
42732.4383657779
61870.0481877479
8891.53535276175
0.0
0.0
0.0
3845.57503470365
8891.53535276175
3847.68762708609
0.0
0.0
0.0
0.0
0.0
0.0
24.7795001400000
0.0
0.0
0.0
0.0
0.0
0.0
24.7795001400000
0.0
0.0
0.0
0.0
0.0
0.0
24.7795001400000
]
.(
5
)
The histogram of the average elasticity matrix for Ash is shown in Figure12.Figure 12
Histogram showing average elasticity matrix of Ash.


#### 5.2.12. MAPLE Evaluations for Beech

(
1
)
The nominal averages of eigenvectors for the single measurements (S. no. 15 of Table5) of Beech are
(84)
[
0.1135176976
0.1764907831
0.9777344911
0.0
0.0
0.0
]
,
[
0.8848520771
-
0.4654963442
-
0.01870707734
0.0
0.0
0.0
]
,
[
-
0.4518302069
-
0.8672739791
0.2090103088
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
2
)
The average eigenvectors for Beech are(85)
[
0.1135176977
0.1764907833
0.9777344921
0.0
0.0
0.0
]
,
[
0.8848520771
-
0.4654963442
-
0.01870707734
0.0
0.0
0.0
]
,
[
-
0.4518302069
-
0.8672739791
0.2090103088
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
3
)
The average eigenvalue for the single measurements (S. no. 15 of Table5) of Beech is 26.63699998.
(
4
)
The average elasticity matrix for Beech is given as(86)
[
26.63700003
0.00000004768023095
0.00000003169803038
0.0
0.0
0.0
0.00000004768023095
26.63699992
0.000000008310743498
0.0
0.0
0.0
0.00000003169803038
0.000000008310743498
26.63700001
0.0
0.0
0.0
0.0
0.0
0.0
26.63699998
26.63699998
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
26.63699998
]
.(
5
)
The histogram of the average elasticity matrix for Beech is shown in Figure13.Figure 13
Histogram showing average elasticity matrix of Beech.The graphics placed in left as well as right positions, showing graphs of I and II eigenvalues of all 15 hardwood species against their apparent densities.
(a)
(b)The graphics placed in left as well as right positions, showing graphs of III and IV eigenvalues of all 15 hardwood species against their apparent densities.
(a)
(b)The graphics placed in left as well as right positions, showing graphs of V and VI eigenvalues of all 15 hardwood species against their apparent densities.
(a)
(b)Figure 17
A combined view of Figures14, 15, and 16.


### 5.3. MAPLE Plotted Graphics for I, II, III, IV, V, and VI Eigenvalues of 15 Hardwood Species against Their Apparent Densities

This subsection is dedicated to the exposition of the extreme quality of MAPLE, which can enable one to simultaneously sketch up graphics for I, II, III, IV, V, and VI eigenvalues of each of the 15 hardwood species (See Tables8 and 9) against the corresponding apparent densities ρ (See Table 5 for apparent density). The MAPLE code regarding this issue has been mentioned in Step 81 of Appendix A.


## 5.1. Computing the Eigenvalues and Eigenvectors for Hardwoods Species, Softwood Species, and Cancellous Bone Using MAPLE

Here we present the outcomes of MAPLE code (which is mentioned in Appendix A) regarding the computation of eigenvalues and eigenvectors of the stiffness matricesC
^ concerning 15 hardwood species (see Table 5).It is well known that the eigenvalues and eigenvectors of a stiffness matrixC
^ or its compliance matrix S
^ are determined from the equations [3]
(44)
(
C
^
-
Λ
I
)
N
^
=
0
or
(
S
^
-
1
Λ
I
)
N
^
=
0
,
where I is the 6 × 6 identity matrix and N
^ represents the normalized eigenvectors of C
^ or S
^. Naturally, C
^ or S
^ being positive definite, therefore there would be six positive eigenvalues and these values are known as Kelvin’s moduli. These Kelvin’s moduli are denoted by Λ
i
,
i
=
1,2
,
…
,
6 and if possible are ordered in such a way that Λ
1
≥
Λ
2
≥
⋯
≥
Λ
6
>
0. Also, the eigenvalues of S
^ can simply de computed by taking the inverse of the eigenvalues of C
^.A MAPLE code mentioned in Step 16 of Appendix A enables us to simultaneously compute the eigenvalues and eigenvectors for all the 15 hardwood species. Also, the results of this MAPLE code are summarized in Tables8 and 9. Further, by simply substituting the elasticity constants from Table 6 in Step 2 to Step 11 of Appendix A, and performing Step 16 again, one can have the eigenvalues and eigenvectors for the 8 softwood species. The computed results for softwood species are summarized in Table 10. Similarly, substitution of elastic constants data for cancellous bone from Table 7 into Step 2 to Step 11 and repeating Step 16 of the MAPLE code yields the desired results tabulated in Table 11.Table 8
The eigenvalues and eigenvectors for 15 hardwoods species calculated using MAPLE.Λ
[
i
] stand for the eigenvalues of ith softwood species, and V
[
i
] stand for the corresponding eigenvectors.
S. no.
Hardwood species
EigenvaluesΛ
[
i
]
EigenvectorsV
[
i
]
1
Quipo
[
1.076866026
0.04067345638
0.2534605191
0.2260000000
0.1180000000
0.07800000000
]
[
0.03276733390
0.9918780499
-
0.1228992897
0.0
0.0
0.0
0.03131140851
-
0.1239237163
-
0.9917976143
0.0
0.0
0.0
0.9989724208
-
0.02865041271
0.03511774899
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
]
2
Quipo
[
3.461963876
0.1399776173
0.4300584948
0.4300000000
0.2800000000
0.1440000000
]
[
0.04079957804
0.9756137601
-
0.2156691559
0.0
0.0
0.0
0.05942480245
-
0.2178361212
-
0.9741745817
0.0
0.0
0.0
0.9973986606
-
0.02692981555
0.06686327745
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
]
3
White
[
10.09116301
0.3556701722
1.333166815
1.442000000
1.344000000
0.02200000000
]
[
0.04028666644
0.9048796003
-
0.4237568827
0.0
0.0
0.0
0.06399313350
-
0.4255671399
-
0.9026613361
0.0
0.0
0.0
0.9971368323
-
0.009247686096
0.07505076253
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
]
4
Khaya
[
10.80102614
0.4606834511
1.475290392
1.800000000
1.196000000
0.4200000000
]
[
0.05368516527
0.9266208315
-
0.3721447810
0.0
0.0
0.0
0.07220768570
-
0.3753089731
-
0.9240829102
0.0
0.0
0.0
0.9959437502
-
0.02273783078
0.08705767064
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
]
5
Mahogany
[
12.10253321
0.6098853915
1.810581412
1.960000000
1.498000000
0.6380000000
]
[
-
0.06484967920
-
0.8666704601
-
0.4946481887
0.0
0.0
0.0
-
0.07817061387
0.4985803653
-
0.8633116316
0.0
0.0
0.0
-
0.9948285648
0.01731853005
0.1000809391
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
]
6
Mahogany
[
13.10854783
0.3895295651
1.814922632
1.218000000
0.9380000000
0.3000000000
]
[
-
0.05490172356
-
0.8720473797
-
0.4863323643
0.0
0.0
0.0
-
0.07547523076
0.4892979419
-
0.8688446422
0.0
0.0
0.0
-
0.9956351187
0.01099502014
0.09268127742
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
]
7
S. Germ
[
12.34053410
0.5212570563
1.922208822
2.318000000
1.582000000
0.5400000000
]
[
0.04967528691
0.9161715971
-
0.3976958243
0.0
0.0
0.0
0.08463937788
-
0.4006166011
-
0.9123280736
0.0
0.0
0.0
0.9951726186
-
0.01165943224
0.09744493938
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
]
8
MAPLE
[
12.05449768
0.6868279555
2.766674361
2.460000000
2.194000000
0.5840000000
]
[
-
0.1387533797
-
0.8476774112
-
0.5120454135
0.0
0.0
0.0
-
0.2031928413
0.5304148448
-
0.8230265872
0.0
0.0
0.0
-
0.9692575344
0.01015375725
0.2458388325
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
]Table 9
In continuation to Table8.
S. no.
Hardwood species
EigenvaluesΛ
[
i
]
EigenvectorsV
[
i
]
9
Walnut
[
12.67869546
0.5204134969
1.919891015
1.922000000
1.400000000
0.4600000000
]
[
0.08621257414
0.8755641188
-
0.4753471023
0.0
0.0
0.0
0.1243595697
-
0.4828494241
-
0.8668282006
0.0
0.0
0.0
0.9884847438
-
0.01561753067
0.1505124700
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
]
10
Birch
[
3.168908680
0.4747157903
0.5946755301
2.346000000
1.816000000
0.3720000000
]
[
0.3961520086
0.9160614623
0.06240981414
0.0
0.0
0.0
0.6338768023
-
0.2236797319
-
0.7403833993
0.0
0.0
0.0
0.6642768888
-
0.3328645006
0.6692812840
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
]
11
Y. Birch
[
15.45414040
0.5549651499
2.059894445
2.120000000
1.942000000
0.4800000000
]
[
-
0.06591789198
-
0.8289381135
-
0.5554425577
0.0
0.0
0.0
-
0.08975775227
0.5593224991
-
0.8240763851
0.0
0.0
0.0
-
0.9937798435
0.004466103673
0.1112729817
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
]
12
Oak
[
17.18666431
0.8648974218
3.239438239
2.380000000
1.532000000
0.7840000000
]
[
0.06975010637
0.9079000793
-
0.4133429180
0.0
0.0
0.0
0.1071019008
-
0.4187725641
-
0.9017531389
0.0
0.0
0.0
0.9917984199
-
0.01862756541
0.1264472547
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
]
13
Ash
[
17.15567967
0.6696042476
2.409716064
2.684000000
1.784000000
0.5400000000
]
[
0.06190313535
0.8746549514
-
0.4807772027
0.0
0.0
0.0
0.09781749768
-
0.4846986500
-
0.8691944279
0.0
0.0
0.0
0.9932772717
-
0.006777439445
0.1155609239
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
]
14
Ash
[
17.42363268
0.7844815431
2.669885744
1.720000000
1.218000000
0.5000000000
]
[
-
0.1004077048
-
0.8554189669
-
0.5081108976
0.0
0.0
0.0
-
0.1363859604
0.5177044741
-
0.8446188165
0.0
0.0
0.0
-
0.9855542417
0.01550704374
0.1686486565
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
]
15
Beech
[
15.99002755
0.9558575345
3.451114905
3.216000000
2.112000000
0.9120000000
]
[
0.1135176976
0.8848520771
-
0.4518302069
0.0
0.0
0.0
0.1764907831
-
0.4654963442
-
0.8672739791
0.0
0.0
0.0
0.9777344911
-
0.01870707734
0.2090103088
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
]Table 10
The eigenvalues and eigenvectors for 8 softwoods species calculated using MAPLE.Λ
[
i
] stand for the eigenvalues of ith softwood species, and  V
[
i
] stand for the corresponding eigenvectors.
S. no.
Softwood species
EigenvaluesΛ
[
i
]
EigenvectorsV
[
i
]
1
Blasa
[
6.385324208
0.09845837744
0.3832174064
0.6240000000
0.4060000000
0.06600000000
]
[
0.01488818113
0.9510211778
-
0.3087669978
0.0
0.0
0.0
0.02575997614
-
0.3090635474
-
0.9506924583
0.0
0.0
0.0
0.9995572851
-
0.006200252135
0.02909967665
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
]
2
Spruce
[
11.98583568
0.4516442734
1.114520065
1.498000000
1.442000000
0.07800000000
]
[
-
0.03300264531
-
0.9144925828
-
0.4032544375
0.0
0.0
0.0
-
0.04689865645
0.4044467539
-
0.9133582749
0.0
0.0
0.0
-
0.9983543167
0.01123114838
0.05623628701
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
]
3
Spruce
[
14.10927768
0.4185382987
1.227183961
1.442000000
1.0
0.06400000000
]
[
-
0.03651783317
-
0.8968065074
-
0.4409132956
0.0
0.0
0.0
-
0.05361649666
0.4423303564
-
0.8952480817
0.0
0.0
0.0
-
0.9978936411
0.009052295016
0.06423657043
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
]
4
Spruce
[
16.30529953
0.3544288592
0.8442715844
1.234000000
1.520000000
0.07200000000
]
[
-
0.02057139660
-
0.9124371483
-
0.4086994866
0.0
0.0
0.0
-
0.02869707087
0.4091564350
-
0.9120128765
0.0
0.0
0.0
-
0.9993764538
0.007032901380
0.03460119282
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
]
5
Spruce
[
17.25845208
0.5092652753
1.171282625
1.250000000
1.706000000
0.07000000000
]
[
-
0.03391880763
-
0.8104111857
-
0.5848788055
0.0
0.0
0.0
-
0.03428302033
0.5858146064
-
0.8097196604
0.0
0.0
0.0
-
0.9988364173
0.007413313465
0.04765347716
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
]
6
Douglas Fir
[
16.13459769
0.6233733771
1.439028943
1.767000000
1.766000000
0.1760000000
]
[
-
0.03639420861
-
0.8057068065
-
0.5911954018
0.0
0.0
0.0
-
0.03697194691
0.5922678885
-
0.8048924305
0.0
0.0
0.0
-
0.9986533619
0.007435777607
0.05134366998
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
]
7
Douglas Fir
[
17.10150156
0.6986859431
2.204812526
2.348000000
1.816000000
0.1600000000
]
[
-
0.04991838192
-
0.8212998538
-
0.5683086323
0.0
0.0
0.0
-
0.06364825251
0.5704773676
-
0.8188433771
0.0
0.0
0.0
-
0.9967231590
0.004703484281
0.08075160717
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
]
8
Pine
[
16.99539133
0.4940019684
1.565606760
3.484000000
1.344000000
0.1320000000
]
[
-
0.03436105770
-
0.8973128643
-
0.4400556096
0.0
0.0
0.0
-
0.05585205149
0.4413516031
-
0.8955943895
0.0
0.0
0.0
-
0.9978476167
0.006195562437
0.06528207193
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
0.0
0.0
0.0
0.0
0.0
0.0
1.0
]Table 11
The eigenvalues and eigenvectors for 3 specimens of cancellous bone calculated using MAPLE.Λ
[
i
] stand for the eigenvalues of ith specimen of cancellous bone, and V
[
i
] stand for the corresponding eigenvectors.
Cancellous bone
EigenvaluesΛ
[
i
]
EigenvectorsV
[
i
]
Specimen 1
[
13.34869880
-
0.2071167473
8.644210768
6.837676373
3.768806080
4.123724855
]
[
0.3316115964
0.8831238953
0.2356441950
0.2242668491
0.05354786914
0.03787816350
0.6475092156
-
0.1099136965
0.1556396475
-
0.6940063479
-
0.1278953723
0.2154647417
0.4800906963
-
0.2832250652
0.07648827487
0.2494278793
0.6410471445
-
0.4585742261
-
0.2737532692
-
0.06354940263
0.4518342532
-
0.06214920576
0.5836952694
0.6101669758
-
0.3706582070
0.3314035780
-
0.1276606259
-
0.6337030234
0.3639432801
-
0.4499474345
0.1671829224
0.1180163925
-
0.8330343502
0.02008124471
0.3109325968
0.4087706408
]
Specimen 2
[
1754.634171
777.3848354
155.0688791
230.1790452
344.5174379
358.9156342
]
[
-
0.3057290719
-
0.3005337771
-
0.9030341401
0.01584500766
-
0.02175914011
-
0.003741933916
-
0.4311831740
-
0.8021570593
0.4130146809
0.0001362645283
-
0.003854432451
0.005396441669
-
0.8488209390
0.5154657137
0.1152135966
-
0.004741957786
0.02229739182
-
0.002068264586
0.0009402467441
-
0.005358584867
0.003987518000
-
0.03989381331
0.03796849613
-
0.9984595017
-
0.01058157817
0.02100470582
0.02092543658
0.006230946613
-
0.9987520802
-
0.03826770492
0.0009823402498
0.006977574167
0.01484146906
0.9990475909
0.008196718814
-
0.03958286493
]
Specimen 3
[
1418.542670
583.8388363
99.59058231
343.3818493
302.4968125
175.7492470
]
[
0.3244960223
-
0.004831643783
-
0.8451713783
0.4062092722
0.1164524794
-
0.04239310198
0.6461152835
-
0.7125334236
0.2668138654
0.04131783689
-
0.02431040234
0.03665187401
0.6621315621
0.7009969173
0.2633438279
0.02836289512
-
0.001711614840
0.005269111592
-
0.1962401199
0.009159859564
0.3740839141
0.8850746598
0.1904577552
0.04284654369
-
0.008597323304
-
0.002526763452
0.01897677365
0.1738799329
-
0.6977939665
-
0.6945566457
0.01533190599
-
0.02803230310
0.06961922985
-
0.1374446216
0.6801869735
-
0.7159517722
]


## 5.2. Computing the Average Elasticity Tensors for Woods and Cancellous Bone Using MAPLE

In order to compute average elasticity tensors for the hardwoods, softwoods, and cancellous bone, let us go through a brief mathematical description regarding this issue [3], as this notion helps us developing an appropriate MAPLE code to have desired results.Suppose, we haveM measurements of the elastic constants data for some particular species or material, C
^
I
,
C
^
I
I
,
…
,
C
^
M and we want to construct a tensor C
AVG representing an average elasticity tensor for that particular material. Then to do so, we need the following key algebra [3].(i)
The eigenvaluesΛ
k
Y, k
=
1,2
,
…
,
6, Y
=
1,2
,
…
,
M and the eigenvectors N
^
k
Y, k
=
1,2
,
…
,
6, Y
=
1,2
,
…
,
M for each of the Yth measurement of the particular species or material.
(ii)
The nominal average (NA),N
^
k
NA of the eigenvectors N
^
k
Y associated with the particular species. The nominal average is given by
(45)
N
^
k
NA
≡
1
M
∑
Y
=
1
M
‍
N
^
k
Y
.
(iii)
The averageN
^
k
AVG, which is obtained by the following formula:
(46)
N
^
k
AVG
≡
J
^
NA
N
^
k
NA
,
whereJ
^
NA stands for the inverse square root of the positive definite tensor |
∑
k
=
1
6
‍
N
^
k
NA
⊗
N
^
k
NA
|; that is,
(47)
J
^
NA
≡
(
∑
k
=
1
6
‍
N
^
k
NA
⊗
N
^
k
NA
)
-
1
/
2
.
(iv)
Now, we proceed to average the eigenvalue. For this, we need to transform each set of eigenvalues corresponding to each eigenvector to the average eigenvectorN
^
k
AVG; that is,
(48)
Λ
q
AVG
≡
1
M
∑
Y
=
1
M
‍
∑
k
=
1
6
‍
Λ
k
Y
(
N
^
k
Y
·
N
^
q
AVG
)
2
.
It is obvious thatΛ
k
AVG
>
0 for all k
=
1,2
,
…
,
6 and Λ
k
AVG
≠
Λ
k
NA, where Λ
k
NA is the nominal average of eigenvalues of material and is defined by
(49)
Λ
k
NA
≡
1
M
∑
Y
=
1
M
‍
Λ
k
Y
.
(v)
Eventually, we can now calculate the average elasticity matrixC
^
AVG in such a way that
(50)
C
^
AVG
=
∑
k
=
1
6
‍
Λ
k
AVG
N
^
k
AVG
⊗
N
^
k
AVG
.Taking care of the above outlined steps, we have developed some MAPLE codes (See Step 17 to Step 81 of Appendix A) that enable us to calculate the nominal averages, average eigenvectors, average eigenvalues, and the average elasticity matrices for each of the 15 hardwood species tabulated in Table 5. In addition to this, Appendix A also retains few steps, for example, Steps 21, 26, 31, 36, 41, 46, 51, 56, 61, 66, 73, and 80, which are particularly developed to plot histograms of average elasticity matrices of each of the hardwood species as well as graphs between I, II, III, IV, V, and VI eigenvalues of all hardwood species and the apparent densities (see Table 5).We summarize the above claimed consequences of each of the hardwood species one by one ss follows.


### 5.2.1. MAPLE Evaluations for Quipo

(
1
)
The nominal averages of eigenvectors for the two measurements (S. nos. 1 and 2 of Table5) of Quipo are
(51)
[
0.0367834559700000
0.0453681054800000
0.998185540700000
0.0
0.0
0.0
]
,
[
0.983745905000000
-
0.170879918750000
-
0.0277901141300000
0.0
0.0
0.0
]
,
[
-
0.169284222800000
-
0.982986098000000
0.0509905132200000
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
2
)
The average eigenvectors for Quipo are(52)
[
0.03679134179
0.04537783172
0.9983995367
0.0
0.0
0.0
]
,
[
0.9859858256
-
0.1712690003
-
0.02785339027
0.0
0.0
0.0
]
,
[
-
0.1697052873
-
0.9854311019
0.05111734310
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
3
)
The averages eigenvalue for the two measurements (S. nos. 1 and 2 of Table5) of Quipo is 3.339500000.
(
4
)
The average elasticity matrix for Quipo is(53)
[
3.34725276837330
0.000112116752981933
0.00198542359441849
0.0
0.0
0.0
0.000112116752981933
3.34773746006714
-
0.000991802322788501
0.0
0.0
0.0
0.00198542359441849
-
0.000991802322788501
3.34013593769726
0.0
0.0
0.0
0.0
0.0
0.0
3.33950000000000
0.0
0.0
0.0
0.0
0.0
0.0
3.33950000000000
0.0
0.0
0.0
0.0
0.0
0.0
3.33950000000000
]
.(
5
)
The histogram of the average elasticity matrix for Quipo is shown in Figure2.Figure 2
Histogram showing average elasticity matrix of Quipo.


### 5.2.2. MAPLE Evaluations for White

(
1
)
The nominal averages of eigenvectors for the single measurements (S. no. 3 of Table5) of White are(54)
[
0.04028666644
0.06399313350
0.9971368323
0.0
0.0
0.0
]
,
[
0.9048796003
-
0.4255671399
-
0.009247686096
0.0
0.0
0.0
]
,
[
-
0.4237568827
-
0.9026613361
0.07505076253
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
2
)
The average eigenvectors for White are(55)
[
0.04028666644
0.06399313350
0.9971368323
0.0
0.0
0.0
]
,
[
0.9048795994
-
0.4255671395
-
0.009247686087
0.0
0.0
0.0
]
,
[
-
0.4237568827
-
0.9026613361
0.07505076253
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
3
)
The average eigenvalue for the single measurement (S. no. 3 of Table5) of White is 14.58799998.
(
4
)
The average elasticity matrix for White is given as(56)
[
14.58799998
0.00000001785571215
-
0.00000001009489597
0.0
0.0
0.0
0.00000001785571215
14.58799997
0.000000002407020197
0.0
0.0
0.0
-
0.00000001009489597
0.000000002407020197
14.58799997
0.0
0.0
0.0
0.0
0.0
0.0
14.58799998
0.0
0.0
0.0
0.0
0.0
0.0
14.58799998
0.0
0.0
0.0
0.0
0.0
0.0
14.58799998
]
.(
5
)
The histogram of the average elasticity matrix for White is shown in Figure3.Figure 3
Histogram showing average elasticity matrix of white.


### 5.2.3. MAPLE Evaluations for Khaya

(
1
)
The nominal averages of eigenvectors for the single measurements (S. no. 4 of Table5) of Khaya are(57)
[
0.05368516527
0.07220768570
0.9959437502
0.0
0.0
0.0
]
,
[
0.9266208315
-
0.3753089731
-
0.02273783078
0.0
0.0
0.0
]
,
[
-
0.3721447810
-
0.9240829102
0.08705767064
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
.
(
2
)
The average eigenvectors for Khaya are(58)
[
0.05368516527
0.07220768570
0.9959437502
0.0
0.0
0.0
]
,
[
0.9266208315
-
0.3753089731
-
0.02273783078
0.0
0.0
0.0
]
,
[
-
0.3721447806
-
0.9240829093
0.08705767055
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
3
)
The average eigenvalue for the single measurement (S. no. 4 of Table5) of Khaya is 16.15299998.
(
4
)
The average elasticity matrix for Khaya is given as(59)
[
16.15299998
0.000000005508173090
-
0.000000008722620038
0.0
0.0
0.0
0.000000005508173090
16.15299996
-
0.000000004345156817
0.0
0.0
0.0
-
0.000000008722620038
-
0.000000004345156817
16.15300000
0.0
0.0
0.0
0.0
0.0
0.0
16.15299998
0.0
0.0
0.0
0.0
0.0
0.0
16.15299998
0.0
0.0
0.0
0.0
0.0
0.0
16.15299998
]
.(
5
)
The histogram of the average elasticity matrix for Khaya is shown in Figure4.Figure 4
Histogram showing average elasticity matrix of Khaya.


### 5.2.4. MAPLE Evaluations for Mahogany

(
1
)
The nominal averages of eigenvectors for the two measurements (S. nos. 5 and 6 of Table5) of Mahogany are(60)
[
-
0.0598757013800000
-
0.0768229223150000
-
0.995231841750000
0.0
0.0
0.0
]
,
[
-
0.869358919900000
0.493939153600000
0.0141567750950000
0.0
0.0
0.0
]
,
[
-
0.490490276500000
-
0.866078136900000
0.0963811082600000
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.(
2
)
The average eigenvectors for Mahogany are(61)
[
-
0.05987730132
-
0.07682497510
-
0.9952584354
0.0
0.0
0.0
]
,
[
-
0.8693926232
0.4939583026
0.01415732393
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
3
)
The average eigenvalue for the two measurements (S. nos. 5 and 6 of Table5) of Mahogany is 18.19400002.
(
4
)
The average elasticity matrix for Mahogany is given as(62)
[
18.19451173
-
0.0001438038661
0.0001358556898
0.0
0.0
0.0
-
0.0001438038661
18.19484018
-
0.0004764690574
0.0
0.0
0.0
0.0001358556898
-
0.0004764690574
18.19454255
0.0
0.0
0.0
0.0
0.0
0.0
18.19400002
18.19400002
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
18.19400002
]
.(
5
)
The histogram of the average elasticity matrix for Mahogany is shown in Figure5.Figure 5
Histogram showing average elasticity matrix of Mahogany.


### 5.2.5. MAPLE Evaluations for S. Germ

(
1
)
The nominal averages of eigenvectors for the single measurement (S. no. 7 of Table5) of S. Gem are
(63)
[
0.04967528691
0.08463937788
0.9951726186
0.0
0.0
0.0
]
,
[
0.9161715971
-
0.4006166011
-
0.01165943224
0.0
0.0
0.0
]
,
[
-
0.3976958243
-
0.9123280736
0.09744493938
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.(
2
)
The average eigenvectors for S. Germ are(64)
[
0.04967528696
0.08463937796
0.9951726196
0.0
0.0
0.0
]
,
[
0.9161715980
-
0.4006166015
-
0.01165943225
0.0
0.0
0.0
]
,
[
-
0.3976958247
-
0.9123280745
0.09744493948
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
3
)
The average eigenvalue for the single measurement (S. no. 7 of Table5) of S. Germ is 19.22399998.
(
4
)
The average elasticity matrix for S. Germ is given as(65)
[
19.22399998
-
0.00000001176508811
-
0.00000001537920006
0.0
0.0
0.0
-
0.00000001176508811
19.22400000
-
0.000000007631928036
0.0
0.0
0.0
-
0.00000001537920006
-
0.000000007631928036
19.22400000
0.0
0.0
0.0
0.0
0.0
0.0
19.22399998
19.22399998
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
19.22399998
]
.(
5
)
The histogram of the average elasticity matrix for S. Germ is shown in Figure6.Figure 6
Histogram showing average elasticity matrix of S. Germ.


### 5.2.6. MAPLE Evaluations for maple

(
1
)
The nominal averages of eigenvectors for the single measurement (S. no. 8 of Table5) of maple are
(66)
[
0.0
0.0
0.0
0.0
0.0
1.0
]
,
[
-
0.1387533797
-
0.2031928413
-
0.9692575344
0.0
0.0
0.0
]
,
[
-
0.8476774112
0.5304148448
0.01015375725
0.0
0.0
0.0
]
,
[
-
0.5120454135
-
0.8230265872
0.2458388325
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
.(
2
)
The average eigenvectors for maple are(67)
[
0.0
0.0
0.0
0.0
0.0
1.0
]
,
[
-
0.1387533798
-
0.2031928415
-
0.9692575354
0.0
0.0
0.0
]
,
[
-
0.8476774112
0.5304148448
0.01015375725
0.0
0.0
0.0
]
,
[
-
0.5120454135
-
0.8230265872
0.2458388325
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
.
(
3
)
The average eigenvalue for the single measurement (S. no. 8 of Table5) of maple is 20.74600000.
(
4
)
The average elasticity matrix for maple is given as(68)
[
20.74599999
0.000000004356660361
0.00000003402343994
0.0
0.0
0.0
0.000000004356660361
20.74600004
0.00000002232269567
0.0
0.0
0.0
0.00000003402343994
0.00000002232269567
20.74600000
0.0
0.0
0.0
0.0
0.0
0.0
20.74600000
20.74600000
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
20.74600000
]
.(
5
)
The histogram of the average elasticity matrix for MAPLE is shown in Figure7.Figure 7
Histogram showing average elasticity matrix of maple.


### 5.2.7. MAPLE Evaluations for Walnut

(
1
)
The nominal averages of eigenvectors for the single measurement (S. no. 9 of Table5) of Walnut are
(69)
[
0.08621257414
0.1243595697
0.9884847438
0.0
0.0
0.0
]
,
[
0.8755641188
-
0.4828494241
-
0.01561753067
0.0
0.0
0.0
]
,
[
-
0.4753471023
-
0.8668282006
0.1505124700
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.(
2
)
The average eigenvectors for Walnut are(70)
[
0.08621257423
0.1243595698
0.9884847448
0.0
0.0
0.0
]
,
[
0.8755641188
-
0.4828494241
-
0.01561753067
0.0
0.0
0.0
]
,
[
-
0.4753471018
-
0.8668281997
0.1505124698
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
3
)
The average eigenvalue for the single measurement (S. no. 9 of Table5) of walnut is 18.90099997.
(
4
)
The average elasticity matrix for Walnut is given as(71)
[
18.90099999
0.00000001209663993
-
0.00000002532733996
0.0
0.0
0.0
0.00000001209663993
18.90099991
0.000000001701090138
0.0
0.0
0.0
-
0.00000002532733996
0.000000001701090138
18.90100001
0.0
0.0
0.0
0.0
0.0
0.0
18.90099997
18.90099997
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
18.90099997
]
.(
5
)
The histogram of the average elasticity matrix for Walnut is shown in Figure8.Figure 8
Histogram showing average elasticity matrix of Walnut.


### 5.2.8. MAPLE Evaluations for Birch

(
1
)
The nominal averages of eigenvectors for the single measurement (S. no. 10 of Table5) of Birch are
(72)
[
0.3961520086
0.6338768023
0.6642768888
0.0
0.0
0.0
]
,
[
0.9160614623
-
0.2236797319
-
0.3328645006
0.0
0.0
0.0
]
,
[
0.06240981414
-
0.7403833993
0.6692812840
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
2
)
The average eigenvectors for Birch are(73)
[
0.3961520086
0.6338768023
0.6642768888
0.0
0.0
0.0
]
,
[
0.9160614614
-
0.2236797317
-
0.3328645003
0.0
0.0
0.0
]
,
[
0.06240981426
-
0.7403834008
0.6692812853
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
3
)
The average eigenvalue for the single measurement (S. no. 10 of Table5) of Birch is 8.772299998.
(
4
)
The average elasticity matrix for Birch is given as
(74)
[
8.772299997
-
0.00000003535236899
0.00000003421196978
0.0
0.0
0.0
-
0.00000003535236899
8.772300024
-
0.00000001649192439
0.0
0.0
0.0
0.00000003421196978
-
0.00000001649192439
8.772299994
0.0
0.0
0.0
0.0
0.0
0.0
8.772299998
8.772299998
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
8.772299998
]
.(
5
)
The histogram of the average elasticity matrix for Birch is shown in Figure9.Figure 9
Histogram showing average elasticity matrix of Birch.


### 5.2.9. MAPLE Evaluations for Y. Birch

(
1
)
The nominal averages of eigenvectors for the single measurement (S. no. 11 of Table5) of Y. Birch are
(75)
[
-
0.06591789198
-
0.08975775227
-
0.9937798435
0.0
0.0
0.0
]
,
[
-
0.8289381135
0.5593224991
0.004466103673
0.0
0.0
0.0
]
,
[
-
0.5554425577
-
0.8240763851
0.1112729817
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
2
)
The average eigenvectors for Y. Birch are(76)
[
-
0.1387533798
-
0.2031928415
-
0.9692575354
0.0
0.0
0.0
]
,
[
-
0.8476774112
0.5304148448
0.01015375725
0.0
0.0
0.0
]
,
[
-
0.5120454135
-
0.8230265872
0.2458388325
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
3
)
The average eigenvalue for the single measurement (S. no. 11 of Table5) of Y. Birch is 22.28057495.
(
4
)
The average elasticity matrix for Y. Birch is given as(77)
[
22.28057494
0.000000004678921127
0.00000003654014285
0.0
0.0
0.0
0.000000004678921127
22.28057499
0.00000002397389829
0.0
0.0
0.0
0.00000003654014285
0.00000002397389829
22.28057495
0.0
0.0
0.0
0.0
0.0
0.0
22.28057495
22.28057495
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
22.28057495
]
.(
5
)
The histogram of the average elasticity matrix for Y. Birch is shown in Figure10.Figure 10
Histogram showing average elasticity matrix of Y. Birch.


### 5.2.10. MAPLE Evaluations for Oak

(
1
)
The nominal averages of eigenvectors for the single measurement (S. no. 12 of Table5) of Oak are(78)
[
0.06975010637
0.1071019008
0.9917984199
0.0
0.0
0.0
]
,
[
0.9079000793
-
0.4187725641
-
0.01862756541
0.0
0.0
0.0
]
,
[
-
0.4133429180
-
0.9017531389
0.1264472547
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.(
2
)
The average eigenvectors for Oak are(79)
[
0.06975010637
0.1071019008
0.9917984199
0.0
0.0
0.0
]
,
[
0.9079000793
-
0.4187725641
-
0.01862756541
0.0
0.0
0.0
]
,
[
-
0.4133429184
-
0.9017531398
0.1264472548
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
3
)
The average eigenvalue for the single measurement (S. no. 12 of Table5) of Oak is 25.98699998.
(
4
)
The average elasticity matrix for Oak is given as(80)
[
25.98699997
-
0.00000001889254939
-
0.000000003638179937
0.0
0.0
0.0
-
0.00000001889254939
25.98700006
0.000000008575709980
0.0
0.0
0.0
-
0.000000003638179937
0.000000008575709980
25.98699998
0.0
0.0
0.0
0.0
0.0
0.0
25.98699998
25.98699998
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
25.98699998
]
.(
5
)
The histogram of the average elasticity matrix for Oak is shown in Figure11.Figure 11
Histogram showing average elasticity matrix of Oak.


### 5.2.11. MAPLE Evaluations for Ash

(
1
)
The nominal averages of eigenvectors for the two measurements (S. nos. 13 and 14 of Table5) of Ash are
(81)
[
-
0.0192522847250000
-
0.0192842313600000
0.00386151499999998
0.0
0.0
0.0
]
,
[
0.00961799224999998
0.0165029120500000
0.00436480214750000
0.0
0.0
0.0
]
,
[
-
0.494444050150000
-
0.856906622200000
0.142104790200000
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
2
)
The average eigenvectors for Ash are(82)
[
-
25.41745843
-
25.45963537
5.098090872
0.0
0.0
0.0
]
,
[
25.05315863
42.98714977
11.36953303
0.0
0.0
0.0
]
,
[
-
0.4949599718
-
0.8578007511
0.1422530677
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
3
)
The average eigenvalue for the two measurements (S. nos. 13 and 14 of Table5) of Ash is 24.77950014.
(
4
)
The average elasticity matrix for Ash is given as(83)
[
31567.9169478799
42732.4383657779
3845.57503470365
0.0
0.0
0.0
42732.4383657779
61870.0481877479
8891.53535276175
0.0
0.0
0.0
3845.57503470365
8891.53535276175
3847.68762708609
0.0
0.0
0.0
0.0
0.0
0.0
24.7795001400000
0.0
0.0
0.0
0.0
0.0
0.0
24.7795001400000
0.0
0.0
0.0
0.0
0.0
0.0
24.7795001400000
]
.(
5
)
The histogram of the average elasticity matrix for Ash is shown in Figure12.Figure 12
Histogram showing average elasticity matrix of Ash.


### 5.2.12. MAPLE Evaluations for Beech

(
1
)
The nominal averages of eigenvectors for the single measurements (S. no. 15 of Table5) of Beech are
(84)
[
0.1135176976
0.1764907831
0.9777344911
0.0
0.0
0.0
]
,
[
0.8848520771
-
0.4654963442
-
0.01870707734
0.0
0.0
0.0
]
,
[
-
0.4518302069
-
0.8672739791
0.2090103088
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
2
)
The average eigenvectors for Beech are(85)
[
0.1135176977
0.1764907833
0.9777344921
0.0
0.0
0.0
]
,
[
0.8848520771
-
0.4654963442
-
0.01870707734
0.0
0.0
0.0
]
,
[
-
0.4518302069
-
0.8672739791
0.2090103088
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
3
)
The average eigenvalue for the single measurements (S. no. 15 of Table5) of Beech is 26.63699998.
(
4
)
The average elasticity matrix for Beech is given as(86)
[
26.63700003
0.00000004768023095
0.00000003169803038
0.0
0.0
0.0
0.00000004768023095
26.63699992
0.000000008310743498
0.0
0.0
0.0
0.00000003169803038
0.000000008310743498
26.63700001
0.0
0.0
0.0
0.0
0.0
0.0
26.63699998
26.63699998
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
26.63699998
]
.(
5
)
The histogram of the average elasticity matrix for Beech is shown in Figure13.Figure 13
Histogram showing average elasticity matrix of Beech.The graphics placed in left as well as right positions, showing graphs of I and II eigenvalues of all 15 hardwood species against their apparent densities.
(a)
(b)The graphics placed in left as well as right positions, showing graphs of III and IV eigenvalues of all 15 hardwood species against their apparent densities.
(a)
(b)The graphics placed in left as well as right positions, showing graphs of V and VI eigenvalues of all 15 hardwood species against their apparent densities.
(a)
(b)Figure 17
A combined view of Figures14, 15, and 16.


## 5.2.1. MAPLE Evaluations for Quipo

(
1
)
The nominal averages of eigenvectors for the two measurements (S. nos. 1 and 2 of Table5) of Quipo are
(51)
[
0.0367834559700000
0.0453681054800000
0.998185540700000
0.0
0.0
0.0
]
,
[
0.983745905000000
-
0.170879918750000
-
0.0277901141300000
0.0
0.0
0.0
]
,
[
-
0.169284222800000
-
0.982986098000000
0.0509905132200000
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
2
)
The average eigenvectors for Quipo are(52)
[
0.03679134179
0.04537783172
0.9983995367
0.0
0.0
0.0
]
,
[
0.9859858256
-
0.1712690003
-
0.02785339027
0.0
0.0
0.0
]
,
[
-
0.1697052873
-
0.9854311019
0.05111734310
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
3
)
The averages eigenvalue for the two measurements (S. nos. 1 and 2 of Table5) of Quipo is 3.339500000.
(
4
)
The average elasticity matrix for Quipo is(53)
[
3.34725276837330
0.000112116752981933
0.00198542359441849
0.0
0.0
0.0
0.000112116752981933
3.34773746006714
-
0.000991802322788501
0.0
0.0
0.0
0.00198542359441849
-
0.000991802322788501
3.34013593769726
0.0
0.0
0.0
0.0
0.0
0.0
3.33950000000000
0.0
0.0
0.0
0.0
0.0
0.0
3.33950000000000
0.0
0.0
0.0
0.0
0.0
0.0
3.33950000000000
]
.(
5
)
The histogram of the average elasticity matrix for Quipo is shown in Figure2.Figure 2
Histogram showing average elasticity matrix of Quipo.


## 5.2.2. MAPLE Evaluations for White

(
1
)
The nominal averages of eigenvectors for the single measurements (S. no. 3 of Table5) of White are(54)
[
0.04028666644
0.06399313350
0.9971368323
0.0
0.0
0.0
]
,
[
0.9048796003
-
0.4255671399
-
0.009247686096
0.0
0.0
0.0
]
,
[
-
0.4237568827
-
0.9026613361
0.07505076253
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
2
)
The average eigenvectors for White are(55)
[
0.04028666644
0.06399313350
0.9971368323
0.0
0.0
0.0
]
,
[
0.9048795994
-
0.4255671395
-
0.009247686087
0.0
0.0
0.0
]
,
[
-
0.4237568827
-
0.9026613361
0.07505076253
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
3
)
The average eigenvalue for the single measurement (S. no. 3 of Table5) of White is 14.58799998.
(
4
)
The average elasticity matrix for White is given as(56)
[
14.58799998
0.00000001785571215
-
0.00000001009489597
0.0
0.0
0.0
0.00000001785571215
14.58799997
0.000000002407020197
0.0
0.0
0.0
-
0.00000001009489597
0.000000002407020197
14.58799997
0.0
0.0
0.0
0.0
0.0
0.0
14.58799998
0.0
0.0
0.0
0.0
0.0
0.0
14.58799998
0.0
0.0
0.0
0.0
0.0
0.0
14.58799998
]
.(
5
)
The histogram of the average elasticity matrix for White is shown in Figure3.Figure 3
Histogram showing average elasticity matrix of white.


## 5.2.3. MAPLE Evaluations for Khaya

(
1
)
The nominal averages of eigenvectors for the single measurements (S. no. 4 of Table5) of Khaya are(57)
[
0.05368516527
0.07220768570
0.9959437502
0.0
0.0
0.0
]
,
[
0.9266208315
-
0.3753089731
-
0.02273783078
0.0
0.0
0.0
]
,
[
-
0.3721447810
-
0.9240829102
0.08705767064
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
.
(
2
)
The average eigenvectors for Khaya are(58)
[
0.05368516527
0.07220768570
0.9959437502
0.0
0.0
0.0
]
,
[
0.9266208315
-
0.3753089731
-
0.02273783078
0.0
0.0
0.0
]
,
[
-
0.3721447806
-
0.9240829093
0.08705767055
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
3
)
The average eigenvalue for the single measurement (S. no. 4 of Table5) of Khaya is 16.15299998.
(
4
)
The average elasticity matrix for Khaya is given as(59)
[
16.15299998
0.000000005508173090
-
0.000000008722620038
0.0
0.0
0.0
0.000000005508173090
16.15299996
-
0.000000004345156817
0.0
0.0
0.0
-
0.000000008722620038
-
0.000000004345156817
16.15300000
0.0
0.0
0.0
0.0
0.0
0.0
16.15299998
0.0
0.0
0.0
0.0
0.0
0.0
16.15299998
0.0
0.0
0.0
0.0
0.0
0.0
16.15299998
]
.(
5
)
The histogram of the average elasticity matrix for Khaya is shown in Figure4.Figure 4
Histogram showing average elasticity matrix of Khaya.


## 5.2.4. MAPLE Evaluations for Mahogany

(
1
)
The nominal averages of eigenvectors for the two measurements (S. nos. 5 and 6 of Table5) of Mahogany are(60)
[
-
0.0598757013800000
-
0.0768229223150000
-
0.995231841750000
0.0
0.0
0.0
]
,
[
-
0.869358919900000
0.493939153600000
0.0141567750950000
0.0
0.0
0.0
]
,
[
-
0.490490276500000
-
0.866078136900000
0.0963811082600000
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.(
2
)
The average eigenvectors for Mahogany are(61)
[
-
0.05987730132
-
0.07682497510
-
0.9952584354
0.0
0.0
0.0
]
,
[
-
0.8693926232
0.4939583026
0.01415732393
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
3
)
The average eigenvalue for the two measurements (S. nos. 5 and 6 of Table5) of Mahogany is 18.19400002.
(
4
)
The average elasticity matrix for Mahogany is given as(62)
[
18.19451173
-
0.0001438038661
0.0001358556898
0.0
0.0
0.0
-
0.0001438038661
18.19484018
-
0.0004764690574
0.0
0.0
0.0
0.0001358556898
-
0.0004764690574
18.19454255
0.0
0.0
0.0
0.0
0.0
0.0
18.19400002
18.19400002
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
18.19400002
]
.(
5
)
The histogram of the average elasticity matrix for Mahogany is shown in Figure5.Figure 5
Histogram showing average elasticity matrix of Mahogany.


## 5.2.5. MAPLE Evaluations for S. Germ

(
1
)
The nominal averages of eigenvectors for the single measurement (S. no. 7 of Table5) of S. Gem are
(63)
[
0.04967528691
0.08463937788
0.9951726186
0.0
0.0
0.0
]
,
[
0.9161715971
-
0.4006166011
-
0.01165943224
0.0
0.0
0.0
]
,
[
-
0.3976958243
-
0.9123280736
0.09744493938
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.(
2
)
The average eigenvectors for S. Germ are(64)
[
0.04967528696
0.08463937796
0.9951726196
0.0
0.0
0.0
]
,
[
0.9161715980
-
0.4006166015
-
0.01165943225
0.0
0.0
0.0
]
,
[
-
0.3976958247
-
0.9123280745
0.09744493948
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
3
)
The average eigenvalue for the single measurement (S. no. 7 of Table5) of S. Germ is 19.22399998.
(
4
)
The average elasticity matrix for S. Germ is given as(65)
[
19.22399998
-
0.00000001176508811
-
0.00000001537920006
0.0
0.0
0.0
-
0.00000001176508811
19.22400000
-
0.000000007631928036
0.0
0.0
0.0
-
0.00000001537920006
-
0.000000007631928036
19.22400000
0.0
0.0
0.0
0.0
0.0
0.0
19.22399998
19.22399998
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
19.22399998
]
.(
5
)
The histogram of the average elasticity matrix for S. Germ is shown in Figure6.Figure 6
Histogram showing average elasticity matrix of S. Germ.


## 5.2.6. MAPLE Evaluations for maple

(
1
)
The nominal averages of eigenvectors for the single measurement (S. no. 8 of Table5) of maple are
(66)
[
0.0
0.0
0.0
0.0
0.0
1.0
]
,
[
-
0.1387533797
-
0.2031928413
-
0.9692575344
0.0
0.0
0.0
]
,
[
-
0.8476774112
0.5304148448
0.01015375725
0.0
0.0
0.0
]
,
[
-
0.5120454135
-
0.8230265872
0.2458388325
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
.(
2
)
The average eigenvectors for maple are(67)
[
0.0
0.0
0.0
0.0
0.0
1.0
]
,
[
-
0.1387533798
-
0.2031928415
-
0.9692575354
0.0
0.0
0.0
]
,
[
-
0.8476774112
0.5304148448
0.01015375725
0.0
0.0
0.0
]
,
[
-
0.5120454135
-
0.8230265872
0.2458388325
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
.
(
3
)
The average eigenvalue for the single measurement (S. no. 8 of Table5) of maple is 20.74600000.
(
4
)
The average elasticity matrix for maple is given as(68)
[
20.74599999
0.000000004356660361
0.00000003402343994
0.0
0.0
0.0
0.000000004356660361
20.74600004
0.00000002232269567
0.0
0.0
0.0
0.00000003402343994
0.00000002232269567
20.74600000
0.0
0.0
0.0
0.0
0.0
0.0
20.74600000
20.74600000
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
20.74600000
]
.(
5
)
The histogram of the average elasticity matrix for MAPLE is shown in Figure7.Figure 7
Histogram showing average elasticity matrix of maple.


## 5.2.7. MAPLE Evaluations for Walnut

(
1
)
The nominal averages of eigenvectors for the single measurement (S. no. 9 of Table5) of Walnut are
(69)
[
0.08621257414
0.1243595697
0.9884847438
0.0
0.0
0.0
]
,
[
0.8755641188
-
0.4828494241
-
0.01561753067
0.0
0.0
0.0
]
,
[
-
0.4753471023
-
0.8668282006
0.1505124700
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.(
2
)
The average eigenvectors for Walnut are(70)
[
0.08621257423
0.1243595698
0.9884847448
0.0
0.0
0.0
]
,
[
0.8755641188
-
0.4828494241
-
0.01561753067
0.0
0.0
0.0
]
,
[
-
0.4753471018
-
0.8668281997
0.1505124698
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
3
)
The average eigenvalue for the single measurement (S. no. 9 of Table5) of walnut is 18.90099997.
(
4
)
The average elasticity matrix for Walnut is given as(71)
[
18.90099999
0.00000001209663993
-
0.00000002532733996
0.0
0.0
0.0
0.00000001209663993
18.90099991
0.000000001701090138
0.0
0.0
0.0
-
0.00000002532733996
0.000000001701090138
18.90100001
0.0
0.0
0.0
0.0
0.0
0.0
18.90099997
18.90099997
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
18.90099997
]
.(
5
)
The histogram of the average elasticity matrix for Walnut is shown in Figure8.Figure 8
Histogram showing average elasticity matrix of Walnut.


## 5.2.8. MAPLE Evaluations for Birch

(
1
)
The nominal averages of eigenvectors for the single measurement (S. no. 10 of Table5) of Birch are
(72)
[
0.3961520086
0.6338768023
0.6642768888
0.0
0.0
0.0
]
,
[
0.9160614623
-
0.2236797319
-
0.3328645006
0.0
0.0
0.0
]
,
[
0.06240981414
-
0.7403833993
0.6692812840
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
2
)
The average eigenvectors for Birch are(73)
[
0.3961520086
0.6338768023
0.6642768888
0.0
0.0
0.0
]
,
[
0.9160614614
-
0.2236797317
-
0.3328645003
0.0
0.0
0.0
]
,
[
0.06240981426
-
0.7403834008
0.6692812853
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
3
)
The average eigenvalue for the single measurement (S. no. 10 of Table5) of Birch is 8.772299998.
(
4
)
The average elasticity matrix for Birch is given as
(74)
[
8.772299997
-
0.00000003535236899
0.00000003421196978
0.0
0.0
0.0
-
0.00000003535236899
8.772300024
-
0.00000001649192439
0.0
0.0
0.0
0.00000003421196978
-
0.00000001649192439
8.772299994
0.0
0.0
0.0
0.0
0.0
0.0
8.772299998
8.772299998
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
8.772299998
]
.(
5
)
The histogram of the average elasticity matrix for Birch is shown in Figure9.Figure 9
Histogram showing average elasticity matrix of Birch.


## 5.2.9. MAPLE Evaluations for Y. Birch

(
1
)
The nominal averages of eigenvectors for the single measurement (S. no. 11 of Table5) of Y. Birch are
(75)
[
-
0.06591789198
-
0.08975775227
-
0.9937798435
0.0
0.0
0.0
]
,
[
-
0.8289381135
0.5593224991
0.004466103673
0.0
0.0
0.0
]
,
[
-
0.5554425577
-
0.8240763851
0.1112729817
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
2
)
The average eigenvectors for Y. Birch are(76)
[
-
0.1387533798
-
0.2031928415
-
0.9692575354
0.0
0.0
0.0
]
,
[
-
0.8476774112
0.5304148448
0.01015375725
0.0
0.0
0.0
]
,
[
-
0.5120454135
-
0.8230265872
0.2458388325
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
3
)
The average eigenvalue for the single measurement (S. no. 11 of Table5) of Y. Birch is 22.28057495.
(
4
)
The average elasticity matrix for Y. Birch is given as(77)
[
22.28057494
0.000000004678921127
0.00000003654014285
0.0
0.0
0.0
0.000000004678921127
22.28057499
0.00000002397389829
0.0
0.0
0.0
0.00000003654014285
0.00000002397389829
22.28057495
0.0
0.0
0.0
0.0
0.0
0.0
22.28057495
22.28057495
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
22.28057495
]
.(
5
)
The histogram of the average elasticity matrix for Y. Birch is shown in Figure10.Figure 10
Histogram showing average elasticity matrix of Y. Birch.


## 5.2.10. MAPLE Evaluations for Oak

(
1
)
The nominal averages of eigenvectors for the single measurement (S. no. 12 of Table5) of Oak are(78)
[
0.06975010637
0.1071019008
0.9917984199
0.0
0.0
0.0
]
,
[
0.9079000793
-
0.4187725641
-
0.01862756541
0.0
0.0
0.0
]
,
[
-
0.4133429180
-
0.9017531389
0.1264472547
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.(
2
)
The average eigenvectors for Oak are(79)
[
0.06975010637
0.1071019008
0.9917984199
0.0
0.0
0.0
]
,
[
0.9079000793
-
0.4187725641
-
0.01862756541
0.0
0.0
0.0
]
,
[
-
0.4133429184
-
0.9017531398
0.1264472548
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
3
)
The average eigenvalue for the single measurement (S. no. 12 of Table5) of Oak is 25.98699998.
(
4
)
The average elasticity matrix for Oak is given as(80)
[
25.98699997
-
0.00000001889254939
-
0.000000003638179937
0.0
0.0
0.0
-
0.00000001889254939
25.98700006
0.000000008575709980
0.0
0.0
0.0
-
0.000000003638179937
0.000000008575709980
25.98699998
0.0
0.0
0.0
0.0
0.0
0.0
25.98699998
25.98699998
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
25.98699998
]
.(
5
)
The histogram of the average elasticity matrix for Oak is shown in Figure11.Figure 11
Histogram showing average elasticity matrix of Oak.


## 5.2.11. MAPLE Evaluations for Ash

(
1
)
The nominal averages of eigenvectors for the two measurements (S. nos. 13 and 14 of Table5) of Ash are
(81)
[
-
0.0192522847250000
-
0.0192842313600000
0.00386151499999998
0.0
0.0
0.0
]
,
[
0.00961799224999998
0.0165029120500000
0.00436480214750000
0.0
0.0
0.0
]
,
[
-
0.494444050150000
-
0.856906622200000
0.142104790200000
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
2
)
The average eigenvectors for Ash are(82)
[
-
25.41745843
-
25.45963537
5.098090872
0.0
0.0
0.0
]
,
[
25.05315863
42.98714977
11.36953303
0.0
0.0
0.0
]
,
[
-
0.4949599718
-
0.8578007511
0.1422530677
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
3
)
The average eigenvalue for the two measurements (S. nos. 13 and 14 of Table5) of Ash is 24.77950014.
(
4
)
The average elasticity matrix for Ash is given as(83)
[
31567.9169478799
42732.4383657779
3845.57503470365
0.0
0.0
0.0
42732.4383657779
61870.0481877479
8891.53535276175
0.0
0.0
0.0
3845.57503470365
8891.53535276175
3847.68762708609
0.0
0.0
0.0
0.0
0.0
0.0
24.7795001400000
0.0
0.0
0.0
0.0
0.0
0.0
24.7795001400000
0.0
0.0
0.0
0.0
0.0
0.0
24.7795001400000
]
.(
5
)
The histogram of the average elasticity matrix for Ash is shown in Figure12.Figure 12
Histogram showing average elasticity matrix of Ash.


## 5.2.12. MAPLE Evaluations for Beech

(
1
)
The nominal averages of eigenvectors for the single measurements (S. no. 15 of Table5) of Beech are
(84)
[
0.1135176976
0.1764907831
0.9777344911
0.0
0.0
0.0
]
,
[
0.8848520771
-
0.4654963442
-
0.01870707734
0.0
0.0
0.0
]
,
[
-
0.4518302069
-
0.8672739791
0.2090103088
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
2
)
The average eigenvectors for Beech are(85)
[
0.1135176977
0.1764907833
0.9777344921
0.0
0.0
0.0
]
,
[
0.8848520771
-
0.4654963442
-
0.01870707734
0.0
0.0
0.0
]
,
[
-
0.4518302069
-
0.8672739791
0.2090103088
0.0
0.0
0.0
]
,
[
0.0
0.0
0.0
1.0
0.0
0.0
]
,
[
0.0
0.0
0.0
0.0
1.0
0.0
]
,
[
0.0
0.0
0.0
0.0
0.0
1.0
]
.
(
3
)
The average eigenvalue for the single measurements (S. no. 15 of Table5) of Beech is 26.63699998.
(
4
)
The average elasticity matrix for Beech is given as(86)
[
26.63700003
0.00000004768023095
0.00000003169803038
0.0
0.0
0.0
0.00000004768023095
26.63699992
0.000000008310743498
0.0
0.0
0.0
0.00000003169803038
0.000000008310743498
26.63700001
0.0
0.0
0.0
0.0
0.0
0.0
26.63699998
26.63699998
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
26.63699998
]
.(
5
)
The histogram of the average elasticity matrix for Beech is shown in Figure13.Figure 13
Histogram showing average elasticity matrix of Beech.The graphics placed in left as well as right positions, showing graphs of I and II eigenvalues of all 15 hardwood species against their apparent densities.
(a)
(b)The graphics placed in left as well as right positions, showing graphs of III and IV eigenvalues of all 15 hardwood species against their apparent densities.
(a)
(b)The graphics placed in left as well as right positions, showing graphs of V and VI eigenvalues of all 15 hardwood species against their apparent densities.
(a)
(b)Figure 17
A combined view of Figures14, 15, and 16.


## 5.3. MAPLE Plotted Graphics for I, II, III, IV, V, and VI Eigenvalues of 15 Hardwood Species against Their Apparent Densities

This subsection is dedicated to the exposition of the extreme quality of MAPLE, which can enable one to simultaneously sketch up graphics for I, II, III, IV, V, and VI eigenvalues of each of the 15 hardwood species (See Tables8 and 9) against the corresponding apparent densities ρ (See Table 5 for apparent density). The MAPLE code regarding this issue has been mentioned in Step 81 of Appendix A.


## 6. Conclusions

In this paper, we have tried to develop a CAS environment that suits best to numerically analyze the theory of anisotropic Hooke’s law for different species of wood and cancellous bone. Particularly, the two fabulous CAS, namely, Mathematica and MAPLE, have been used in relation to our proposed research work. More precisely, a Mathematica package “Tensor2Analysis.m” has been employed to rectify the issue regarding unfolding the tensorial form of anisotropic Hooke’s law, whereas a MAPLE code consisting of almost 81 steps has been developed to dig out the following numerical analysis.(i)
Computation of eigenvalues and eigenvectors for all the 15 hardwood species, 8 softwood species, and 3 specimens of cancellous bone using MAPLE.
(ii)
Computation of the nominal averages of eigenvectors and the average eigenvectors for all the 15 hardwood species.
(iii)
Computation of the average eigenvalues for all 15 hardwood species.
(iv)
Computation of the average elasticity matrices for all the 15 hardwood species.
(v)
Plotting the histograms for elasticity matrices of all the 15 hardwood species.
(vi)
Plotting the graphics for I, II, III, IV, V, and VI eigenvalues of 15 hardwood species against their apparent densities.It is also claimed that the developed MAPLE code cannot only simultaneously tackle the 6 × 6 matrices relevant to the different wood species and cancellous bone specimen, but also bears the capacity of handling 100 × 100 matrix whose entries vary with respect to some parameter. Thus, from the viewpoint of author, the MAPLE code developed by him is of great interest in those fields where higher order matrix algebra is required.



---

*Source: 487314-2014-03-23.xml*