Next Article in Journal
Total Roman {3}-domination in Graphs
Next Article in Special Issue
A GL Model on Thermo-Elastic Interaction in a Poroelastic Material Using Finite Element Method
Previous Article in Journal
Blind Image Watermarking in Canonical and Cepstrum Domains Based on 4-Connected t-o’clock Scrambling
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Study of Deformations in a Thermoelastic Dipolar Body with Voids

1
Department of Mathematics and Computer Science, Transilvania University of Brasov, 500093 Brasov, Romania
2
Department of Mathematics, Faculty of Science, Sohag University, Sohag 82524, Egypt
3
Department of Mechanical Engineering, Transilvania University of Brasov, 500093 Brasov, Romania
4
Faculty of Mechanical, Industrial and Maritime Engineering, Ovidius University of Constanta, 900527 Constanta, Romania
*
Author to whom correspondence should be addressed.
Symmetry 2020, 12(2), 267; https://doi.org/10.3390/sym12020267
Submission received: 5 January 2020 / Revised: 17 January 2020 / Accepted: 31 January 2020 / Published: 9 February 2020
(This article belongs to the Special Issue Composite Structures with Symmetry)

Abstract

:
In this paper, we consider the mixed initial boundary value problem in the context of a thermoelastic porous body having a dipolar structure. We intend to analyze the rate of decay of solutions to this problem to ensure that in a finite time, they become null. In our main result, we find that the combined contribution of the dipolar constitution of the body together with voids dissipation and thermal behavior cannot cause vanishing of the deformations in a finite time.

1. Introduction

We must outline that our study is dedicated to the dipolar structure for a porous thermoelastic body. The unanimous opinion of the specialists is that the porous media with a dipolar structure very accurately shapes the structure and behavior of the bones of humans and, obviously, of animals.
This kind of material is an integral part of a more general theory, namely the microstructure, whose initiator was Eringen (see, for instance, [1,2]). The utility of the dipolar bodies can be observed from the great number of appeared papers in this subject, see for instance [3,4,5,6,7]. Therefore, our present study is a continuation of research in this domain.
Another kind of the microstructure is also obtained by considering the pores (voids) in the materials, and this was starting by the paper [8] of Nunziato and Cowin. Then, the number of papers with this subject has become impressive. Here, we mention some of these [9,10,11,12,13,14,15,16]. In the opinion of many researchers, the backward in time problem and the problem of localization in a finite time of solutions have many points in common. It is considered that the initiator of the backward in time problems is Serrin which approached this kind of problem by taking into account the Navier–Stokes Equations (see [17]). Some results of uniqueness for the forward in time problem can be found in the study [17]. Recently, the number of articles having as a subject the problem of backward in time has increased greatly. Therefore, we can recommend the work in [18,19,20,21,22,23,24,25,26,27]. As a particular observation, the results of Quintanilla in [25] improves the studies of Ciarletta [23] and Ciarletta and Chirita [24]. Other studies of Quintanilla solve some questions regarding the location in time considering the solutions for problems back in time even in the case of theory of thermoelastic bodies with voids and of the theory of Green and Naghdi for thermoelasticity [26,27]. Iovane and Passarella in [28] approached also the elastic bodies with voids. Passarella and Tibullo consider in [29] the problem of backward in time for the thermo-microstretch elastic bodies. Similar results can be found in [30,31].
We should point out that our idea to consider the problem of localization in a finite time of solutions in the case of the theory of dipolar bodies with voids was inspired by the paper [32] of Quintanilla and Straughan.
The structure of our paper is the following.
In Section 2, we write the system of the differential equations that govern the evolution of our body. After that, as usual, we add the initial conditions and, to complete the mixed problem, we consider the boundary data. Then, we continue with few auxiliary results, which offer some estimations regarding the displacements, the displacements of dipolar type and regarding the function of voids. The main result of our work is regarding the localization in time of solution four our mixed problem. In fact, we prove that is not possible to locate a finite time for vanishing of the deformations.

2. Preliminaries

