# A Result regarding the Seismic Dislocations in Microstretch Thermoelastic Bodies

Authors: M. Marin; O. Florea; S. R. Mahmoud
Journal: Mathematical Problems in Engineering
 (2015)

---

## Abstract

The aim of our study is to derive a relation of De Hoop-Knopoff type for displacement fields within context of thermoelastic microstretch bodies. Then, as a consequence, an explicit expression of the body loadings equivalent to a seismic dislocation is obtained. The results are extensions of those from the classical theory of elastic bodies.


---

## Body

## 1. Introduction

The theory of thermomicrostretch elastic solids was first elaborated by Eringen [1], and, in short, this is a theory of thermoelasticity with microstructure that includes intrinsic rotations and microstructural expansion and contractions.The purpose of this theory is to eliminate discrepancies between classical elasticity and experiments, since the classical elasticity failed to present acceptable results when the effects of material microstructure were known to contribute significantly to the body’s overall deformations, for example, in the case of granular bodies with large molecules (e.g., polymers), graphite, or human bones.These cases are becoming increasingly important in the design and manufacture of modern day advanced materials, as small-scale effects become paramount in the prediction of the overall mechanical behaviour of these materials.Other intended applications of this theory are to composite materials reinforced with chopped fibers and various porous materials.This theory can be useful in the applications which deal with porous materials as geological materials, solid packed granular materials, and many others.On the other hand, materials which operate at elevated temperatures will invariably be subjected to heat flow at some time during normal use. Such heat flow will involve a nonlinear temperature distribution which will inevitably give rise to thermal stresses. For these reasons, the development, design, and selection of materials for high temperature applications require a great deal of care. The role of the pertinent material properties and other variables which can affect the magnitude of thermal stress must be considered.The main difficulty of the thermomicrostretch materials is the large number of the thermoelastic coefficients and, as such, the problem of their determination in the laboratory. Yet many authors consider that this problem will be solved in the future.Already, in the isotropic case, when the number of coefficients decreases a lot, they are calculated as can be seen in many works due to Eringen or Iesan.The present paper must be considered a first step to a better understanding of microstretch and thermal stress in the study of above enumerated materials.The reciprocity and representation relations that appear in our study constitute powerful theoretical tools in the assessment of the theory of seismic-sources mechanism, in the studies connected with seismic wave propagation.Also, we think that this paper is a good help to understand the application of microstretch mechanism to earthquake problems.There are many results regarding the mechanism of earthquake, as in [2–8].For instance, in [6], the authors establish a reciprocity relation which forms the basis for a uniqueness result, a continuous dependence of solutions upon initial data and body loads and a variational characterization of solutions. The effect of a concentrated heat supply and of a concentrated heat volume charge density in an unbounded homogeneous and isotropic electromagnetic body is investigated.We find other results regarding thermoelasticity of nonclassical materials, as in [1, 9–12]. The main results of our study are extensions of some similar results in the classical elasticity in order to cover the thermoelasticity of microstretch bodies (see the results due to Saccomandi in [13]).


## 2. Basic Equations

For convenience the notations and terminology chosen are almost identical to those of our study [11]. Our paper is concerned with an anisotropic and homogeneous material.Let the body occupy, at timet
=
0, a properly regular region B of the three-dimensional Euclidian space, bounded by the piece-wise smooth surface ∂
B and we denote the closure of B by B
-. We refer the motion of the body to a fix system of rectangular Cartesian axes O
x
i, i
=
1,2
,
3 and adopt the Cartesian tensor notation. Points in B are denoted by x
j and t
∈
[
0
,
∞
) is temporal variable. Throughout this work the Einstein summation convention over repeated indices is used. The subscript j after comma indicates partial differentiation with respect to the spatial argument x
j. All Latin subscripts are understood to range over the integers (
1,2
,
3
), while the Greek indices have the range (
1,2
). A superposed dot denotes the derivatives with respect to the t—time variable. Also, the spatial argument and the time argument of a function will be omitted when there is no likelihood of confusion.Let us denote byu
i the components of the displacement vector and by φ
i the components of the microrotation vector. Also, we denote by ω the microstretch function and by θ the temperature measured from the constant absolute temperature T
0 of the body in its reference state.As usual, we denote byt
i
j the components of the stress tensor and by m
i
j the components of the couple stress tensor over B. Also, we denote by λ
i the components of the microstress vector.For clarity and simplicity in presentation, the regularity hypotheses on the considered functions will be ommited.On these grounds, the field equations in the dynamic theory of thermoelasticity of microstretch bodies are as follows (see [1, 11, 14]):(i)
the equations of motion(1)
t
j
i
,
j
+
ϱ
F
i
=
ϱ
u
¨
i
,
m
j
i
,
j
+
ε
i
j
k
t
j
k
+
ϱ
G
i
=
I
i
j
φ
¨
j
;
(ii)
the balance of the equilibrated forces(2)
λ
i
,
i
-
s
+
ϱ
L
=
J
ω
¨
.The equation of energy is given by
(3)
ϱ
T
0
η
˙
=
q
i
,
i
+
ϱ
r
.
In the above equations we have used the following notations:(i)
F
i are the components of body force;
(ii)
G
i are the components of body couple;
(iii)
L is the generalized external body load;
(iv)
s is the generalized internal body load;
(v)
ϱ is the reference constant mass density;
(vi)
J and I
i
j
=
I
j
i are the coefficients of microinertia;
(vii)
η is the entropy per unit mass;
(viii)
r is the heat supply per unit mass;
(ix)
q
i are the components of heat flux vector.For an anisotropic and homogeneous microstretch thermoelastic material, the constitutive equations have the form:(4)
t
i
j
=
A
i
j
r
s
ε
r
s
+
B
i
j
r
s
μ
r
s
+
D
i
j
r
γ
r
+
a
i
j
ω
-
E
i
j
θ
,
m
i
j
=
B
r
s
i
j
ε
r
s
+
C
i
j
r
s
μ
r
s
+
E
i
j
r
γ
r
+
b
i
j
ω
-
D
i
j
θ
,
λ
i
=
D
r
s
i
ε
r
s
+
E
r
s
i
μ
r
s
+
C
i
j
γ
j
+
d
i
ω
-
L
i
θ
,
s
=
a
i
j
ε
i
j
+
b
i
j
μ
i
j
+
d
i
γ
i
+
m
ω
-
α
θ
,
η
=
η
0
+
E
i
j
ε
i
j
+
D
i
j
μ
i
j
+
L
i
γ
i
+
α
ω
+
a
T
0
θ
,
q
i
=
k
i
j
θ
,
j
,
where A
i
j
r
s, B
i
j
r
s, C
i
j
r
s, D
i
j
r, E
i
j
r, a
i
j, b
i
j, c
i
j, d
i, E
i
j, D
i
j, L
i, m, α, a, and k
i
j are the characteristic constitutive coefficients.The components of the strain tensorsε
i
j, μ
i
j, and γ
i are defined by means of the geometric equations:
(5)
ε
i
j
=
u
j
,
i
+
ε
i
j
k
φ
k
,
μ
i
j
=
φ
j
,
i
,
γ
i
=
ω
,
i
,
where ε
i
j
k is the alternating symbol.The constitutive coefficients obey the following symmetry relations:(6)
A
i
j
r
s
=
A
r
s
i
j
,
C
i
j
r
s
=
C
r
s
i
j
,
C
i
j
=
C
j
i
,
k
i
j
=
k
j
i
.
One can assume that a positive constant λ
0 exists such that
(7)
I
i
j
ξ
i
ξ
j
≥
λ
0
ξ
i
ξ
i
,
∀
ξ
i
.
Also, the Second Law of Thermodynamics implies that
(8)
k
i
j
ξ
i
ξ
j
≥
0
,
∀
ξ
i
.
We denote by t
i the components of surface traction, m
i the components of surface couple, p the microsurface traction, and q the heat flux. These quantities are defined by
(9)
t
i
=
t
j
i
n
j
,
m
i
=
m
j
i
n
j
,
p
=
λ
i
n
i
,
q
=
q
i
n
i
,
at regular points of the surface ∂
B.Here,n
i are the components of the outward unit normal of the surface ∂
B.Along with the system of field equations (1)–(5) we consider the following initial conditions:
(10)
u
i
x
,
0
=
u
i
0
x
,
u
˙
i
x
,
0
=
u
i
1
x
,
φ
i
x
,
0
=
φ
i
0
x
,
φ
˙
i
x
,
0
=
φ
i
1
x
,
ω
x
,
0
=
ω
0
x
,
ω
˙
x
,
0
=
ω
1
x
,
θ
x
,
0
=
θ
0
x
,
x
∈
B
-
and the following prescribed boundary conditions:
(11)
u
i
=
u
-
i
on
⁡
∂
B
1
×
0
,
t
0
,
t
i
=
t
-
i
on
⁡
∂
B
1
c
×
0
,
t
0
,
φ
i
=
φ
-
i
on
⁡
∂
B
2
×
0
,
t
0
,
m
i
=
m
-
i
on
⁡
∂
B
2
c
×
0
,
t
0
,
ω
=
ω
-
on
⁡
∂
B
3
×
0
,
t
0
,
p
=
p
-
on
⁡
∂
B
3
c
×
0
,
t
0
,
θ
=
θ
-
on
⁡
∂
B
4
×
0
,
t
0
,
q
=
q
-
on
⁡
∂
B
4
c
×
0
,
t
0
,
where t
0 is some instant that may be infinite.Also,∂
B
1, ∂
B
2, ∂
B
3, and ∂
B
4 with respective complements ∂
B
1
c, ∂
B
2
c, ∂
B
3
c, and ∂
B
4
c are subsets of the surface ∂
B such that
(12)
∂
B
1
∩
∂
B
1
c
=
∂
B
2
∩
∂
B
2
c
=
∂
B
3
∩
∂
B
3
c
=
∂
B
4
∩
∂
B
4
c
=
∅
,
∂
B
1
∪
∂
B
1
c
=
∂
B
2
∪
∂
B
2
c
=
∂
B
3
∪
∂
B
3
c
=
∂
B
4
∪
∂
B
4
c
=
∂
B
.Also,u
i
0, u
i
1, φ
i
0, φ
i
1, ω
0, ω
1, θ
0, u
-
i, t
-
i, φ
-
i, m
-
i, ω
-, p
-, θ
-, and q
- are prescribed functions in their domains.By a solution of the mixed initial boundary value problem of the theory thermoelasticity of dipolar bodies with voids in the cylinderΩ
0
=
B
×
[
0
,
t
0
) we mean an ordered array (u
i, φ
i, ω, and θ) which satisfies (1), (2), and (3) for all (
x
,
t
)
∈
Ω
0, the boundary conditions (11) and the initial conditions (10).