As already seen, our study is dedicated to a thermoelastic dipolar body with voids. We consider the general of an anisotropic and non-homogeneous body, which occupies a regular region D, included in the Euclidean space R 3 . The border of the domain D is a regular surface, denoted by D . As usual, we use the notation D ¯ for the closure of the domain D, and we have D ¯ = D D . We will suppose that the vectors and tensors have components with indices over 1, 2, and 3. The convention for summation in the case of repeated index is used and the derivation operations are corresponding defined. Therefore, a partial derivative with respect to a spatial coordinate is designates by a subscript preceded by a comma. A superposed dot is used to designate the derivative with regards to t, the time time variable. All functions used in our study are supposed be sufficient regular as necessary. Sometimes, the dependence of functions regarding the time or spatial independent variables is omitted, of course, when is no possibility of confusion.
To describe the evolution of a dipolar thermoelastic body, we will use the following variables,
u i ( t , x ) , ϕ i j ( t , x ) , θ ( t , x ) , φ ( t , x ) ,
where we denoted by u = u i the vector field for displacement, and by ϕ = ϕ i j the tensor field for dipolar displacement. Also, here, θ is the temperature variation and φ is a function for the volume distribution.
All these functions are defined for ( t , x ) [ 0 , t 0 ) × D .
Using the above variables u i ( x , t ) , ϕ i j ( x , t ) , we will introduce the components of the tensors of strain, namely, e i j , ν i j and χ i j k , as follows,
2 e i j = u j , i + u i , j , ν i j = u j , i ϕ i j , χ i j k = ϕ j k , i .
We must specify that our further considerations refers to linear theory as such we can consider the Helmholtz’s free energy as a quadratic form, regarding to all its independent variables. Let us denote by W the Helmholtz’s free density of energy in the initial configuration. As a consequence, in accordance with the principle of conservation of energy, we will develop the function W in series and we retain only the terms of first order and of second order. If we take into account that the initial state (the reference state) is supposed free of charges, we are led to the conclusion that the energy density of Helmholtz has the following expression,
W = 1 2 A i j m n e i j e m n + D i j m n e i j ν m n + F i j m n r e i j χ m n r + 1 2 B i j m n ν i j ν m n     + G i j m n r ν i j χ m n r + 1 2 C i j k m n r χ i j k χ m n r + a i j e i j θ + b i j ν i j θ + c i j k χ i j k θ     1 2 c θ 2 m θ φ + 1 2 A i j φ , i φ , j + α i j e i j φ + β i j ν i j φ + δ i j k χ i j k φ + 1 2 ξ φ 2 .
If we substitute the Helmholtz’s energy in the entropy production inequality, we obtain a relation from where the equations of motion are deduced. Also, the same relation can be used to deduce the constitutive equations. As is known, the constitutive equations give the expression of stress tensors as functions of the strain tensors and some constants of the material. In what follows, we will use the notations t i j , σ i j and η i j k for the stress tensors. As such, the connections between the tensors t i j , σ i j , η i j k and the tensors e i j , ν i j , χ i j k are the constitutive equations
Inspired by the procedure used by Green and Rivlin in the paper [6], we adopt a similar technique, so that by taking into account the Helmholtz’s energy (3) we are led to the following constitutive equations,
t i j = W e i j = A i j m n e m n + D m n i j ν m n + F m n r i j χ m n r + a i j θ + α i j φ , σ i j = W ν i j = D i j m n e m n + B i j m n ν m n + G i j m n r χ m n r + b i j θ + β i j φ , η i j k = W χ i j k = F i j k m n e m n + G m n i j k ν m n + C i j k m n r χ m n r + c i j k θ + δ i j φ , h i = W φ , j = A i j φ , i , S = W θ = a i j e i j b i j ν i j c i j k χ i j k + c θ + m φ , g = W φ = α i j e i j β i j ν i j δ i j k χ i j k + m θ ξ φ ,
where we use the notation S to designate the entropy (per unit mass).
We will consider that the above equations take place in the cylinder [ 0 , t 0 ) × D .
Denoting by q the heat flux vector of components q i , we can deduce a constitutive relation, similar to that from the classical theory, namely,
q j = K i j θ , i ,
where we denoted by K i j the tensor of the thermal conductivity and assume that K i j is a symmetric tensor.
In the absence of supply terms, the differential equations that govern the thermoelasticity of dipolar bodies with voids, obtained as in [5,6], are as follows.
The equations of motion:
t i j + σ i j , j = ρ u ¨ i , η i j k , i + σ j k = I k r ϕ ¨ j r ;
The balance of the equilibrated forces:
h i , i + g = J φ ¨ ;
The equation of energy:
ρ T 0 S ˙ = q i , i + ρ r .
In the preceding equations, we used some notions with following signification; ρ -the constant mass density, and I i j -the tensor of microinertia, which is a symmetric tensor.
A i j m n , B i j m n , , a i j are the coefficients that characterize the properties of the material with regards to the elasticity. Clearly, from (3), we can deduce the following symmetry relations,
A i j m n = A j i m n = A m n i j , B i j m n = B m n i j , A i j = A j i , a i j = a j i , C i j k m n r = C m n r i j k , F i j k m n = F i j k n m , D i j m n = D i j n m , α i j = α j i .
If we consider Equations (6)–(8), and then involve the kinematic equations (2) and the constitutive relations (4) and (5), from the equations of motion and the energy equations, we obtain the following system of equations,
ρ u ¨ i = C i j m n + G i j m n e m n + G m n i j + B i j m n ν m n + + F m n r i j + D i j m n r χ m n r + a i j + b i j θ + α i j + β i j φ , j , I k r ϕ ¨ j r = F i j k m n e m n + D m n i j k ν m n + A i j k m n r χ m n r + c i j k θ + δ i j k φ , i + + G j k m n e m n + B j k m n ν m n + D j k m n r χ m n r + b j k θ + β j k φ , c θ ˙ = K i j θ , i , j + a i j e ˙ i j + b i j ν i j + c i j k χ ˙ i j k m φ ˙ , J φ ¨ = A i j φ , i , j α i j e ˙ i j β i j ν i j δ i j k χ ˙ i j k + m θ ξ φ ,
which take place for ( t , x ) [ 0 , ) × D .
We also suppose that the Equations (2) and (4) take place for ( t , x ) [ 0 , ) × D .
We now want to define a surface tractions vector of components t i , a surface couple tensor of components η j k , the flux of heat denoted by q, and the traction of the equilibrated stress vector h. For this we need the unit normal to the border D , of components n i , outward oriented. These have the following expressions,
t i = t i j + σ i j n j , η j k = η i j k n i , q = q i n i , h = h i n i ,
and are defined at regular points of the border D .
Having these tractions, we can now consider the following boundary data, in their nonhomogeneous form,
u i ( t , x ) = u ˜ i , ( t , x ) [ 0 , ) × D u , t i ( t , x ) = t ˜ i , ( t , x ) [ 0 , ) × D u c , ϕ i j ( t , x ) = ϕ ˜ i j , ( t , x ) [ 0 , ) × D ϕ , η j k ( t , x ) = η ˜ j k , ( t , x ) [ 0 , ) × D ϕ c , θ ( t , x ) = θ ˜ , ( t , x ) [ 0 , ) × D θ , q ( t , x ) = q ˜ , ( t , x ) [ 0 , ) × D θ c , φ ( t , x ) = φ ˜ , ( t , x ) [ 0 , ) × D φ , h ( t , x ) = h ˜ , ( t , x ) [ 0 , ) × D φ c ,
where the subsets D u , D ϕ , D θ , and D φ and its complements D u c , D ϕ c , D θ c , D φ c are surfaces from the boundary D , having the following properties,
D ¯ u D u c = D ¯ ϕ D ϕ c = D ¯ θ D θ c = D ¯ φ D φ c = D , D u D u c = D ϕ D ϕ c = D θ D θ c = D φ D φ c = .
The mixed problem will be complete if we add the initial restrictions:
u i ( x , 0 ) = u i 0 ( x ) , u ˙ i ( x , 0 ) = u i 1 ( x ) , θ ( x , 0 ) = θ 0 ( x ) , ϕ i j ( x , 0 ) = ϕ i j 0 ( x ) , ϕ ˙ i j ( x , 0 ) = ϕ i j 1 ( x ) , φ ( x , 0 ) = φ 0 ( x ) ,
which take place on the closed domain D ¯ .
Here, u i 0 ( x ) , u i 1 ( x ) , ϕ i j 0 ( x ) , ϕ i j 1 ( x ) , θ 0 ( x ) , and φ 0 ( x ) are given functions and are assume be continuous. Furthermore, we must suppose that the above functions are adequacy with the restrictions (11) on the respective subsurfaces of D .
By Ψ we denote the internal energy and consider it as a quadratic form, as follows,
Ψ = 1 2 A i j m n e i j e m n + D i j m n e i j ν m n + F i j m n r e i j χ m n r + + 1 2 B i j m n ν i j ν m n + G i j m n r ν i j χ m n r + 1 2 C i j k m n r χ i j k χ m n r .
Let us denote by P our initial boundary value problem, which includes the differential equations from (10), the initial restrictions from (12), and the data to the limit (11).
The results we want to expose will only be possible if certain restrictions are met, imposed on all the functions we will deal with.
First, we denote by J m ( x ) the smallest eigenvalue of the inertia tensor I i j ( x ) , and we must suppose that J m is a continuous function. Also, the density ρ and the constitutive coefficients are assumed be functions of class C 1 ( D ) . Also, we must suppose that
(a)
ρ ( x ) a 1 , J m ( x ) a 2 , c ( x ) a 3 , where the real constants a 1 , a 2 , a 3 are positive;
(b)
K i j is a positive definite tensor;
(c)
the quadratic form Ψ is positive definite.
As a consequence of the hypothesis (b), we can deduce the inequalities:
K m θ , i θ , j K i j θ , i θ , j K M θ , i θ , j ,
where K m and K M are two positive constants.
If we take into account the hypothesis (c), then we can determine two positive numbers M 1 and M 2 to satisfy the double inequality that follows,
M 1 2 e i j e i j + ν i j ν i j + χ i j k χ i j k Ψ M 2 2 e i j e i j + ν i j ν i j + χ i j k χ i j k .
The above restrictions are commonly imposed in mechanics of solids, therefore they cannot be considered as very restrictive.
As such, our mixed problem P consists of the following restrictions and equations.
-
the equations of motion (10) 1 and (10) 2 , which take place in cylinder D × [ 0 , ) ;
-
the equation of energy:
K i j θ , j , i = T 0 a i j u ˙ i , j + b i j u ˙ j , i ϕ ˙ i j + c i j k ϕ ˙ i j , k + c θ ˙ , in [ 0 , ) × D ;
-
the kinematic equations (2), which take place in cylinder [ 0 , ) × D ;
-
the constitutive conditions (4), which take place in cylinder [ 0 , ) × D ;
-
the initial restrictions (11), which take place in D ¯ ;
-
the conditions to the limit:
u i ( x , t ) = 0 , ( x , t ) [ 0 , ) × D u , t i = 0 , ( x , t ) [ 0 , ) × D u c , ϕ i j ( x , t ) = 0 , ( x , t ) [ 0 , ) × , D ϕ m j k = 0 , ( x , t ) [ 0 , ) × D ϕ c , θ ( x , t ) = 0 , ( x , t ) [ 0 , ) × D θ , q = 0 , ( x , t ) [ 0 , ) × D θ c , φ ( x , t ) = 0 , ( x , t ) [ 0 , ) × D φ , h = 0 , ( x , t ) [ 0 , ) × D φ c .