## 3. Main Results

Letu and v be functions defined on B
-
×
[
0
,
∞
) and continuous on [
0
,
∞
) with respect to the time variable t for each spatial variable x
∈
B
-.We denote byu
*
v
the convolution of u and v; that is,
(13)
u
*
v
x
,
t
=
∫
0
t
‍
u
x
,
t
-
τ
v
x
,
τ
d
τ
.
Let us introduce the notations
(14)
g
t
=
t
,
h
t
=
1
,
f
i
=
ϱ
g
*
F
i
+
ϱ
t
u
i
1
x
+
u
i
0
x
,
g
i
=
ϱ
g
*
G
i
+
I
i
j
t
φ
j
1
(
x
)
+
φ
j
0
x
,
l
=
ϱ
g
*
L
+
J
t
ω
1
x
+
ω
0
x
,
w
=
ϱ
h
*
r
+
ϱ
T
0
η
0
.Following the same procedure used by Ieşan in [5], it is easy to prove the following result that enables us to give an alternative formulation of the initial boundary value problem in which the initial data are incorporated into the field of equations.Theorem 1.
The functionsu
i, φ
j
k, σ, θ, τ
i
j, η
i
j, μ
i
j
k, and q
i satisfy (1), (2), (3), and the initial conditions (10) if and only if they satisfy the following system of equations:
(15)
g
*
t
j
i
,
j
+
f
i
=
ϱ
u
i
,
g
*
m
j
i
,
j
+
ε
i
j
k
t
j
k
+
g
i
=
I
i
j
φ
j
,
g
*
λ
i
,
i
-
s
+
l
=
J
ω
,
h
*
q
i
,
i
+
w
=
ϱ
T
0
η
.In our following estimations, we will use formulation (15) of the mixed problem. We wish to find the behavior of the considered medium when embedded in B there is a discontinuity surface Σ for the displacements, the microration vector, the microstretch function, and the temperature. The sides of Σ are denoted by Σ
- and Σ
+.Letν
i be the components of the unit normal vector of Σ, directed from the side (
-
) to the side (
+
).Then on surfaceΣ we have the conditions
(16)
u
i
+
-
u
i
-
=
U
i
,
t
j
i
+
ν
j
=
t
j
i
-
ν
j
,
φ
i
+
-
φ
i
-
=
Φ
i
,
m
i
j
+
ν
j
=
m
i
j
-
ν
j
,
ω
+
-
ω
-
=
Ψ
,
λ
j
+
ν
j
=
λ
j
-
ν
j
,
θ
+
-
θ
-
=
Θ
,
q
j
+
ν
j
=
q
j
-
ν
j
,
where f
+ and f
- are the limits of the function f
(
x
) as x approaches a point on the side (
+
) or (
-
) of the surface Σ, respectively, and U
i, Φ
i, Ψ, and Θ are prescribed functions. In this way we can consider (15) in the domain B
∖
Σ.Let us consider two different systems of loadings for the body(17)
G
α
=
F
i
α
,
G
i
α
,
L
α
,
r
α
,
u
-
i
α
,
φ
-
i
α
,
ω
-
i
α
,
θ
-
α
,
G
α
=
a
a
t
-
i
α
,
m
-
i
α
,
h
-
α
,
q
-
α
,
U
i
,
Φ
i
,
Ψ
,
Θ
,
α
=
1,2
and two corresponding solutions
(18)
S
(
α
)
=
u
i
(
α
)
,
φ
i
(
α
)
,
ω
(
α
)
,
θ
(
α
)
,
ε
i
j
(
α
)
,
μ
i
j
(
α
)
,
t
i
j
(
α
)
,
m
i
j
(
α
)
,
S
α
=
a
a
λ
i
α
,
s
α
,
q
i
α
,
α
=
1,2
.
For the sake of simplicity, we now introduce the notations
(19)
t
i
=
t
i
j
n
j
,
T
i
=
t
i
j
+
ν
j
,
m
i
=
m
i
j
n
j
,
M
i
=
m
i
j
+
ν
j
,
λ
=
λ
i
n
i
,
Λ
=
λ
i
+
ν
i
,
q
=
q
i
n
i
,
Q
=
q
i
+
ν
i
.
In the following theorem, we prove a reciprocity relation of Betti type.Theorem 2.
If a microstretch thermoelastic body is subjected to two systems of loadingsG
(
α
) then between the corresponding solutions S
(
α
) there is the following reciprocity relation:
(20)
∫
B
‍
f
i
(
1
)
*
u
i
(
2
)
+
g
i
(
1
)
*
φ
i
(
2
)
+
l
(
1
)
*
ω
(
2
)
1
T
0
a
a
a
a
-
1
T
0
g
*
w
(
1
)
*
θ
(
2
)
d
V
+
∫
∂
B
‍
g
*
t
i
1
*
u
i
2
+
m
i
1
*
φ
i
2
+
λ
1
*
ω
2
1
T
0
+
∫
∂
B
‍
g
*
aaa
-
1
T
0
h
*
q
1
*
θ
2
d
A
-
∫
Σ
‍
g
*
T
i
1
*
U
i
2
+
M
i
1
*
Φ
i
2
+
Λ
1
*
Ψ
2
1
T
0
+
∫
∂
B
‍
g
*
a
-
1
T
0
h
*
Q
1
*
Θ
2
d
A
=
∫
B
‍
f
i
2
*
u
i
1
+
g
i
2
*
φ
i
1
+
l
2
*
ω
1
1
T
0
+
∫
∂
B
‍
-
1
T
0
g
*
w
2
*
θ
1
d
V
+
∫
∂
B
‍
g
*
t
i
(
2
)
*
u
i
(
1
)
+
m
i
(
2
)
*
φ
i
(
1
)
+
λ
(
2
)
*
ω
(
1
)
1
T
0
+
∫
∂
B
‍
g
*
aaa
-
1
T
0
h
*
q
2
*
θ
1
d
A
-
∫
Σ
‍
g
*
T
i
2
*
U
i
1
+
M
i
2
*
Φ
i
1
+
Λ
2
*
Ψ
1
1
T
0
+
∫
∂
B
‍
g
*
aa
-
1
T
0
h
*
Q
2
*
Θ
1
d
A
.Proof.
In view of symmetry relations (6) and with the aid of the constitutive relations (4), by direct calculations it is easy to obtain
(21)
t
i
j
(
1
)
*
ε
i
j
(
2
)
+
m
i
j
(
1
)
*
μ
i
j
(
2
)
+
λ
i
(
1
)
*
γ
i
(
2
)
+
s
(
1
)
*
ω
(
2
)
-
ϱ
θ
1
*
η
2
=
t
i
j
2
*
ε
i
j
1
+
m
i
j
2
-
ϱ
θ
1
*
η
2
=
aaaaaa
*
μ
i
j
1
+
λ
i
2
*
γ
i
1
+
s
2
-
ϱ
θ
1
*
η
2
=
aaaaaa
*
ω
1
-
ϱ
θ
2
*
η
1
.
Let us introduce the notation
(22)
I
α
β
=
∫
B
‍
g
*
t
i
j
α
*
ε
i
j
β
+
m
i
j
α
*
μ
i
j
β
+
λ
i
α
*
γ
i
β
I
α
β
=
∫
B
‍
g
*
aa
+
s
(
α
)
*
ω
(
β
)
-
ϱ
θ
(
α
)
*
η
(
β
)
d
V
for α
,
β
=
1,2
Based on identity (21) and notation (22), it is easy to see that
(23)
I
α
β
=
I
β
α
.
With the aid of the equations of motion and (15), we can write
(24)
g
*
t
i
j
(
α
)
*
ε
i
j
(
β
)
+
m
i
j
(
α
)
*
μ
i
j
(
β
)
+
λ
i
(
α
)
*
γ
i
(
β
)
aaaa
+
s
α
*
ω
β
-
ϱ
θ
α
*
η
β
=
g
*
t
j
i
α
*
u
j
β
+
m
j
i
α
*
φ
j
β
+
λ
i
α
*
ω
β
1
T
0
=
g
a
a
a
a
-
1
T
0
h
*
q
i
α
*
θ
β
,
i
+
f
i
α
*
u
i
β
+
g
i
α
*
φ
i
β
+
l
α
*
ω
β
-
1
T
0
g
*
w
α
*
θ
β
-
ϱ
u
i
α
*
u
i
β
+
I
i
j
φ
i
α
*
φ
j
β
+
J
ω
α
*
ω
β
+
1
T
0
g
*
h
*
k
i
j
θ
,
i
(
α
)
*
θ
,
j
(
β
)
.
By integrating in (24) and using the divergence theorem, we are lead to
(25)
I
α
β
=
∫
B
‍
f
i
α
*
u
i
β
+
g
i
α
*
φ
i
β
+
l
α
*
ω
β
1
T
0
I
α
β
=
∫
B
‍
a
a
-
1
T
0
g
*
w
α
*
θ
β
d
V
+
∫
∂
B
‍
g
*
t
i
α
*
u
i
β
+
m
i
α
*
φ
i
β
+
λ
α
*
ω
β
1
T
0
+
∫
∂
B
‍
g
*
aa
-
1
T
0
h
*
q
α
*
θ
β
d
A
-
∫
Σ
‍
g
*
T
i
α
*
U
i
β
+
M
i
α
*
Φ
i
β
+
Λ
α
*
Ψ
β
1
T
0
+
∫
∂
B
‍
g
*
aa
-
1
T
0
h
*
Q
α
*
Θ
β
d
A
.
Finally, introducing (25) into (22), we arrive at the desired result (20).It is easy to see that in the absence of the discontinuities we obtain the generalization, in the context of the thermoelasticity of thermoelastic microstretch bodies, of the previous results established in the classical thermoelastodynamics.Based on relation (20), we now calculate the thermomechanical body loadings equivalent to a given dislocation. To this aim, we assume that u
i
(
2
), φ
i
(
2
), ω
(
2
), and θ
(
2
), as functions of (
t
,
x
), where x
=
(
x
i
), are of class C
∞
(
B
×
[
0
,
∞
)
). Of course, if the functions u
i
(
2
), φ
i
(
2
), ω
(
2
), and θ
(
2
) are given, then by means of (15), we can determine the functions F
i
(
2
), G
i
(
2
), L
(
2
), and r
(
2
).Now, we restrict our considerations only to the case whenU
i
(
2
)
=
Φ
i
(
2
)
=
0, i
=
1,2
,
3, Ψ
(
2
)
=
Θ
(
2
)
=
0, and η
(
2
) correspond to the faulted medium. Then, in the case of identical boundary conditions, we obtain
(26)
∫
B
‍
ϱ
F
i
1
*
u
i
2
+
G
i
1
*
φ
i
2
+
L
1
*
ω
2
1
T
0
a
a
a
a
a
-
1
T
0
h
*
r
1
*
θ
2
d
V
=
∫
B
‍
ϱ
F
i
2
*
u
i
1
+
G
i
2
*
φ
i
1
+
L
2
*
ω
1
1
T
0
a
a
a
a
a
a
a
a
a
-
1
T
0
h
*
r
2
*
θ
1
d
V
-
∫
Σ
‍
g
*
T
i
2
*
U
i
1
+
M
i
2
*
Φ
i
1
+
Λ
2
*
Ψ
1
1
T
0
a
a
a
a
a
a
a
a
a
a
a
a
-
1
T
0
h
*
Q
2
*
Θ
1
d
A
.
In view of (19), we have
(27)
T
i
(
2
)
=
A
i
j
m
n
ε
m
n
+
B
i
j
m
n
μ
m
n
+
D
i
j
m
γ
m
+
a
i
j
ω
-
E
i
j
θ
ν
j
,
M
i
(
2
)
=
B
i
j
m
n
ε
m
n
+
C
i
j
m
n
μ
m
n
+
E
i
j
m
γ
m
+
b
i
j
ω
-
D
i
j
θ
ν
j
,
Λ
(
2
)
=
D
m
n
i
ε
m
n
+
E
m
n
i
μ
m
n
+
C
i
j
γ
j
+
d
i
ω
-
L
i
θ
ν
i
,
Θ
(
2
)
=
k
i
j
θ
,
i
ν
j
.
Taking into account the definition of the Dirac translated measure, δ, we can prove the relation of the following type:
(28)
ψ
i
ξ
,
t
=
∫
B
‍
ψ
i
x
,
t
δ
x
-
ξ
d
V
,
ψ
i
,
j
ξ
,
t
=
∫
B
‍
ψ
i
x
,
t
δ
,
j
(
x
-
ξ
)
d
V
and then the relation (26) can be rewritten as follows:
(29)
∫
B
‍
ϱ
F
i
1
+
F
i
*
u
i
2
+
G
i
1
+
G
i
*
φ
i
2
1
T
0
∫
B
‍
ϱ
a
a
+
L
1
+
L
*
ω
2
-
1
T
0
*
r
1
+
R
*
θ
2
d
V
=
∫
B
‍
ϱ
F
i
2
*
u
i
1
+
G
i
2
*
φ
i
1
+
L
2
*
ω
1
1
T
0
=
∫
B
‍
a
a
a
-
1
T
0
h
*
r
(
2
)
*
θ
(
1
)
d
V
.
In the above relations we have used the notations
(30)
F
k
=
-
1
ϱ
∫
Σ
‍
A
j
i
r
k
U
i
1
+
B
j
i
r
k
Φ
i
1
+
D
j
k
r
Ψ
1
a
a
a
a
a
a
a
i
·
δ
,
r
x
-
ξ
ν
j
d
A
ξ
,
G
k
=
-
1
ϱ
∫
Σ
‍
B
j
i
r
k
U
i
1

+
C
j
i
r
k
Φ
i
1

+
E
j
k
r
Ψ
1
δ
,
r
x
-
ξ
G
k
=
-
1
ϱ
∫
Σ
‍
a
+
ε
m
n
k
A
j
i
m
n
U
i
1
+
B
j
i
m
n
Φ
i
1
+
D
j
m
n
Ψ
1
G
k
=
-
1
ϱ
∫
Σ
‍
a
B
j
i
r
k
U
i
1
+
δ
,
r
·
δ
x
-
ξ
ν
j
d
A
ξ
,
L
=
-
1
ϱ
∫
Σ
‍
D
j
i
m
U
i
1
+
E
j
i
m
Φ
i
1
+
C
j
m
Ψ
1
δ
,
m
L
=
-
1
ρ
∫
Σ
‍
aa
·
x
-
ξ
+
a
j
i
U
i
1
+
b
j
i
Φ
i
1
+
d
j
Ψ
1
L
=
-
1
ρ
∫
Σ
‍
ai
U
i
1
·
δ
x
-
ξ
ν
j
d
A
ξ
,
R
=
1
ϱ
∫
Σ
‍
T
0
E
j
i
U
i
(
1
)
+
D
j
i
Φ
i
(
1
)
+
L
j
Ψ
(
1
)
δ
(
x
-
ξ
)
R
=
1
ρ
∫
Σ
‍
T
0
a
i
E
j
i
U
i
1
+
δ
-
k
i
j
δ
,
i
x
-
ξ
ν
j
d
A
ξ
.In order to illustrate theoretical results derived in the previous sections, we present some numerical results taking into account the initial conditions given in relations (10).Figure1 shows the 3D variation of displacement with respect to time and space. Because we are interested only in the positive values of time, it is observed that the displacement increases sharply with x. The 3D variation of microstretch function with respect to time and space is shown in Figure 2. It is observed that for an increasing of time the microstretch function tends to become constant and the microstretch surface becomes a parallel plane with x
O
t.Figure 1
The variation ofu
(
x
,
t
).Figure 2
The variation ofω
(
x
,
t
).In Figure3 the 3D variation of microrotation function is represented and Figure 4 shows the variation of the temperature with increasing x under the effects of thermoelastic body.Figure 3
The variation ofϕ
(
x
,
t
).Figure 4
The variation ofθ
(
x
,
t
).Figure5 presents the displacement profiles u
(
x
,
t
) against x at different times. The displacement achieves a maximum for all chosen times for x
∈
(
-
2.5
;
-
2.35
) The displacement is an increasing function with respect to space fot x
∈
[
-
5
;
-
2.5
]
∪
[
3.4
;
6.3
], it is a decreasing function with respect to space for x
∈
(
-
2.5
;
3.4
) and tends to a steady-state.Figure 5
The behavior ofu
(
x
,
t
) at different values of time.Figure6 presents the microstretch profiles ω
(
x
,
t
) against x at different times. The microstretch is an increasing function with respect to time. It is clear that the microstretch function is not stable.Figure 6
The behavior ofω
(
x
,
t
) at different values of time.The profiles of microrotation functionϕ
(
x
,
t
) over space are shown in Figure 7 for four values of time. The microrotation function, as it results from Figure 4, is a decreasing function with respect to time.Figure 7
The behavior ofϕ
(
x
,
t
) at different values of time.Figure8 presents the temperature profiles θ
(
x
,
t
) against x at different times. The medium temperature is an increasing function with respect to time for x
>
0.Figure 8
The behavior ofθ
(
x
,
t
) at different values of time.Based on papers by Eringen [1, 14] and Kumar et al. [12] and based on the simulations performed in the laboratory of virtual engineering, we have used the following characteristic constitutive coefficients for the microstretch thermoelastic materials (in the isotropic case):
(31)
λ
=
9.4
·
1
0
10
N
m
-
2
,
μ
=
4
·
1
0
10
N
M
-
2
,
ρ
=
1.74
·
1
0
3
Kg
m
-
3
,
κ
=
1
0
10
N
m
-
2
,
α
=
0.779
·
1
0
-
9
N
,
η
=
1.22
kJ
/
mol
,
a
1
=
0.779
·
1
0
-
9
N
,
a
2
=
0.19
·
1
0
-
19
,
b
1
=
0.72
·
1
0
-
17
,
c
1
=
0.227
·
1
0
-
11
,
d
1
=
0.723
·
1
0
-
19
,
e
1
=
0.315
·
1
0
-
17
,
a
3
=
1.7
·
1
0
-
4
,
b
2
=
1.26
·
1
0
-
6
,
c
2
=
0.97
·
1
0
-
5
;
c
3
=
0.779
·
1
0
-
9
;
e
2
=
1.48
·
1
0
-
5
;
d
2
=
0.84
·
1
0
-
3
,
l
1
=
0.64
·
1
0
-
7
.It is observed that for the displacement vector greater values thanx
=
3.2 the solution is stabilized. The microstretch function increases asymptotically to x
=
5. The microrotation vector is stabilized for all chosen moment starting with x in the range (
3.6,3.85
). The temperature increases for x
∈
(
-
5
,
-
3.5
), decreases for x
∈
(
-
3.5,0
), and then increases asymptotically to x
=
5.


## 4. Conclusions

In the absence of the discontinuities we obtain the generalization, in the context of the thermoelasticity of thermoelastic microstretch bodies, of the previous results established in the classical thermoelastodynamics.In the classical elasticity in the system of equations of motion occur only three unknown functions, while in our case the number of unknowns is seven.Also, while in the classical elasticity the equations of motion are parabolic, hyperbolic, or elliptic only, in our case we have a mixed system of equations.Taking into account relation (29) we deduce that the effect of the discontinuities across the surface Σ can be represented by extra external body loads and heat supply.Although these are supposed to act in an unfaulted medium and cannot in any sense represent real forces acting in the real medium, they may nevertheless provide, as pointed in [3, 4, 7], a useful theoretical tool and this because if two dislocations have the same equivalent force, they also emit the same radiation.



---

*Source: 850261-2015-03-02.xml*