3. Auxiliary Results

At the beginning of this section, we will establish some identities of integral type, with regards to the solutions u = u i , ϕ i j , θ , φ of the above mixed problem P . The important results of our study will be based on these identities.
Proposition 1.
If u = u i , ϕ i j , θ is a solution of our above problem P , then the following identity is satisfied
1 2 B ρ u ˙ i ( t ) u ˙ i ( t ) + I j k ϕ ˙ j m ( t ) ϕ ˙ k m ( t ) + J φ ˙ 2 ( t ) + c θ 2 ( t ) d V + + B 1 2 A i j m n e i j ( t ) e m n ( t ) + D i j m n e i j ( t ) ν m n ( t ) + F i j m n r e i j ( t ) χ m n r ( t ) + + 1 2 B i j m n ν i j ( t ) ν m n ( t ) + G i j m n r ν i j ( t ) χ m n r ( t ) + 1 2 C i j k m n r χ i j k ( t ) χ m n r ( t ) d V + + B A i j φ , i φ , j + α i j e i j φ + β i j ν i j φ + δ i j k χ i j k φ + ξ φ 2 d V = + 0 t D 1 T 0 K i j θ , i ( s ) θ , j ( s ) d V d s , t 0 .
Proof. 
We multiply the equations of motion (10) 1 by u ˙ i and consider the kinematic Equation (2), after simple calculations, the following identity is obtained,
ρ u ˙ i u ˙ i = A i j m n + G i j m n e m n + G m n i j + B i j m n ν m n + + F m n r i j + D i j m n r χ m n r + a i j + b i j θ + α i j + β i j φ u ˙ i , j A i j m n + G i j m n e m n + G m n i j + B i j m n ν m n + + F m n r i j + D i j m n r χ m n r + a i j + b i j θ + α i j + β i j φ u ˙ i , j .
If we multiply the equations of motion (10) 2 by ϕ ˙ i j and consider the kinematic Equation (2), after simple calculations, the following identity is obtained,
I j k ϕ ˙ j m ϕ ˙ k m = F i j k m n e m n + D m n i j k ν m n + C i j k m n r χ m n r + c i j k θ + δ i j k φ ϕ ˙ j k , i F i j k m n e m n + D m n i j k ν m n + C i j k m n r χ m n r + c i j k θ + δ i j k φ ϕ ˙ j k , i + + G j k m n e m n + B j k m n ν m n + D j k m n r χ m n r + b j k θ + β j k φ ϕ j k .
Now we multiply the equations of motion (10) 3 by θ . After simple calculations the following identity is obtained,
c θ θ ˙ = K i j θ , i θ , j + K i j θ , i θ , j + + a i j e ˙ i j θ + b i j ν ˙ i j θ + c i j k χ ˙ i j k θ m θ φ ˙ .
Finally, we multiply the equations of motion (10) 4 by φ , so that after simple calculations, the following identity is obtained,
J φ ˙ φ ¨ = A i j φ , i φ ˙ , j A i j φ , i φ ˙ , j α i j e i j φ ˙ β i j ν i j φ ˙ δ i j k χ i j k φ ˙ ξ φ φ ˙ + m θ φ ˙ .
Now we sum up, term with term, the equalities (19)–(22) and take into account the geometric relations (5) so that we reach the equality
ρ u ˙ i u ¨ i + I j k ϕ ˙ j m ϕ ¨ k m + c θ θ ˙ + J φ ˙ φ ¨ = = A i j m n + G i j m n e m n + G m n i j + B i j m n ν m n + + F m n r i j + D i j m n r χ m n r + a i j + b i j θ + α i j + β i j φ u ˙ i , j + + F i j k m n e m n + D m n i j k ν m n + C i j k m n r χ m n r + c i j k θ + δ i j k φ ϕ ˙ j k , i A i j m n + G i j m n e m n + G m n i j + B i j m n ν m n + + F m n r i j + D i j m n r χ m n r + a i j + b i j θ + α i j + β i j φ u ˙ i , j F i j k m n e m n + D m n i j k ν m n + C i j k m n r χ m n r + c i j k θ + δ i j k φ ϕ ˙ j k , i + + G j k m n e m n + B j k m n ν m n + D j k m n r χ m n r + b j k θ + β j k φ ϕ j k K i j θ , i θ , j + K i j θ , i θ , j + a i j e ˙ i j θ + b i j ν ˙ i j θ + c i j k χ ˙ i j k θ + + A i j φ , i φ ˙ , j A i j φ , i φ ˙ , j + α i j e i j φ ˙ + β i j ν i j φ ˙ + δ i j k χ i j k φ ˙ ξ φ φ ˙ ,
which can be restated in the following form
ρ u ˙ i u ¨ i + I j k ϕ ˙ j m ϕ ¨ k m + c θ θ ˙ + J φ ˙ φ ¨ + + A i j m n e m n e ˙ i j + G m n i j ν m n e ˙ i j + ν ˙ m n e i j + F m n r i j χ m n r e ˙ i j + χ ˙ m n r e i j + B i j m n ν m n ν ˙ i j + D m n i j k ν m n ν m n χ ˙ i j k + ν ˙ m n χ i j k + C i j k m n r χ m n r χ ˙ i j k + + α i j e i j φ ˙ + e ˙ i j φ + β i j ν i j φ ˙ + ν ˙ i j φ + δ i j k χ i j k φ ˙ + χ ˙ i j k φ + + A i j φ , i φ ˙ , j + ξ φ φ ˙ = = A i j m n + G i j m n e m n + G m n i j + B i j m n ν m n + + F m n r i j + D i j m n r χ m n r + a i j + b i j θ + α i j + β i j φ u ˙ i , j + + F i j k m n e m n + D m n i j k ν m n + C i j k m n r χ m n r + c i j k θ + δ i j k φ ϕ ˙ j k , i + + A i j φ , i φ ˙ , j K i j θ , i θ , j + K i j θ , i θ , j .
We can now integrate this identity over domain [ 0 , t ] × D . Taking into account the divergence theorem and considering that we have null boundary data, we obtain the identity (18), as such the proof of Proposition 1 is finished. □
A new energy relation is obtained in next proposition. We consider the case of null initial data.
Proposition 2.
Let u i , ϕ i j , θ , φ be a solution of the mixed problem, P , corresponding to zero initial conditions. Then, the following equality is satisfied:
1 2 B ρ u ˙ i u ˙ i + I j k ϕ ˙ j m ϕ ˙ k m + J φ ˙ 2 c θ 2 d V + + B 1 2 A i j m n e i j e m n + D i j m n e i j ν m n + F i j m n r e i j χ m n r + + 1 2 B i j m n ν i j ν m n + G i j m n r ν i j χ m n r + 1 2 C i j k m n r χ i j k χ m n r d V + + B 1 2 A i j φ , i φ , j + α i j e i j φ + β i j ν i j φ + δ i j k χ i j k φ + 1 2 ξ φ 2 d V = = 0 t D K i j θ , i θ , j a i j e ˙ i j θ b i j ν ˙ i j θ c i j k χ ˙ i j k θ m θ φ ˙ d V d τ ,
for any t [ 0 , ) .
Proof. 
If we multiply the equation of motion (10) 1 by u ˙ i , then Equation (19) is obtained. If we multiply the equation of motion (10) 2 by ϕ ˙ i j , then Equation (20) is obtained. If we multiply the equation of motion (10) 4 by φ ˙ , then Equation (22) is obtained. Now, we multiply the equation of motion (10) 3 by θ so that the following identity is obtained,
c θ θ ˙ = K i j θ , i θ , j K i j θ , i θ , j a i j e ˙ i j θ b i j ν ˙ i j θ c i j k χ ˙ i j k θ + m θ φ ˙ .
From Equations (19), (20), (22), and (26) we deduce
ρ u ˙ i u ¨ i + I j k ϕ ˙ j m ϕ ¨ k m c θ θ ˙ + J φ ˙ φ ¨ + + A i j m n e m n e ˙ i j + G m n i j ν m n e ˙ i j + ν ˙ m n e i j + F m n r i j χ m n r e ˙ i j + χ ˙ m n r e i j + B i j m n ν m n ν ˙ i j + D m n i j k ν m n ν m n χ ˙ i j k + ν ˙ m n χ i j k + C i j k m n r χ m n r χ ˙ i j k + + α i j e i j φ ˙ + e ˙ i j φ + β i j ν i j φ ˙ + ν ˙ i j φ + δ i j k χ i j k φ ˙ + χ ˙ i j k φ + + A i j φ , i φ ˙ , j + ξ φ φ ˙ = = A i j m n + G i j m n e m n + G m n i j + B i j m n ν m n + + F m n r i j + D i j m n r χ m n r + a i j + b i j θ + α i j + β i j φ u ˙ i , j + + F i j k m n e m n + D m n i j k ν m n + C i j k m n r χ m n r + c i j k θ + δ i j k φ ϕ ˙ j k , i + + A i j φ , i φ ˙ , j + K i j θ , i θ , j K i j θ , i θ , j + + a i j e ˙ i j θ + b i j ν ˙ i j θ + c i j k χ ˙ i j k θ + m θ φ ˙ .
Finally, we integrate the last identity over domain [ 0 , t ] × D , so that considering the divergence theorem and taking into account that we have null initial data, we obtain the identity (25), as such the proof of Proposition 2 is finished. □
In the following proposition, we deduce another energy relation, obtained by a method of the Lagrange identity type.
Proposition 3.
Let u i , ϕ i j , θ , φ be a solution of the mixed problem P , corresponding to zero initial conditions and null boundary data. Then, the following equality is satisfied,
B ρ u ˙ i u ˙ i + I j k ϕ ˙ j m ϕ ˙ k m + J φ ˙ 2 c θ 2 d V = = B A i j m n e i j e m n + 2 D i j m n e i j ν m n + 2 F i j m n r e i j χ m n r + + B i j m n ν i j ν m n + 2 G i j m n r ν i j χ m n r + C i j k m n r χ i j k χ m n r d V + + B A i j φ , i φ , j + 2 α i j e i j φ + 2 β i j ν i j φ + 2 δ i j k χ i j k φ + ξ φ 2 d V ,
for any t [ 0 , ) .
Proof. 
Let us consider the Equation (10) 1 and fix an arbitrary t ( 0 , ) . By simple calculations, we reach the identity
τ ρ u ˙ i ( τ ) u ˙ i ( 2 t τ ) = ρ u ¨ i ( τ ) u ˙ i ( 2 t τ ) ρ u ˙ i ( τ ) u ¨ i ( 2 t τ ) = = A i j m n + G i j m n e m n ( τ ) + G m n i j + B i j m n ν m n ( τ ) + + F m n r i j + D i j m n r χ m n r ( τ ) + a i j + b i j θ ( τ ) + α i j + β i j φ ( τ ) u ˙ i ( 2 t τ ) , j A i j m n + G i j m n e m n ( τ ) + G m n i j + B i j m n ν m n ( τ ) + + F m n r i j + D i j m n r χ m n r ( τ ) + a i j + b i j θ ( τ ) + α i j + β i j φ ( τ ) u ˙ i , j ( 2 t τ ) A i j m n + G i j m n e m n ( 2 t τ ) + G m n i j + B i j m n ν m n ( 2 t τ ) + + F m n r i j + D i j m n r χ m n r ( 2 t τ ) + a i j + b i j θ ( 2 t τ ) + α i j + β i j φ ( 2 t τ ) u ˙ i ( τ ) , j + + A i j m n + G i j m n e m n ( 2 t τ ) + G m n i j + B i j m n ν m n ( 2 t τ ) + + F m n r i j + D i j m n r χ m n r ( 2 t τ ) + a i j + b i j θ ( 2 t τ ) + α i j + β i j φ ( 2 t τ ) u ˙ i , j ( τ ) .
Now we take into account the Equation (10) 2 , for an arbitrary fixed t ( 0 , ) . After some calculations, we reach the identity
τ I j k ϕ ˙ j m ( τ ) ϕ ˙ k m ( 2 t τ ) = I j k ϕ ¨ j m ( τ ) ϕ ˙ k m ( 2 t τ ) I j k ϕ ˙ j m ( τ ) ϕ ¨ k m ( 2 t τ ) = = F i j k m n e m n ( τ ) + D m n i j k ν m n ( τ ) + C i j k m n r χ m n r ( τ ) + c i j k θ ( τ ) + δ i j k φ ( τ ) ϕ ˙ j k ( α ) , i F i j k m n e m n ( τ ) + D m n i j k ν m n ( τ ) + C i j k m n r χ m n r ( τ ) + c i j k θ ( τ ) + δ i j k φ ( τ ) ϕ ˙ j k , i ( α ) + + G j k m n e m n ( τ ) + B j k m n ν m n ( τ ) + D j k m n r χ m n r ( τ ) + b j k θ ( τ ) + β j k φ ( τ ) ϕ j k ( α ) F i j k m n e m n ( α ) + D m n i j k ν m n ( α ) + C i j k m n r χ m n r ( α ) + c i j k θ ( α ) + δ i j k φ ( α ) ϕ ˙ j k ( τ ) , i + + F i j k m n e m n ( α ) + D m n i j k ν m n ( α ) + C i j k m n r χ m n r ( α ) + c i j k θ ( α ) + δ i j k φ ( α ) ϕ ˙ j k , i ( τ ) G j k m n e m n ( α ) + B j k m n ν m n ( α ) + D j k m n r χ m n r ( α ) + b j k θ ( α ) + β j k φ ( α ) ϕ j k ( τ ) ,
where, to simplify the writing, we used the notation 2 t τ = α .
Now, similarly, for an arbitrarily fixed t ( 0 , ) we consider the Equation (10) 3 to deduce
τ c θ ( τ ) θ ( 2 t τ ) = c θ ˙ ( τ ) θ ( 2 t τ ) c θ ( τ ) θ ˙ ( 2 t τ ) = = K i j θ , i ( τ ) θ ( 2 t τ ) , j + a i j e ˙ i j ( τ ) + b i j ν ˙ i j ( τ ) + c i j k χ ˙ i j k ( τ ) m φ ˙ ( τ ) θ ( 2 t τ ) + + K i j θ , i ( 2 t τ ) θ ( τ ) , j a i j e ˙ i j ( 2 t τ ) + b i j ν ˙ i j ( 2 t τ ) + c i j k χ ˙ i j k ( 2 t τ ) m φ ˙ ( 2 t τ ) θ ( τ ) .
Finally, for an arbitrarily fixed t ( 0 , ) we will use the Equation (10) 4 to obtain
τ J φ ˙ ( τ ) φ ˙ ( 2 t τ ) = J φ ¨ ( τ ) φ ˙ ( 2 t τ ) J φ ˙ ( τ ) φ ¨ ( 2 t τ ) = = A i j φ , i ( τ ) φ ˙ ( 2 t τ ) , j A i j φ , i ( τ ) φ ˙ , j ( 2 t τ ) α i j e i j ( τ ) φ ˙ ( 2 t τ ) β i j ν i j ( τ ) φ ˙ ( 2 t τ ) δ i j k χ i j k ( τ ) φ ˙ ( 2 t τ ) ξ φ ( τ ) φ ˙ ( 2 t τ ) + m θ ( τ ) φ ˙ ( 2 t τ ) A i j φ , i ( 2 t τ ) φ ˙ ( τ ) , j + A i j φ , i ( 2 t τ ) φ ˙ , j ( τ ) + α i j e i j ( 2 t τ ) φ ˙ ( τ ) + + β i j ν i j ( 2 t τ ) φ ˙ ( τ ) + δ i j k χ i j k ( 2 t τ ) φ ˙ ( τ ) + ξ φ ( 2 t τ ) φ ˙ ( τ ) m θ ( 2 t τ ) φ ˙ ( τ ) .
Now we use the geometric Equation (2) so that from Equations (29)–(32) we are led to the identity
τ ρ u ˙ i ( τ ) u ˙ i ( 2 t τ ) + I j k ϕ ˙ j m ( τ ) ϕ ˙ k m ( 2 t τ ) c θ ( τ ) θ ( 2 t τ ) + J φ ˙ ( τ ) φ ˙ ( 2 t τ ) = + A i j m n + G i j m n e m n ( τ ) + G m n i j + B i j m n ν m n ( τ ) + + F m n r i j + D i j m n r χ m n r ( τ ) + a i j + b i j θ ( τ ) + α i j + β i j φ ( τ ) u ˙ i ( 2 t τ ) , j A i j m n + G i j m n e m n ( 2 t τ ) + G m n i j + B i j m n ν m n ( 2 t τ ) + + F m n r i j + D i j m n r χ m n r ( 2 t τ ) + a i j + b i j θ ( 2 t τ ) + α i j + β i j φ ( 2 t τ ) u ˙ i ( τ ) , j + + F i j k m n e m n ( τ ) + D m n i j k ν m n ( τ ) + C i j k m n r χ m n r ( τ ) + c i j k θ ( τ ) + δ i j k φ ( τ ) ϕ ˙ j k ( 2 t τ ) , i F i j k m n e m n + D m n i j k ν m n + C i j k m n r χ m n r + c i j k θ + δ i j k φ ( 2 t τ ) ϕ ˙ j k ( τ ) , i + + K i j θ , j ( τ ) θ ( 2 t τ ) , i K i j θ , i ( 2 t τ ) θ ( τ ) , j + A i j φ , i ( τ ) φ ˙ ( 2 t τ ) , j A i j φ , i ( 2 t τ ) φ ˙ ( τ ) , j + τ A i j m n e i j ( τ ) e m n ( 2 t τ ) + D i j m n e i j ( τ ) ν m n ( 2 t τ ) + e i j ( 2 t τ ) ν m n ( τ ) + + F i j k m n e m n ( τ ) χ i j k ( 2 t τ ) + e m n ( 2 t τ ) χ i j k ( τ ) + B i j m n ν i j ( τ ) ν m n ( 2 t τ ) + + G i j k m n ν m n ( τ ) χ i j k ( 2 t τ ) + ν m n ( 2 t τ ) χ i j k ( τ ) + C i j k m n r χ i j k ( τ ) χ m n r ( 2 t τ ) + + τ α i j e i j ( 2 t τ ) φ ˙ ( τ ) + e i j ( τ ) φ ˙ ( 2 t τ ) + β i j ν i j ( 2 t τ ) φ ˙ ( τ ) + ν i j ( τ ) φ ˙ ( 2 t τ ) + + τ δ i j k χ i j k ( 2 t τ ) φ ˙ ( τ ) + χ i j k ( τ ) φ ˙ ( 2 t τ ) + ξ φ ( 2 t τ ) φ ˙ ( τ ) .
This last equality is integrated over [ 0 , t ] × D . Applying the divergence theorem and taking into account the initial conditions (12) and the boundary data (17), we arrive at the desired identity (28). As such, the proof of proposition is finished. □

4. On Localization in Time of Solutions

The following theorem is essential to obtain the main result of our study.
Theorem 1.
Suppose the main hypotheses (a)–(c) are satisfied. If u i , ϕ i j , θ , φ is a solution of the mixed problem P , which corresponds to zero initial conditions (12) and null boundary data (17), then this solution can only be null, that is,
u i ( t , x ) = 0 , ϕ i j ( t , x ) = 0 , θ ( t , x ) = 0 , φ ( t , x ) = 0 , ( t , x ) [ 0 , ) × D .
Proof. 
By combining the results from Equations (25) and (28), we obtain
B A i j m n e i j e m n + 2 D i j m n e i j ν m n + 2 F i j m n r e i j χ m n r + + B i j m n ν i j ν m n + 2 G i j m n r ν i j χ m n r + C i j k m n r χ i j k χ m n r d V + + B A i j φ , i φ , j + 2 α i j e i j φ + 2 β i j ν i j φ + 2 δ i j k χ i j k φ + ξ φ 2 d V = = 0 t D K i j θ , i θ , j a i j e ˙ i j θ b i j ν ˙ i j θ c i j k χ ˙ i j k θ m θ φ ˙ d V d τ ,
for any t [ 0 , ) .
Inspired by the identities (18) and (35), we introduce the functions G, E and F, defined by
E ( t ) = 1 2 B ρ u ˙ i u ˙ i + I j k ϕ ˙ j m ϕ ˙ k m + J φ ˙ 2 + c θ 2 d V + + B 1 2 A i j m n e i j e m n + D i j m n e i j ν m n + F i j m n r e i j χ m n r + + 1 2 B i j m n ν i j ν m n + G i j m n r ν i j χ m n r + 1 2 C i j k m n r χ i j k χ m n r d V + + B 1 2 A i j φ , i φ , j + α i j e i j φ + β i j ν i j φ + δ i j k χ i j k φ + 1 2 ξ φ 2 d V , G ( t ) = B 1 2 A i j m n e i j e m n + D i j m n e i j ν m n + F i j m n r e i j χ m n r + + 1 2 B i j m n ν i j ν m n + G i j m n r ν i j χ m n r + 1 2 C i j k m n r χ i j k χ m n r d V + + B 1 2 A i j φ , i φ , j + α i j e i j φ + β i j ν i j φ + δ i j k χ i j k φ + 1 2 ξ φ 2 d V , F ( t ) = G ( t ) + α E ( t ) ,
for any t [ 0 , ) .
In (36) 3 α is a positive constant, conveniently of small chosen.
It is easy to see that F ( t ) can be rewritten in the form
F ( t ) = 1 2 B α ρ u ˙ i u ˙ i + I j k ϕ ˙ j m ϕ ˙ k m + J φ ˙ 2 + c θ 2 + + ( 2 + α ) A i j m n e i j e m n + 2 D i j m n e i j ν m n + 2 F i j m n r e i j χ m n r + + B i j m n ν i j ν m n + 2 G i j m n r ν i j χ m n r + C i j k m n r χ i j k χ m n r + + A i j φ , i φ , j + 2 α i j e i j φ + 2 β i j ν i j φ + 2 δ i j k χ i j k φ + ξ φ 2 d V .
Based on Equations (35) and (37), the function F receives the following from
F ( t ) = ( α 1 ) 0 t D K i j θ , i θ , j d V d τ 2 0 t D a i j e ˙ i j θ + b i j ν ˙ i j θ + c i j k χ ˙ i j k θ d V d τ 2 0 t D m θ φ ˙ d V d τ ,
so we can immediately deduce that
d F ( t ) d t = ( α 1 ) D K i j θ , i θ , j d V 2 D a i j e ˙ i j θ + b i j ν ˙ i j θ + c i j k χ ˙ i j k θ d V 2 D m θ φ ˙ d V .
Now, we can choose a sufficiently small positive constant β and a positive constant M 1 , computed in terms of β and constitutive coefficients, so that
D a i j e ˙ i j θ + b i j ν ˙ i j θ + c i j k χ ˙ i j k θ d V β D K i j θ , i θ , j d V + + M 1 D ρ u ˙ i u ˙ i + I j k ϕ ˙ j m ϕ ˙ k m + J φ ˙ 2 + c θ 2 d V .
On the other hand, can be determined the positive constant M 2 so that the following inequality takes place,
D m θ φ ˙ d V M 2 B J φ ˙ 2 + c θ 2 d V .
If we choose β 1 α and take into account the inequalities (39)–(41), we can find a positive constant M 3 so that
d F ( t ) d t M 3 D ρ u ˙ i u ˙ i + I j k ϕ ˙ j m ϕ ˙ k m + J φ ˙ 2 + c θ 2 d V .
Now, considering the inequalities (40)–(42) we can find a positive constant c so that for every t 0 , inequality occurs
d F ( t ) d t c F ( t ) ,
from where by direct integration, we are led to the inequality
F ( t ) F ( 0 ) e c t , t 0 ,
so that we can deduce that F ( t ) = 0 , for any t 0 , taking into account the null initial data and the expression of F ( t ) from (36).
Also, from (36), we deduce u ˙ i ( x , t ) = 0 , ϕ ˙ i j ( x , t ) = 0 , θ ˙ ( x , t ) = 0 , φ ˙ ( x , t ) = 0 , for any t 0 , x D .
Now the conclusion of Theorem 1 is easy to obtain, again considering that the initial data is null. □
Finally, we want to formulate and prove the main theorem of our work. We will show that it is impossible to localizate in a finite time any solutions of our mixed problem P consisting of Equations (2), (4), and (6)–(8), the null initial conditions (12) and homogeneous boundary conditions (17). In fact, we show that the only solution to this mixed problem that disappears in a finite time is the null solution.
Theorem 2.
Assume that the main hypotheses (a)–(c) are satisfied. If t 1 0 is a finite time and u i , ϕ i j , θ , φ is a solution of the mixed problem P , so that u i ( t ) = 0 , ϕ i j ( t ) = 0 , θ ( t ) = 0 , φ ( t ) = 0 , for any t t 1 , then we have u i ( t ) = 0 , ϕ i j ( t ) = 0 , θ ( t ) = 0 , φ ( t ) = 0 , for any t t 1 .
Proof. 
To obtain the conclusion of Theorem 2, we consider the mixed problem P in the case of homogeneous boundary conditions and considering the initial null data
u i ( t 1 , x ) = 0 , u ˙ i ( t 1 , x ) = 0 , ϕ i j ( t 1 , x ) = 0 , ϕ ˙ i j ( t 1 , x ) = 0 , θ ( t 1 , x ) = 0 , φ ( t 1 , x ) = 0 , φ ˙ ( t 1 , x ) = 0 , t 1 0 , x D .
If we take into account these null initial data, according to Theorem 1, the problem P admits only the null solution. □

5. Conclusions

This paper is dedicated to the linear mixed initial-boundary value problem in the context of the theory of thermoelasticity for bodies with voids that have a dipolar structure. After we put down some auxiliary identities, we approach the issue of the uniqueness for the solution for the mixed problem. Then, we approach the question of the possibility of locating in time of solution for this problem, namely, we have shown that it is impossible to locate in a finite time of the solution of mixed problem. Even if we combined the effect of the dipolar structure with the effect of the voids and with the thermal one, not were ensured the sufficient conditions for deformations vanish in a finite time.
If we take into account the fact that the porous media with a dipolar structure very accurately models the structure and behavior of the bones, we can deduce that a concrete application of our calculations is the following. If a bone is hit in one part accidentally, the time interval after which the pain caused by the blow completely disappears cannot be established exactly. However, the porous bodies with a dipolar structure also model other practical media, such as polystyrene, artificial porous materials, rocks, soils, and so on. A number of concrete applications of the theoretical results from our study can be made for these particular media. To this aim we must consider the isotropic case because in this situation the number of coefficients decreases significantly. Then, we will look for a specialized laboratory, of high performance, that can provide us with concrete values of the coefficients, for specific materials.

Author Contributions

All four authors conceived the framework and structured the whole manuscript, checked the results, and completed the revision of the paper. The authors have equally contributed to the elaboration of this manuscript. All authors have read and approved the final form of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Eringen, A.C. Theory of thermo-microstretch elastic solids. Int. J. Eng. Sci. 1990, 28, 1291–1301. [Google Scholar] [CrossRef]
  2. Eringen, A.C. Microcontinuum Field Theories; Springer: New York, NY, USA, 1999. [Google Scholar]
  3. Marin, M.; Nicaise, S. Existence and stability results for thermoelastic dipolar bodies with double porosity. Contin. Mech. Therm. 2016, 28, 1645–1657. [Google Scholar] [CrossRef]
  4. Marin, M.; Ellahi, R.; Chirila, A. On solutions of Saint-Venant’s problem for elastic dipolar bodies with voids. Carpathian J. Math. 2017, 33, 219–232. [Google Scholar]
  5. Mindlin, R.D. Micro-structure in linear elasticity. Arch. Ration. Mech. Anal. 1964, 16, 51–78. [Google Scholar] [CrossRef]
  6. Green, A.E.; Rivlin, R.S. Multipolar continuum mechanics. Arch. Ration. Mech. Anal. 1964, 17, 113–147. [Google Scholar] [CrossRef]
  7. Fried, E.; Gurtin, M.E. Thermomechanics of the interface between a body and its environment. Contin. Mech. Therm. 2007, 19, 253–271. [Google Scholar] [CrossRef]
  8. Nunziato, J.W.; Cowin, S.C. A nonlinear theory of materials with voids. Arch. Ration. Mech. Anal. 1979, 72, 175–201. [Google Scholar] [CrossRef]
  9. Cowin, S.C.; Nunziato, J.W. Linear elastic materials with voids. J. Elast. 1983, 13, 125–147. [Google Scholar] [CrossRef]
  10. Goodman, M.A.; Cowin, S.C. A continuum theory of granular material. Arch. Ration. Mech. Anal. 1971, 44, 249–266. [Google Scholar] [CrossRef]
  11. Marin, M.; Baleanu, D.; Vlase, S. Effect of microtemperatures for micropolar thermoelastic bodies. Struct. Eng. Mech. 2017, 61, 381–387. [Google Scholar] [CrossRef]
  12. Abbas, I.A. A GN model based upon two-temperature generalized thermoelastic theory in an unbounded medium with a spherical cavity. Appl. Math. Comput. 2014, 245, 108–115. [Google Scholar] [CrossRef]
  13. Marin, M.; Florea, O. On temporal behaviour of solutions in thermoelasticity of porous micropolar bodies. An. St. Univ. Ovidius Constanta 2014, 22, 169–188. [Google Scholar] [CrossRef] [Green Version]
  14. Othman, M.I.A.; Hasona, W.M.; Abd-Elaziz, E.M. Effect of Rotation on Micropolar Generalized Thermoelasticity with Two-Temperatures using a Dual-Phase-Lag Model. Can. J. Phys. 2014, 92, 149–158. [Google Scholar] [CrossRef]
  15. Marin, M. Lagrange identity method for microstretch thermoelastic materials. J. Math. Anal. Appl. 2010, 363, 275–286. [Google Scholar] [CrossRef] [Green Version]
  16. Marin, M.; Craciun, E.M. Uniqueness results for a boundary value problem in dipolar thermoelasticity to model composite materials. Compos. Part B Eng. 2017, 126, 27–37. [Google Scholar] [CrossRef]
  17. Serrin, J. The Initial Value Problem for the Navier–Stokes Equations. In Nonlinear Problems; University of Wisconsin Press: Madison, WI, USA, 1963; pp. 69–98. [Google Scholar]
  18. Ciarletta, M.; Scalia, A. Some Results in Linear Theory of Thermomicrostretch Elastic Solids. Meccanica 2004, 39, 191–206. [Google Scholar] [CrossRef]
  19. Knops, R.J.; Payne, L.E. On the Stability of Solutions of the Navier–Stokes Equations Backward in Time. Arch. Ration. Mech. Anal. 1968, 29, 331–335. [Google Scholar] [CrossRef]
  20. Galdi, G.P.; Straughan, B. Stability of Solutions of the Navier–Stokes Equations Backward in Time. Arch. Ration. Mech. Anal. 1988, 101, 107–114. [Google Scholar] [CrossRef]
  21. Payne, L.; Straughan, B. Improperly Posed and Nonstandard Problems for Parabolic Partial Differential Equations. In Elasticity: Mathematical Methods and Applications; Eason, R.W.O.G., Ed.; Ellis Horwood: Devon, UK, 1990; pp. 273–300. [Google Scholar]
  22. Ames, K.A.; Payne, L.E. Stabilizing Solutions of the Equations of Dynamical Linear Thermoelasticity Backward in Time. Stab. Appl. Anal. Contin. Media 1991, 1, 243–260. [Google Scholar]
  23. Ciarletta, M. On the Uniqueness and Continuous Dependence of Solutions in Dyna-mical Thermoelasticity Backward in Time. J. Therm. Streses 2002, 25, 969–984. [Google Scholar] [CrossRef]
  24. Ciarletta, M.; Chirita, S. Spatial Behavior in Linear Thermoelasticity Backward in Time. In Proceedings of the Fourth International Congress on Thermal Stresses, Osaka, Japan, 8–11 June 2001; Chao, C.K., Lin, C.Y., Eds.; pp. 485–488. [Google Scholar]
  25. Quintanilla, R. Impossibility of Localization in Linear Thermoelasticity with Voids. Mech. Res. Commun. 2007, 34, 522–527. [Google Scholar] [CrossRef]
  26. Green, A.E.; Naghdi, P.M. On Undamped Heat Waves in an Elastic Solid. J. Therm. Stress. 1992, 15, 253–264. [Google Scholar] [CrossRef]
  27. Green, A.E.; Naghdi, P.M. Thermoelasticity without Energy Dissipation. J. Elast. 1993, 31, 189–208. [Google Scholar] [CrossRef]
  28. Iovane, G.; Passarella, F. Saint-Venant’s Principle in Dynamic Porous Thermoelastic Media with Memory for Heat Flux. J. Therm Stress. 2004, 27, 983–999. [Google Scholar] [CrossRef]
  29. Passarella, F.; Tibullo, V. Some Results in Linear Theory of Thermoelasticity Backward in Time for Microstretch Materials. J. Therm. Stress. 2010, 33, 559–576. [Google Scholar] [CrossRef]
  30. Chirila, A.; Marin, M. The theory of generalized thermoelasticity with fractional order strain for dipolar materials with double porosity. J. Mater. Sci. 2018, 53, 3470–3482. [Google Scholar] [CrossRef]
  31. Marin, M.; Agarwal, R.P.; Codarcea, L. A mathematical model for three-phase-lag dipolar thermoelastic bodies. J. Ineq. Appl. 2017, 2017, 109. [Google Scholar] [CrossRef] [Green Version]
  32. Quintanilla, R.; Straughan, B. Energy Bounds for Some Non-standard Problems in Thermoelasticity. Proc. R. Soc. Lond. A 2005, 461, 1–15. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Marin, M.; Abbas, I.; Vlase, S.; Craciun, E.M. A Study of Deformations in a Thermoelastic Dipolar Body with Voids. Symmetry 2020, 12, 267. https://doi.org/10.3390/sym12020267

AMA Style

Marin M, Abbas I, Vlase S, Craciun EM. A Study of Deformations in a Thermoelastic Dipolar Body with Voids. Symmetry. 2020; 12(2):267. https://doi.org/10.3390/sym12020267

Chicago/Turabian Style

Marin, Marin, Ibrahim Abbas, Sorin Vlase, and Eduard M. Craciun. 2020. "A Study of Deformations in a Thermoelastic Dipolar Body with Voids" Symmetry 12, no. 2: 267. https://doi.org/10.3390/sym12020267

APA Style

Marin, M., Abbas, I., Vlase, S., & Craciun, E. M. (2020). A Study of Deformations in a Thermoelastic Dipolar Body with Voids. Symmetry, 12(2), 267. https://doi.org/10.3390/sym12020267

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop